【已解决】month_to_season函数报错

潇潇

想求80-99年冬季的月平均值，用到month_to_season函数，程序如下：f=addfile("E:/data/ncep/hgt/hgt.mon.mean.nc","r")

hgt=f->hgt(:,5,:,:)
;printVarSummary(hgt)

month_abbr = (/"","Jan","Feb","Mar","Apr","May","Jun","Jul","Aug","Sep","Oct","Nov","Dec"/)
time =f->time
    time@units = "hours since 1800-1-1 00:00:00"
utc_date = cd_calendar(time, 0)
year   = tointeger(utc_date(:,0))   
month  = tointeger(utc_date(:,1))   
MMYYYY=month_abbr(month)+sprinti("%0.4i", year)
;print(MMYYYY)

ymStrt="Mar1979"
ymLast="Feb2000"
iStrt=ind(MMYYYY.eq.ymStrt)
iLast=ind(MMYYYY.eq.ymLast)
;print(iStrt)
;print(iLast)
hgtx=hgt(iStrt(0):iLast(0),:,:)
printVarSummary(hgtx)
hgtx!0="time"
hgtx!1="lat"
hgtx!2="lon"

hgt_DJF=month_to_season(hgtx,"DJF")
;printVarSummary(hgt_DJF)



提示错误：

file:///C:\Users\Administrator.WIN-S3THKM90QAH\AppData\Roaming\Tencent\Users\29678350\QQ\WinTemp\RichOle\$HDJ@Y3AK%BUJ{_$JY%@$E1.png

zongqing123

我说的month——to——season怎么用不了，原来是数据类型的问题

潇潇

没人吗？ 菜鸟一枚，求大神指点！

井中月

  temp1&time = cd_calendar(temp1&time,1)，用这个函数可以把时间坐标转化为yyyymm格式，然后就可以用{}截取数据，比如你的hgt({198001:199912},:,:),就可以直接month_to_season(hgt({198001:199912},:,:),"DJF")。使用时最好经常pringVarSummary，确保是你想要的数据

潇潇

井中月 发表于 2016-7-14 08:54
temp1&time = cd_calendar(temp1&time,1)，用这个函数可以把时间坐标转化为yyyymm格式，然后就可以用{}截 ...

我已经将时间格式转换为MMYYYY了，最后用month_to_season函数,
hgt_DJF=month_to_season(hgt({May1979:Feb2000,:,:),"DJF")
但是还是报错
错误是：Variable (hgt_DJF) is undefined
但是我在官网上看，month_to_season函数左边的变量是根据右边的变量走的，不需要定义啊。。。我就不懂了

潇潇

井中月 发表于 2016-7-14 08:54
temp1&time = cd_calendar(temp1&time,1)，用这个函数可以把时间坐标转化为yyyymm格式，然后就可以用{}截 ...

问题已解决 其实是hgt的文件是short类型，转换成float就可以了~
不过还是要谢谢你！

nelflower

很有用谢谢