;$Id: pdecay.pro,v 1.1 2010/02/19 22:18:52 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=22,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  mv idl.ps ~/tex/tina/Kazantsev/fig/pdecay_256d_noforce2.ps
;  mv idl.ps ~/tex/tina/Kazantsev/fig/pdecay_u256d2_noforce2.ps
;
cwd,run
@parameters
!p.charsize=1.6
!p.multi=[0,1,2]
pc_read_ts,o=ts
tt=ts.t-t1
good=where(tt gt t1good)
slope=deriv(alog(ts.t-t1),alog(ts.brms^2))
slopem=mean(slope(good))
slopek=deriv(alog(ts.t-t1),smooth(alog(ts.urms^2),5))
slopekm=mean(slopek(good))
tgood=tt(good)
;
plot_oo,tt,ts.brms^2,tit=run,xtit='t-t0',yr=[.0001,.4]*.01
oplot,tt,ts.urms^2,li=2
plot,tt,slope,yr=[-1.5,-.6],xtit='t-t0'
oplot,tt,slopek,li=1
oplot,tgood,0.*tgood+slopem,col=122
print,slopem
END