0.5度降水数据求极值

老男孩

各位大虾：我才接触grads40小时,安装的是官方网站下载的opengrads，现紧急需要处理气象局下载的0.5度日降水数据。我下载了1960-2009年的数据，同时根据说明书修改了ctl文件。如下：
dset ^SURF_CLI_CHN_PRE_DAY_GRID_0.50-%y4
*
title  Daily Rain Analysis at 0.5
*
options little_endian template
*
undef -999.0  
*
xdef  140 linear  70.25  0.5
*
ydef   90 linear  15.25  0.5
*
zdef    1 linear  1 1
*
tdef  18263 linear  01jan1960 1dy
*
vars 2
rain   1 0 Total Rain Analysis (0.1mm/day)
nobs   1 0 number of gauge
ENDVARS
现在我需要的是每年的累计1天，2天，3天的最大值格点数据（后期导出到GIS软件中进一步分析）。憋了一下午使用搜索功能多次，尝试自己写GS也失败，自己不知怎么修改。但又要的很急。感觉短期无法解决，可能大家也都处理过这个问题，请给予帮助，以实现目标，谢谢。
'reinit'
'open /cygdrive/d/precipation/SURF_CLI_CHN_PRE_DAY_GRID_0.50.ctl'
'set x 1 140'
'set y 1 90'
'set z 1'
'set t 1 18263'
'set gxout fwrite'
'set fwrite annualmax'
n=1
while(n<=60)
'set t 1'
starttime=(n-1)*365+1
endtime=n*365
'define pj=max(rain,t='starttime',t='endtime')'
'd pj'
n=n+j
endwhile
'disable fwrite'
;

river

每年的累计1天，2天，3天的最大值格点数据是什么意思啊？而且做这个不一定非要grads，fortran也可以处理

老男孩

嗯，就是连续几天的降水量和，1天就是当天，2天就是当天+第二天，以此类推。
感觉grads还能快点上手。fortran目前更不敢想了。
先考虑1天的年最大值，这个相对难度低一些，如果只求几年年的1天最大值，我在opengrads gui里面结合一些手工命令就直接max(rain,t=1,t=365)这样一个命令一个命令敲出来了。但是50年，还是考虑要循环，同时还得考虑闰年天数不一样。
再要考虑2天累计，3天累计。我一开始考虑可以先计算2天或3天的均值，然后乘3，就是想要的2，3天的累计值。然后再求每年的累计最大值。

老男孩

请有兴趣大侠帮忙，可付辛苦费。

jiangli

你好，请问你上边的问题解决了，我现在也面临类似的数据处理难题