pas(al, ap, bl, bp)
[SOFA] Returns the position angle of b with respect to a.
al = longitude of point a in radians,
latitude of point a in radians,
bl = longitude of point b
bp = latitude of point b in radians. All
arguments must have the same dimensions, which are also the dimensions
of the result.