麻烦高手，计算湿位涡gs 文件不不知哪儿错了，画出来图是空白，谢谢

siqin

*水汽通量：(g*s)/kg=g/(s*hpa*cm)（常用）量级1---10 两种都对g/(cm*hpa*s)  正确河套地区经纬度'set lon 98 118''set lat 35 48'，东部地区'set lon 110 130''set lat 38 52'
'reinit'
'open D:\by\fnl_20140707.ctl'
*'set mpdset hires cnworld'
*'set mpdset neimenggu neimengguq'
*'set map 1 1 1'
'set grads off'  
'set xlopts 1 4 0.16'
'set ylopts 1 4 0.16'
'set clopts -1 -1 0.13'
'set gxout shaded'
'set cterp on'
'set csmooth on'

i=1
while (i<=4)
'set t ' i
'set lon 105 130'
'set lat 35 52'
'set lev 1000 100'
'set zlog on'
'define t=TMPprs'
'define rh=RHprs'
'define prs=lev'
'define es=(6.112*exp((17.67*(t-273.16))/(t-29.65)))'
'define q=rh*(0.62197*es/(prs-es))/100.'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(t)-log(e)-4.805)'
'undefine e'
'define theta=t*pow((1000./prs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'
*'d theta'
'define u=UGRDprs'
'define v=VGRDprs'
'define vor= 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=v(z-1)-v(z+1)'
'define dx=cdiff(lon,x)*3.14159/180.0'
'define dy=cdiff(lat,y)*3.14159/180.0'
'define dtx=cdiff(eqt,x)'
'define dty=cdiff(eqt,y)'
'define pv1=-g*(vor+f)*deqt/dp'
'define pv2=g*((dv/dp)*(dtx/dx)-(du/dp)*(dty/dy))'
'define pv=pv1+pv2'
*'d hcurl(UGRDprs,VGRDprs)'
*'d uGRDprs;vGRDprs'
*'q time'
*    res = subwrd(result,3)
*    year= substr(res,1,2)
     'q time'
     res = subwrd(result,3)
     hour = substr(res,1,2)
     say hour
     'q time'
     res = subwrd(result,3)
     day=substr(res,4,2)
     say day
     'q time'
     res = subwrd(result,3)
     month=substr(res,6,3)
     say month
'set lon 121'
'set lat 35 52'
'set grid off'
'set grads off'
'set gxout contour'
'set csmooth on'
'set cthick 7'
*'set clopts -1 -1 0.10'
'd pv1*10e5'
*'d pv2*10e6'
*'d pv*10e5'
*'print'
*'disable print'
*'q time'
*res=subwrd(result,3)
'draw title 2014.'month'.'day'.'hour'  vorp'
'printim  d:\pv1.201447.'month'.'day'.'hour'.png white'
'c'
i=i+1
endwhile
;

kzpyezi

多逛帖子，才有收获呀

xiaohundanpipi

只贴这一部分程序，没有任何说明和错误信息，无法判断哪里的问题

就要唯一

要详细点，把报错的地方发一下