Q矢量及锋生函数量级问题

木果果

自己用grads算了q矢量和锋生函数，可是量级总是对不上啊。一般量级为10的-16次方，可是我必须乘以10的19次方才行，可否有人能帮我看看问题出在哪儿呢？万分感谢~~~~~~~~~下面是我的计算q矢量的gs：

'reinit'
'open E:\973\fnldata\2013070506.ctl'
'set lon 100 130'
'set lat 20 45'
*'set lev 1000 100'
'set lev 850'
'set grads off'
'set grid off'
'set mpdset cnworld '


****计算假相当位温***********
'define prs=lev'
'define es=(6.1078*exp(17.269*(TMPprs-273.16)/(TMPprs-35.86)))'
'define qs=RHprs*(0.62197*es/(prs-es))/100.'
'define e=prs*qs/(0.62197+qs)+1e-10'
********上式加上1e-10'，是因为下式需要对e去log，而log（x）中x必须大于0，为防止算出的e不满足大于0而加********
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.805)'
'define theta=TMPprs*pow((1000/prs),(0.2854*(1.0-0.28*qs)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*qs*(1.0+0.81*qs))'

******计算地转Q矢量（梅雨锋用假相当位温tidu代替温度梯度）**********
'define f=2*7.292e-5*sin(lat*3.14159/180)'
'define ug= -9.8/f*cdiff(HGTprs,y)/(6.37e6*cdiff(lat,y))'
'define vg= 9.8/f*cdiff(HGTprs,x)/(6.37e6*cdiff(lon,x)*cos(lat*3.14159/180))'
'define dugdx=cdiff(ug,x)/(6.37e6*cdiff(lon,x)*cos(lat*3.14159/180))'
'define dvgdx=cdiff(vg,x)/(6.37e6*cdiff(lon,x)*cos(lat*3.14159/180))'
'define dtdx=cdiff(eqt,x)/(6.37e6*cdiff(lon,x)*cos(lat*3.14159/180))'
'define dugdy= cdiff(ug,y)/(6.37e6*cdiff(lat,y))'
'define dvgdy= cdiff(vg,y)/(6.37e6*cdiff(lat,y))'
'define dtdy= cdiff(eqt,y)/(6.37e6*cdiff(lat,y))'

'define qx=(-287.05/prs)*(dugdx*dtdx+dvgdx*dtdy)'
'define qy=(-287.05/prs)*(dugdy*dtdx+dvgdy*dtdy)'
'define div=hdivg(qx,qy)*1e19'

'set mpdset cnworld'
'set map 15 4 6'
'set lon 110 125'
'set lat 25 40'
'set gxout contour'
'set ccolor 2'
'set cthick 10'
'set clab forced'
'set clskip 1'
'set cmax 0'
'set black 0 0'
'set cint 0.1'
'd div'

*'set gxout barb'
*'set ccolor 1'
*'d 2.5*u;2.5*v'

'set mpdset cnriver'
'set map  3  1 8'
'draw map'
'q time'
xx=subwrd(result,3)
'draw title 'xx' 850hPa_Q'
*'printim E:\973\qshiliang\'xx'_850hpa.png white x1200 y800'
'printim E:\973\qshiliang\'xx'_850hpa.png white '

;