x,y,s,r= np.loadtxt('remd_all.stat',usecols=(7,3,0,9),unpack=True)
x0,s0,r0= np.loadtxt('remd_all0.stat',usecols=(7,0,9),unpack=True)
x,y,s,r= np.loadtxt('remd_all.stat',usecols=(7,3,0,9),unpack=True)
x0,s0,r0= np.loadtxt('remd_all0.stat',usecols=(7,0,9),unpack=True)