fnl资料画整层水汽通量

Crystal_PP

如题，附件为fnl资料画整层水汽通量及散度的gs，但个人觉得这个gs有问题，我画出来的图和文献的图对不上。请问高人们问题出在哪了？

'reinit'
'open D:\lianxi\4\fnl_20080828_00_00.ctl'

'set gxout fwrite'
'set fwrite D:\lianxi\4\q.dat'  
'set lon 0 359'
'set lat -90 90'

iz=1
while(iz<=21)
'set z 'iz''
'define prs=lev'
'define es=6.11*exp(17.27*(tmpprs-273.16)/(tmpprs-35.86))'
'define qs=0.622*es/(prs-0.378*es)'
'd qs*rhprs/100'
iz=iz+1
endwhile
'disable fwrite'
'reinit'


'reinit'
'open D:\lianxi\4\fnl_20080828_00_00.ctl'
'set gxout fwrite'
'set fwrite D:\lianxi\4\u.dat'
'set lon 0 359'
'set lat -90 90'

iz=1
while(iz<=21)
'set z 'iz''
'd ugrdprs'

iz=iz+1
endwhile
'disable fwrite'
'reinit'


'reinit'
'open D:\lianxi\4\fnl_20080828_00_00.ctl'
'set gxout fwrite'
'set fwrite D:\lianxi\4\v.dat'
'set lon 0 359'
'set lat -90 90'

iz=1
while(iz<=21)
'set z 'iz''
'd vgrdprs'
iz=iz+1
endwhile
'disable fwrite'
'reinit'


'reinit'
'open D:\lianxi\4\fnl_20080828_00_00.ctl'
'set gxout fwrite'
'set fwrite D:\lianxi\4\ps.dat'
'set lon 0 359'
'set lat -90 90'
'd PRESsfc'
'disable fwrite'
'reinit'


'reinit'
'enable print D:\lianxi\4\2008082808quv.gmf'
'open D:\lianxi\4\q.ctl'
'open D:\lianxi\4\u.ctl'
'open D:\lianxi\4\v.ctl'
'open D:\lianxi\4\ps.ctl'

'set lon 89 121'
'set lat 14 41'
'set map 1 1 6'
'set mpdset cnworld'

'set lev 1000 300'
'define qx=u.2*q.1'
'define qy=v.3*q.1'

'set z 1'
'define a=vint(Ps.4/100,u.2*q.1/9.8,300)'
'define b=vint(Ps.4/100,v.3*q.1/9.8,300)'
'define div=hdivg(a,b)*100000'


'set gxout vector'
'set arrscl 0.3 60'
'd a;b'

'print'
'disable print'

river

和文献对不上是怎么对不上？拿图来做对比一下······

thukd35

是不是vint函数的vint(expr,a,b)里的参数位置写乱了啊，expr是被积函数，a是积分下限，b是积分上限

may6117

请教：'open D:\lianxi\4\q.ctl'
'open D:\lianxi\4\u.ctl'
'open D:\lianxi\4\v.ctl'
'open D:\lianxi\4\ps.ctl'
目前没有这四个ctl，是要自己写吗？

may6117

请教：'open D:\lianxi\4\q.ctl'
'open D:\lianxi\4\u.ctl'
'open D:\lianxi\4\v.ctl'
'open D:\lianxi\4\ps.ctl'
目前没有这四个ctl，是要自己写吗？