WRF后处理数据用fortran或者NCL读取

robin_ouc

由于投影方式以及垂直坐标类型的问题，想读取wrfout用ARWpost处理之后的二进制dat文件。按照grads存储数据的顺序，分别用fortran和NCL读了前几个数据，结果数量级明显不对。
还有，ctl里面有：
pdef  150 180 lcc  33.712  121.190   75.500   90.500  60.00000  30.00000  123.00000  10000.000  10000.000
xdef  431 linear  111.24580   0.04504505
ydef  391 linear   24.80210   0.04504505
pdef我知道是数据差值用的，不知道会不会改变数据的储存格式？一个时次一层水平数据是150*180个还是431*391个？
求大神指教！

Mid_Farmer

pdef是投影数据的说明 一个时次一层水平数据是150*180个 xydef定义的投影背景范围 与dat保存格式无关 与grads作图有关 我是这么理解的

robin_ouc

Mid_Farmer 发表于 2015-12-14 12:31
pdef是投影数据的说明 一个时次一层水平数据是150*180个 xydef定义的投影背景范围 与dat保存格式无关 与gra ...

非常感谢！但是我读x方向风速，读出来的值都很离谱：

                               
登录/注册后可看大图

，读其他变量也都差不多。不知您知不知道是什么原因？

robin_ouc

Mid_Farmer 发表于 2015-12-14 12:31
pdef是投影数据的说明 一个时次一层水平数据是150*180个 xydef定义的投影背景范围 与dat保存格式无关 与gra ...

截的图在这里：http://bbs.06climate.com/forum.p ... mp;page=1#pid563568

Mid_Farmer

robin_ouc 发表于 2015-12-14 15:13
截的图在这里：http://bbs.06climate.com/forum.php?mod=viewthread&tid=40581&page=1#pid563568

不知道你的哪个截图是风速的截图 而且没有看到你的程序 不清楚是否正确读取了 不过需要注意的是：grads写数据是按经度 纬度 高度 时间来循环的 所以读取的时候要按逆序来 底层高度的所有xy范围内的所有变量读完了才能读高层的其他变量 不知道你是不是这样的读取顺序 有问题多交流哈~~