用isen脚本画等熵位涡哪错了？

.雾霭流岚

我用isen脚本画等熵位涡，总是画不出图，不知是哪错了，请各位高手指点！
画的是两个等熵面350K和310K
grads错误：cannot contourgrid - all undefined values

ctl是fnl批处理得到的

gs：(function后的是脚本，直接粘贴的)
'reinit'
'open D:\date\bolaven\fnl_20120825_12_00_c.ctl'
'set lon 100 160'
'set lat 10 70'
'set t 16'
'set lev 900 100'
'set mpdset mres'
'set grid off'
'set display color white'
'c'
'set grads off'

"define tt=TMPprs"
"define u=ugrdprs"
"define v=vgrdprs"
"define pp=lev"
"define coriol=2*7.29e-5*sin(lat*3.1415/180)"
"define dudy=cdiff(u,y)/(111177*cdiff(lat,y))"
"define dvdx=cdiff(v,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))"
"define dt=tt(z-1)*pow(1000/pp(z-1),0.286)-tt(z+1)*pow(1000/pp(z+1),0.286)"
"define dp=100*(pp(z-1)-pp(z+1))"
"define dtdp=dt/dp"

'set gxout shaded'
'set csmooth on'
"define part1="isen(dvdx,tt,PP,350)
"define part2="isen(dudy,tt,PP,350)
"define part3="isen(dtdp,tt,PP,350)
"define pv350=-9.8*(coriol+part1-part2)*part3"
"set z 1"
'set lon 110 150'
'set lat 20 60'
"d pv350*1000000"

'set gxout contour'
'set csmooth on'
"define part1="isen(dvdx,tt,PP,310)
"define part2="isen(dudy,tt,PP,310)
"define part3="isen(dtdp,tt,PP,310)
"define pv310=-9.8*(coriol+part1-part2)*part3"
"set z 1"
'set lon 110 150'
'set lat 20 60'
'set cstyle 3'
'set cthick 4'
"d pv310*1000000"
'q w2xy 129.4 46.3'
x1=subwrd(result,3)
y1=subwrd(result,6)
'draw wxsym 41 'x1' 'y1' 0.3<1<2>>'
'enable print d:\date\bpvu16.gmf'
'print'
'disable'
;

.雾霭流岚

出图区是一片空白

.雾霭流岚

把isen脚本那部分去掉了，重点问题在cannot contourgrid - all undefined values
ctl是正确的，已经做出了其他物理量的图

river

你运行了你的gs有什么错误提示，发上来让别人看看那啊。这么长的程序没有几个人会从头看到尾检查错误的

1417920146

试试把ipv350、ipv310的各四行计算移动到作图的语句上面，作图时对经纬度重新设置了会不不会有影响呢

Joucly

想问问楼主画出来没有

Joucly

我和你一样的错误

Joucly

我和你一样的错误

clare

lz尝试一下把第10行的‘c'去掉~

zhukefu

问题怎么解决了呀？