;$Id: pcomp_diss.pro,v 1.2 2026/10/05 07:22:10 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=19,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
!p.charsize=1.7
!x.margin=[10.8,1.5]
!y.margin=[ 3.2,0.5]
!y.range=[2e-2,2.3]
dot='!s!9!A .!n!6!r'
;
;  sample to compare time series form different directories
;
!y.range=[0,1]*1.5
dir0='OT_cs50_1024_8em2_x2'
dir1='OT_cs50_1024_4em2_x2b'
dir2='OT_cs50_1024_2em2_x2'
dir3='OT_cs50_1024_1em2_x2'
dir4='OT_cs50_1024_5em3_x2'
dir5='OT_cs50_1024_5em3'
;
pc_read_ts,o=ts0,datadir='../'+dir0+'/data',/double
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
pc_read_ts,o=ts4,datadir='../'+dir4+'/data',/double
pc_read_ts,o=ts5,datadir='../'+dir5+'/data',/double
;
!p.multi=[0,1,2]
!x.range=[0,2.0]
!y.range=[0,1.7]
;
!x.title='!8t!6'
!y.title='!7e!6!dK!n+!7e!6!dM!n'
;
plot,ts0.t,ts0.epsK+ts0.epsM
oplot,ts1.t,ts1.epsK+ts1.epsM,col=55 & loadct,6
oplot,ts2.t,ts2.epsK+ts2.epsM,col=144 & loadct,5
oplot,ts3.t,ts3.epsK+ts3.epsM,col=155
oplot,ts4.t,ts4.epsK+ts4.epsM,col=122
;
!x.range=[0,14.]
!y.range=[0,.16]
plot,ts5.t,ts5.epsK+ts5.epsM,/nodata
oplot,ts5.t,ts5.epsK+ts5.epsM,col=122
oplot,2.*ts4.t,.25*(ts4.epsK+ts4.epsM),col=122,li=2
;
!p.multi=0
print,"$mv idl.ps ~/Overleaf/Liubov/Dissipation-code-comparison/fig/pcomp_diss.eps"
END