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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3263|回复: 2

[分享资料] 请各位大神帮我看看利用grib2这样求比湿有没问题?

[复制链接]

新浪微博达人勋

发表于 2018-12-9 18:06:02 | 显示全部楼层 |阅读模式

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

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

x
比如说我是想求850hPa上的比湿,gs如下:
'set lev 850'
'define prs=lev'
'define tc=TMPprs-273.16'
'define es=(6.112*exp((17.67*tc)/(tc+243.5)))'
'define u=(0.62197*es/(prs-es))'
'define tv=TMPprs*(1+0.61*u)'
'define q=tvrh2q(tv,RHprs)*1000'
'd q'
grib2中数据描述如下:
TMPprs    26,100  0,0,0 ** (1000 975 950 925 900.. 70 50 30 20 10) Temperature [K]  (每层的温度,单位是开尔文)
RHprs    21,100  0,1,1 ** (1000 975 950 925 900.. 300 250 200 150 100) Relative Humidity [%](每层的相对湿度)

tvrh2q计算的结果是g/g,我后面乘以1000,换成g/kg.
这是借鉴http://bbs.06climate.com/forum.p ... hlight=%B1%C8%CA%AA帖子的计算方法写的。
后面发现http://bbs.06climate.com/forum.p ... id=16776&page=2帖子也是用同样的方法。

有人告诉我这么算,对不对?我始终没找到对应计算虚温的公式。
运行结果:
45.png
出图结果
00Z18APR2012850hPa-q.png
希望大神们说说这算法对还是错哈?请不吝赐教。


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

新浪微博达人勋

 楼主| 发表于 2018-12-9 18:52:57 | 显示全部楼层
自己顶顶!!
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2018-12-10 09:19:04 | 显示全部楼层
你可以下个ncep2.5°的资料,算一下对比看看
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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