关于水汽通量散度计算的gs问题

fanyuxiang

最近用grads绘水汽通量散度，结果与文献的图有出入，水汽通量散度好像偏小。请各位大咖帮忙看一下问题出在哪里。下面是gs脚本：'reinit'
'open d:/grads/fnl_20110717_12_00.ctl'
'set mpdset cnworld'
'set lon 104 136'
'set lat 10 40'
'set lev 925'
'set gxout shaded'
'set cmax 0'
'define p=lev'
'define e=(6.112*exp((17.67*(TMPprs-273.16))/(TMPprs-29.65)))'
'define qs=(0.622*e/(p-0.378*e))'
'define q=qs*rhprs/100'
'define qdiv=hdivg(q*ugrdprs/9.8,q*vgrdprs/9.8)'
'set csmooth on'
'd qdiv*1e7'
'cbarn'
'printim d:/grads/print/qdiv.png white x1000 y800'
;

出图效果

yuchen

e的计算公式换个试试。

fanyuxiang

yuchen 发表于 2017-12-31 13:08
e的计算公式换个试试。

换了公式，效果不佳

fanyuxiang

yuchen 发表于 2017-12-31 13:08
e的计算公式换个试试。

换了公式，效果不佳

yuchen

文献结果贴出来看看啊。

fanyuxiang

yuchen 发表于 2018-1-5 14:42
文献结果贴出来看看啊。

D:\Documents\Pictures\无标题.png

fanyuxiang

yuchen 发表于 2018-1-5 14:42
文献结果贴出来看看啊。

D:\Documents\Pictures\无标题.png

yuchen

还是单位换算的问题，为了使水汽通量散度的单位为g/(cm^2*hPa*s)，qdiv要乘以10。

gaoyuan215

各位大神辛苦了。

四叶草

'qs=622*es/(lev-0.378*es)'本人平时计算比湿的公式是这个，与楼主的存在差别，主要来自单位换算的差别。