Description: Position-angle from p-vectors.
The source code for the Fortran routine PAP is available: