;$Id: pcomp_for.pro,v 1.6 2026/07/25 09:08:41 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=22,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
!p.charsize=1.7
!x.margin=[9.0,.5]
!y.margin=[3.2,.5]
print,"$mv idl.ps ~/Overleaf/Oksana/Schwinger-effect-in-axion-inflation/fig/pbzmmax.eps"
;
;  sample to compare time series form different directories
;
dir1='for512c'
dir1='for512e2'
dir1='for512d3'
dir1='for512d5'
dir1='for512d6c'
dir2='for512b'
dir2='for512d'
dir2='for512e3'
dir2='for512d4'
dir2='for512d6b'
dir2='for512d6e'
dir3='for512a'
dir3='for512e'
dir3='for512d'
dir3='for512d5'
dir3='for512d6g'
dir4='Hy512f_noshock' & nu4=1e-2
dir4='Hy512g_noshock' & nu4=4e-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'
;
!p.multi=[0,1,3]
!p.multi=[0,1,4]
!x.range=[0,20.5e-2]
!x.range=[.140,4.8e-1]
;
yr=[1e-6,1e-5]
plot_io,ts1.t,ts1.dt,yr=yr,ytit='dt'
oplot,ts2.t,ts2.dt,col=122,li=0,th=3
oplot,ts3.t,ts3.dt,col=155,li=0,th=4
;
yr=[1e-2,5e3]
yr=[1e-1,5e3]
plot_io,ts1.t,ts1.rhomax,yr=yr,ytit='rho'
oplot,ts2.t,ts2.rhomax,col=122,li=0,th=3
oplot,ts3.t,ts3.rhomax,col=155,li=0,th=4
;
oplot,ts1.t,ts1.rhomin
oplot,ts2.t,ts2.rhomin,col=122,li=0,th=3
oplot,ts3.t,ts3.rhomin,col=155,li=0,th=4
;
;plot,ts1.t,ts1.orms,yr=yr
;oplot,ts2.t,ts2.orms,col=122,li=2,th=3
;oplot,ts3.t,ts3.orms,col=155,li=1,th=4
;oplot,ts4.t,ts4.orms,col=55,th=4
;
yr=[0.,.55]
yr=[1e-2,1]
yr=[1e-1,2]
nu1=2e-2
nu2=3e-2
nu3=nu1
;plot,ts1.t,ts1.urms,yr=yr
plot_io,ts1.t,ts1.urms,yr=yr,ytit='urms'
oplot,ts2.t,ts2.urms,col=122,li=0,th=3
oplot,ts3.t,ts3.urms,col=155,li=0,th=4
;
yr=[1e-17,1e-16]
plot_io,ts1.t,ts1.brms,yr=yr,ytit='brms'
oplot,ts2.t,ts2.brms,col=122,li=0,th=3
oplot,ts3.t,ts3.brms,col=155,li=0,th=4
;oplot,ts4.t,ts4.orms/ts4.urms,col=55,th=4
;
!p.multi=0
print,"$mv idl.ps ~/Overleaf/Axel/2phaseDynamo/fig/pcomp_for.eps"
END