;$Id: pcomp_mesh_Re_urms.pro,v 1.1 2026/07/01 19:21:58 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=[6.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' & 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
!x.range=[0,3.3]
!x.title='!8t!6'
!y.title='!6Re!dmesh!n'
;
cs=1.
k1=1.
k0=5.*k1
;
length=4*!pi
yr=[1e-1,1.6]
dx=length/nx1 & s=dx/nu1 & plot_io,ts1.t,s*ts1.urms,yr=yr,li=2
dx=length/nx0 & s=dx/nu0 & oplot,ts0.t,s*ts0.urms
dx=length/nx2 & s=dx/nu2 & oplot,ts2.t,s*ts2.urms,col=55
loadct,6
dx=length/nx3 & s=dx/nu3 & oplot,ts3.t,s*ts3.urms,col=144
dx=length/nx4 & s=dx/nu4 & oplot,ts4.t,s*ts4.urms,col=144,li=2
loadct,5
dx=length/nx5 & s=dx/nu5 & oplot,ts5.t,s*ts5.urms,col=155
dx=length/nx6 & s=dx/nu6 & oplot,ts6.t,s*ts6.urms,col=122
;
xr=!x.range
oplot,yr*0+2.8,yr,li=3
oplot,xr,xr*0+.7,li=3
;
!p.multi=0
print,"$mv idl.ps ~/Overleaf/Eva/BE_collapse/fig/pcomp_mesh_Re.eps"
END