import numpy as np phminsYSTD = np.loadtxt("phmins_Esm75_YSTD.dat") phmaxsYSTD = np.loadtxt("phmaxs_Esm75_YSTD.dat") print(np.mean(phminsYSTD),np.std(phminsYSTD)) print(np.mean(phmaxsYSTD),np.std(phmaxsYSTD))