求助水汽通量散度的量级问题

269984288

这是自己编的fnl1*1的求水汽通量散度的gs文件:

'reinit'
'open d:\fnl\fnl_2011051100.ctl'
'set lev 850'
'set lon 99.99 120.01'
'set lat 17.99 30.01'
'set mpdset cnworld'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set grads off'
'set timelab off'
'set grid off'
'enable print d:\fnl\051100_850_qfdiv1.gmf'

'define g=9.8'
'define es=6.112*exp(17.67*(TMPprs-273.16)/(TMPprs-29.65))'
'define qs=0.62197*es/(lev-0.378*es)'
'define q=qs*RHprs/100'
'define qhdivg=hdivg(q*ugrdprs/g,q*vgrdprs/g)*1000'

'set clskip 2'
'set clopts 1 1 0.12'
'set gxout shaded'
'd qhdivg'
'set gxout contour'
'd qhdivg'

'print'
'disable print'
;


得出来的图为（截图）：

发现量级是10-5，但看到很多文献及帖子都说量级是10-6~10-7之间，单位为g/(s•cm2•hPa)，
所以不知道我编的gs是不是有问题，最后要乘以多少才是标准单位，求助求助~~~

269984288

bluexlm 发表于 2011-8-24 13:01

                               
登录/注册后可看大图

100要除，因为fnl资料里rh是几十的量级；1000要乘，因为比湿要求g/kg，最后量级是e-5没错，有些文献是这个量 ...

哇哇~~~谢谢啦~~~热心人士啊~~~{:3_48:}{:3_62:}，还总结了那么多先，膜拜……

风行浪尖

色彩丰富，我喜欢

269984288

风行浪尖 发表于 2011-8-4 22:38

                               
登录/注册后可看大图

色彩丰富，我喜欢

不是色不色彩的问题，是量级对不对的问题，汗啊~~~~{:3_62:}

尽头的尽头

'define qhdivg=hdivg(q*ugrdprs/g,q*vgrdprs/g)*1000'
这句后面的1000是干什么的？

269984288

尽头的尽头 发表于 2011-8-5 01:26

                               
登录/注册后可看大图

'define qhdivg=hdivg(q*ugrdprs/g,q*vgrdprs/g)*1000'
这句后面的1000是干什么的？

水汽通量散度中的q一般取克/千克，而这里q的单位为：克/克，所以乘以1000，我是这么认为的，不懂对不对

晴天雨天

刚,把资料发给我,让我试试.
luojyemail@163.com

269984288

晴天雨天 发表于 2011-8-5 10:26

                               
登录/注册后可看大图

刚,把资料发给我,让我试试.

哇哈哈~~~是罗副吧~~~ 我是小戴，文件传过去了~~~哇哈哈，麻烦了~~~

mofangbao

269984288 发表于 2011-8-5 11:03

                               
登录/注册后可看大图

哇哈哈~~~是罗副吧~~~ 我是小戴，文件传过去了~~~哇哈哈，麻烦了~~~

你们竟然认识哈
ps:送你的高亮卡这么快就用上了哈

269984288

mofangbao 发表于 2011-8-5 12:14

                               
登录/注册后可看大图

你们竟然认识哈
ps:送你的高亮卡这么快就用上了哈

非常感谢哈~~~我只想尽快找到答案而已~~~{:3_51:}

尽头的尽头

269984288 发表于 2011-8-5 09:44

                               
登录/注册后可看大图

水汽通量散度中的q一般取克/千克，而这里q的单位为：克/克，所以乘以1000，我是这么认为的，不懂对不对

q*ugrdprs/g这个求出来不就是水汽通量么，水汽通量的单位是g/（s·hpa·cm），单位已经是克了么，求出的水汽通量散度应该不用乘以1000了吧
define q=qs*RHprs/100；这句中RHprs是相对湿度，在FNL的CTL文件中，已经有%号了，应该不用在除以100了。
我是这么认为的，我用我手头的资料试了下，按我的做法量级是10-6，不知道对不对。