Description: Solve for tangent point, spherical.
The source code for the Fortran routine TPORS is available: