【已解决】用Grads求 月数据的冬季平均

学而不厌

'reinit'
'sdfopen D:\shuju\1_nc\hgt.mon.mean.nc'
'set lon 0 360'
'set lat 0 90'
'set lev 500'

a=12
while(a<=100)
   'define om=ave(hgt,t='a',t='a+2')'
   'd om'

   'set t 'a''
   time = subwrd(result,4)
   yearmon=substr(time,1,4)

   'gxprint D:\'yearmon'.png white'
   a=a+12
   'c'
endwhile
;


这个程序只能输出每一年的冬季平均，如果想输出多年平均的冬季该如何做呢？希望就在一个程序里解决，过去都是输出每年的冬季，再编写一个程序求平均。太麻烦了!

经过大神指导，解决了。

'reinit'
'sdfopen D:\shuju\1_nc\hgt.mon.mean.nc'
'set lon 0 360'
'set lat 0 90'
'set lev 500'
'set t 1 '
'define om=0'
bb=0
a=12
while(a<=100)
   bb=bb+1
   'om=om+ave(hgt,t='a',t='a+2')'
   a=a+12
endwhile
'set t 1'
'd om/'bb''
'gxprint D:\ave.png white'

学而不厌

想了许久，没有解决掉。求高手指导！

学而不厌

'reinit'
'sdfopen D:\shuju\1_nc\hgt.mon.mean.nc'
'set lon 0 360'
'set lat 0 90'
'set lev 500'
'set t 1 8'
'define om=0'
bb=0
a=12
while(a<=100)
   bb=bb+1
   'set t 'bb''
   'om=ave(hgt,t='a',t='a'+2)'
   'd om'
   'gxprint D:\'bb'.png white'
   a=a+12
   'c'
endwhile
'set t 1'
'd ave(om,t=1,t=8)'
'gxprint D:\ave.png white'
;

尝试这样子做，但是不对，AVE.png只是和最后一张图一样。

学而不厌

求Grads高手指导，主要问题是Grads好像没有数组的概念，不能将计算结果暂时保存。

男紫汗

om=om+ave(hgt,t='a',t='a'+2)/8.0

学而不厌

男紫汗 发表于 2018-5-24 09:39
om=om+ave(hgt,t='a',t='a'+2)/8.0

感谢指导，问题解决了。

学而不厌

男紫汗 发表于 2018-5-24 09:39
om=om+ave(hgt,t='a',t='a'+2)/8.0

感谢指导，问题解决了。

学而不厌

男紫汗 发表于 2018-5-24 09:39
om=om+ave(hgt,t='a',t='a'+2)/8.0

感谢指导，问题解决了。

学而不厌

男紫汗 发表于 2018-5-24 09:39
om=om+ave(hgt,t='a',t='a'+2)/8.0

感谢指导，问题解决了。