; print the output of the fits
; assume the fitiar has been loaded from the save file
;
;restore,'fitiarAll.sav',/verb
; look at some of the radii
;
; try printing out the fits info
;
;' x0       y0       z0         Radius
;' xx.xxxx  xx.xxxx  xx.xxxx   xxx.xxxx  fit
;                                       errors
printFit:
;
openw,lunout,"fitResults.dat",/get_lun
printf,lunout,"Fit scan19 data to a sphere."
printf,lunout,"Iterate 10 times clipping at 3sig each iteration"
printf,lunout," "
lab=$
'  X0    sigX    Y0      sigY   Z0       sigZ   Radius   sigRadius fitErr Npnts'
labT=strarr(2)
labT[0]='Fit center and radius (all values in meters)'
labT[1]='Fit center. Radius fixed at 265.176 meters'

;' 1.2406 0.0005  1.3673 0.0005 262.8882 0.0021  263.2915 0.0020
for j=0,1 do begin &$
	printf,lunout,labT[j] &$
	printf,lunout,lab &$
	for i=0,toloop-1 do begin &$
		print,"j,i:",j,i &$
		lin1=string(format=$
'(f7.4,1x,f6.4,1x,f7.4,1x,f6.4,f9.4,1x,f6.4,1x,f9.4,1x,f6.4,4x,f6.4,1x,i8)',$
		fitiar[i,j].fitCoef[0],fitiar[i,j].sigmas[0],$
		fitiar[i,j].fitCoef[1],fitiar[i,j].sigmas[1],$
		fitiar[i,j].fitCoef[2],fitiar[i,j].sigmas[2],$
		fitiar[i,j].fitCoef[3],fitiar[i,j].sigmas[3],$ 
		fitiar[i,j].fitErr,fitiar[i,j].npnts) &$
		printf,lunout,lin1 &$
	endfor  &$
	printf,lunout," " &$
endfor
free_lun,lunout
;
end