求助大神计算湿位涡问题！！！！在线等！！！

当年明月

第一张图是我在运行下面的脚本程序之后出现的，第二张图是在我试图打开nc自带描述程序时出现的，求大神解答啊！！！！！！！！！！！！！！！！！在线等！！！！
'reinit'
'sdfopen f:\nc\nc1\air.2015.nc'
'sdfopen f:\nc\nc1\rhum.2015.nc'
'sdfopen f:\nc\nc1\uwnd.2015.nc'
'sdfopen f:\nc\nc1\vwnd.2015.nc'


'set mpdset cnworld'
'draw map'
'set lat 0 50'
'set lon 100 120'


'set ylevs 995 900 800 700 600 500 400 300'
'set lev 1000 100'
'set zlog on'


'set t 880'

*****计算假相当位温*****
'set grads off'
'define prs=lev'
'define es=6.112*exp(17.67*(air-273.15)/(air-29.65))'
'define qs=0.62197*es/(prs-0.378*es)'
'define qv=rhum.1*qs/100'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(air)-log(e)-4.775)'

'define theta=air*pow((1000./prs),(0.2854*(1.0-0.28*qv)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'

*****计算湿位涡3*****
'define vo=hcurl(uwnd.3,vwnd.4)'
'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=uwnd.3(z-1)-uwnd.3(z+1)'
'define dv=vwnd.4(z-1)-vwnd.4(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 107'
'set lat 10 40'
'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 'res'    'pv

'printim f:\0606.jpg white'

当年明月

求助大神，，，在线等啊！自顶！！

river

建议加强grads的基础知识······

爽歪歪

。。。。rhum.2吧

当年明月

爽歪歪 发表于 2015-11-22 20:36
。。。。rhum.2吧

是呀，我改成rhum.2啦

547680911

{:eb513:}{:eb513:}{:eb513:}{:eb513:}{:eb513:}

547680911

{:eb502:}{:eb502:}{:eb502:}{:eb502:}{:eb502:}

547680911

{:5_214:}{:5_214:}{:5_214:}{:5_214:}{:5_214:}

547680911

{:5_214:}{:5_214:}{:5_214:}{:5_214:}{:5_214:}

547680911

{:eb513:}{:eb513:}{:eb513:}{:eb513:}{:eb513:}{:eb513:}