;$Id: pcomp_ts_ampl.pro,v 1.7 2026/07/29 07:19:05 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=18,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
!p.charsize=1.7
!x.margin=[7.2,.5]
!y.margin=[3.2,.5]
;
;  sample to compare time series form different directories
;
dir1='Hy512d_noshock'         & nu1=2e-2 & ampl1=.2
dir2='Hy512d_noshock_ampl01'  & nu2=2e-2 & ampl2=.1
dir3='Hy512d_noshock_ampl005' & nu3=2e-2 & ampl3=.05
;
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=[0,3.3]
!x.title='!6'
!y.title='!6'
;
cs=1.
k1=1.
k0=5.*k1
;
s=1./(cs)
yr=[2e-3,3.]
yr=[2e-4,3.]
plot_io,ts1.t,s*ts1.urms,yr=yr
oplot,ts2.t,s*ts2.urms,col=55
oplot,ts3.t,s*ts3.urms,col=122
;
li=2
s=1./(cs^3*k0)
yr=[1e-5,1.]
oplot,ts1.t,s*ts1.epsK,li=li
oplot,ts2.t,s*ts2.epsK,col=55,li=li
oplot,ts3.t,s*ts3.epsK,col=122,li=li
;
siz=1.7
xyouts,siz=siz,0.1,1.2,'!6(a)'
xyouts,siz=siz,2.0,0.23,'!6Ma=!8u!6!drms!n/!8c!6!ds!n'
xyouts,siz=siz,1.5,0.012,'!7e!6!dK!n/!8c!6!s!ds!n!r!u3!n!8k!6!d0!n'
Ma='!8u!6!drms!n/!8c!6!ds!n'
;
;-----------------------------------------------------------------------------
!x.title='!8t!6'
s=1./(cs*k0)^2
yr=[3e-5,1.]
yr=[3e-7,1.]
plot_io,ts1.t,s*ts1.divu2m,yr=yr,li=li
oplot,ts2.t,s*ts2.divu2m,col=55,li=li
oplot,ts3.t,s*ts3.divu2m,col=122,li=li
;
s=1./(cs*k0)^2
yr=[1e-4,1e-1]
oplot,ts1.t,s*ts1.orms^2
oplot,ts2.t,s*ts2.orms^2,col=55
oplot,ts3.t,s*ts3.orms^2,col=122
;
tcrit=.5
print,'tcrit=',tcrit
print,ampl1,ts1.divu2m(findex(tcrit,ts1.t)),ts1.orms(findex(tcrit,ts1.t))^2,ts1.epsK(findex(tcrit,ts1.t))
print,ampl2,ts2.divu2m(findex(tcrit,ts2.t)),ts2.orms(findex(tcrit,ts2.t))^2,ts2.epsK(findex(tcrit,ts2.t))
print,ampl3,ts3.divu2m(findex(tcrit,ts3.t)),ts3.orms(findex(tcrit,ts3.t))^2,ts3.epsK(findex(tcrit,ts3.t))
;
xyouts,siz=siz,2.9,0.2,'!6(b)'
xyouts,siz=siz,1.5,0.10,'!6<(!9G.!8u!6)!u2!n>/!8c!6!s!ds!n!r!u2!n!8k!6!s!d0!n!r!u2!n'
xyouts,siz=siz,1.8,2e-5,'!6<(!9GX!8u!6)!u2!n>/!8c!6!s!ds!n!r!u2!n!8k!6!s!d0!n!r!u2!n'
;
!p.multi=0
print,"$mv idl.ps ~/Overleaf/Eva/BE_collapse/fig/pcomp_ts_ampl.eps"
END