Description: Project tangent plane to celestial, spherical.
The source code for the Fortran routine TPSTS is available: