Routine to convert the 3D position vector (x,y,z) of point into its position angles $(\theta,\phi)$ on the sphere with $x = \sin\theta\cos\phi$ , $y=\sin\theta\sin\phi$ , $z=\cos\theta$ .

Location in HEALPix directory tree: src/f90/mod/pix_tools.f90

FORMAT

call vec2ang( vector, theta, phi )

name&dimensionality	kind	in/out	description

vector(3)	DP	IN	three dimensional cartesian position vector (x,y,z). The north pole is (0,0,1)
theta	DP	OUT	colatitude in radians measured southward from north pole (in [0, $\pi$ ]).
phi	DP	OUT	longitude in radians measured eastward (in [0, $2\pi$ ]).