如何将grd格式的多年月平均海温数据输出成多年冬季平均的grd数据呢

supermann

各位前辈好！我的原始grd格式的海温数据是1950年1月到2016年12月的，现在我需要求出这些年的每年的冬季平均数据，也就是每年的12月到次年2月的，一共66年的，可是我的脚本运行不了，请各位大神们帮我看看哪里出错了，附上ctl文件和脚本，感谢大家了！

ctl:

dset d:/change/sst_195001_201612.grd
title mean monthly NCEP Reanalysis （1950-2016）
undef -999000000.000000
xdef 61 linear -150 1
ydef  11 linear -5 1
zdef 1 levels 1000
tdef 804 linear Jan1950 1mo
vars 1
ps 0 99 surface pressure
endvars


脚本：
'reinit'
'open D:/change/winner.ctl'
'set fwrite D:/change/sst_ave_winner.grd'
'set gxout fwrite'
tt=12
while(tt<=794)
'set t 'tt'  'tt+2''
'define a=ave(sst,t='tt',t='tt+2')'
'd a'
tt=tt+12
endwhile
'disable fwrite'

扫眉才子

你这几个维度都没固定，最不该固定的就是时间维你却写了。按照标准的写法，将lon lat都加上，将set t 'tt' 'tt+2' 改为set t 1 或者set t ‘tt’

mozhiting

你的ctl里面变量是ps不是sst啊