if ncolumns==14:
x,y= np.loadtxt('md.stat',usecols=(0,10),skiprows=10,unpack=True)
- x1,e,r,gy= np.loadtxt('md.stat',usecols=(0,3,5,12),skiprows=0,unpack=True)
+ x1,e,r,gy,nc= np.loadtxt('md.stat',usecols=(0,3,5,12,6),skiprows=0,unpack=True)
else:
x,y= np.loadtxt('md.stat',usecols=(0,6),skiprows=10,unpack=True)
x1,e,gy= np.loadtxt('md.stat',usecols=(0,3,8),skiprows=0,unpack=True)
plt.ylabel('RMSD')
plt.plot(x1,r,'.')
plt.savefig('md_rms.png')
+
+ plt.clf()
+ plt.xlabel('step')
+ plt.ylabel('fraction of native side-chain concacts')
+ plt.plot(x1,nc,'.')
+ plt.savefig('md_fracn.png')
\ No newline at end of file