[gradsusr] Terrain data of WRF and ARWpost for GrADS

[Mark Wang]

I am comparing the terrain input data for WRF (e.g. topo_10m), and WRF output terrain height (hgt).
The fortran code seems working correctly to read those geography data,e.g. 00001-00120.00001-00120, ... and ... 02041-02160.00961-01080,total 162 files for 10m resolution, and 648 files for 30s resolution globally,with "halo" (tile_bdr=3 in the index file) in each file (or tile).
And, the fortran code wrote out the data in binary format as:=======================open(23,file=file(3),status='unknown',form='unformatted'+ ,convert='little_endian')=======================
The ctl file for 10m globally is:=======================xdef 2160 linear 180.083333 0.16666667ydef 1080 linear -89.916667 0.16666667=======================
The ctl file for 30s regionally:=======================xdef 4800 linear 100.00416 0.00833333ydef 3600 linear 10.00416 0.00833333=======================
However, the pattern between WRF hgt and topo_10m looks okay,but the WRF hgt shifted eastward of the topo_30s data.
If the Fortran code is wrong in I/O, or format of ctl is wrong,both 10m and 30s should be wrong simultaneously.
In fact, I did make a mistake of starting lat/lon of x/y def in ctl.But after I corrected the error, the shifted pattern of 30s data remained.
Have anyone dealt with this problem?Thanks for reading this email.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20171020/a0b0db06/attachment.html