;$Id: pcomp3.pro,v 1.1 2026/09/20 04:21:19 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=[10.8,.5]
!y.margin=[3.2,.5]
print,"$mv idl.ps ~/Overleaf/Oksana/Schwinger-effect-in-axion-inflation/fig/pbzmmax.eps"
!y.range=[2e-2,2.3]
;
;  sample to compare time series form different directories
;
!y.range=[-1,1]*1e-5
dir3='D4pi_mu01a_dt3em4'
dir2='D4pi_mu01a_10th_dt3em4'
;
!y.range=[-1,1]*2e-3
dir3='D4pi_mu01b_dt3em4'
dir2='D4pi_mu01b_10th_dt3em4'
;
!y.range=[-1,1]*2e-3
dir3='D4pi_mu01c_dt3em4'
dir2='D4pi_mu01c_10th_dt3em4'
dir1='1_mu01a'
;
!y.range=[-1,1]*2e-7
dir3='OT_cs20_512_5em3_dt2em4b'
dir2='OT10th_cs20_512_5em3_dt2em4b'
dir1='1_mu01a'
;
!y.range=[-1,1]*2.5e-6
dir3='OT_cs20_512_2em3_dt2em4'
dir2='OT10th_cs20_512_2em3_dt2em4'
dir1='1_mu01a'
;
pc_read_ts,o=ts1,datadir='../'+dir1+'/data',/double
pc_read_ts,o=ts2,datadir='../'+dir2+'/data',/double
pc_read_ts,o=ts3,datadir='../'+dir3+'/data',/double
;
;!p.multi=[0,1,3]
!x.range=[0,4.1]
;
!x.title='!8t!6'
!y.title='!7de!6!dK!n  and  !7de!6!dM!n'
;
;plot_io,ts1.t,ts1.epsK,/nodata
plot,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,ts3.t,ts3.epsK-ts2.epsK,col=55
print,ts3.epsK-ts2.epsK
;
;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
oplot,ts3.t,ts3.epsM-ts2.epsM,col=122
oplot,ts3.t,(ts3.epsM-ts2.epsM)+(ts3.epsK-ts2.epsK)
;
run=dir3
print,"$mv idl.ps ../tex/fig/pcomp3_"+run+".eps"
END