- 积分
- 3168
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-10-24
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
求的是10年内季节平均水汽通量代码如下:
'reinit'
'set grid off'
'set grads off'
'set poli on'
'sdfopen D:\Grads\shum.mon.mean.nc'
'sdfopen D:\Grads\uwnd.mon.mean.nc'
'sdfopen D:\Grads\vwnd.mon.mean.nc'
'sdfopen D:\Grads\pres.mon.mean.nc'
'set fwrite d:\aut_shuiqi.grd'
'set gxout fwrite'
'set lat -10 55'
'set lon 10 140'
'set xlpos -20'
'set ylpos -20'
'set mpdset cnworld'
'set map 1 1 1'
'set cthick 5'
'set t 601 721'
'set z 1 8'
'define uq=shum.1*uwnd.2'
'define vq=shum.1*vwnd.3'
'define v=vwnd.3'
'define u=uwnd.2'
'set z 1'
'define uqs=vint(pres.4(z=1),uq,300)/1000'
'define vqs=vint(pres.4(z=1),vq,300)/1000'
'define sq=mag(uqs,vqs)'
'set t 1 12'
'define w1=ave(sq,t=601,t=721,12)'
'define w2=ave(uqs,t=601,t=721,12)'
'define w3=ave(vqs,t=601,t=721,12)'
'set t 1'
'define w1ave=sum(w1,t=9,t=11)'
'define w2ave=sum(w2,t=9,t=11)'
'define w3ave=sum(w3,t=9,t=11)'
*发现修改t=9,t=11、t=6,t=8或者t=3,t=5,出来的图都一样,求助如何修改,谢谢
'set grads off'
'set gxout shaded'
'set ccolor 1'
'set gxout vector'
'set ccolor 1'
'd w2ave;w3ave'
'cbarn 1 0 5.5 7.8'
'printim e:\aut.png x1000 y800 white'
'disable fwrite'
|
|