;$Id: pcomp_ts_Grho.pro,v 1.1 2026/07/29 06:18:09 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
;
;  nu     ratm   ratm2    ratdm    genmax    dismax    dynmax     N      epsK    divurms     orms    knu   kNy   ampl  ratdx   umax  run
;2.0e-02  0.632  0.969   1.4e-01   1.2e-01   8.8e-02   2.8e-03   512   4.5e-02   1.0e+00   6.9e-02   8.4   128   0.20   0.36   2.83  Hy512d_noshock
;2.0e-02  0.332  0.893   2.9e-02   2.3e-01   1.7e-01   4.6e-03   512   2.7e-02   8.0e-01   3.2e-02   7.3   128   0.20   0.58   2.74  Hy512d_noshock_k10
;1.0e-02  0.119  0.519   4.3e-02   8.2e-02   3.2e-01   7.9e-03   512   4.3e-02   1.4e+00   1.3e-01  13.9   128   0.20   0.31   2.93  Hy512h_noshock2
;1.0e-02  0.269  0.506   1.1e-02   1.6e-01   3.7e-01   1.8e-02   512   2.4e-02   1.0e+00   6.2e-02  12.0   128   0.20   0.44   2.81  Hy512h_noshock_k10
;
dir1='Hy512d_noshock_k10' & nu1=1e-2 & k0_1=5.
dir2='Hy512d_noshock_k2' & nu2=1.4e-2 & k0_2=1.
dir1='Hy512d_noshock' & nu4=2e-2 & k0_4=2.5
;
;dir0='Hy512h_noshock_k2'  & nu0=1e-2 & k0_0=2.  ;(black dashed)
dir2='Hy512h_noshock_k10' & nu2=1e-2 & k0_2=10. ;(blue)
dir3='Hy512h_noshock_k20' & nu3=1e-2 & k0_3=20. ;(red)
;
dir1='Hy512h_noshock_G02b' & nu2=2e-2 & k0_3=5.
dir1='Hy512h_noshock_G1p5' & nu1=2e-2 & k0_2=5.
;
dir1='Hy512h_noshock_G02b' & nu1=2e-2 & k0_1=5.
dir1='Hy512h_noshock_G03b' & nu2=2e-2 & k0_2=5.
dir1='Hy512h_noshock_G05b' & nu3=2e-2 & k0_3=5.
dir1='Hy512h_noshock_G07b' & nu3=2e-2 & k0_3=5.
dir1='Hy512h_noshock2'    & nu1=1e-2 & k0_1=5.  ;(black)
dir1='Hy512h_noshock_G1p2' & nu1=2e-2 & k0_2=5.
dir2='Hy512h_noshock_G1p5' & nu1=2e-2 & k0_2=5.
dir3='Hy512h_noshock_G2b' & nu1=2e-2 & k0_2=5.
;
;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=ts4,datadir='../'+dir4+'/data'
;pc_read_ts,o=ts5,datadir='../'+dir5+'/data'
;pc_read_ts,o=ts6,datadir='../'+dir6+'/data'
;
!p.multi=[0,1,2]
!x.title='!6'
!y.title='!6'
!x.range=[0,40.]
!x.range=[0,8.]
!x.range=[0,3.2]
;
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,ts0.t,s*ts0.urms,li=2
  oplot,ts2.t,s*ts2.urms,col=55
  oplot,ts3.t,s*ts3.urms,col=122
 ;oplot,ts4.t,s*ts4.urms,col=122
 ;oplot,ts5.t,s*ts5.urms,col=155
 ;oplot,ts6.t,s*ts6.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_0) & oplot,ts0.t,s*ts0.epsK,li=2
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_0) & oplot,ts4.t,s*ts4.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=[9e-7,1.]
s=1./(cs*k0_1)^2 & plot_io,ts1.t,s*ts1.divu2m,yr=yr,li=li
;s=1./(cs*k0_0)^2 & oplot,ts0.t,s*ts0.divu2m,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_4)^2 & oplot,ts4.t,s*ts4.divu2m,col=122,li=li
;oplot,ts5.t,s*ts5.divu2m,col=155
;oplot,ts6.t,s*ts6.divu2m,col=122
;

yr=[1e-4,1e-1]
s=1./(cs*k0_1)^2 & oplot,ts1.t,s*ts1.orms^2
;s=1./(cs*k0_0)^2 & oplot,ts0.t,s*ts0.orms^2,li=1
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_4)^2 & oplot,ts4.t,s*ts4.orms^2,col=122
;oplot,ts5.t,s*ts5.orms^2,col=155
;oplot,ts6.t,s*ts6.orms^2,col=122
;
tcrit=2.0
print,k0_1,ts1.divu2m(findex(tcrit,ts1.t)),ts1.orms(findex(tcrit,ts1.t))^2
print,k0_2,ts2.divu2m(findex(tcrit,ts2.t)),ts2.orms(findex(tcrit,ts2.t))^2
print,k0_3,ts3.divu2m(findex(tcrit,ts3.t)),ts3.orms(findex(tcrit,ts3.t))^2
;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