GRADS和FORTRAN将NCEP数据写成MICAPS第四类数据

yanglei_nuist · 发表于 2013-11-20 21:59:18

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 yanglei_nuist 于 2013-11-20 22:01 编辑

最近学习GRADS处理NCEP数据，现在想和大家分享一点成果。目的就是同时使用GRADS和FORTRAN将NCEP数据写成MICAPS第四类数据。（同时处理多个时次多个高度的水汽通量散度）1、GRADS软件读取NCEP1°×1° FNL数据，输出为FORTRAN可以读取的grd文件，这里需要注意该如何输出文件，建议每个FNL对应一个GRD，之所以为什么，我也没搞懂！set fwrite ....路径、set gxout fwrite和disable fwrite放在哪个循环处都很讲究，具体请看write.qdiv.day文件！水汽通量的计算也是参考之前气象家园的帖子计算的。
2、FORTRAN通过循环设置可以读取多个grd格式数据，然后按照MICAPS第四类数据的要求写成相应的文件，这里需要注意的当然就是open（file）时该如何设置，因为读取的是grd文件，所以
open(1,file='qdiv.grd',form='binary',recl=NX*NY*NZ*4)中的recl设置也要和1中输出的时候对应。这里recl应为NX*NY*NZ*4。详见write.qdiv.day.f90程序
3、最后就可以把生成的文件拖到MICAPS中，显现图形了。附件里有使用GRADS画的图和MICAPS呈现最后结果得到的图！
QQ图片20131120215519.jpg

从图中对比可以看到两幅图是很接近的！因为初学写这个东西花了几个晚上，想和大家分享下，如果大家有在用GRADS和MICAPS、NCEP数据的同仁，大家可以互相交流下，请大家多多指正！

write.qdiv.day.f90 (4.84 KB, 下载次数: 176)

二斤幸福 · 发表于 2013-11-20 22:03:45

雋雨 · 发表于 2013-11-20 22:10:33

虽然我看不懂，但我知道楼主说的很高端！

FLYING_YGY · 发表于 2013-11-20 22:12:55

1024

捷然不同 · 发表于 2013-11-20 22:15:11

1024

mofangbao · 发表于 2013-11-20 23:10:51

提问的那个后来发现是什么原因吗

river · 发表于 2013-11-21 08:21:25

感谢分享，很不错

善人/jw · 发表于 2013-11-21 08:49:18

很好的总结啊，感谢分享
话说，现在看不了谁加了多少分嘛？@mofangbao

mofangbao · 发表于 2013-11-21 09:17:34

善人/jw 发表于 2013-11-21 08:49

登录/注册后可看大图

很好的总结啊，感谢分享
话说，现在看不了谁加了多少分嘛？@mofangbao

看的了啊。。。你啥眼神。。。

yanglei_nuist · 发表于 2013-11-21 10:11:10

mofangbao 发表于 2013-11-20 23:10 提问的那个后来发现是什么原因吗

之前的那个还是没找到原因，后来合计这么每个fnl对应一个输出一个grd文件应该可以，试了一下还真可以。我估计之前的那个fortran里读语句中rec还是应该设置为一个定值。

		自动登录	找回密码
密码			立即注册

GRADS和FORTRAN将NCEP数据写成MICAPS第四类数据

登录后查看更多精彩内容~

评分

点评

点评