grads 画风场的地面到300hpa的垂直积分出错

nunu18

小白，弱弱的求帮助。用fnl资料画水汽通量积分和风场
本人用的部分程序如下，水汽通量积分是对的，风场我大概和水汽通量矢量比了一下，也没有太大的出入，但是量级明显的错了，我是不是应该再除以重力加速度？tu=vint ( PRESsfc/100,UGRDprs/9.8,300)'，还有输出的时候为什么要除以1000？求大神能指点

'set z 1 17'
'tc = (TMPprs-273.16)'
'td = tc-((14.55+0.114*tc)*(1-0.01*RHprs)+pow((2.5+0.007*tc)*(1-0.01*RHprs),3)+(15.9+0.117*tc)*pow((1-0.01*RHprs),14))'
'vapr = 6.112*exp((17.67*td)/(td+243.5))'
'e = vapr*1.001+(lev-100)/900*0.0034'
'q = 622*e/(lev-0.378*e)'
'define qx=UGRDprs*q/9.8'
'define qy=VGRDprs*q/9.8'
'set z 1'
'define tqx=vint(PRESsfc/100,qx,300)'
'define tqy=vint(PRESsfc/100,qy,300)'
'define   tu=vint ( PRESsfc/100,UGRDprs,300)'
'define   tv=vint ( PRESsfc/100,VGRDprs,300)'

'set gxout shaded'
'd mag(tqx/1000,tqy/1000)'
'cbarn'
'set gxout vector'
'set ccolor 1'
'set arrscl 0.5  10 '
'd tu/1000;tv/1000'

'printim e:\2016.08.rain\grads\sqjf2\300hpa2_'month'.'day'.'hour'.png  white'
'c'
i=i+1
endwhile
;

男紫汗

按照公式推一遍单位就好了