;$Id: pcomp.pro,v 1.1 2026/06/26 09:06:24 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]
print,"$mv idl.ps ~/Overleaf/Oksana/Schwinger-effect-in-axion-inflation/fig/pbzmmax.eps"
;
;  sample to compare time series form different directories
;
dir1='Hy512d_noshock'
dir2='Hy512e_noshock'
dir3='Hy1024a_noshock'
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,2]
;!x.range=[0,300]
;
yr=[0.,.13]
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,th=4
oplot,ts4.t,ts4.orms,col=55,th=4
;
yr=[0.,.4]
nu1=2e-2
nu2=3e-2
nu3=nu1
plot,ts1.t*nu1,ts1.orms/ts1.urms,yr=yr
oplot,ts2.t*nu2,ts2.orms/ts2.urms,col=122,li=2,th=3
oplot,ts3.t*nu3,ts3.orms/ts3.urms,col=155,th=4
oplot,ts4.t*nu4,ts4.orms/ts4.urms,col=55,th=4
;
!p.multi=0
print,"$mv idl.ps ~/Overleaf/Eva/BE_collapse/fig/pcomp.eps"
END