涡度方程中涡度水平平流项(A项)总是提示entire grid undefined?

bear_dancing

'set t 12'
'set z 2 20'
'set lon 108 111'
'set lat 22 24'
**************
'define f=2*7.292*sin(lat*3.14159/180.0)*1e-5'
'define vor=hcurl(UGRDprs,VGRDprs)'
'define dx=cdiff(lon,x)*3.1416/180*6370949'
'define dy=cdiff(lat,y)*3.1416/180*6370949'
'define dp=(lev(z+1)-lev(z-1))'
'define dvorx=cdiff(vor,x)'   
'define dvory=cdiff(vor,y)'
'define dfy=cdiff(f,y)'
'define A=-(UGRDprs*(dvorx/(cos(lat*3.1416/180)*dx))+VGRDprs*(dfy/dy

+dvory/dy))*1e9'（这一项怎么改都提示entire grid undefined！！！实在不知道怎么办了）
'define dvorp=hcurl(UGRDprs(z+1),VGRDprs(z+1))*10000-hcurl(UGRDprs(z-

1),VGRDprs(z-1))'
'define dvorpdp=dvorp/dp'
'define B=-10*VVELprs*dvorpdp*1e9'
'define dux=cdiff(UGRDprs,x)'
'define dvy=cdiff(VGRDprs,y)'
'define C=(f+vor)*(dux/(cos(lat*3.1416/180)*dx)+dvy/dy)*1e9'

'define dup=UGRDprs(z+1)-UGRDprs(z-1)'
'define dvp=VGRDprs(z+1)-VGRDprs(z-1)'
'define amiga=-10*VVELprs'
'define dVVELprsx=cdiff(amiga,x)'
'define dVVELprsy=cdiff(amiga,y)'
'define D=-((dVVELprsx/(cos(lat*3.1416/180)*dx))*(dvp/dp)-(dVVELprsy/dy)*

(dup/dp))*1e9'
'define E=A+B+C+D'

********************************************
'define aveA=aave(A,lon=108.8,lon=110.4,lat=22.5,lat=23.6)'
'define aveB=aave(B,lon=108.8,lon=110.4,lat=22.5,lat=23.6)'
'define aveC=aave(C,lon=108.8,lon=110.4,lat=22.5,lat=23.6)'
'define aveD=aave(D,lon=108.8,lon=110.4,lat=22.5,lat=23.6)'
'define aveE=aave(E,lon=108.8,lon=110.4,lat=22.5,lat=23.6)'
'set zlog on'
'set ylab on'
'set grads off'
'set grid off'
'set t 11'
'set z 2 20'
'set x 1'
'set y 1'
*'set vrange -12 8'

'set gxout line'
'set ccolor 1'
'set cthick 8'
'set cstyle 1'
'set cmark 1'
'd aveE'
'cbar_line -x 9.00 -y 6.77 -c 1 -m 1 -l 1 -t "E" -p'

'set gxout line'
'set ccolor 1'
'set cthick 8'
'set cstyle 2'
'set cmark 2'
'd aveA'
'cbar_line -x 9.00 -y 6.77 -c 1 -m 2 -l 2 -t "A" -p'

'set gxout line'

'set ccolor 1'
'set cthick 8'
'set cstyle 4'
'set cmark 3'
'd aveB'
'cbar_line -x 9.00 -y 6.77 -c 1 -m 3 -l 4 -t "B" -p'

'set gxout line'

'set ccolor 1'
'set cthick 8'
'set cstyle 5'
'set cmark 4'
'd aveC'
'cbar_line -x 9.00 -y 6.77 -c 1 -m 4 -l 5 -t "C" -p'

'set gxout line'

'set ccolor 1'
'set cthick 8'
'set cstyle 6'
'set cmark 5'
'd aveD'
'cbar_line -x 9.00 -y 6.77 -c 1 -m 5 -l 6 -t "D" -p'

'q time'
title=sublin(result,1)
ret=subwrd(title,3)
'draw title 'ret''
'enable print g:\df1\buget-'ret'.gmf'
'print'
'disable print'

dcx

在你求平均之后，set zlog之前加上'c' 'reset'两条语句试一下，不知道对不对

dcx

请问你的问题解决了么

bear_dancing

dcx 发表于 2016-1-15 18:31
请问你的问题解决了么

木有想不明白，你有办法么

仗剑海天

请问你的问题解决了么