X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;ds=sidebyside;f=files%2Fmatplotlib_hist.py;h=9521e267f68913a1c3135c1159b23b393d26171b;hb=2fca6b6df0ebd2df200cddb1e51049afa13399dd;hp=7bbf4e9af19bcd323d817d1fc989ca720f7ed977;hpb=e0e7c8852118f1e9fc1b66a5a45c86e21abef7b2;p=django_unres.git

diff --git a/files/matplotlib_hist.py b/files/matplotlib_hist.py
index 7bbf4e9..9521e26 100755
--- a/files/matplotlib_hist.py
+++ b/files/matplotlib_hist.py
@@ -14,8 +14,10 @@ with open('remd_all.stat','r') as f:
 
 if ncolumns==14:  
  x,y,s,r,ek,rms= np.loadtxt('remd_all.stat',usecols=(11,3,0,13,2,5),unpack=True)
+ x0,s0,r0= np.loadtxt('remd_all0.stat',usecols=(11,0,13),unpack=True) 
 else:
  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) 
 
 hall,binall=np.histogram(y,bins=40,density=False)
 
@@ -62,8 +64,8 @@ plt.xlabel('step*replica')
 replica=range(int(sys.argv[2]))
 colors = cm.rainbow(np.linspace(0, 1, len(replica)))
 for i,c in zip(replica,colors):
- yt=x[r==i]
- xt=(s+r*max(s))[r==i]
+ yt=x0[r0==i]
+ xt=(s0+r0*max(s0))[r0==i]
  plt.plot(xt,yt,'-',color=c)
    
 plt.savefig('remd_ex.png')