;$Id: pcomp2.pro,v 1.2 2026/09/20 04:21:12 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
;
dir3='D4pi_mu01a'
dir3='4pi_mu04pia'
dir3='D4pi_mu01a_10th'
dir2='4pi_mu01a'
dir2='D4pi_mu01a'
dir1='1_mu01a'
;
pc_read_ts,o=ts1,datadir='../'+dir1+'/data'
pc_read_ts,o=ts2,datadir='../'+dir2+'/data'
pc_read_ts,o=ts3,datadir='../'+dir3+'/data'
;
;!p.multi=[0,1,3]
!x.range=[0,.9]
!y.range=[2e-2,1.3]
;
!x.title='!8t!6'
!y.title='!7e!6!dK!n  and  !7e!6!dM!n'
;
plot_io,ts1.t,ts1.epsK,/nodata
;oplot,ts1.t,ts1.epsK,col=55
oplot,ts2.t,ts2.epsK,col=55,li=2
oplot,ts3.t,ts3.epsK,col=55,li=1,th=3
;
;oplot,ts1.t,ts1.epsM,col=122
oplot,ts2.t,ts2.epsM,col=122,li=2
oplot,ts3.t,ts3.epsM,col=122,li=1,th=3
;
print,"$mv idl.ps ../tex/fig/pcomp.eps"
END