ave语句定义变量问题

酒国

程序是实现沿109.5E，时间25时次到43时次，水汽通量散度经向-高度分布的平均值。下面是程序：
'reinit'
'open e:/wrf/lingd02/lingd02.ctl'
'set grid off'
'set grads off'
'set zlog on'
'set parea 1 9.5 1 8'
'set lon 109.5'
'set lat 18.2 20.2'
'set lev 900 300'
'set t 25 43'
'set gxout contour'
'set xlopts 1 6 0.15'
'set ylopts 1 6 0.15'
'set clopts 1 6 0.10'
'set cthick 6'
'set ccolor 1'
'set xlint 0.2'
'set csmooth on'
'define pre=lev'
'define tc=(tk-273.16)'
'define a1=(14.55+0.114*tc)*(1-0.01*rh)'
'define aa=(2.5+0.007*tc)*(1-0.01*rh)'
'define a2=pow(aa,3)'
'define a3=(15.9+0.37*tc)*pow((1-0.01*rh),14)'
'define td=tc-a1-a2-a3'
'define vapr=6.112*exp((17.67*td)/(td+243.5))'
'define e=vapr*1.001+(pre-100)/900*0.0034'
'define mixr=0.62137*(e/(pre-e))*1000'
'define qx=u*mixr'
'define qy=v*mixr'
'define hd=hdivg(qx,qy)*1e4/9.8'
'define sdave=ave(hd,t=25,t=43)'
'set lon 109.5'
'set lat 18.2 20.2'
'set lev 900 300'
'd sdave'
运行后的错误提示为：

男紫汗

是不是应该先全部算出来再画经向-高度图

huyuhengwan

不能这样直接拿定义的变量用在ave里面吧

LYM

你这个问题解决了吗