请教高手，我这用1*1fnl资料画850hpa比湿的gs文件正确不？

bluoshe

'reinit'
'clear'
'open E:\grads\fnl\fnl.ctl'

'set mpdset cnworld'
'set map 1 1 1'
'set csmooth on'
'set lon 105  125'
'set lat 18  30'
'set xlopts 1 6 0.16'
'set ylopts 1 6 0.16'
tt=1
while(tt<=12)
'set grads off'
'set grid off'
'set t 'tt''
'set lev 850'
'define prs=lev'
'define tc=TMPprs-273.16'
'define es=(6.112*exp((17.67*tc)/(tc+243.5)))'
'define q=(0.62197*es/(prs-es))'
'define tv=TMPprs*(1+0.61*q)'
'define q=tvrh2q(tv,RHprs)'
'set ccolor 1'
'set cint 3'
'set cthick 8'
'd q'
'q dims'
rec=sublin(result,5)recs=subwrd(rec,6)
'draw title  'recs'''printim E:\grads\tu\q\q'tt'850hpa.jpg x800 y600 white'
'c'
tt=tt+1
endwhile
;
运行以上gs以后，grads提示以下问题：



画出来的第一张图是这样的：



请教各位高手，gs是不是有问题？问题在哪？为什么我画出来的图，比湿的量级这么小？

bluoshe

river 发表于 2013-8-24 08:20

                               
登录/注册后可看大图

公式是没有问题的，用混合比的值近似的等于饱和比湿的值，然后计算出虚温的值，最后带入函数计算比湿。那你 ...

恩，有道理，不过这个计算公式我是在网上找的，然后自己编的gs，至于单位，比湿的一般都是g/g或g/kg，这套公式计算的单位应该是g/Kg吧

bluoshe

有人给我解答一下吗？我这计算比湿的公式正确吗？

river

公式是没有问题的，用混合比的值近似的等于饱和比湿的值，然后计算出虚温的值，最后带入函数计算比湿。那你知道你计算出来的比湿单位是什么吗？还有你和哪的结果对比说你的小，那别人计算出来的结果单位是什么？如果没有单位，说什么都是没有意义的

river

bluoshe 发表于 2013-8-26 09:12

                               
登录/注册后可看大图

恩，有道理，不过这个计算公式我是在网上找的，然后自己编的gs，至于单位，比湿的一般都是g/g或g/kg，这套 ...

那要看你带进去的各个物理量的单位是什么。没有应该一说，你得彻底搞清楚是g/g还是g/kg，要是确定不了，那你这个结果就是没有意义的

yaya83

帮大忙了 谢谢

jiajian0512

你的这个比湿算出来是kg/kg吗？

zaozedi

单位是g/g，tvrh2q(tvexpr,rhexpr)：

根据给出的虚温以及相对湿度, tvrh2q将返回比湿q,单位g/g.

tvexpr     表示虚温的GrADS表达式, 单位:K.
rhexpr     表示相对湿度的GrADS表达式, 单位: %(从0到100).

此函数仅适用于格点数据.

nicocle

我用这个GS做。。图是空白的，只有行政边界，GRADS显示Using fixed pressure level 850 mb
是怎么回事？

jiaohe1990

牛牛牛牛牛牛牛牛牛牛牛牛牛牛牛牛牛牛牛牛牛牛牛牛