;$Id: pfit.pro,v 1.3 2026/05/13 15:22:18 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=24,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
i1=21 & i2=29  ;(Orang2)
i1=12 & i2=17  ;(White)
i1=29 & i2=33  ;(Red)
i1=2 & i2=4  ;(Blue)
i1=15 & i2=22  ;(Green)
i1=32 & i2=35  ;(Orange)
!p.charsize=1.7
!x.margin=[9.4,.5]
!y.margin=[3.2,.3]
;
yr=[1e2,1e5]
!x.range=[1e-1,2e6]
!p.multi=[0,1,3]
plot_oo,t,tA,ps=-1,yr=yr
nt=n_elements(t)
pp1=fltarr(nt)
pp0=fltarr(nt)
for it=0,nt-1 do begin
  tt=t[it:*]
  ttA=tA[it:*]
  p=linfit(tt,ttA)
  pp1[it]=p[1]
  pp0[it]=p[0]
  oplot,tt,tt*p[1]+p[0],col=122
  print,t[it],p[1],p[0]
stop
endfor
;
i=indgen(nt)
good=where(i ge i1 and i le i2)
yr=[18,24]
yr=[9,9.5]
yr=[0,30]
plot_oi,t,1./pp1,ps=-1,yr=yr
pp1m=mean(pp1(good)) & oplot,t(good),t(good)*0+1./pp1m,col=122
pp2m=mean(1./pp1(good)) & oplot,t(good),t(good)*0+pp2m,col=155
print,'CM=',1./pp1m,pp2m
;
yr=[3e3,4e3] ;(Orang2)
yr=[5e3,1e4] ;(Orang2)
yr=[5e-1,1e5] ;(Orang2)
yr=[5e3,1e4] ;(Red)
yr=[5e0,1e4] ;(Blue)
yr=[5e-1,1e4] ;(Blue)
;plot_oi,t,pp0,ps=-1,yr=yr
plot_oo,t,pp0,ps=-1,yr=yr
pp0m=mean(pp0(good)) & oplot,t(good),t(good)*0+pp0m,col=122
print,'t0=',pp0m
;
print,"$mv idl.ps ~/Overleaf/Mattia/Isochrones/fig/pfit_"+dir+".eps"
END