求助

qxtjls

'reinit'
'open c:\grads20\gs1\fnl_20120828_12.ctl'
'enable c:\grads20\gs1\weiwo5.jpg'
'set z 2 20'
'set zlog on'
'set lon 70 150'
'set lat 20 70'
'define tc=tmpprs-273.15'
'define rh=rhprs'
'define prs=lev'  
'define es=(6.112*exp((17.67*tc)/(tc+243.5)))'
'define qs=0.62197*es/(prs-0.378*es)'
'define q=rh*qs/100'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(tc+273.16)-log(e)-4.805)'
'define theta=(tc+273.16)*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'

'define u=UGRDprs'
'define v=VGRDprs'
'define vo=hcurl(u,v)'
'define f=2*7.292*sin(lat*3.14159/180.0)*0.00001'
'define g=9.8'
'define dp=100*(lev(z-1)-lev(z+1))'
'define deqt=eqt(z-1)-eqt(z+1)'
'define du=u(z-1)-u(z+1)'
'define dv=u(z-1)-u(z+1)'
'define dx=2.0*6370949.0*cos(lat*3.14159/180.0)*3.14159/180.0'
'define dy=2.0*6370949.0*3.14159/180.0'
'define dtx=cdiff(eqt,x)'
'define dty=cdiff(eqt,y)'
'define pv1=-g*(vo+f)*deqt/dp'  
'define pv2=g*((dv/dp)*(dtx/dx)-(du/dp)*(dty/dy))'
'define pv=pv1+pv2'

'set lon 75 125'
'set lat 20 60'
'set grid off'
'set grads off'
'set gxout contour'
'set csmooth on'
'set cthick 7'
'set clopts -1 -1 0.10'
'd pv*10e5'
*'d pv2*10e6'
*'d pv*10e5'
'print'
'printim c:\grads20\gs1\weiwo5.jpg white'
'disable print'

ctl如下：
dset C:\GrADS20\gs\fnl_20120828_12_00_c
index C:\GrADS20\gs\fnl_20120828_12.idx
undef 9.999E+20
title fnl_20120828_12_00_c
dtype grib 3
options yrev
ydef 181 linear -90.000000 1
xdef 360 linear 0.000000 1.000000
tdef 1 linear 12Z28aug2012 1mo
zdef 26 levels
1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10
vars 23
ABSVprs 26 41,100,0 ** Absolute vorticity
ACPCPsfc  0 63,1,0 ** surface Convective precipitation [kg/m^2]
APCPsfc  0 61,1,0 ** surface Total precipitation [kg/m^2]
CAPEsfc  0 157,1,0 ** surface Convective Avail. Pot. Energy [J/kg]
CINsfc  0 156,1,0 ** surface Convective inhibition [J/kg]
CLWMRprs 21 153,100,0 ** Cloud water [kg/kg]
HGTprs 26 7,100,0 ** Geopotential height [gpm]
PRESsfc  0 1,1,0 ** surface Pressure [Pa]
PRESmsl  0 1,102,0  ** Pressure [Pa]
PRES 2 1,117,0 ** Pressure [Pa]
RHprs 21 52,100,0 ** Relative humidity [%]
RH2m  0 52,105,2 ** 2 m Relative humidity [%]
SPFH2m  0 51,105,2 ** 2 m Specific humidity [kg/kg]
TMAX2m  0 15,105,2 ** 2 m Max. temp. [K]
TMIN2m  0 16,105,2 ** 2 m Min. temp. [K]
TMPsfc  0 11,1,0 ** surface Temp. [K]
TMPprs 26 11,100,0 ** Temp. [K]
TMP2m  0 11,105,2 ** 2 m Temp. [K]
UGRDprs 26 33,100,0 ** u wind [m/s]
UGRD10m  0 33,105,10 ** 10 m u wind [m/s]
VGRDprs 26 34,100,0 ** v wind [m/s]
VGRD10m  0 34,105,10 ** 10 m v wind [m/s]
VVELprs 21 39,100,0 ** Pressure vertical velocity [Pa/s]
ENDVARS

就是不出图，请教各路大侠帮帮忙，看那里需要改，不胜感激！！！

river

在gs文件最后加上一个  ;   试试

qxtjls

多谢你的帮助，但还是不出图，错误信息如下
a->swwgs1.gs
No hardcopy metafile open
All grads attributes have been reinitialized
cannot contour grid-all undefined values
Eps file written to grads.esp
No hardcopy metafile open