;$Id: pbalance.pro,v 1.3 2014/04/06 14:01:16 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=20,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  mv idl.ps ~/tex/tina/Decay/fig/pbalance_short2048pm1b4_noforce.ps
;  mv idl.ps ~/tex/tina/Decay/fig/pbalance_short2048pm100b4_noforce.ps
;  mv idl.ps ~/tex/tina/Decay/fig/pbalance_short2304pm1_kf60b_noforce.ps
;
!p.charsize=1.0
!x.margin=[8.8,.5]
!y.margin=[3.2,1.8]
;
pc_read_ts,o=ts
tt=smooth(ts.t,9)
EMdot=deriv(tt,smooth(ts.brms^2/2.,9))
EKdot=deriv(tt,smooth(ts.urms^2/2.,19))
 ujxbm=ts.ujxbm
epsM=ts.epsM
epsK=ts.epsK
epsT=epsM+epsK
;
cwd,run
!p.multi=[0,2,1]
;
!p.title='!6'+run
plot,tt,EMdot/epsT,yr=[-1.5,1]
oplot,tt,epsM/epsT,li=1
 oplot,tt,-ujxbm/epsT,li=2,col=55
 oplot,tt,-epsM/epsT-ujxbm/epsT,li=3,col=122
;
!p.title='!6'
plot,tt,EKdot/epsT,yr=[-1.5,1]
oplot,tt,epsK/epsT,li=1
 oplot,tt,+ujxbm/epsT,li=2,col=55
 oplot,tt,-epsK/epsT+ujxbm/epsT,li=3,col=122
;
;oplot,tt,EKdot/epsT,li=2
;oplot,tt,epsK/epsT,li=2
END