<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_144907">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_144908"><br id="yui_3_16_0_ym19_1_1508231998965_144909"></div><div id="yui_3_16_0_ym19_1_1508231998965_144908"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_144910">The fortran code seems working correctly to read those geography data,</div><div id="yui_3_16_0_ym19_1_1508231998965_144911"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_144911">e.g. 00001-00120.00001-00120, ... and ... 02041-02160.00961-01080,</div><div id="yui_3_16_0_ym19_1_1508231998965_144912"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_144912">total 162 files for 10m resolution, and 648 files for 30s resolution globally,</div><div id="yui_3_16_0_ym19_1_1508231998965_144913"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_144913">with "halo" (tile_bdr=3 in the index file) in each file (or tile).</div><div id="yui_3_16_0_ym19_1_1508231998965_144914"><br id="yui_3_16_0_ym19_1_1508231998965_144915"></div><div id="yui_3_16_0_ym19_1_1508231998965_144914"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_144916">And, the fortran code wrote out the data in binary format as:</div><div id="yui_3_16_0_ym19_1_1508231998965_144917"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_144917">=======================</div><div id="yui_3_16_0_ym19_1_1508231998965_144918"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_144918">open(23,file=file(3),status='unknown',form='unformatted'</div><div id="yui_3_16_0_ym19_1_1508231998965_144919">+ ,convert='little_endian')<br></div><div id="yui_3_16_0_ym19_1_1508231998965_144920"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_144920">=======================</div><div id="yui_3_16_0_ym19_1_1508231998965_144921"><br id="yui_3_16_0_ym19_1_1508231998965_144922"></div><div id="yui_3_16_0_ym19_1_1508231998965_144923"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_144923">The ctl file for 10m globally is:</div><div id="yui_3_16_0_ym19_1_1508231998965_144924"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_144924">=======================</div><div id="yui_3_16_0_ym19_1_1508231998965_144925"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_144925">xdef 2160 linear 180.083333 0.16666667</div><div id="yui_3_16_0_ym19_1_1508231998965_144926"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_144926">ydef 1080 linear -89.916667 0.16666667</div><div id="yui_3_16_0_ym19_1_1508231998965_144927"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_144927">=======================</div><div id="yui_3_16_0_ym19_1_1508231998965_144928"><br id="yui_3_16_0_ym19_1_1508231998965_144929"></div><div id="yui_3_16_0_ym19_1_1508231998965_144930"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_144930">The ctl file for 30s regionally:</div><div id="yui_3_16_0_ym19_1_1508231998965_144931"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_144931">=======================</div><div id="yui_3_16_0_ym19_1_1508231998965_144932"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_144932">xdef 4800 linear 100.00416 0.00833333</div><div id="yui_3_16_0_ym19_1_1508231998965_144933"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_144933">ydef 3600 linear 10.00416 0.00833333</div><div id="yui_3_16_0_ym19_1_1508231998965_144934"><br></div><div id="yui_3_16_0_ym19_1_1508231998965_144934">=======================</div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_144935"><br id="yui_3_16_0_ym19_1_1508231998965_144936"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_144935"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_144937">However, the pattern between WRF hgt and topo_10m looks okay,</div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_144938"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_144938">but the WRF hgt shifted eastward of the topo_30s data.</div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_144939"><br id="yui_3_16_0_ym19_1_1508231998965_144940"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_144939"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_144941">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_144942"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_144942">both 10m and 30s should be wrong simultaneously.</div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_144943"><br id="yui_3_16_0_ym19_1_1508231998965_144944"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_144943"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_144945">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_144946"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_144946">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_144947"><br id="yui_3_16_0_ym19_1_1508231998965_144948"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_144947"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_144949">Have anyone dealt with this problem?</div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_144950"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_144950">Thanks for reading this email.</div><div dir="ltr" id="yui_3_16_0_ym19_1_1508231998965_144951"><br id="yui_3_16_0_ym19_1_1508231998965_144952"></div></div></body></html>