; look at the results of the spher fitting.
; assume fitiar had been restored from the .sav file
; 
;restore,'fitiarAll.sav',/verb
;ao9todish=.7   ; the starting offset p50 to dish
;         x0 y0  z0                radius
;a=size(xyz)
;ntot=a[2]
;
; start with all

elMax=34
maxDif=.3
iloop=2
ifit=1
coef=fitiar[iloop,ifit].fitcoef
rtp=blmxyztosph(xyzS,center=coef[0:2],/deg)
ii=where((abs(rtp[0,*] - coef[3]) lt maxdif) and (rtp[2,*] lt elMax),cnt)
;
; for ao9 x,y center .. doesn't work..
;
;cen=[0.,0.,265.176-.34]
;rtp=blmxyztosph(xyzS,center=cen,/deg)
;r0=265.176
;ii=where((abs(rtp[0,*] - r0) lt maxdif) and (rtp[2,*] lt elMax),cnt)

print,cnt
;
!p.multi=0
hor
ver
plot,xyzS[0,ii],xyzS[1,ii],psym=3
plot,xyzS[0,ii],xyzS[2,ii],psym=3
; histogram
;
elbin=.162 ; ..75 meters
nh=blmhistscan(rtp[2,ii],rtp[0,ii],elbin,hI,hrevI)
;
r0=265.176
cs=2.0
eps=.3
!p.multi=[0,1,3]
ver
plot,hI.binX,hI.binC,chars=cs
ver,r0-eps,r0+eps
plot,hI.binX,hI.avgB,chars=cs
oplot,hI.binX,hI.medB,col=colph[2]
oplot,[0,90],[r0,r0],col=colph[3]
;
ver
plot,hI.binx,hI.rmsB,chars=cs
;
; now split el bins into az steps
;
; bin daz= r0*sin(el)
r0=265.176
daz=1.
for i=1,32 do begin &$
	daz=.75/(r0*sin(i*!dtor)*!dtor) &$
	print,i,daz &$
endfor



;
; look at intensity
nh=blmhistscan(rtp[2,ii],in[ii],bin,hI,hrevI)
;
;
r0=265.176
cs=2.0
eps=.3
!p.multi=[0,1,3]
ver
plot,hI.binX,hI.binC,chars=cs
ver
plot,hI.binX,hI.avgB,chars=cs
oplot,hI.binX,hI.medB,col=colph[2]
;
ver
plot,hI.binx,hI.rmsB,chars=cs
;
end