circularize() subtract pi/2 from theta coordinate, so correction is applied to the correct quadrant
circularize() subtract pi/2 from theta coordinate, so correction is applied to the correct quadrant

stggh/PyAbel,PyAbel/PyAbel,DhrubajyotiDas/PyAbel