用GrADS求夏季海温平均温度

千里草

原始数据为HadISST环流中心的海温月平均数据，时间维1870年1月到2013年，精度：1度x1度。想用GrADS求1950年到2005年夏季（6-8月）平均SST，请教大家。
gs文件：'reinit'

'set fwrite F:\nc\HadSST\sst_JJA.grd'
'set gxout fwrite'

'sdfopen F:\nc\HadSST\HadISST_sst.nc'
'set lon -180 180'
'set lat -90 90'
'set z 1'

i=978----------------------------------------------------------------------------->从1950年6月起
while(i<1641)
'set t 'i''
'd ave(sst,t='i',t='i+2')'-------------------------------------------------------->夏季平均

i=i+12----------------------------------------------------------------------------->下一年
endwhile

'disable fwrite'

gs文件没报错能写出.grd文件，但配上.ctl文件后画出图不正确，请教了

ctl文件:
dset F:\nc\HadSST\sst_JJA.grd
undef    -999
title  HadISST
xdef  360  linear  0  1
ydef  180   linear  -89.5  1
zdef  1  levels  0
tdef  55  linear  12Z16jul1951   1yr
vars 1
sst 0 -999 Monthly 1 degree resolution SST
endvars

open F:\nc\HadSST\sst_JJA.ctl
set lon -180 180
set lat -90 90
set t 2
d sst
               

徐徐向前

谢谢，谢谢

千里草

已解决，'reinit'

'set fwrite F:\nc\HadSST\sst_JJA.grd'
'set gxout fwrite'



'sdfopen F:\nc\HadSST\HadISST_sst.nc'
*'set lon -180 180'
*'set lat -90 90'
'set z 1'

i=978

while(i<1641)
'set x 1 360'
'set y 1 180'
'set t 'i''
'd ave(sst,t='i',t='i+2')'

i=i+12
endwhile

'disable fwrite'
'reinit'

zglnew

谢谢，很好

听雨轩

顶一个                                

小傅之追梦

能自己找出错误，不错~~

千里草

小傅之追梦 发表于 2014-3-25 21:16
能自己找出错误，不错~~

也是请教别人的

xiaocaoqiqiao

弱弱地问一句，1870年1月到1950年6月不应该是966个月吗？为什么是978呢。1870年1月到1950年1月是80年，一共12×80=960个月，从1950年1月到6月是6个月，所以一共是966个月啊

千里草

xiaocaoqiqiao 发表于 2014-4-20 19:37
弱弱地问一句，1870年1月到1950年6月不应该是966个月吗？为什么是978呢。1870年1月到1950年1月是80年，一共 ...

观察的很仔细，计算的也很准确，但我用到的是1951年到2005年的资料，所以再加12就是978了

xiaocaoqiqiao

还有一个问题想请教一下，你增加的两句：
'set x 1 360'
'set y 1 180'  
这里的360和180是怎么确定的呢？