pc_read_ts,obj=ts
;
t=ts.t
eheatcool=-integr(ts.Gamm-ts.rhoLm,x=ts.t)
ekin=ts.ekintot
eth=ts.ethtot
etot=eheatcool+ekin+eth
;
;!p.multi=[0,1,2]
plot,t,etot,yr=[-1,1]*400.
oplot,t,eheatcool,col=122
oplot,t,eth,li=2
oplot,t,eth,li=2
;
;plot,t,ts.Gamm-ts.rhoLm
;
END