求助！计算地形强迫抬升造成的垂直速度

pyt007

'reinit'
'open E:\grads\cx\20180608\h\2.ctl'
'set grads off'
'set grid off'
'set vpage 0 11 0 8.5'
'set parea 0.6 10.2 0.5 7.8'
'set map 1 1 10'
'set gxout line'

*------calulation data u----------------
tt=1
while(tt<=3)
'set xlopts 1 6 0.15'
'set ylopts 1 6 0.15'
'set clopts 1 10 0.16'
'set ylevs -3.0 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5'
*'set vrange -0.06 0.04'
'set cthick 8'
'set grads off'
'set gxout shaded'
'set csmooth on'
'set t 'tt
'set x 1 361'
'set y 1 181'
'define hg=HGTsfc'
'define u=UGRDsig995'
'define v=VGRDsig995'
'define r=6371000'
'define p=3.1416/180.0'
*------h---------
'define A=lat*p'
'define g=9.7803*(1+0.0053024*sin(A)*sin(A)-0.000005*sin(2*A)*sin(2*A))'
'define h=hg*9.8/g'
*--------dh----------
'define dhx=cdiff(h,x)'
'define dhy=cdiff(h,y)'
*----dx dy---du zhuanwei hudu------
'define dx=cdiff(lon,x)*p'
'define dy=cdiff(lat,y)*p'
*-----uh=u*dh/dx--------------
'define uh=u*dhx/(cos(lat*p)*dx*r)'  
*-----vh=v*dh/dy--------------
'define vh=v*dhy/(dy*r)'
'define w=uh+vh'
'set lon 98 108'
'set lat 20 30'
'set lat 25'
*'set lon 102.6'
*'d uh*10+1'
'd w*100'
'q time'
x=subwrd(result,3)
'draw title w\'x''
*'print'
tt=tt+1
endwhile
'printim E:\grads\pc\20180608\h\2015\4ww_102_'x'.png white x800 y600'
结果与文献中相差100倍，不知是哪里出错，求助各位帮忙，万分感谢

付亚男

有一个我不太明白的地方，最后w乘100，是为了调整单位吗？

pyt007

付亚男 发表于 2020-5-6 10:54
有一个我不太明白的地方，最后w乘100，是为了调整单位吗？

是的，算出来的值太小

19921219

好像计算有问题，显示outside file ，你用的公式是什么呢？

pyt007

付亚男 发表于 2020-5-6 10:54
有一个我不太明白的地方，最后w乘100，是为了调整单位吗？

是的，就是为了调整单位

zhoudog

计算中使用define hg=HGTsfc是地面位势高度场，但是公式中应该是使用地形高度计算吧，亲亲，这是错了吗？

canghaidie

define hg=HGTsfc用位势高度对不对呢？