ncl批量读取nc文件

dmkk

我需要用1979到2017年的各月数据，但是日本气象中心下载的数据是多个nc文件，想请问一下各位前辈如何进行批量读取，我自己写了一个do循环读取的程序，但是显示 ，我的程序如下：begin
ys=1979
ye=2017
year=ispan(ys,ye,1)
ny=dimsizes(year)
mon = (/"01","02","03","04","05","06","07","08","09","10","11","12"/)
vs=new((/(ye-ys+1)*12,288,145/),float)
nt=ny*12
do i=0,ny-1
do j=0,11
  hcdc=addfile("/media/usb/Reanalysis/JRA/JRA55_1.25/Daily/fcst_surf125/High_cloud_cover/"+hcdc+"."+year(i)+"_"+mon(j)+".nc","r")
  v:=nc_vs->vs(:,{lev0s:lev0e},:)
  vmon=calculate_monthly_values(v, "avg", 0, False)
  vs(i*12+j,:,:)=vmon
end do
end do

这是我的数据和单个文件的内容：
。
希望各位前辈多多帮助！

dmkk

请各位前辈多多帮忙呀

Sebastian

Reading multiple supported files using 'addfiles'
https://www.ncl.ucar.edu/Applications/addfiles.shtml

BilibiliHoon

diri="E:/201607/201607/grib2/"
files=systemfunc("ls " + diri + "*.grib2")
这个是我用的，文件和脚本在不同文件夹里

白衣风月倾天下

files     = systemfunc("ls 路径/EIHR*.nc")
f = addfiles(files,"r")
我的文件名都是eihr开头的 所以这么写了 你可以根据自己实际文件名修改