- 积分
- 122
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-6-14
- 最后登录
- 1970-1-1
|
发表于 2015-9-12 16:49:39
|
显示全部楼层
'reinit'
'open D:\2014\fnl.ctl'
'set grid off'
'set grads off'
'set t 9 '
'set lat 10 35'
'set lon 100 140'
'set lev 1000 100'
*'set display color white'
*'c'
'set zlog on'
'set ylint 100'
'set xlint 2'
'set clopts -1 -1 0.14'
*'set xlopts 1 4 0.14'
*'set ylopts 1 4 0.14'
'set csmooth on'
'define div=hdivg(ugrdprs,vgrdprs)*100000'
'define vor=hcurl(UGRDprs,VGRDprs)*100000'
'define tmp=TMPprs'
'define rh=RHprs'
'define prs=lev'
'define es=(6.112*exp(17.67*(tmp-273.15)/(tmp-29.65)))'
'define q=rh*(0.62197*es/(prs-es))/100.'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(tmp)-log(e)-4.805)'
'define theta=tmp*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 lat 15'
'set lon 100 140'
'set t 9'
'set grads off'
*'set gxout shaded'
*'set cmax 0'
*'d vor'
'set gxout contour'
'set cint 8'
'set clopts -1 -1 0.12'
'd theta'
*'cbarn'
* 'printim d:\2014\theta_1620.png white'
; |
|