假相当位温的图画不出来

shuling216

假相当位温的图画不出来假相当位温的图画不出来假相当位温的图画不出来

独奏古筝到天明

'reinit'
'open e:\nc\fnl_20160415_00_00.ctl'
'set lon 100 140'
'set lat 15 50'
'set zlog on'

i=1
while(i<=27)
'set lev 1000 100'
'set lon 120'
'set lat 15 45'
'set mpdset hires'
'set grid on'
'set grads off'

'set t 'i
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set clopts 1 4 0.08'
'define t=tmpprs'
'define rh=rhprs'
'define prs=lev'
'define es=(6.112*exp(17.67*(t-273.15)/(t-29.65)))'
'define q=rh*(0.62197*es/(prs-0.378*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))'

'd eqt'
'q time'
t=subwrd(result,3)
'draw title  120`3.`1E TKETA`bSE`n_'i'_'t''
'printim e:\nc\120Etketa_'i'_'t'.png white'
'c'
i=i+1
endwhile
'reinit'