- 积分
- 941
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-4-27
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2016-12-25 13:38:49
|
显示全部楼层
'reinit'
'open e:\datafiles\fnl_20161030_18_00.ctl'
'set lat 30 50'
'set lon 100 125'
'set lev 1000 100'
'set display color white'
'c'
'set gxout shaded'
'set mpdset cnworld'
'set csmooth on'
'set grads off'
'define w=vvelprs'
'set lon 123'
'set lat 35 45'
'set lev 1000 400'
'set clevs -2 -1.5 -1 -0.5 0'
'set ccols 14 9 4 5 3 0'
'd w'
'cbarn'
'define prs=lev'
'define es=(6.112*exp(17.67*(tmpprs-273.15)/(tmpprs-29.65)))'
'define q=RHprs*(0.62197*es/(prs-es))/100.'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(tmpprs)-log(e)-4.805)'
'define theta=tmpprs*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 cthick 4'
'set clskip 1'
'set gxout contour'
'set cthick 4'
'set cint 1'
'set gxout contour'
'd eqt'
'draw title 30-18-eqt-vvelprs'
'printim e:\datafiles\30-18-eqt-vvelprs.jpg'
'print'
’reinit' |
|