ncl批量读取nc文件

dmkk · 发表于 2018-4-22 16:19:43

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 dmkk 于 2018-4-22 16:19 编辑

我需要用1979到2017年的各月数据，但是日本气象中心下载的数据是多个nc文件，想请问一下各位前辈如何进行批量读取，我自己写了一个do循环读取的程序，但是显示 QQ截图20180422161203.png

，我的程序如下：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

这是我的数据和单个文件的内容：
QQ截图20180422160608.png

。
希望各位前辈多多帮助！

dmkk · 发表于 2018-4-22 16:20:42

请各位前辈多多帮忙呀

Sebastian · 发表于 2018-4-22 20:16:07

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

BilibiliHoon · 发表于 2018-4-22 21:14:58

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

白衣风月倾天下 · 发表于 2018-4-23 16:46:23

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

		自动登录	找回密码
密码			立即注册

[作图] ncl批量读取nc文件

登录后查看更多精彩内容~