yt=rms[r==i]
xt=(s+r*max(s))[r==i]
tt=x[r==i]
- plt.scatter(xt,yt,c=tt,edgecolors='face',s=0.1,cmap=cm.rainbow)
+ plt.scatter(xt,yt,c=tt,edgecolors='face',s=0.1,cmap=cm.rainbow,vmin=Tremd[0],vmax=Tremd[-1])
plt.xlim(0,max(s)+max(s)*max(r))
plt.savefig('remd_step_rms.png')
plt.clf()
plt.xlabel('bath temperature [K]')
plt.ylabel('heat capacity')
-plt.xlim(Tremd[1]-10, Tremd[-1]+10)
+plt.xlim(Tremd[0], Tremd[-1])
plt.plot(x,y,'-',color=c)
plt.savefig('remd_cv.png')