cnames = ['CN', 'CM', 'CA', 'CY', 'CLN'] * 1000
nalpha, nmach, nbeta, nalt = 10, 20, 30, 40
numofc = len(cnames)
final_data = {}
for i in range(numofc):
    final_data[cnames[i]] = np.zeros((nalpha, nmach, nbeta, nalt))
Словарь массива 4-D numpy: как оптимизировать?