;$Id: pcomp_ts.pro,v 1.5 2014/11/23 13:49:59 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=10,yoffset=3
  !p.charthick=2.3 & !p.thick=2.3 & !x.thick=2.3 & !y.thick=2.3
endif
thick0=6
;
;  mv idl.ps ~/tex/tina/Decay/fig/pcomp_ts.ps
;
ihel=1 ;(helical)
;
!p.charsize=1.6
!x.margin=[5.8,.7]
!y.margin=[3.1,.7]
!x.title='!8t v!6!dA0!n!8k!dM!60'
!p.multi=0
siz=1.2
;
;  sample to compare time series form different directories
;
dir1='short2304pm1_kf60b_noforce'
dir2='short2304pm100_kf60b_noforce'
;
pc_read_ts,o=ts1,datadir='../'+dir1+'/data'
pc_read_ts,o=ts2,datadir='../'+dir2+'/data'
;
!p.multi=0
!x.range=[1e-0,2.3e3]
!y.range=[3e-3,1]
dt=.1
;
kM0=60. & vA0=.15 & s=vA0*kM0
plot_oo,s*(ts1.t-ts1.t(0)+dt),ts1.brms,/nodata
oplot,s*(ts1.t-ts1.t(0)+dt),ts1.brms,thick=thick0
oplot,s*(ts2.t-ts2.t(0)+dt),ts2.brms,li=1
;
oplot,s*(ts1.t-ts1.t(0)+dt),ts1.urms,li=0,col=55,thick=thick0
oplot,s*(ts2.t-ts2.t(0)+dt),ts2.urms,li=1,col=55
;
oplot,s*(ts1.t-ts1.t(0)+dt),ts1.urms/ts1.brms,li=0,col=122,thick=thick0
oplot,s*(ts2.t-ts2.t(0)+dt),ts2.urms/ts2.brms,li=1,col=122
;
;  slopes
;
xx=[120,2000] & oplot,xx,.68/xx^.5
xx=[120,2000] & oplot,xx,.24/xx^.5,col=55
xyouts,1000,.026,'!8t!6!u-1/2!n',siz=siz
xyouts,1000,.009,'!8t!6!u-1/2!n',siz=siz,col=55
;
xyouts,30.,.45,'!8u!6!drms!n/!8v!6!dA!n',siz=siz,col=122
xyouts,30.,.10,'!8v!6!dA!n/!8c!6!ds!n',siz=siz
xyouts,20.,.042,'!8u!6!drms!n/!8c!6!ds!n',siz=siz,col=55
;
END