求一个过程的假相当位温场，大神们帮忙看看哪里出错了

右手年华

'reinit'
'sdfopen E:\rhum.2015.nc'
'sdfopen E:\uwnd.2015.nc'
'sdfopen E:\vwnd.2015.nc'
'sdfopen E:\air.2015.nc'

*绘图区域
'set vpage 0 11 0 8.5'
'set parea 1.2 10.5 1.5 8'
'set grads off'
'set grid off'
'set mpdset cnworld'
'set xlopts 1 4 0.25'
'set ylopts 1 4 0.25'
'set ylevs 10N 20N 30N 40N 50N 60N 70N 80N'
'set xlevs 40E 50E 60E 70E 80E 90E 100E 110E 120E 130E 140E 150E 160E'
'set clopts 1 1 0.2'
'set ccolor 1'

'set lev 850'
'set lon 90 140'
'set lat 0 60'
'set t 155'
'define l1=ave(air.4,t=155,t=170)'
if(l1>263.0);
'define a1=0.622*6.11*EXP(17.26*(l1-273.16)/(l1-35.86))'
'define b1=850-0.278*EXP(17.26*(l1-273.16)/(l1-35.86))'
'define  qs1=a1/b1'
endif
if(l1<=263.0);
'define a1=0.622*6.11*EXP(21.87*(l1-273.16)/(l1-7.66))'
'define b1=850-0.278*EXP(21.87*(l1-273.16)/(l1-7.66))'
'define  qs1=a1/b1'
endif
'define  q1=qs1*ave(rhum.1,t=155,t=170)'
'define e1=850*q1/100./(0.62197+q1/100.0)'
'define tk1=55.0+2840.0/(3.5*log(l1)-log(e1)-4.805)'
'define pot1=l1*pow((1000/850),(0.2854*(1.0-0.28*q1/100.0)))'
'define ept1=pot1*exp(((3376./tk1)-2.54)*q1/100.0*(1.0+0.81*q1/100.0))'

'set strsiz 0.25'
'set string 1 tl 1 90'
*'draw string 0.1 4.5 P/hPa'

'set gxout contour'
'set clab forced'
'set cterp on'
'set csmooth on'
'set cint 2'
'set ccolor 1'
*'set lev 1000 100'
*'set lon 118'
*'set lat 10 50'
*'set ylevs 1000 925 850 700 600 500 400 300 200 100'

*'d ave(ave(ept1,lon=116,lon=120),t=155,t=170)'
'd ept1'
'enable print E:\850se\se.gmf 1024 768 white'
'print'
'disable print'
'c'

勿忘草

根据提示，首先你的坐标设置是无效的，其次gs文件中新变量的定义方法好像有问题，是一些细节的问题，课本上应该有这方面的例子

右手年华

勿忘草 发表于 2016-1-4 21:23
根据提示，首先你的坐标设置是无效的，其次gs文件中新变量的定义方法好像有问题，是一些细节的问题，课本上 ...

恩，发给我朋友，他那边没改居然能出图

超级大坏蛋

嘿嘿，你要做这个？

右手年华

超级大坏蛋 发表于 2016-1-9 15:40
嘿嘿，你要做这个？

这不是伟哥嘛