0

我有 NETCDF 格式的每日数据，每个“.nc”文件都包含一年的数据，而我有 1948 年到 2017 年的数据。文件名如下：tmax.2m.gauss.1948.nc tmax.2m.gauss .1949.nc .. tmax.2m.gauss.2017.nc

我想将它们全部合并到一个“.nc”文件中。我不知道如何准确地使用“CDO”或“ncrcat”。我写了这个：

for year in range(1948,2017):
    year_files = glob.glob('D:/Dissertation/Data/NCEP-NCAR/tmax.2m.gauss.}.nc'.format (year))
    subprocess.call(['ncrcat'] + year_files + ['-O'])
    ncfile = netCDF4.Dataset(year_files, 'r')

但我收到错误，它不起作用。有谁知道我如何将所有 70 个 NETCDF 文件合并到一个文件中？

4

0 回答 0