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