给大家分享做剖面

蒹葭

function b( args  )
dat = subwrd( args,1)
'open f:\20160811ncep\81112.ctl'
'enable print f:\20160811ncep\thetse-t.gmf'
'set lon 90 140 1'
'set lat 10 60 1'
'set lev 1000 100'
'set t 1 19'
'define prs=lev'
'set ccolor 1'
'set grid off'
'define  l1=TMPprs'
if(l1>263.0);
'define a1=0.622*6.11*EXP(17.26*(TMPprs-273.16)/(TMPprs-35.86))'
'define b1=prs-0.278*EXP(17.26*(TMPprs-273.16)/(TMPprs-35.86))'
'define  qs1=a1/b1'
endif
if(l1<=263.0);
'define a1=0.622*6.11*EXP(21.87*(TMPprs-273.16)/(TMPprs-7.66))'
'define b1=lev-0.278*EXP(21.87*(TMPprs-273.16)/(TMPprs-7.66))'
'define  qs1=a1/b1'
endif
'define  q1=qs1*RHprs'
'define e1=lev*q1/100./(0.62197+q1/100.0)'
'define tk1=55.0+2840.0/(3.5*log(TMPprs)-log(e1)-4.805)'
'define pot1=TMPprs*pow((1000/lev),(0.2854*(1.0-0.28*q1/100.0)))'
'define ept1=pot1*exp(((3376./tk1)-2.54)*q1/100.0*(1.0+0.81*q1/100.0))'
'set lon 110'
'set lat 38'
'set t 1 19'
'set lev 850 300'
'set xlopts 1 4 0.14'
'set ylopts 1 4 0.14'
*'set mpdset cnworld'
'set grads off'
*'set yaxis 20  50 1'
*'set xaxis 105 120 1'
'set cint 1'
'd ept1-273.16'
'print'
'c'
'disable print'
*'reinit'

numa

谢谢分享！上个图看看效果呗！

rainbow007

numa 发表于 2017-9-25 14:51
谢谢分享！上个图看看效果呗！

，谢谢，我试试先