合并一个目录中的所有NETCDF文件

Question

我有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文件合并到一个文件中？