kepde Function

private pure function kepde(de, dm, sigma0, sqrta, a, r)

Elliptic Kepler's equation written in terms of the eccentric anomaly difference. See Battin, eqn 4.43.

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: de

eccentric anomaly difference
real(kind=wp), intent(in) :: dm

mean anomaly difference
real(kind=wp), intent(in) :: sigma0
real(kind=wp), intent(in) :: sqrta
real(kind=wp), intent(in) :: a
real(kind=wp), intent(in) :: r

Return Value real(kind=wp)


Source Code

    pure function kepde(de, dm, sigma0, sqrta, a, r)

    implicit none

    real(wp),intent(in) :: de      !! eccentric anomaly difference
    real(wp),intent(in) :: dm      !! mean anomaly difference
    real(wp),intent(in) :: sigma0
    real(wp),intent(in) :: sqrta
    real(wp),intent(in) :: a
    real(wp),intent(in) :: r
    real(wp) :: kepde

    kepde = -dm + de + sigma0 / sqrta * (one - cos(de)) - &
            (one - r / a) * sin(de)

    end function kepde