季节平均 区域平均 区域的时间序列

lsmm

现有数据54年的全国逐月降水资料，现求东北各季节的时间序列，（四张图 每个季节一张，横坐标要年份 如：1961 1962 以此类推）个人想法 先求季节平均 再求区域平均，但出不来。求助，本人新手，急需帮助！感谢大家！gs如下
'reinit'
'sdfopen E:\lsm\CN05.1_pre_1961_2014_month_05x05.nc'
'set fwrite E:\lsm\1.grd'
'set gxout fwrite'
a=3
while(a<=639)
'define spr=ave(pre,t='a',t='a+2')'
'set t 'a''
'define bb=aave(spr,lon=115,lon=135,lat=35,lat=55)'
'd bb'
a=a+1
endwhile
'disable fwrite'
;

lsmm

不好意思 不知道为什么不能选择脚本编辑
真的急需帮助 感谢大家了！

halexxxxx

楼主解决了吗

ldwloveyyz

set t 1 648
define aaa=tloop(aave(xxxxxx))
d aaa
将结果复制到电子表格中
查找替换非数字部分(将字母就剔除了)
直接在电子表格中求季度平均就可以了
不需要先平均再求时间序列，换个思路也可以解决

ldwloveyyz

我也是遇到这个问题了，但是我换了种方法

lsmm

ldwloveyyz 发表于 2019-4-6 11:14
set t 1 648
define aaa=tloop(aave(xxxxxx))
d aaa

谢谢，我最后也换了方法，已解决。谢谢你

LNYZ

遇见同样的问题，想问下楼主解决了吗？

lsmm

LNYZ 发表于 2019-4-8 15:46
遇见同样的问题，想问下楼主解决了吗？

单独算季节平均再求区域平均出来了，分开写gs

ldwloveyyz

lsmm 发表于 2019-4-8 20:02
单独算季节平均再求区域平均出来了，分开写gs

求分开写的GS学习下，感谢

lsmm

ldwloveyyz 发表于 2019-4-8 20:12
求分开写的GS学习下，感谢

'reinit'
'sdfopen E:\lsm\CN05.1_pre_1961_2014_month_05x05.nc'
'set fwrite E:\lsm\CN05.1_pre_1961_2014_month_NE.grd'
'set gxout fwrite'
'set lon 115 135'
'set lat 35 55'
'set lev 0 0'
a=3
while(a<=639)
a0=a+2
'set t 'a0 ''
'd ave(pre,t='a',t='a0')'
a=a+12
endwhile
'disable fwrite'
;

'reinit'
'open E:\lsm\pre.ctl'
'set fwrite E:\lsm\CN05.1.grd'
'set gxout fwrite'
'set lon 115 135'
'set lat 35 55'
'set lev 0 0'
'set t 1 54'
'define spr=aave(pre,lon=115,lon=135,lat=35,lat=55)'
'd spr'
;