interface angleBetweenCelCoords
function angleBetweenCelCoordsSingle(vectorA, vectorB, isRadian) result(angle)
real(single) :: angle
real(single), intent(in) :: vectorA(2), vectorB(2)
logical, intent(in), optional :: isRadian
end function angleBetweenCelCoordsSingle
function angleBetweenCelCoordsDouble(vectorA, vectorB, isRadian) result(angle)
real(double) :: angle
real(double), intent(in) :: vectorA(2), vectorB(2)
logical, intent(in), optional :: isRadian
end function angleBetweenCelCoordsDouble
end interface angleBetweenCelCoords
Returns the angle between the two pair of input celestial coordinates. The default unit is radian (isRadian is T), but can be degree (isRadian is F).