+ plt.clf()
+ plt.xlabel('step*replica')
+ plt.ylabel('rmsd')
+ for i in replica:
+ yt=rms0[r0==i]
+ xt=(s0+r0*max(s0))[r0==i]
+ tt=x0[r0==i]
+ 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')