画假相当位温的垂直剖面出错

zlh

求大神，改了几遍还是这个样，不知道哪里的原因'reinit''sdfopen D:\20081228\NC\air.2008.nc'
'sdfopen D:\20081228\NC\rhum.2008.nc'


'set grid off'
'set grads off'
'set t 1447'
'set xlint 5'
'set ylint 100'
'set ylopts 1 6 0.15'
'set xlopts 1 6 0.15'
'set lat 35'
'set lon -110 -80'
'set lev 1000 300'
'set gxout contour'




'define t=air'
'define rh=rhum'
'define prs=lev'
'define es=(6.112*exp(17.67*(t-273.15)/(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)'
'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))'
'set cint 4'


*'set gxout contour'


'set cthick 5'
'd eqt'
;

zlh

看了这个贴子，http://bbs.06climate.com/forum.p ... 0%B5%B1%CE%BB%CE%C2
问题已解决