Description: Position-angle from spherical coordinates.
The source code for the ANSI C routine Pas is available: