grads批量处理nc文件求区域平均

stefan_scofield

海温资料是1965-2013年的逐月资料（每年的4,5,6月），我需要求每年4,5,6月nino3.4区域的平均海温，然后三个月再平均，就是每年得一个值。报错：Data request warning: request is completely outside file limits

gs如下：
'reinit'
'set gxout fwrite'
'set fwrite d:/sst/4-6nino.grd'
yy=1965
while(yy<=2013)
'sdfopen d:/sst/ersst.'yy'04.nc'
'sdfopen d:/sst/ersst.'yy'05.nc'
'sdfopen d:/sst/ersst.'yy'06.nc'
;'set z 1'
;'set t 1'
;'set lat -5 5'
;'set lon 170 240'
'a=aave(sst.1,lon=170,lon=240,lat=-5,lat=5)'
'b=aave(sst.2,lon=170,lon=240,lat=-5,lat=5)'
'c=aave(sst.3,lon=170,lon=240,lat=-5,lat=5)'
's=(a+b+c)/3.0'
'd s'
'close 3'
'close 2'
'close 1'
yy=yy+1
endwhile
'disable fwrite'
;

river

计算的那几句前面都加上define试试

stefan_scofield

river 发表于 2014-11-12 16:10
计算的那几句前面都加上define试试

没用，还是一样的错误

river

stefan_scofield 发表于 2014-11-12 18:45
没用，还是一样的错误

前面那几句加上
'set z 1'
'set t 1'
'set lat 1'
'set lon 1'

stefan_scofield

river 发表于 2014-11-12 19:30
前面那几句加上
'set z 1'
'set t 1'

一样的错误

river

stefan_scofield 发表于 2014-11-12 19:43
一样的错误

你先试试少几年的，不要一上来就好几十年的

stefan_scofield

river 发表于 2014-11-12 20:56
你先试试少几年的，不要一上来就好几十年的

没有用。只做一年也会出现两行warning：Data request warning: request is completely outside file limits
每增加一年就多两行相同的warning

river

stefan_scofield 发表于 2014-11-13 09:32
没有用。只做一年也会出现两行warning：Data request warning: request is completely outside file limi ...

先少做几年，把三个close注释掉试一下

stefan_scofield

river 发表于 2014-11-13 11:41
先少做几年，把三个close注释掉试一下

啥意思？怎么改close注释？

river

stefan_scofield 发表于 2014-11-13 14:59
啥意思？怎么改close注释？

三个close语句前面加 * 号，相当于注释行