求助 ncl读取特定的批量数据

云一林

ncl代码如下
begin
  ;f     = addfile ("uwnd.mon.mean.nc","r")  
  ;time= f->time(372:374)                    

  Din = "D:/123/MONTH/*/*/"
  ;Dout = "D:/123/MONTH/"

  files=systemfunc("ls " + Din + "CRA40_VVP_*_GLB_MONTH_V1_0_0.grib2")
  ;printVarSummary(files)
  fin = addfiles(files,"r")
  ListSetType (fin,"join")
  vv = fin[:]->VVEL_P0_L100_GGA0

  ;vv!0="time"                              
  ;vv&time=time                              
  printVarSummary(vv)

  ;system("rm -f " + Dout + "test.nc")  
  ;fout = addfile( Dout + "test.nc","c")

  ;fout->vv = vv(:,0:16,:,:)                 
end


文件路径如下：

代码是批量读取全部的数据的，请问如果想提取197901-201812共480个月的数据，有什么办法修改读取数据的代码（systemfunc）又不影响后续的ListSetType等的代码？在数据文件无法删除的情况下？

学气象的小弱鸡

楼主我也遇到了类似的问题，您知道怎么从一堆nc文件中只提取某一月份的数据吗?errst.v5.201001.nc 一直到202010.nc  怎么取其中的10月份，编程小白

对方的凤飞飞

请问楼主解决了吗

云一林

对方的凤飞飞 发表于 2020-11-6 16:33
请问楼主解决了吗

解决了，后来发现，读的时候把文件全部读进去了，输出的时候可以;fout->vv = vv(:,0:16,:,:)改成;fout->vv = vv(0:2,0:16,:,:)就可以输出前三个时间不输出所有时间了，或者读取数据处理数据时就能设置，但因为我ncl太烂了还不会弄，你如果数据量并不是特别大的话，并且是生成一个新的文件的话可以跟我一样在输出文件那里改

云一林

学气象的小弱鸡 发表于 2020-11-6 16:10
楼主我也遇到了类似的问题，您知道怎么从一堆nc文件中只提取某一月份的数据吗?errst.v5.201001.nc 一直到20 ...

不好意思啊，我不太清楚这个怎么弄，我才刚自学ncl，帮不了你。。。你去论坛搜吧，感觉应该挺多人做过提取月份的，或者你去ncl官网看下有没例子可以对照着写http://www.ncl.ucar.edu/Applications/