;$Id: pcomp_ts.pro,v 1.5 2026/06/30 11:18:14 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
;
dir0='Hy2048a_noshock' & nu0=1e-2 & nx0=2048
dir1='Hy1024b_noshock' & nu1=1e-2 & nx1=1024
dir2='Hy1024c_noshock' & nu2=1.4e-2 & nx2=1024
dir3='Hy1024a_noshock' & nu3=2e-2 & nx3=1024
dir4='Hy512d_noshock' & nu4=2e-2 & nx4=512
dir5='Hy512e_noshock' & nu5=3e-2 & nx5=512
dir6='Hy512g_noshock' & nu6=4e-2 & nx6=512
;
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,3.3]
!x.title='!6'
!y.title='!6'
;
cs=1.
k1=0.5
k0=5.*k1
;
s=1./(cs)
yr=[2e-3,3.]
plot_io,ts1.t,s*ts1.urms,yr=yr,li=2
oplot,ts0.t,s*ts0.urms
oplot,ts2.t,s*ts2.urms,col=55
loadct,6
oplot,ts3.t,s*ts3.urms,col=144
oplot,ts4.t,s*ts4.urms,col=144,li=2
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,li=2
oplot,ts2.t,s*ts2.epsK,col=55
loadct,6
oplot,ts3.t,s*ts3.epsK,col=144
oplot,ts4.t,s*ts4.epsK,col=144,li=2
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.1,1.23,'!6(a)'
xyouts,siz=siz,2.0,0.3,'!6Ma=!8u!6!drms!n/!8c!6!ds0!n'
xyouts,siz=siz,1.5,0.0134*2,'!7e!6!dK!n/!8c!6!s!u3!n!r!ds0!n!8k!6!d0!n'
Ma='!8u!6!drms!n/!8c!6!ds0!n'
;
yy=1.8
yyy=1.7
si2=1.2
arrow,/data,1.72,yy,1.92,yy
arrow,/data,1.68,yy,1.48,yy
xyouts,2.0,yyy,siz=si2,'developed collapse'
xyouts,0.68,yyy,siz=si2,'early collapse'
;-----------------------------------------------------------------------------
!x.title='!8t!6'
s=1./(cs*k0)^2
yr=[3e-5,1.]*4.
plot_io,ts1.t,s*ts1.divu2m,yr=yr,li=2
oplot,ts0.t,s*ts0.divu2m
oplot,ts2.t,s*ts2.divu2m,col=55
loadct,6
oplot,ts3.t,s*ts3.divu2m,col=144
oplot,ts4.t,s*ts4.divu2m,col=144,li=2
loadct,5
oplot,ts5.t,s*ts5.divu2m,col=155
oplot,ts6.t,s*ts6.divu2m,col=122
;
Remeshcrit=.7
length=4.*!pi
s=1./(cs*k0)^2
yr=[1e-4,1e-1]
oplot,ts0.t,s*ts0.orms^2
dx=length/nx0 & ss=dx/nu0 & it=findex(Remeshcrit,ss*ts0.urms)
circ_sym,1.3,1 & oplot,[1,1]*ts0.t[it],[1,1]*s*ts0.orms[it]^2,ps=8
;
oplot,ts1.t,s*ts1.orms^2,li=2
dx=length/nx1 & ss=dx/nu1 & it=findex(Remeshcrit,ss*ts1.urms)
circ_sym,1.3,0 & oplot,[1,1]*ts1.t[it],[1,1]*s*ts1.orms[it]^2,ps=8
;
oplot,ts2.t,s*ts2.orms^2,col=55
dx=length/nx2 & ss=dx/nu2 & it=findex(Remeshcrit,ss*ts2.urms)
circ_sym,1.3,1 & oplot,[1,1]*ts2.t[it],[1,1]*s*ts2.orms[it]^2,ps=8,col=55
;
loadct,6
oplot,ts3.t,s*ts3.orms^2,col=144
dx=length/nx3 & ss=dx/nu3 & it=findex(Remeshcrit,ss*ts3.urms)
circ_sym,1.3,1 & oplot,[1,1]*ts3.t[it],[1,1]*s*ts3.orms[it]^2,ps=8,col=144
;
oplot,ts4.t,s*ts4.orms^2,col=144,li=2
dx=length/nx4 & ss=dx/nu4 & it=findex(Remeshcrit,ss*ts4.urms)
circ_sym,1.3,0 & oplot,[1,1]*ts4.t[it],[1,1]*s*ts4.orms[it]^2,ps=8,col=144
;
loadct,5
oplot,ts5.t,s*ts5.orms^2,col=155
dx=length/nx5 & ss=dx/nu5 & it=findex(Remeshcrit,ss*ts5.urms)
circ_sym,1.3,1 & oplot,[1,1]*ts5.t[it],[1,1]*s*ts5.orms[it]^2,ps=8,col=155
;
oplot,ts6.t,s*ts6.orms^2,col=122
dx=length/nx6 & ss=dx/nu6 & it=findex(Remeshcrit,ss*ts6.urms)
circ_sym,1.3,1 & oplot,[1,1]*ts6.t[it],[1,1]*s*ts6.orms[it]^2,ps=8,col=122
;
xyouts,siz=siz,2.8,0.3*4.,'!6(b)'
xyouts,siz=siz,1.5,0.18*4.,'!6<(!9G.!8u!6)!u2!n>/!8c!6!s!u2!n!r!ds0!n!8k!6!s!d0!n!r!u2!n'
xyouts,siz=siz,1.8,0.0014*4.,'!6<(!9GX!8u!6)!u2!n>/!8c!6!s!u2!n!r!ds0!n!8k!6!s!d0!n!r!u2!n'
xyouts,siz=si2,0.7,0.00008*4.,'!6Run G',col=122
xyouts,siz=si2,0.7,0.0022*4.,'!6Run A'
;
!p.multi=0
print,"$mv idl.ps ~/Overleaf/Eva/BE_collapse/fig/pcomp_ts.eps"
END