[SOFA] Returns the position angle of p-vector
ref. It is in the range from minus pi to plus
pi radians. The vectors
not be of unit length.
have the same dimensions, and their first dimension must have 3
elements. The result has the same dimensions as
except that the first dimension (equal to 3) is omitted.