;$Id: pcomp_ts_nopretend.pro,v 1.2 2026/06/29 04:16:43 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
;
;  nu     ratm    ratdm     genmax    dismax    dynmax   run
;1.0e-02  0.628   1.7e-01   9.2e-02   8.6e-02   1.4e-02  Hy1024b_noshock
;1.4e-02  0.691   1.6e-01   8.8e-02   6.4e-02   5.4e-03  Hy1024c_noshock
;2.0e-02  0.657   1.2e-01   1.2e-01   8.6e-02   2.7e-03  Hy1024a_noshock
;2.0e-02  0.629   1.4e-01   1.2e-01   8.8e-02   5.9e-03  Hy512d_noshock
;3.0e-02  0.604   9.6e-02   1.6e-01   1.2e-01   2.4e-03  Hy512e_noshock
;4.0e-02  0.596   7.0e-02   1.8e-01   1.4e-01   2.0e-03  Hy512g_noshock
;
;dir0='Hy2048a_noshock' & nu1=1e-2
;dir1='Hy512d_noshock_ampl01' & nu1=1e-2
;dir2='Hy512d_noshock_ampl005' & nu2=1.4e-2
;dir3='Hy1024a_noshock' & nu3=2e-2
dir2='Hy512d_noshock' & nu2=2e-2
dir1='Hy512d_noshock_nopretend2' & nu1=2e-2
;dir5='Hy512e_noshock' & nu5=3e-2
;dir6='Hy512g_noshock' & nu6=4e-2
;
;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.range=[0,7.42]
!x.title='!6'
!y.title='!6'
;
cs=1.
k1=1.
k0=5.*k1
;
s=1./(cs)
yr=[2e-3,3.]
yr=[1.2e-4,3.]
plot_io,ts1.t,s*ts1.urms,yr=yr
;oplot,ts0.t,s*ts0.urms
oplot,ts2.t,s*ts2.urms,li=1
;loadct,6
;oplot,ts3.t,s*ts3.urms,col=144
;oplot,ts4.t,s*ts4.urms,col=122
;loadct,5
;oplot,ts5.t,s*ts5.urms,col=155
;oplot,ts6.t,s*ts6.urms,col=122
;
s=1./(cs^3*k0)
yr=[1e-5,1.]
;oplot,ts0.t,s*ts0.epsK
oplot,ts1.t,s*ts1.epsK
oplot,ts2.t,s*ts2.epsK,li=1
;loadct,6
;oplot,ts3.t,s*ts3.epsK,col=144
;oplot,ts4.t,s*ts4.epsK,col=122
;loadct,5
;oplot,ts5.t,s*ts5.epsK,col=155
;oplot,ts6.t,s*ts6.epsK,col=122
;
siz=1.7
xyouts,siz=siz,0.2,1.0,'!6(a)'
xyouts,siz=siz,5.0,0.20,'!6Ma=!8u!6!drms!n/!8c!6!ds!n'
xyouts,siz=siz,2.5,0.002,'!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=[1.2e-4,1.]
plot_io,ts1.t,s*ts1.divu2m,yr=yr
;oplot,ts0.t,s*ts0.divu2m
oplot,ts2.t,s*ts2.divu2m,li=1
;loadct,6
;oplot,ts3.t,s*ts3.divu2m,col=144
;oplot,ts4.t,s*ts4.divu2m,col=122
;loadct,5
;oplot,ts5.t,s*ts5.divu2m,col=155
;oplot,ts6.t,s*ts6.divu2m,col=122
;
s=1./(cs*k0)^2
yr=[1e-4,1e-1]
;oplot,ts0.t,s*ts0.orms^2
oplot,ts1.t,s*ts1.orms^2
oplot,ts2.t,s*ts2.orms^2,li=1
;loadct,6
;oplot,ts3.t,s*ts3.orms^2,col=144
;oplot,ts4.t,s*ts4.orms^2,col=122
;loadct,5
;oplot,ts5.t,s*ts5.orms^2,col=155
;oplot,ts6.t,s*ts6.orms^2,col=122
;
xyouts,siz=siz,6.8,0.34,'!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,4.0,9e-4,'!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_nopretend.eps"
END