function skyToCartesian(ra, dec) result(vector)
real(double), intent(in) :: ra, dec
real(double) :: vector(3)
end function skyToCartesian
Returns a vector of direction cosines of the celestial coordinates
and
. NOTE
and
must be in radians.