;$Id: ppower_comp.pro,v 1.9 2026/06/07 09:13:18 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.8
!x.margin=[8.0,2.1]
!y.margin=[3.2,0.3]
;
;  sample to compare time series form different directories
;
dir1='8192_k1em2_k1em0_vA1em1_nu1em5' & istart4=2+6 ;green, running, CB=2.9-2.95
dir2='8192_k1em3_k1em1_vA9em3_nu1em5u2' & istart5=0+6*0 & i1_dir2=21 & i2_dir2=29  ;Orange
;
siz=1.7
si2=1.3
xout=3e-6
yout=2e-23
fo="(e8.1)"
;
;  times must come after parameters, because it overwrites t1 and t2.
;
power,'_mag','_kin',k=k1,spec1=magn1,spec2=kine1,i=n,tt=t1,/noplot,/lks,datatopdir='../'+dir1+'/data'
power,'_mag','_kin',k=k2,spec1=magn2,spec2=kine2,i=n,tt=t2,/noplot,/lks,datatopdir='../'+dir2+'/data'
nk=n_elements(k)
nt=n_elements(t1)
xr=[6e-4,60.]
yr=[2d-16,3d-6]
!p.multi=0
;
;-----------------------------------------------------------------------------
it0=0
it1=30
it2=38
it3=32
it4=43
s=10.
!x.title='!8k!6'
!y.title='!8E!6(!8k!6)'
plot_oo,k1[1:*],magn1[1:*,it0],xr=xr,yr=yr
oplot,k1[1:*],kine1[1:*,it0],li=1
print,'black: ',t1[it0],t1[it1],t1[it2]
;
oplot,k1[1:*],magn1[1:*,it1],li=0
oplot,k1[1:*],kine1[1:*,it1],li=1
;
oplot,k1[1:*],magn1[1:*,it2],li=0
oplot,k1[1:*],kine1[1:*,it2],li=1
xyouts,siz=si2,.20,6e-7,'!64600'
xyouts,siz=siz,7.,1e-7,'!8t!6=0'
;
oplot,k2[1:*],s*magn2[1:*,it3],li=0,col=122
oplot,k2[1:*],s*kine2[1:*,it3],li=1,col=122
xyouts,siz=si2,.06,1.7e-7,'!8t!6=0',col=122
xyouts,siz=si2,.05,5.8e-7,'!8t!6=10!u5!n'
print,t2[it3]
;
oplot,k2[1:*],s*magn2[1:*,it4],li=0,col=55
oplot,k2[1:*],s*kine2[1:*,it4],li=1,col=55
xyouts,siz=si2,.005,4e-8,'!8t!6=1.5!9X!610!u6!n',col=55
;
loadct,0
xx=[.01,1.5] & oplot,xx,1.4e-6*xx,col=188,li=2
loadct,5
;
siz=1.2
si2=1.0
xyouts,siz=siz,5.4,1.4e-10,'!6A4'
xyouts,siz=siz,1.7,1.2e-11,'!6C4',col=122
xyouts,siz=siz,0.45,1.0e-13,'!6Run C4',col=55
;
siz=1.7
xx=[.002,.012] & oplot,li=3,xx,40.*xx^5 & xyouts,.0010,5e-12,siz=siz,'!9A!8k!6!u5!n'
xx=[.001,.005] & oplot,li=3,xx,2e-4*xx^2 & xyouts,.0012,3e-9,siz=siz,'!9A!8k!6!u2!n'
;
print,"$mv -i idl.ps ~/Overleaf/Mattia/Isochrones/fig/ppower_comp.eps"
END