
<html><head></head><body><div style="color:#000; background-color:#fff; font-family:Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:16px"><div id="yui_3_16_0_ym19_1_1508231998965_138551">I am comparing the terrain input data for WRF (e.g. topo_10m), and WRF output terrain height (hgt).</div><div id="yui_3_16_0_ym19_1_1508231998965_138551"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_138645">The fortran code seems working correctly to read those geography data,</div><div id="yui_3_16_0_ym19_1_1508231998965_138646">e.g. 00001-00120.00001-00120, ... and ... 02041-02160.00961-01080,</div><div id="yui_3_16_0_ym19_1_1508231998965_138647">total 162 files for 10m resolution, and 648 files for 30s resolution globally,</div><div id="yui_3_16_0_ym19_1_1508231998965_138648">with "halo" (tile_bdr=3 in the index file) in each file (or tile).</div><div id="yui_3_16_0_ym19_1_1508231998965_138649"><br id="yui_3_16_0_ym19_1_1508231998965_138650"></div><div id="yui_3_16_0_ym19_1_1508231998965_138651">And, the fortran code wrote out the data in binary format as:</div><div id="yui_3_16_0_ym19_1_1508231998965_138652">=======================</div><div id="yui_3_16_0_ym19_1_1508231998965_138653">open(23,file=file(3),status='unknown',form='unformatted'</div><div id="yui_3_16_0_ym19_1_1508231998965_138654">+ ,convert='little_endian')</div><div id="yui_3_16_0_ym19_1_1508231998965_138655">=======================</div><div id="yui_3_16_0_ym19_1_1508231998965_138656"><br id="yui_3_16_0_ym19_1_1508231998965_138657"></div><div id="yui_3_16_0_ym19_1_1508231998965_138658">The ctl file for 10m globally is:</div><div id="yui_3_16_0_ym19_1_1508231998965_138659">=======================</div><div id="yui_3_16_0_ym19_1_1508231998965_138660">xdef 2160 linear 180.083333 0.16666667</div><div id="yui_3_16_0_ym19_1_1508231998965_138661">ydef 1080 linear -89.916667 0.16666667</div><div id="yui_3_16_0_ym19_1_1508231998965_138662">=======================</div><div id="yui_3_16_0_ym19_1_1508231998965_138663"><br id="yui_3_16_0_ym19_1_1508231998965_138664"></div><div id="yui_3_16_0_ym19_1_1508231998965_138665">The ctl file for 30s regionally:</div><div id="yui_3_16_0_ym19_1_1508231998965_138666">=======================</div><div id="yui_3_16_0_ym19_1_1508231998965_138667">xdef 4800 linear 100.00416 0.00833333</div><div id="yui_3_16_0_ym19_1_1508231998965_138668">ydef 3600 linear 10.00416 0.00833333</div><div id="yui_3_16_0_ym19_1_1508231998965_138669">=======================</div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_138670"><br id="yui_3_16_0_ym19_1_1508231998965_138671"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_138670">However, the pattern between WRF hgt and topo_10m looks okay,</div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_138670">but the WRF hgt shifted eastward of the topo_30s data.</div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_138670"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_138670">If the Fortran code is wrong in I/O, or format of ctl is wrong,</div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_138670">both 10m and 30s should be wrong simultaneously.</div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_138670"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_138670">In fact, I did make a mistake of starting lat/lon of x/y def in ctl.</div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_138670">But after I corrected the error, the shifted pattern of 30s data remained.</div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_138670"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_138670">Have anyone dealt with this problem?</div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_138670">Thanks for reading this email.</div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_138670"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_138670"><br></div></div></body></html>