;$Id: panastrophy.pro,v 1.5 2026/05/17 09:13:45 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
!p.charsize=1.7
!x.margin=[9.4,.5]
!y.margin=[3.2,.3]
;
!x.title='!6Lu'
!y.title='!8s!6(Lu)'
file='anastrophy.txt'
a=rtable(file,2,head=1)
plot_oo,xr=[30,3000],a[0,*],a[1,*],yr=[.01,.3],ps=4
p=linfit(alog(a[0,*]),alog(a[1,*]))
print,exp(p[0]),-p[1]
xx=[60.,1400.] & oplot,xx,exp(p[0])/xx^(-p[1]),li=3
xyouts,400.,.07,siz=siz,'4.9 !6Lu!u0.74'
;
print,"$mv idl.ps ~/Overleaf/Mattia/Isochrones/fig/panastrophy.eps"
END