请教关于“比湿”的问题

abd

比湿与相对湿度有什么区别啊？有了相对湿度后比湿怎么算啊？也就是确切计算公式怎么写？
目前我知道的算法是：
'define prs=pressure'    'define tc=tk-273.16'
'define es=(6.112*exp((17.67*tc)/(tc+243.5)))'
'define q=(0.62197*es/(prs-es))'
'define tv=tk*(1+0.61*q)'
'define q=tvrh2q(tv,rh)'
'd q'
其中prs=pressure(气压)，tk,tc为温度（单位分别k,和摄氏度），tv为虚温，q为比湿，tvrh2q为Grads自带的函数但是我用WRF模式模拟出来的气压，温度，相对湿度等以上的物理量来计算出比湿时gs文件如下：'reinit'
'open /gpfsTMP/ywang/abu/aar/ARWpost/t20c-19-d4-rh.ctl'
****************************** bishi
'set gxout fwrite'
'set fwrite /gpfsTMP/ywang/abu/aar/ARWpost/t20c-19/cross/Bishi.dat'
t=1
while(t<=18)
'set t 't
nz=1
while(nz<=66)
'set lev 'nz''
'define prs=pressure'
'define tc=tk-273.16'
'define es=(6.112*exp((17.67*tc)/(tc+243.5)))'
'define q=(0.62197*es/(prs-es))'
'define tv=tk*(1+0.61*q)'
'define q=tvrh2q(tv,rh)'
'd q'
nz=nz+1
endwhile
t=t+1
endwhile
'disable fwrite'

但是输出数据时，刚开始貌似挺好的，但是从第5层开始出现这样的问题
Using fixed pressure level 1 mb
Using fixed pressure level 2 mb
Using fixed pressure level 3 mb
Using fixed pressure level 4 mb
Using fixed pressure level 5 mb
Warning from TVRH2Q: Convergence failed for 32 grid points
Using fixed pressure level 6 mb
Warning from TVRH2Q: Convergence failed for 8033 grid points
Using fixed pressure level 7 mb
Warning from TVRH2Q: Convergence failed for 25934 grid points

。。。。。
Using fixed pressure level 24 mb
Warning from TVRH2Q: Convergence failed for 5611 grid points
Using fixed pressure level 25 mb
Warning from TVRH2Q: Convergence failed for 390 grid points
Using fixed pressure level 26 mb
Using fixed pressure level 27 mb
Using fixed pressure level 28 mb
Using fixed pressure level 29 mb
Using fixed pressure level 30 mb
Using fixed pressure level 31 mb

。。。。。

到了26层后又貌似正常，真不明白是什么问题，请大家指点一下!非常感谢！

鬼王

我不知道怎么弄，哈哈，我是新人

river

是不是有缺测值什么的造成的，你的算法没有什么问题

abd

river 发表于 2012-10-14 18:22

                               
登录/注册后可看大图

是不是有缺测值什么的造成的，你的算法没有什么问题

缺测值？那个是WRF模拟出来的数据啊！也会有缺测值吗？

river

abd 发表于 2012-10-14 18:47

                               
登录/注册后可看大图

缺测值？那个是WRF模拟出来的数据啊！也会有缺测值吗？

啊哦，我没弄过这个，只是猜测啊。你问问版主清风吧

abd

river 发表于 2012-10-14 18:53

                               
登录/注册后可看大图

啊哦，我没弄过这个，只是猜测啊。你问问版主清风吧

嗯，好的，谢谢！

river

比湿计算好像不太对啊，应该是q=相对湿度*(0.62197*es/(prs-0.378*es))/100.  。

abd

river 发表于 2012-10-14 19:04

                               
登录/注册后可看大图

比湿计算好像不太对啊，应该是q=相对湿度*(0.62197*es/(prs-0.378*es))/100.  。

你确定是这个公式吗？
那我试着这样画出来看看，谢谢你的回复哦！

abd

river 发表于 2012-10-14 19:04

                               
登录/注册后可看大图

比湿计算好像不太对啊，应该是q=相对湿度*(0.62197*es/(prs-0.378*es))/100.  。

你告诉我的这个算法好像不太对哦，用你那个公式画出来的比湿图为：

而用我原来的公式算出来的为：

而且，我用的那个tvrh2q(tv,rh) 这个函数式grads自带的用相对湿度和虚温算比湿的函数，但核心算法就不知道了，更郁闷的事出现那种诡异的错误，。。。。。

river

abd 发表于 2012-10-14 19:45

                               
登录/注册后可看大图

你告诉我的这个算法好像不太对哦，用你那个公式画出来的比湿图为：

而用我原来的公式算出来的为：

不太懂你的gs，下面的那个q是比湿我知道，你第一次计算的q是什么？是不是这个变量定义成一样的名字的缘故