根据逐日的nc资料截取夏季时间

Sky221

各位前辈们，我想要截取300hPa1981-2010年6,7,8月份uwnd的nc数据，时间这该怎么截取呢？以下是所写的程序,请大佬帮忙看一下，谢谢！
;设置经度和纬度的范围
latstart = -90
latend = 90
lonstart = 0
lonend = 360
;高度为300hPa
level = 300

;读取文件
  list_of_files = systemfunc("ls D:/uwnd1981-2010/uwnd.*.nc")
  f = addfiles(list_of_files, "r")
  ListSetType(f, "cat")
  time = f[:]->time
  uwnd = f[:]->uwnd(::,{level},{latstart:latend},{lonstart:lonend})

Sky221

  ;设置经度和纬度的范围
  latstart = -90
  latend = 90
  lonstart = 0
  lonend = 360
  ;高度为300hPa
  level = 300
程序如下：
  ;读取u风场
  list_of_files1 = systemfunc("ls D:/uwnd1981-2010/uwnd.*.nc")
  f_u = addfiles(list_of_files1, "r")
  ListSetType(f_u, "cat")
  time1 = f_u[:]->time
  uwnd1 = f_u[:]->uwnd(:,{level},{latstart:latend},{lonstart:lonend})
   ;提取1981-2010年6/7/8月份的数据
    YYYYMMDD1 = cd_calendar(time1, -2)
    itClm_s1 = 19810101
    itClm_e1 = 20101231
    recClm_s1 = ind(YYYYMMDD1 .eq. itClm_s1)
    recClm_e1 = ind(YYYYMMDD1 .eq. itClm_e1)
    time1_1 = time1(recClm_s1:recClm_e1)
    ymdhms1 = ut_calendar(time1_1, 0)
    mm1 = floattoint(ymdhms1(:,1))
    iSummer1 = ind(mm1.ge.6 .and. mm1.le.8) ;提取夏季的数据

FDU.EDU

可以利用CDO来截取

Sky221

好的谢谢 已经解决了

myj0219

111111111111111111111111111

qiqicai

Sky221 发表于 2022-4-16 16:12
好的谢谢 已经解决了

请问怎么解决的