- 积分
- 1594
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-11-14
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 VanisherV 于 2013-3-23 23:14 编辑
求多年季节平均水汽通量,下面的程序我求月平均时没问题,求季节平均opengrads的即时显示窗口有图显示(不完整,少了水汽通量填图),但输出的却是空白图,不知道为什么啊。。
'reinit'
'sdfopen f:\shuiqi\uwnd.mon.1981-2010.ltm.nc'
'sdfopen f:\shuiqi\vwnd.mon.1981-2010.ltm.nc'
'sdfopen f:\shuiqi\shum.mon.1981-2010.ltm.nc'
'sdfopen f:\shuiqi\pres.mon.ltm.nc'
'set t 1 12'
'set lon 10 180'
'set lat -20 60'
'set z 1 8'
'set grads off'
'set grid off'
'set xlint 10'
'set ylint 5'
'set map 20 1 7'
'set rgb 20 27 222 27'
'set rgb 21 255 255 255'
'set rgb 22 102 153 255'
'set rgb 23 50 101 253'
'set rgb 24 51 50 254'
'set rgb 25 102 0 254'
'define uwnd=ave(uwnd,t=3,t=5,1)'
'define vwnd=ave(vwnd.2,t=3,t=5,1)'
'define shum=ave(shum.3,t=3,t=5,1)'
'define pres=ave(pres.4(z=1),t=3,t=5,1)'
'define qu=uwnd*shum/9.8'
'define qv=vwnd*shum/9.8'
'define qus=vint(pres,qu,300)*9.8'
'define qvs=vint(pres,qv,300)*9.8'
'set lon 39.999 160.001'
'set lat -14.999 45.001'
'set z 1'
'set clevs 100 200 300 400 '
'set ccols 21 22 23 24 25'
'set gxout shaded'
'd mag(qus/1000,qvs/1000)'
'cbarn'
'd qus/1000;qvs/1000'
'printim f:\shuiqi\spring.png white'
;
|
|