NCL如何从全年的气候态处理中截取3-7月的数据组成时间序列？

sysululu

  各位大神，在下小白一枚，想选取南海地区进行气候态区域平均，选取1979-2012年每年3-7月每一天的气候态组成一个时间序列

    f = addfile("/home/lewlee/OLR/olr.day.mean.f25.nc","r")
    time = f->time
    YYYY = cd_calendar(time, -2)
    t_index_start = ind(YYYY.eq.19790101)
    t_index_end   = ind(YYYY.eq.20121231)
    time_olr = f->time(time|t_index_start:t_index_end)
    lon_olr  = f->lon({105:118})
    lat_olr  = f->lat({4:21})
    hgt = short2flt(f->olr(t_index_start:t_index_end,{4:21},{105:118}))
    printVarSummary(hgt)
;------------------------------------------------------------------
    TIME    = cd_calendar(time, 0)         
    year    = toint( TIME(:,0) )
    month   = toint( TIME(:,1) )
    day     = toint( TIME(:,2) )
    ddd     = day_of_year(year, month, day)
    yyyyddd = year*1000 + ddd                           
    printVarSummary(yyyyddd)
   
    hClmDay = clmDayTLL(hgt, yyyyddd)
    printVarSummary(hClmDay)

现在写出这些，想先算全年的，可是总是报错

而且想做个循环截取3-7月成时间序列也没有头目，求各位大神拯救指点！

赤炼修

可以取34年每年3-7月逐日资料组成一个（34，153）的数组再来做你要的处理
函数不熟悉的话，最简单的思路，做个循环.假设每年365天为例
do i=0,33
   x2(i,0:152)=x1(59:211)
end do
x2第0维是年，第1维是3-7月的逐日

linhaiqiu123

真是个好帖子！

杨大只

作用很大，经常用到