ax1.plot(t,Fs0,"-o",color="red",markersize=8,lw=3.5,label=r"$T=1.0$",zorder=2) t,Fs0,Fs1 = np.loadtxt("./glass/T0.72/Fs_T0.72_rho0.80.dat",comments='!', unpack=True ...