Description: Position-angle from spherical coordinates.
The source code for the Fortran routine PAS is available: