Position Angle (East of North) between two points on a sphere.
Longitude and latitude of the two points. Quantities should be in
angular units; floats in radians.
The (positive) position angle of the vector pointing from position 1 to
position 2. If any of the angles are arrays, this will contain an array
following the appropriate numpy broadcasting rules.