This IDL facility returns the number and indices of the topological immediate neighbours
of a central pixel. The pixels
are ordered in a clockwise sense (when watching the sphere from the outside)
about the central pixel with the
southernmost pixel in first element. For the four pixels in the southern
corners of the
equatorial faces which have two equally southern neighbours the
routine returns the southwestern pixel first and proceeds clockwise.