- 积分
- 36
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-12-20
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
'reinit'
'open D:\lianxi\130111_721\T721\T721.ctl'
'open D:\lianxi\130111_721\RH721\RH721.ctl'
'enable print D:\lianxi\130111_721\THETA\THETA_1.gmf'
i=69
while(i<=109)
'set grid off'
'set t 'i
'set lon 45 165'
'set lat 30 50'
'set lev 500'
'define m=Tprs'
'define rh=Rprs'
'define prs=lev'
'define es=(6.112*exp(17.67*(m-273.15)/(m-29.65)))'
'define q=rh*(0.62197*es/(prs-0.378es))/100.'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(m)-log(e)-4.805)'
'define theta=m*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 gxout contour'
'set cthick 5'
'd eqt'
'print'
'c'
i=i+1
endwhile
'disable print'
'reinit'
;
|
|