- 积分
- 5297
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-6-23
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2014-12-12 22:10:03
|
显示全部楼层
*计算比湿
'reinit'
'open e:/kt/7/20130716/071506.ctl'
'set mpdset cnworld'
'set map 1 1 8'
'set lon 72 97'
'set lat 34 50'
'set lev 700'
tt=1
while(tt<=4)
'set grads off'
'set grid off'
'set t 'tt
'define prs=lev'
'define tc=TMPprs-273.16'
'define es=(6.112*exp((17.67*tc)/(tc+243.5)))'
'define q=(0.62197*es/(prs-es))'
'define tv=TMPprs*(1+0.61*q)'
'define q=tvrh2q(tv,RHprs)'
'set ccolor 1'
*'set cint 1'
'set cthick 6'
'd q'
'q w2xy 87.9 43.5'
x=subwrd(result,3)
y=subwrd(result,6)
'set line 1'
'draw mark 3 'x' 'y' 0.05'
'printim e:/kt/7/20130716/'tt'.png white x1000 y800'
'clear'
tt=tt+1
endwhile
; |
|