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).