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