yr=[.1,5]
plot_io,ts.t,sqrt(ts.rho2m/ts.rhom^2-1.),yr=yr
END