;$Id: pcomp_decay.pro,v 1.4 2012/01/27 21:33:34 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=24,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  mv idl.ps ~/tex/tina/Decay/fig/pcomp_decay.ps
;  mv idl.ps ~/tex/tina/Decay/fig/pcomp_decay_hel.ps
;
!p.charsize=1.7
!x.margin=[8.8,.5]
!y.margin=[3.2,.2]
;
;  sample to compare time series form different directories
;
dir1='128c5_noforce3'
dir2='128c5_noforce3_pm10'
dir3='128c5_noforce3_pm100'
;
dir1='short256pm1000b_noforce'
dir2='short256pm1000a_noforce'
dir3='128c5_noforce3_pm100'
dir3='short256pm1b_noforce'
;
;pc_read_ts,o=ts1,datadir='../'+dir1+'/data'
;pc_read_ts,o=ts2,datadir='../'+dir2+'/data'
;pc_read_ts,o=ts3,datadir='../'+dir3+'/data'
;
!p.multi=[0,1,2]
!x.range=[1.,500]
;
print,'start times:'
print,ts1.t(0)
print,ts2.t(0)
print,ts3.t(0)
;
eps=.1
!x.title=''
!y.title='!8E!6!dM!n'
 plot_oo,ts1.t-ts1.t(0)+eps,ts1.brms^2,yr=[1e-4,1.]
oplot   ,ts2.t-ts2.t(0)+eps,ts2.brms^2,col=122
oplot   ,ts3.t-ts3.t(0)+eps,ts3.brms^2,col=55
xx=[6,100] & oplot,xx,.04/xx^1.
;xx=[12,120] & oplot,xx,.4/xx^1.5,li=2
;
!x.title='!8t!6'
!y.title='!8E!6!dK!n'
 plot_oo,ts1.t-ts1.t(0)+eps,ts1.urms^2,yr=[1e-7,1e-1]
oplot   ,ts2.t-ts2.t(0)+eps,ts2.urms^2,col=122
oplot   ,ts3.t-ts3.t(0)+eps,ts3.urms^2,col=55
xx=[6,100] & oplot,xx,.002/xx^1.
;xx=[1.2,12] & oplot,xx,.04/xx^1.5,li=2
;
END