[gradsusr] Terrain data of WRF and ARWpost for GrADS (correct email format)

Mark Wang markymwang at yahoo.com
Fri Oct 20 12:31:17 EDT 2017


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.16666667
ydef 1080 linear -89.916667 0.16666667
=======================

The ctl file for 30s regionally:
=======================
xdef 4800 linear 100.00416 0.00833333
ydef 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/0ba981aa/attachment.html 


More information about the gradsusr mailing list