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))
Dictionary of 4-D numpy array: how to optimize?