- 积分
- 10
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-12-1
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
'reinit'
'sdfopen d:/monthly/uwnd.mon.mean.nc'
'sdfopen d:/monthly/vwnd.mon.mean.nc'
'sdfopen d:/monthly/hgt.mon.mean.nc'
'set mproj scaled'
'set vpage 0 11 0 8.5'
'set parea 2 8 1 7'
'set mproj scaled'
'draw rec 0 0 10.9 8.49'
'draw rec 0 0 10.9 8.49'
'draw rec 2 1 8 7 '
'set lev 200'
'set lon 100 150'
'set lat 15 60'
i=433
m=1984
while (i<=793)
'set t 'i
m=1984+((i-433)/12)
'define a=uwnd(lev-50)-uwnd(lev+50)'
'define b=vwnd.2(lev-50)-vwnd.2(lev+50)'
'define c=mag(a,b)'
'define d=hgt.3(lev-50)-hgt.3(lev+50)'
'define e=(c/(0.001*d))'
'define dux=cdiff(uwnd,x)'
'define dvy=cdiff(vwnd.2,y)'
'define dx=cdiff(lon,x)*3.1416/180'
'define dy=cdiff(lat,y)*3.1416/180'
'define p=dux/dx/637800'
'define q=dvy/dy/637800'
'define r=mag(p,q)*100000'
'define Dutton=(1.25*r+0.25*e*e+10.5)'
'set gxout shaded'
'set cint 5'
'set cmin 25'
'd Dutton'
'draw title 200.Dutton.'i'.'m''
'enable print d:/gradsst/200.Dutton.'i'.'m''
'print'
'disable print'
'c'
i=i+12
endwhile
;
|
|