Перебор файлов и сшивание трехмерных массивов

Question

У меня есть несколько файлов netCDF, которые имеют массив в форме (365,585,1386), и я пытаюсь прочитать в каждом новом массиве и соединить их вместе вдоль оси = 0; или добавить все дни года (365). Двумя другими измерениями являются широта и долгота, поэтому в идеале у меня должно быть несколько лет данных для каждой широты / долготы (каждый файл netCDF представляет собой один календарный год данных).

import glob
from netCDF4 import Dataset
import numpy as np

data = '/Users/sjakober/Documents/ResearchSpring2020/erc_1979.nc'
files = sorted(glob.glob('erc*'))

for x, f in enumerate(files):

    nc = Dataset(data, mode='r')   

    print(f)

    if x == 0:
        a = nc.variables['energy_release_component-g'][:]
    else:
        b = nc.variables['energy_release_component-g'][:]
        np.hstack((a, b))

    nc.close()