;$Id: pcomp_ts_k0_G0.pro,v 1.1 2026/07/31 10:45:53 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=[6.0,.5]
!y.margin=[3.2,.5]
;
;  sample to compare time series form different directories
;
dir1='Hy512h_noshock_G0_k5_kcut100'  & nu1=1e-2 & k0_1=5.  ;(black)
dir2='Hy512h_noshock_G0_k10_kcut100' & nu2=1e-2 & k0_2=10. ;(blue)
dir3='Hy512h_noshock_G0_k20_kcut100' & nu3=1e-2 & k0_3=20. ;(red)
dir0='Hy512h_noshock_G0_k50_kcut100' & nu0=1e-2 & k0_0=50.  ;(black)
;
dir01='Hy512h_noshock_G0_k5'  & nu1=2e-2 & k0_1=5.
dir02='Hy512h_noshock_G0_k10' & nu1=2e-2 & k0_2=10.
dir03='Hy512h_noshock_G0_k20' & nu2=2e-2 & k0_3=20.
;
pc_read_ts,o=ts0,datadir='../'+dir0+'/data'
pc_read_ts,o=ts1,datadir='../'+dir1+'/data'
pc_read_ts,o=ts2,datadir='../'+dir2+'/data'
pc_read_ts,o=ts3,datadir='../'+dir3+'/data'
;
pc_read_ts,o=ts01,datadir='../'+dir01+'/data'
pc_read_ts,o=ts02,datadir='../'+dir02+'/data'
pc_read_ts,o=ts03,datadir='../'+dir03+'/data'
;
!p.multi=[0,1,2]
!x.range=[0,32.]
!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
  oplot,ts0.t,s*ts0.urms,col=155
;
  oplot,ts01.t,s*ts01.urms
  oplot,ts02.t,s*ts02.urms,col=55
  oplot,ts03.t,s*ts03.urms,col=122
;
li=2
yr=[1e-5,1.]
s=1./(cs^3*k0_1) & oplot,ts1.t,s*ts1.epsK,li=li
s=1./(cs^3*k0_2) & oplot,ts2.t,s*ts2.epsK,col=55,li=li
s=1./(cs^3*k0_3) & oplot,ts3.t,s*ts3.epsK,col=122,li=li
;
s=1./(cs^3*k0_1) & oplot,ts01.t,s*ts01.epsK,li=li
s=1./(cs^3*k0_2) & oplot,ts02.t,s*ts02.epsK,col=55,li=li
s=1./(cs^3*k0_3) & oplot,ts03.t,s*ts03.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.014,'!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'
li=2
yr=[3e-5,1.]
yr=[2e-9,1.]
s=1./(cs*k0_1)^2 & plot_io,ts1.t,s*ts1.divu2m,yr=yr,li=li
s=1./(cs*k0_2)^2 & oplot,ts2.t,s*ts2.divu2m,col=55,li=li
s=1./(cs*k0_3)^2 & oplot,ts3.t,s*ts3.divu2m,col=122,li=li
;
s=1./(cs*k0_1)^2 & oplot,ts01.t,s*ts01.divu2m,li=li
s=1./(cs*k0_2)^2 & oplot,ts02.t,s*ts02.divu2m,col=55,li=li
s=1./(cs*k0_3)^2 & oplot,ts03.t,s*ts03.divu2m,col=122,li=li
;
s=1./(cs*k0_1)^2 & oplot,ts1.t,s*ts1.orms^2
s=1./(cs*k0_2)^2 & oplot,ts2.t,s*ts2.orms^2,col=55
s=1./(cs*k0_3)^2 & oplot,ts3.t,s*ts3.orms^2,col=122
s=1./(cs*k0_0)^2 & oplot,ts0.t,s*ts0.orms^2,col=155
;
s=1./(cs*k0_1)^2 & oplot,ts01.t,s*ts01.orms^2
s=1./(cs*k0_2)^2 & oplot,ts02.t,s*ts02.orms^2,col=55
s=1./(cs*k0_3)^2 & oplot,ts03.t,s*ts03.orms^2,col=122
;
tcrit=2.0
print,'tcrit=',tcrit
fo="(4f10.6,2x,a)"
print,k0_1,ts1.divu2m(findex(tcrit,ts1.t)),ts1.orms(findex(tcrit,ts1.t))^2,ts1.epsK(findex(tcrit,ts1.t)),dir1,fo=fo
print,k0_2,ts2.divu2m(findex(tcrit,ts2.t)),ts2.orms(findex(tcrit,ts2.t))^2,ts2.epsK(findex(tcrit,ts2.t)),dir2,fo=fo
print,k0_3,ts3.divu2m(findex(tcrit,ts3.t)),ts3.orms(findex(tcrit,ts3.t))^2,ts3.epsK(findex(tcrit,ts3.t)),dir3,fo=fo
;print,k0_4,ts4.divu2m(findex(tcrit,ts4.t)),ts4.orms(findex(tcrit,ts4.t))^2
;
xyouts,siz=siz,2.8,0.2,'!6(b)'
xyouts,siz=siz,1.5,0.15,'!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,7e-6,'!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_k0.eps"
END