;$Id: pkpm4096a_u.pro,v 1.5 2026/05/21 11:53:37 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,2.5]
;
xr=[5.,110.] & yr=[7.0e-4,3e-2]
xr=[1.,40.] & yr=[2.0e-5,3e-2]
!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='4096_k1em3_k1em1_vA1em1_nu2em4'        ;0
;dir2='4096_k1em3_k1em1_vA1em1_nu2em4_u5em2'  ;5e-2
;dir3='4096_k1em3_k1em1_vA1em1_nu2em5_u5em2'  ;5e-2
dir2='4096_k1em3_k1em1_vA1em3_nu2em6_u5em2'  ;7e-3
dir3='4096_k1em3_k1em1_vA1em1_nu2em6_u5em2'  ;5e-2
dir4='4096_k1em3_k1em1_vA1em2_nu2em6_u5em2'  ;5e-2
;
;  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
xx1=xiM[0] & yy1=vA[0] & yy2=max(vA) & xx2=xiM[findex(yy2,vA)]
oplot,[xx1,xx2],[yy1,yy1],li=1
oplot,[xx2,xx2],[yy1,yy2],li=1
arrow,/data,xx2,yy1,xx2,yy2
;
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
xx1=xiM[0] & yy1=vA[0] & yy2=max(vA) & xx2=xiM[findex(yy2,vA)]
oplot,[xx1,xx2],[yy1,yy1],li=1
oplot,[xx2,xx2],[yy1,yy2],li=1
arrow,/data,xx2,yy1,xx2,yy2
;
;a=rtable('../'+dir5+file,head=1,6)
;t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
;oplot,xiM,vA,col=155
;
xx=9.4 & dx=.5
;legend,xx,dx,10.^(-2.95),0,col=0,'0'
;legend,xx,dx,10.^(-2.97),0,col=55,'0.002'
;legend,xx,dx,10.^(-2.99),0,col=122,'0.003'
;
;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/pkpm4096a_u.eps"
END