Объединить все файлы NETCDF из одного каталога

У меня есть ежедневные данные в формате NETCDF, и каждый файл «.nc» включает данные за один год, в то время как у меня есть данные за период с 1948 по 2017 год. Имя файлов выглядит следующим образом: tmax.2m.gauss.1948.nc tmax.2m.gauss .1949.nc .. tмакс.2м.гаусс.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 в один файл?