Description: Angular separation from spherical coordinates.
The source code for the Fortran routine SEPS is available: