;$Id: peps.pro,v 1.3 2014/04/08 10:08:37 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=12,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  mv idl.ps ~/tex/tina/Decay/fig/peps.ps
;
!p.charsize=1.7
!x.margin=[8.8,.5]
!y.margin=[3.2,.5]
;
pc_read_ts,o=ts
dt=ts.t-ts.t(0)+1.
!p.multi=[0,1,2]
;!p.multi=0
;
!x.title='!8t!6'
plot_oo,dt,ts.epsM+ts.epsK,yr=[.00005,.01],xr=[1,100]
;oplot,dt,ts.epsK
oplot,dt,-.5*deriv(dt,ts.brms^2+ts.urms^2),col=122
;
xyouts,20,.002,'!6-d!8E!6/d!8t!6',col=122
xyouts,15,.0004,'!7e!6'
;plot,dt,ts.epsK/ts.epsM
xx=[20,70] & oplot,xx,.13/xx^1.667,li=3
;
plot,dt,ts.epsK/ts.epsM
;plot_oo,((ts.epsM+ts.epsK)/1e-15)^.25,xr=[1,100],yr=[100,2000]
;
END