- 积分
- 24
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-5-18
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2012-9-28 17:35:08
|
显示全部楼层
'reinit'
'sdfopen c:\pcgrads\air.2011-1.nc'
'sdfopen c:\pcgrads\shum.2011-1.nc'
'sdfopen c:\pcgrads\rhum.2011-1.nc'
'sdfopen c:\pcgrads\slp.2011-1.nc'
'set lon 110 130'
'set lat 25 45'
'set t 873'
'd slp.4'
'set lev 500'
'set gxout shaded'
'd rhum.3'
'cbarn 1 1'
'set gxout contour'
'set ccolor 0'
'define prs=lev'
'define e1=17.67*(air.1-273.16)/(air.1-29.65)'
'define es=6.112*exp(e1)'
'define qs=0.62197*es/(prs-0.378*es)'
'define a=prs*shum.2/(0.62197+shum.2)+1e-10'
'define tk=55.0+2840.0/(3.5*log(air.1)-log(a)-4.805)'
'define pot=air.1*pow((1000/prs),(0.2854*(1.0-0.28*shum.2)))'
'define eqt=pot*exp(((3376./tk)-2.54)*shum.2*(1.0+0.81*shum.2))'
'set cint 4'
'd eqt'
; |
|