;
; if you exit idl, this will get you to after the fits.
; this assumes that the .sav files have already been created.
;
ftom=.3048D
designRadiusF=870
s2002RadiusF=869.883
designRadiusM=designRadiusF*ftom
s2002RadiusM=s2002RadiusF*ftom
designRadiusM=designRadiusF*ftom

@../p50init
@../dishinit
addpath,'/share/megs/phil/x101/p50/200311/pro'
.run inpscan19
;
; get fit results
;
restore,'fitsphere/fitiarAll.sav',/verb
;
restore,'fitsphere/iikeep_scan19.sav',/verb
;
restore,'img.sav',/verb
;
; at this point
; xyz  has 17million points, no fit clipping
; fitiar[10,2] has the 4param,3 param fits itererated 9 times
; get the best fit radius and centers
rad4=fitiar[9,0].fitcoef[3]
center4=fitiar[9,0].fitcoef[0:2]
rad3=fitiar[9,1].fitcoef[3]
center3=fitiar[9,1].fitcoef[0:2]