;$Id: pkpm1024a_u.pro,v 1.10 2026/05/24 20:14:57 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
;
siz=1.7
si2=1.4
!p.charsize=1.7
!x.margin=[9.4,.5]
!y.margin=[3.2,2.5]
;
xr=[7.,8.] & yr=[1.0e-3,1.2e-3]
xr=[5.8,41.] & yr=[1.0e-4,1.4e-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_nu1em5b'   ;0
dir2='1024_k1em3_k1em1_vA7em3_nu1em5u6'  ;2e-3
dir3='1024_k1em3_k1em1_vA7em3_nu1em5u4'  ;3e-3
dir4='1024_k1em3_k1em1_vA7em3_nu1em5u7'  ;4e-3
dir5='1024_k1em3_k1em1_vA7em3_nu1em5u8'  ;7e-3
dir6='1024_k1em3_k1em1_vA7em3_nu1em5u9'  ;7e-3
dir7='1024_k1em3_k1em1_vA7em3_nu1em5u2'  ;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
;
loadct,6
a=rtable('../'+dir5+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,xiM,vA,col=144
;
a=rtable('../'+dir6+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,xiM,vA,col=144,li=2
;
a=rtable('../'+dir7+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,xiM,vA,col=144,li=1
;
xx=22. & dx=6. & siz=1.2 & y0=-2.95 & dy=.07
legend,siz=siz,xx,dx,10.^(y0-4*dy),0,col=144,'0.71'
legend,siz=siz,xx,dx,10.^(y0-5*dy),2,col=144,'0.86'
legend,siz=siz,xx,dx,10.^(y0-6*dy),1,col=144,'1'
loadct,5
;
;  plot start point
;
circ_sym,1.0,1
x=7.51 & y=1.166e-3
oplot,[1,1]*x,[1,1]*y,ps=8
;
print,[.002,.003,.004,.005,.006,.007]/.007
legend,siz=siz,xx,dx,10.^(y0-0*dy),0,'!8u!6!d0!n/!8B!6!d0!n=0'
legend,siz=siz,xx,dx,10.^(y0-1*dy),0,col=55,'0.29'
legend,siz=siz,xx,dx,10.^(y0-2*dy),0,col=122,'0.43'
legend,siz=siz,xx,dx,10.^(y0-3*dy),0,col=155,'0.57'
;
yy=7.8e-5
siz=1.6
si2=1.5
xyouts,siz=siz,5.90,yy,'6'
xyouts,siz=siz,7.90,yy,'8'
xyouts,siz=siz,19.00,yy,'20'
xyouts,siz=siz,28.60,yy,'30'
;xyouts,siz=si2,6.55,.00078,'0.0008'
;
print,"$mv idl.ps ~/Overleaf/Mattia/Isochrones/fig/pkpm1024a_u.eps"
END