;
; recompute errors setting x,y offsets to 0.. use same z offset
; see how this affects the errors.
;
; assume xyz and iikeepf3 have already been input (see initrestart or doit)
; center3 is the 3parameter fit center. just use z
;
mtocm=100.
axisF=1
th=-129.25996D
rotateGrid=th
cent=center3
cent[0:1]=0d ; for x,y offsets to 0
radius=designRadiusM
rtp=blmxyztosph(xyz,cent=cent,/deg)
radErrCm=(rtp[0,*]-radius)*mtocm
; max Err
maxErr=5.
ii=where(abs(radErrCm) lt maxErr,cnt)
useTr=0
verb=1
; rotating before gridding causes triangulate to fail try doing it after gridding
istat=griddatap50(xyz[*,ii],radErrCm[ii],img,xout,yout,useTr=useTr,tr=tr,gridLen=gridLen,$
                 axisF=axisF,rotateGrid=rotateGrid,verb=verb)