- 积分
- 92
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-3-1
- 最后登录
- 1970-1-1
|
发表于 2012-3-23 22:03:55
|
显示全部楼层
'reinit'
'open d:/data/data09/fnl_20090801_06_00_c.ctl'
'enable print d:/data/vel_cita_0106.gmf'
'set grads off'
'set lat 15 55'
'set lon 98 118'
'set lev 1000 100'
'set clab forced'
'define t0=TMPprs'
'define es=(6.112*exp((17.67*(t0-273.15))/(t0-29.65)))'
'define qs=0.62197*es/(hgtprs-0.378*es)'
'define q=qs*RHprs/100'
'define e=hgtprs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(t0)-log(e)-4.805)'
'define theta=t0*pow((1000./hgtprs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.28*q))'
'set csmooth on'
'set lat 31'
'set gxout contour'
'd smth9(eqt)'
'run cbarn.gs'
'print'
'disable print' |
|