NCL将月资料转化为年资料

nuist.atm

想用month_to_annual函数把月资料转化为年资料，总报错说有一个变量没有定义，跪求各位大神帮忙看一下。
脚本如下
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin

  f=addfile("E:/eof/HadISST_sst.nc","r")

  x   =   f->sst(361:1669:1,:,:)

  x_year   =   month_to_annual(x,1)

  x_year&year = ispan(1900,2009,1)

  b=addfile("1900-2009.nc","c")

  b->sst_year=x_year

end

deemo7

  x   =   f->sst(361:1669:1,:,:)这句写为  x   =   f->sst(361:1669,:,:)就可以了吧，如果你的一维是时间的话

x_year   =   month_to_annual(x,1)

时间是第一维x_year   =   month_to_annual(x,0)应该是这样吧

nuist.atm

deemo7 发表于 2017-11-15 09:02
x   =   f->sst(361:1669:1,:,:)这句写为  x   =   f->sst(361:1669,:,:)就可以了吧，如果你的一维是时间 ...

试了下，还是提示 Variable (nyrs) is undefined，但是我不太懂nyrs这个变量是哪里出现的
错误命令是
Execute: Error occurred at or near line 10215 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl
fatal:["Execute.c":8640]:Execute: Error occurred at or near line 17 in file /cygdrive/e/eof/f+eof/year.ncl

x_year   =   month_to_annual(x,1)  我0和1都试了一下，还是一直有错{:lxm_21:}{:lxm_21:}

刘勇的下铺

可以参考下这个
假设有月平均数据 hgt(time, level, lat, lon)
可以用var= reshape(hgt,(/nyear,nmon,level,lat,lon/))  变成一个变为5维数组。

nuist.atm

刘勇的下铺 发表于 2017-11-15 11:36
可以参考下这个
假设有月平均数据 hgt(time, level, lat, lon)
可以用var= reshape(hgt,(/nyear,nmon,lev ...

嗯嗯，非常感谢，是我的年份算错了

nuist.atm

已经解决了，谢谢大家