;$Id: pkpm1024a_k.pro,v 1.4 2026/06/06 19:11:33 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=[9.4,.5]
!y.margin=[3.2,.5]
;
xr=[7.,11.] & yr=[7.0e-4,1.3e-3]
xr=[7.,110.] & yr=[7.0e-5,1.3e-3]
!x.title='!7n!6!dM!n(!8t!6)'
!y.title='!8v!6!dA!n(!8t!6)'
;
file='/kpm.txt'
;
;-----------------------------------------------------------------------------
; 1024_k1em3_k1em1_vA7em3_nu1em5u1b ampluu=2e-2
; 1024_k1em3_k1em1_vA7em3_nu1em5u3 ampluu=14e-3
;dir3='1024_k1em3_k1em1_vA7em3_nu1em5u5'  ;1e-3
;
dir1='1024_k1em3_k1em1_vA7em3_nu1em5a'  ;k^3
dir2='1024_k1em3_k1em1_vA7em3_nu1em5b'  ;k^4
dir3='1024_k1em3_k1em1_vA7em3_nu1em5c'  ;k^5
dir4=''  ;7e-3
;
;  i   t      EMag       xiMag        CM          Lu
;  & CM=a[4,*] & Lu=a[5,*]
;
a=rtable('../'+dir1+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
plot_oo,xiM,vA,xr=xr,yr=yr,xtickv=[10,20]
;
a=rtable('../'+dir2+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,xiM,vA,col=55
;
a=rtable('../'+dir3+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,xiM,vA,col=122
;
;a=rtable('../'+dir4+file,head=1,6)
;t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
;oplot,xiM,vA,col=155
;
xx=35. & dx=5.
;legend,xx,dx,10.^(-3.10),0,col=0,'!8k!6!u3!n'
;legend,xx,dx,10.^(-3.20),0,col=55,'!8k!6!u4!n'
;legend,xx,dx,10.^(-3.30),0,col=122,'!8k!6!u5!n'
;
xx=35. & dx=8.
legend,xx,dx,10.^(-3.10),0,col=0,'!7a!s!8!dB!n!r!u!6sub!n=3'
legend,xx,dx,10.^(-3.25),0,col=55,'!7a!s!8!dB!n!r!u!6sub!n=4'
legend,xx,dx,10.^(-3.40),0,col=122,'!7a!s!8!dB!n!r!u!6sub!n=5'
;
siz=1.6
si2=1.5
;xyouts,siz=siz,7.96,.000659,'8'
;xyouts,siz=siz,8.95,.000659,'9'
;xyouts,siz=si2,6.55,.00078,'0.0008'
;
print,"$mv idl.ps ~/Overleaf/Mattia/Isochrones/fig/pkpm1024a_k.eps"
END