;$Id: pkpm8192a_u.pro,v 1.8 2026/06/08 03:24:41 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=[8.3,1.3]
!y.margin=[3.2,.4]
;
xr=[.5,10.] & yr=[1.0e-3,2.5e-2]
!x.title='!7n!6!dM!n(!8t!6)'
!y.title='!8v!6!dA!n(!8t!6)'
;
file='/kpm.txt'
;
;-----------------------------------------------------------------------------
dir1='8192_k1em2_k1em0_vA1em1_nu1em5'        ;0
dir2='8192_k1em2_k1em0_vA1em1_nu1em5_alp20'        ;0
dir3='8192_k1em2_k1em0_vA1em1_nu1em5_alp20_cookup5' & s3=1.   ;0
;
;  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)]
;
a=rtable('../'+dir3+file,head=1,6)
t=a[1,*] & vA=sqrt(2.*a[2,*]) & xiM=a[3,*]
oplot,xiM,s3*vA,col=122,th=4
xx1=xiM[0] & yy1=vA[0] & yy2=max(vA) & xx2=xiM[findex(yy2,vA)]
;
xx=.6 & dx=.2 & siz=1.4
legend,siz=siz,xx,dx,10.^(-2.55),0,'!7a!8!s!dB!n!r!u!6ine!n=-2, !8k!6!dp!n=1'
legend,siz=siz,xx,dx,10.^(-2.70),0,'!7a!8!s!dB!n!r!u!6ine!n=-20, !8k!6!dp!n=1',col=55
legend,siz=siz,xx,dx,10.^(-2.85),0,'!7a!8!s!dB!n!r!u!6ine!n=-20, !8k!6!dp!n=2',col=122
;
;-----------------------------------------------------------------------------
!p.position=[.65,.65,.95,.95]
dir3='8192_k1em2_k1em0_vA1em1_nu1em5_alp20_cookup4' & s3=1.32   ;0
dir2='8192_k1em2_k1em0_vA1em1_nu1em5_alp20'        ;0
power,'_mag','_kin',k=k,spec1=magn1,spec2=kine1,i=n,tt=t1,/noplot,/lks,datatopdir='../'+dir2+'/data'
!p.charsize=1.4
!x.title='!8k!6'
!y.title='!8E!6!dM!n(!8k,t!6)'
plot_oo,k[1:*],magn1[1:*,0],/noerase,yr=[1e-15,1e-5],li=0
;oplot,k[1:*],magn1[1:*,14],li=1
;oplot,k[1:*],magn1[1:*,17],li=1
oplot,k[1:*],magn1[1:*,18],li=1
oplot,k[1:*],magn1[1:*,21],li=1
oplot,k[1:*],magn1[1:*,24],li=2
good=[0,18,21,24]
print,t1(good)
;
!p.position=0
print,"$mv idl.ps ~/Overleaf/Mattia/Isochrones/fig/pkpm8192a_u.eps"
END