500hpa水汽通量计算结果为负值

cafflybird

我计算的水汽通量结果为何是负值
贴上其中一年的图


所用的数据
air.mon.mean.nc
rhum.mon.mean.nc
uwnd.mon.mean.nc
vwnd.mon.mean.nc

计算过程
'reinit'
'sdfopen e:\data0220\air.mon.mean.nc'
'sdfopen e:\data0220\rhum.mon.mean.nc'
'sdfopen e:\data0220\uwnd.mon.mean.nc'
'sdfopen e:\data0220\vwnd.mon.mean.nc'
'enable print e:\data0220\test.gmf '
nt=145
while(nt<=768)
'set t 'nt
'set lev 500'
'define es=6.112*exp(17.269*(air-273.15)/(air-29.65))'
'define qs=0.622*es/(lev-0.378*es)'
'define q=rhum.2*qs/100'
'define v=mag(uwnd.3,vwnd.4)'
'define uvq=v*q/9.8'
'set grads off'
'set grid off'
'set mpdset cnworld'
'set lon 40 170'
'set lat 10 70'
'set gxout shaded'
'd uvq'
'cbarn'
'print'
'c'
nt=nt+1
endwhile
'disable print'

Spider-Woman

'define uvq=v*qs/9.8'
应该是'v*q/9.8' 吧 定义错了

river

水汽通量没有那么算的吧，要先分开计算，r然后再合成或者画图的时候再合成。
define qu=(uwnd.3*qs/9.8)
define qv=(vwnd.4*qs/9.8)
define uqv=mag(qu,qv)

cafflybird

Spider-Woman 发表于 2013-3-8 11:05

                               
登录/注册后可看大图

'define uvq=v*qs/9.8'
应该是'v*q/9.8' 吧 定义错了

恩，表达错了，我是用的v*q/9.8,结果还是负值，，，我这就改下帖子

cafflybird

river 发表于 2013-3-8 12:13

                               
登录/注册后可看大图

水汽通量没有那么算的吧，要先分开计算，r然后再合成或者画图的时候再合成。
define qu=(uwnd.3*qs/9.8)
...

我马上试试

xiayu1990426

river 发表于 2013-3-8 12:13

                               
登录/注册后可看大图

水汽通量没有那么算的吧，要先分开计算，r然后再合成或者画图的时候再合成。
define qu=(uwnd.3*qs/9.8)
...

水汽通量不可以用vint来算么

river

xiayu1990426 发表于 2013-8-5 21:07

                               
登录/注册后可看大图

水汽通量不可以用vint来算么

那是计算整层水汽通量采用vint积分

Paranoia

整层水汽通量用VINT吧