爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索

[经验总结] NCL画水汽通量+调用fortran程序

  [复制链接]
发表于 2014-10-13 21:30:34 | 显示全部楼层
给群主个赞,感觉ncl配上fortran就谁也挡不住了呢~~~
密码修改失败请联系微信:mofangbao
发表于 2014-10-27 09:57:36 | 显示全部楼层
你好,最近在做WRF初始场修改工作,想通过NCL调用外部fortran程序对输入的fnl资料本身进行滤波,不另外输出文件。看到网上说NCL下grib格式只可读,不可写,那是不是这种想法就不可行了?能不能先对fnl处理,生成netcdf格式文件,再把它用NCL转成grib格式的?
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-10-27 10:42:04 | 显示全部楼层
schliezer 发表于 2014-10-27 09:57
你好,最近在做WRF初始场修改工作,想通过NCL调用外部fortran程序对输入的fnl资料本身进行滤波,不另外输出 ...

NETCDF转GRIB?我倒是没见过。WRF初始场不能是NC格式的数据么?
密码修改失败请联系微信:mofangbao
发表于 2014-10-27 14:48:06 | 显示全部楼层
海豚∩_∩浩 发表于 2014-10-27 10:42
NETCDF转GRIB?我倒是没见过。WRF初始场不能是NC格式的数据么?

我用的fnl资料,是grib格式,但是NCL貌似只能读这种格式,不知道你有没有什么办法可以直接读写这种格式的数据,即使不是NCL,fortran不行吧?
密码修改失败请联系微信:mofangbao
发表于 2014-10-28 15:29:18 | 显示全部楼层
感谢分享,认真学习下~
密码修改失败请联系微信:mofangbao
发表于 2014-11-16 23:47:40 | 显示全部楼层
         我想请问一下楼主,程序里面做这样的判断是什么意思啊?
    if(t(ix,iy,it)>=9.96921e+36)then
            t(ix,iy,it)=0.0
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-11-17 09:05:09 | 显示全部楼层
甜果果 发表于 2014-11-16 23:47
我想请问一下楼主,程序里面做这样的判断是什么意思啊?
    if(t(ix,iy,it)>=9.96921e+36)then
       ...

嗯,这个就是将缺测值改为0.
密码修改失败请联系微信:mofangbao
发表于 2014-11-17 10:11:18 | 显示全部楼层
向大神学习,谢谢
密码修改失败请联系微信:mofangbao
发表于 2014-11-18 08:54:43 | 显示全部楼层
感谢分享!
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

发表于 2014-11-18 20:44:34 | 显示全部楼层

同求NCL水汽通量散度函数,请问大神是哪一个呢
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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