【已解决】计算冬季均值的gs文件，求改正，可以用

开始起飞1

先说下哈，这些数据都是每天为一个文件的，下面出现的ctl文件是一年一个，因为部分天数缺测，没想到别的好办法，请多多给出出主意吧，小弟在此谢过~
'reinit'
'set gxout fwrite'
'set fwrite  c:\GrADS19\data\winter.grd'
'open 1978data.ctl'
'set x 1 288'
'set y 1 180'
'define aaa3=ave(a,t=29,t=55)'***** 这是求得78年12月份的均值，部分时间缺测，根据数据的实际情况计算的
'close 1'
'open 1979data.ctl'
'set x 1 288'
'set y 1 180'
'define a791=ave(a,t=1,t=27)'*****这是79年1月份的均值
'define a792=ave(a,t=28,t=52)'*****这是79年2月份的均值
'define a793=ave(a,t=306,t=336)'*******这是79年12月份的均值，不知道可以不可以这样写（后面也打算根据这个样子写）
'd (aaa3+a791+a792)/3'*****3个月求平均
'close 1'
因为这最后是要提取每年的冬季均值的数据的，所以就先写了一部分，但不确定写的对不对，所以还没敢往下写，求各位大神们指导啊

开始起飞1

麻烦各位多多帮助啊

做个霸气的木头

这类问题，你先自己试，确实遇到问题再上来，不然你每写一段岂不是要发个帖子

开始起飞1

我试了下，如果*掉提取数据的那几步（gs如下），是可以出图的，但我不知道怎么去检验对不对
'reinit'
*'set gxout fwrite'
*'set fwrite  c:\GrADS19\data\winter2.grd'
'open 1978data.ctl'
*'set x 1 288'
*'set y 1 180'
'define aaa3=ave(a,t=29,t=55)'
'close 1'
'open 1979data.ctl'
*'set x 1 288'
*'set y 1 180'
'define a791=ave(a,t=1,t=27)'
'define a792=ave(a,t=28,t=52)'
'define a793=ave(a,t=306,t=336)'
'd (aaa3+a791+a792)/3'
*'close 1'
*'disable print'
;

baggio

没看懂什么意思！