求助！关于计算2011年7月整层水汽通量与常年的差值场

from1214

求助！我写了一个gs文件计算2011年7月整层水汽通量与常年的差值场，我觉得量级和方向都有问题（图片所示），但是我不知道是哪里出错，望高手指导，急


'reinit'
'sdfopen d:\nc\shum.mon.mean.nc'
'sdfopen d:\nc\uwnd.mon.mean.nc'
'sdfopen d:\nc\vwnd.mon.mean.nc'
'sdfopen d:\nc\pres.mon.mean.nc'
'enable print d:\gs\gmf\qtl11.gmf'
'set parea 1 9 1 8'
'set xlopts 1 4 0.2'
'set ylopts 1 4 0.2'
'set mpdset cnworld'
'set lon 60 140'
'set lat 0 40'
'set lev 1000 300'
'set grads off'
'set grid off'
'set t 763'
'define uuq=uwnd.2*shum'
'define vvq=vwnd.3*shum'
'pp=pres.4(z=1)'
'define ffuq=-vint(pp,uuq,300)/9.8'
'define ffvq=-vint(pp,vvq,300)/9.8'
'define shtll=mag(ffuq,ffvq)'

'define u=ave(uwnd.2,t=283,t=751,12)'
'define v=ave(vwnd.3,t=283,t=751,12)'
'define s=ave(shum,t=283,t=751,12)'
'define uq=u*s'
'define vq=v*s'
'set z 1'
'define p=ave(pres.4(z=1),t=283,t=751,12)'
'define fuq=-vint(p,uq,300)/9.8'
'define fvq=-vint(p,vq,300)/9.8'
'define shtl=mag(fuq,fvq)'  
'set gxout shaded'
'run cbar.gs'

'set cthick 6'
'd shtll-shtl'                              
'd ffuq;ffvq'
'cbarn 0.7 1 9.5 4.5'
'print'
c
'reinit'

曾记否

'define ffuq=-vint(pp,uuq,300)/9.8'
'define ffvq=-vint(pp,vvq,300)/9.8'
去掉负号，9.8改成1000

from1214

曾记否 发表于 2013-4-15 22:03

                               
登录/注册后可看大图

'define ffuq=-vint(pp,uuq,300)/9.8'
'define ffvq=-vint(pp,vvq,300)/9.8'
去掉负号，9.8改成1000

不好意思，请问不应该除以g么？公式中有一个负号，删掉可以吗？这样得出的单位是g/cm·s，为什么要改成1000呢？

river

整层水汽通量积分前面问为什么要加负号呢，不需要的啊

曾记否

from1214 发表于 2013-4-15 22:14

                               
登录/注册后可看大图

不好意思，请问不应该除以g么？公式中有一个负号，删掉可以吗？这样得出的单位是g/cm·s，为什么要改成10 ...

以前做过整层水汽通量，具体的公式中单位换算过程忘了，只保留了程序！就按我说的那样！

from1214

river 发表于 2013-4-15 22:19

                               
登录/注册后可看大图

整层水汽通量积分前面问为什么要加负号呢，不需要的啊

我是看文献积分公式前面有个负号

from1214

曾记否 发表于 2013-4-15 23:05

                               
登录/注册后可看大图

以前做过整层水汽通量，具体的公式中单位换算过程忘了，只保留了程序！就按我说的那样！

好的，感谢你的帮助

from1214

曾记否 发表于 2013-4-15 23:05

                               
登录/注册后可看大图

以前做过整层水汽通量，具体的公式中单位换算过程忘了，只保留了程序！就按我说的那样！

不好意思，再问一个，这样得出来的单位是？

曾记否

from1214 发表于 2013-4-16 08:38

                               
登录/注册后可看大图

不好意思，再问一个，这样得出来的单位是？

应该是kg/m/s，记不太清了

river

from1214 发表于 2013-4-16 08:35

                               
登录/注册后可看大图

我是看文献积分公式前面有个负号

你用的又不是书上的那个公式，你使用的一个函数啊，里面自然已经包括了这个公式，所以没必要自己加负号了