;$Id: pcomp_divu.pro,v 1.2 2026/06/26 11:36:14 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=12,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
!p.charsize=1.7
!x.margin=[7.8,.5]
!y.margin=[3.2,.5]
;
;  sample to compare time series form different directories
;
dir1='Hy512d_noshock'
dir2='Hy512e_noshock'
dir3='Hy1024a_noshock'
dir4='Hy512g_noshock' & nu4=4e-2
dir5='Hy512f_noshock' & nu5=1e-2
;
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'
;
!p.multi=[0,1,2]
;!x.range=[0,300]
;
yr=[0.,5.]
plot,ts1.t,ts1.divu2m^.5,yr=yr
oplot,ts2.t,ts2.divu2m^.5,col=122,li=2,th=3
oplot,ts3.t,ts3.divu2m^.5,col=155,th=4
oplot,ts4.t,ts4.divu2m^.5,col=55,th=4
loadct,6
oplot,ts5.t,ts5.divu2m^.5,col=144,th=5
loadct,5
;
yr=[0.,8.]
nu1=2e-2
nu2=3e-2
nu3=nu1
plot,ts1.t*nu1,ts1.divu2m^.5/ts1.urms,yr=yr
oplot,ts2.t*nu2,ts2.divu2m^.5/ts2.urms,col=122,li=2,th=3
oplot,ts3.t*nu3,ts3.divu2m^.5/ts3.urms,col=155,th=4
oplot,ts4.t*nu4,ts4.divu2m^.5/ts4.urms,col=55,th=4
loadct,6
oplot,ts5.t*nu5,ts5.divu2m^.5/ts5.urms,col=144,th=5
loadct,5
;
!p.multi=0
print,"$mv idl.ps ~/Overleaf/Eva/BE_collapse/fig/pcomp_divu.eps"
END