- 积分
- 30
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-9-23
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2013-11-7 21:49:08
|
显示全部楼层
'reinit'
'set grid off'
'set grads off'
'set poli on'
'set parea 1 10 1 8'
'sdfopen D:\GrADS\NC\shum.mon.mean.nc'
'sdfopen D:\GrADS\NC\uwnd.mon.mean.nc'
'sdfopen D:\GrADS\NC\vwnd.mon.mean.nc'
'sdfopen D:\GrADS\NC\pres.mon.mean.nc'
'enable print D:\GrADS\summerflux.gmf'
'set lat 4 53'
'set lon 73 135'
'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6'
'set mpdset cnworld hires'
'set map 1 1 1'
'draw map'
'set ylopts 1 4 0.12'
'set xlopts 1 4 0.12'
'set cthick 4'
'set z 1 8'
'define qu=uwnd.2*shum.1/9.8'
'define qv=vwnd.3*shum.1/9.8'
'set z 1'
'set t 456 660'
'define uqs=vint(pres.4(z=1),qu,300)/1000'
'define vqs=vint(pres.4(z=1),qv,300)/1000'
'define sq=mag(uqs,vqs)'
'define uqs=vint(pres(z=1),qu,300)/1000'
'define vqs=vint(pres(z=1),qv,300)/1000'
'define sq=mag(uqs,vqs)'
'set t 6 8'
'define sq1=ave(sq,t=456,t=660,12)'
'modify sq1 seasonal'
'define sq2=sum(sq1,t=462,t=464)'
'set gxout shaded'
'd sq2'
'run cbarn'
'print'
'draw title 1986-2003 year summer average vapor flux china'
'printim D:\GrADS.gmf x1000 y800 white'
'disable print' |
|