SEPS Subroutine

  1. astro_module.f90
  2. astro_module
  3. SEPS

public subroutine SEPS(al, ap, bl, bp, s)

Angular separation between two sets of spherical coordinates.

Status: vector/matrix support routine.

History

  • IAU SOFA revision: 2006 November 13

Arguments

TypeIntentOptionalAttributesName
real(kind=wp), intent(in) :: al

first longitude (radians)
real(kind=wp), intent(in) :: ap

first latitude (radians)
real(kind=wp), intent(in) :: bl

second longitude (radians)
real(kind=wp), intent(in) :: bp

second latitude (radians)
real(kind=wp), intent(out) :: s

angular separation (radians)

Calls

proc~~seps~~CallsGraph proc~seps SEPS proc~s2c S2C proc~seps->proc~s2c proc~sepp SEPP proc~seps->proc~sepp proc~pxp PXP proc~sepp->proc~pxp proc~pdp PDP proc~sepp->proc~pdp

Called by

proc~~seps~~CalledByGraph proc~seps SEPS proc~pmsafe PMSAFE proc~pmsafe->proc~seps

Contents

Source Code

SEPS

Source Code

    subroutine SEPS ( al, ap, bl, bp, s )

    implicit none

    real(wp),intent(in) :: al !! first longitude (radians)
    real(wp),intent(in) :: ap !! first latitude (radians)
    real(wp),intent(in) :: bl !! second longitude (radians)
    real(wp),intent(in) :: bp !! second latitude (radians)
    real(wp),intent(out) :: s !! angular separation (radians)

    real(wp) :: ac(3), bc(3)

    !  Spherical to Cartesian.
    call S2C ( al, ap, ac )
    call S2C ( bl, bp, bc )

    !  Angle between the vectors.
    call SEPP ( ac, bc, s )

    end subroutine SEPS