爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 64004|回复: 157

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

  [复制链接]

新浪微博达人勋

发表于 2013-11-7 12:02:32 | 显示全部楼层 |阅读模式

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

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 海豚∩_∩浩 于 2014-5-20 16:23 编辑

       学校运动会,无心学术,看到NCL调用FROTRAN,小尝试了一下,不多说,如下:  

                               
登录/注册后可看大图
       FNL(girb)资料用NCL提取变量后,有两种方法可以计算水汽通量:        一个是可以直接在ncl脚本中通过算术公式(或者使用ncl内部部分函数,该函数能求得比湿)计算。这个之前NCL论坛里有相关的一个帖子;
         另一个是通过在ncl脚本中调用fortran程序计算(之前论坛里有过一个调用fortran的帖子,受其启发,故自己尝试了一遍,个人觉得对于比较大的程序还是使用该方法好,因为NCL仅仅是解释性语言,而Fortran和c是编译语言)

        下面来尝试一下,先上张图:
850hpa vapour_flux(0 18).png

                               
登录/注册后可看大图
;*******************************************************************************************
;计算水汽通量,提供脚本中的部分语句
;*********************************************************************************************
xx=41
yy=41
zz=1
tt=13
    MY::vaporflux(p,t,rhprs,u,v,qu,qv,mm,xx,yy,zz,tt,g)

                               
登录/注册后可看大图

;=====================================================================
本帖主要是尝试了第二种方法。附近包括:NCL脚本,f90程序,相对于的stub文件、.so文件(可能对于不同版本的NCL,需要自己通过WRAPIT重新生成.so文件,我的是6.1.0版本的)。


当然FORTRAN程序编得不咋滴,缺测值处理有待提高,学弱弱只是想尝试一下想法,大家见谅哈!

                               
登录/注册后可看大图

(ps:只要一点贡献得了,毕竟也花了我一段时间来摸索的,付出需要一点回报哒) vaporflux.stub (218 Bytes, 下载次数: 213)

评分

参与人数 1金钱 +4 贡献 +1 收起 理由
river + 4 + 1

查看全部评分

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

新浪微博达人勋

发表于 2013-11-7 12:39:42 | 显示全部楼层
顶,浩哥厉害!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-7 13:35:34 | 显示全部楼层
浩班长真大神也,顶一个~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-7 14:36:02 | 显示全部楼层
有点小贵,不过为了鼓励你,我来购买了,希望你分享更多成果。谢谢!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-7 15:46:13 | 显示全部楼层
浩班长无敌
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-7 16:20:56 | 显示全部楼层
顶,我收了!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-7 16:40:03 | 显示全部楼层
感谢分享            
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-7 16:46:44 | 显示全部楼层
顶一个!!!!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-7 17:06:07 | 显示全部楼层
感谢分享!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-7 19:10:43 | 显示全部楼层
呵呵,支持一下,虽然我不太懂NCL
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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