关于ncl批量读取文件的问题

yiersan

想请教一下大家是怎样读取格式为数据时间的文件吗
我是这样写了：
path="anl_surf125."
fils=systemfunc("ls " + path + "*")
f=addfiles(fils,"r")
但ncl显示读取不了这些文件，请问大家是怎么处理的

wsx

begin
   diri  = "/mnt/e/data/hgt_daily/"
   files = systemfunc("ls " + diri + "hgt.*.nc" ) ;"文件名相同的部分，不一样的地方用*代替")
    f     = addfiles(files, "r")
    ListSetType(f, "cat")  
    ;var = new((/10,14610,37,73/), float)
    ;hh_100 = f[:]->hgt(:,{100},{0:90},{0:180});要合并的数据变量  
    ;hh_150 = f[:]->hgt(:,{150},{0:90},{0:180});要合并的数据变量
    hh_200 = f[:]->hgt(:,{250},:,:);要合并的数据变量
    ;hh_250 = f[:]->hgt(:,{250},:,:);要合并的数据变量
    ;hh_300 = f[:]->hgt(:,{300},{0:90},{0:180});要合并的数据变量
    ;hh_400 = f[:]->hgt(:,{400},{0:90},{0:180});要合并的数据变量
    ;hh_500 = f[:]->hgt(:,{500},{0:90},{0:180});要合并的数据变量
    ;hh_600 = f[:]->hgt(:,{600},{0:90},{0:180});要合并的数据变量
    ;hh_700 = f[:]->hgt(:,{700},{0:90},{0:180});要合并的数据变量
    ;hh_850 = f[:]->hgt(:,{850},{0:90},{0:180});要合并的数据变量
    ; hh_850 = f[:]->hgt(:,{850},:,:);要合并的数据变量
    ;hh = f[:]->$var(i)$(:,:,:,:)                  
    ;print(files)
;================================================================数据筛选
;================================================================数据输出
fout = addfile("/mnt/e/h250.nc", "c") ; create 文件用以输出
fout@title = "250hpa_hgt"      ; 给文件加上一个全域属性
fout->hgt = hh_200       ; 将 hh 存储到 fout 里 ，变量名命名为 h
;var = new((/14610,10,37,73/), float)
;var(:,0,:,:) = hh_100
;var(:,1,:,:) = hh_150
;var(:,2,:,:) = hh_200
;var(:,3,:,:) = hh_250
;var(:,4,:,:) = hh_300
;var(:,5,:,:) = hh_400
;var(:,6,:,:) = hh_500
;var(:,7,:,:) = hh_600
;var(:,8,:,:) = hh_700
;var(:,9,:,:) = hh_850



;fname = "/mnt/d/hgt_250.nc"
;system("rm -f "+fname)
;fout := addfile(fname,"c")
;fout->hgt=hh_250

end

LIqy

没办法识别文件类型吗？你看你的文件类型是.195802之类的，是不是不是nc文件呀？

yiersan

LIqy 发表于 2022-3-6 10:29
没办法识别文件类型吗？你看你的文件类型是.195802之类的，是不是不是nc文件呀？

确实不是，好像什么格式都不是，只是二进制文件

sugazz

我也读过类似数字结尾的，你先看一下能否用记事本查看文件，然后用asciiread函数读应该是可以的

yiersan

谢谢大家，安装了最新版的ncl后就可以读入了

myj0219

666666666666666666666666