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))
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))