请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 21230|回复: 37

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

  [复制链接]

新浪微博达人勋

发表于 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, 下载次数: 173)

评分

参与人数 5金钱 +70 贡献 +16 收起 理由
kongfeng0824 + 10
善人/jw + 20 + 2 赞一个!
传说中的谁 + 10 + 4
river + 10 + 2 赞一个!
mofangbao + 20 + 8

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-20 22:03:45 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-20 22:10:33 | 显示全部楼层
虽然我看不懂,但我知道楼主说的很高端!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-20 22:12:55 | 显示全部楼层
      1024  

点评

咳咳。。。这里不是某农业论坛,正经回复哈  发表于 2013-11-21 08:41
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-20 22:15:11 | 显示全部楼层
      1024

点评

咳咳。。。这里不是某农业论坛,正经回复哈  发表于 2013-11-21 08:41
密码修改失败请联系微信:mofangbao

新浪微博达人勋

0
早起挑战累计收入
发表于 2013-11-20 23:10:51 | 显示全部楼层
提问的那个后来发现是什么原因吗
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-21 08:21:25 | 显示全部楼层
感谢分享,很不错
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-21 08:49:18 | 显示全部楼层
很好的总结啊,感谢分享
话说,现在看不了谁加了多少分嘛?@mofangbao
密码修改失败请联系微信:mofangbao

新浪微博达人勋

0
早起挑战累计收入
发表于 2013-11-21 09:17:34 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-11-21 10:11:10 来自手机 | 显示全部楼层
mofangbao 发表于 2013-11-20 23:10 提问的那个后来发现是什么原因吗

之前的那个还是没找到原因,后来合计这么每个fnl对应一个输出一个grd文件应该可以,试了一下还真可以。我估计之前的那个fortran里读语句中rec还是应该设置为一个定值。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表