;$Id: pcomp_mesh_Re_umax.pro,v 1.3 2026/07/30 07:48:00 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=11,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='Hy2048b_noshock2' & nu0=5e-3 & nx0=2048
dir1='Hy2048a_noshock' & nu1=1e-2 & nx1=2048
dir2='Hy1024b_noshock' & nu2=1e-2 & nx2=1024
;dir2='Hy1024c_noshock' & nu2=1.4e-2 & nx2=1024 ;(now suppressed)
dir3='Hy1024a_noshock' & nu3=2e-2 & nx3=1024
dir4='Hy512d_noshock' & nu4=2e-2 & nx4=512
;dir5='Hy512e_noshock' & nu5=3e-2 & nx5=512 ;(now suppressed)
dir5='Hy512g_noshock' & nu5=4e-2 & nx5=512
dir6='Hy1024d_noshock' & nu6=4e-2 & nx6=1024
;
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
!x.range=[0,3.3]
!x.range=[0,4.0]
!x.title='!8t!6'
!y.title='!6Re!dmesh!n'
;
cs=1.
k1=1.
k0=5.*k1
;
length=4*!pi
yr=[1e-1,10.]
dx=length/nx1 & s=dx/nu1 & plot_io,ts1.t,s*ts1.umax,yr=yr,li=2,/nodata
dx=length/nx0 & s=dx/nu0 & oplot,ts0.t,s*ts0.umax
dx=length/nx1 & s=dx/nu1 & oplot,ts1.t,s*ts1.umax,col=55
dx=length/nx2 & s=dx/nu2 & oplot,ts2.t,s*ts2.umax,col=55,li=2
loadct,6
dx=length/nx3 & s=dx/nu3 & oplot,ts3.t,s*ts3.umax,col=144
dx=length/nx4 & s=dx/nu4 & oplot,ts4.t,s*ts4.umax,col=144,li=2
loadct,5
dx=length/nx5 & s=dx/nu5 & oplot,ts5.t,s*ts5.umax,col=122,li=2
dx=length/nx6 & s=dx/nu6 & oplot,ts6.t,s*ts6.umax,col=122
;
;plot_io,ts1.t,s*ts1.divu2m,yr=yr,li=2,/nodata
;oplot,ts0.t,s*ts0.divu2m
;oplot,ts1.t,s*ts1.divu2m,col=55
;oplot,ts2.t,s*ts2.divu2m,col=55,li=2
;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=122,li=2
;oplot,ts6.t,s*ts6.divu2m,col=122
;
xr=!x.range
oplot,yr*0+2.7,yr,li=3
oplot,xr,xr*0+2.,li=3
;
!p.multi=0
print,"$mv idl.ps ~/Overleaf/Eva/BE_collapse/fig/pcomp_mesh_Re.eps"
END