东西南北边界整层水汽通量

1092995269

下面的gs是计算各边界整层水汽通量的，空心圆是东边界、实心圆是西边界、南边界是空心方框、北边界是实心方框，文献和我使用的资料都是FNL1*1的时间一样，结果不一样，帮看下是什么地方出错，谢了！！！
'reinit'
'open d:\521\2007\fnl.ctl'
'set grads off'
'set grid off'
'set t 49 68'
***比湿计算***
'define a=17.67'
'define b=29.65'
**定义饱和  水气压es(TMPprs-温度,单位:K)**
'define es=6.1078*exp(a*(TMPprs-273.15)/(TMPprs-b))'
**定义饱和比湿qs(p-气压,单位:mb=Pa)**
'define qs=0.622*es/(lev-0.378*es)'
**定义比湿q(RHprs是相对湿度)**
'q=RHprs*qs/100'
***U方向水汽通量***
'define a1=q*UGRDprs/9.8*1000'
***垂直积分***
'define aa=vint(PRESsfc/100,a1,100)'
'set vrange -40 100'
'set gxout line'
'set ccolor 1'
'set cstyle 1'
'set cmark 2'
'set lon 105'
***经向积分***
'define F1=sum(aa,lat=22,lat=28)'
'd F1/10000'
'set gxout line'
'set ccolor 1'
'set cstyle 2'
'set cmark 3'
'set lon 98'
'define F2=sum(aa,lat=22,lat=28)'
'd F2/10000'
'define a2=q*VGRDprs/9.8*1000'
'define bb=vint(PRESsfc/100,a2,100)'
'define F3=sum(bb,lon=98,lon=105)'
'set vrange -40 100'
'set gxout line'
'set ccolor 1'
'set cstyle 3'
'set cmark 4'
'set lat 20'
'd F3/10000'
'set gxout line'
'set ccolor 1'
'set cstyle 4'
'set cmark 5'
'set lat 28'
'define F4=sum(bb,lon=98,lon=105)'
'd F4/10000'
'q time'
tt=subwrd(result,5)
'printim d:\521\2007\500\500'%tt%'.png  white'
;

hshse'fe'f

请问这个图 用什么软件加载出来的

hshse'fe'f

请问这个图 用什么软件加载出来的

tulalang

楼主你的v的分量呢。。为何你水汽通量只用了u方向的。。

看云淡风清

楼主，vint函数似乎用的有问题，这个函数默认的除了9.8，还有你边界的水汽通量是按F1=sum(aa,lat=22,lat=28)计算的，有没有问题？

xuwei

楼主找到问题所在了么？