grads画水汽通量距平合成图

自由行

最近在做整层水汽通量的距平合成图，选取的是多雨年作为异常年，看了相关帖子，还是不是太确定gs里面的求多雨年和所有年的水汽通量求平均的语句有问题，不知道错在了哪里，请版主和大神求救。非常急需要用到整层水汽通量的分析。'reinit'
'open H:\shuiqitongliang\strongyear\2\strongshum.ctl'
'open H:\shuiqitongliang\strongyear\2\strongu.ctl'
'open H:\shuiqitongliang\strongyear\2\strongpres.ctl'

'set lon 60 160'
'set lat 10 50'
'set t 6'
'd qstr=ave(vint(pres.3(z=1),u.2*shum.1/1000,300)/9.8,t=1,t=6)'
*'d qu=ave(vint(pres.3(z=1),u.2*shum.1/1000,300)/9.8,t=1,t=6)'
'close 3'
'close 2'
'close 1'

'sdfopen H:\shuiqitongliang\strongyear\2\shum.mon.mean.nc'
'sdfopen H:\shuiqitongliang\strongyear\2\uwnd.mon.mean.nc'
'sdfopen H:\shuiqitongliang\strongyear\2\pres.mon.mean.nc'
'set lon 60 160'
'set lat 10 50'
'd qa=ave(vint(pres.6(z=1),u.5*shum.4/1000,300)/9.8,t=271,t=760,12)'
*'d quu=ave(vint(pres.6(z=1),u.5*shum.4/1000,300)/9.8,t=271,t=760,12)'

'set gxout vector'
'set arrscl 0.8 80'
'd  (qstr-qa)'






'enable print H:\shuiqitongliang\strongyear\2\shuiqi.gmf'
'print'
'disable print'

zhouling

这下载还要积分，恐怕没人会帮你，你还不如把代码贴出来，也许有人会看呢