【求助】用addfiles读取多个nc文件并求逐日平均出错

sunny21

目标：读取1991-2020年850hPa vwnd，将时间维串联起来，用clmDayTLL()实现逐日气候态。代码：   
begin

latstart = 0
latend = 50
lonstart = 110
lonend = 120
level = 850

list_of_files = systemfunc("ls D:/lwsj/vwnd/vwnd.*.nc")
f = addfiles(list_of_files, "r")
ListSetType(f, "cat")
time = f[:]->time
vwnd = f[:]->vwnd(:,{level},{latstart:latend},{lonstart:lonend})

    YYYYMMDD = cd_calendar(time, -2)
    itClm_s = 19910101
    itClm_e = 20201231
    recClm_s = ind(YYYYMMDD .eq. itClm_s)
    recClm_e = ind(YYYYMMDD .eq. itClm_e)
    vwnd_Clm = vwnd(recClm_s:recClm_e,:,:) ;


    TIME    = cd_calendar(time, 0)          ; type float
    year    = toint( TIME(:,0) )
    month   = toint( TIME(:,1) )
    day     = toint( TIME(:,2) )
    ddd     = day_of_year(year, month, day)
    yyyyddd = year*1000 + ddd                                  ; needed for input
    vwndClmDay = clmDayTLL(vwnd_Clm,yyyyddd)
    vwndClmDay!0 = "year_day"
    vwndClmDay!1 = "lat"
    vwndClmDay!2 = "lon"

    delete(time)
    delete(year)
    delete(month)
    delete(day)
    delete(ddd)
    delete(yyyyddd)

报错：
warning:Aggregated dimension coordinate values are non-monotonic; check aggregated file ordering
warning:Aggregated dimension coordinate values are non-monotonic; check aggregated file ordering
fatal:Illegal Subscript. Only scalar values are allowed in subscript ranges.

vwnd_Clm = vwnd(recClm_s:recClm_e,:,:) 错误

ncl小白感谢各位大佬的帮助！

645345

我没有用过这个函数，感觉问题出在recClm_s和recClm_e，你可以试试笨方法，加print，一行一行地排查。希望能帮上你

    YYYYMMDD = cd_calendar(time, -2)
    print(YYYYMMDD)
    itClm_s = 19910101
    itClm_e = 20201231
    recClm_s = ind(YYYYMMDD .eq. itClm_s)
    recClm_e = ind(YYYYMMDD .eq. itClm_e)
    print(recClm_s)
    print(recClm_e)
    vwnd_Clm = vwnd(recClm_s:recClm_e,:,:) ;

你是最美七月天

问题应该是出在recClm_s和recClm_e，你试试从这之前print一下

回想之后

没引入包哦。

coliane

fatal:Illegal Subscript. Only scalar values are allowed in subscript ranges.
这句的意思是，你的索引并不是标量，也就是说recClm_s和recClm_e这两个参数可能是数列，建议print出来看看

mmyin_nju

你这个数据是不是每一个日期有四个时次？
    recClm_s = ind(YYYYMMDD .eq. itClm_s)
    recClm_e = ind(YYYYMMDD .eq. itClm_e)
    vwnd_Clm = vwnd(recClm_s:recClm_e,:,:) ;
所以recClm_s/e都不是标量

是冉冉升起的冉

    date = ind((YYYYMMDD .ge. itClm_s).and.(YYYYMMDD .le. itClm_e))
    vwnd_Clm = vwnd(date,:,:) ;
试试这样。

木昜人蕉

回想之后 发表于 2022-2-25 17:33
没引入包哦。

你好，请问是缺少了哪个包？