rdot Function

private pure function rdot(rv) result(rd)

Compute radial velocity magnitude

Arguments

TypeIntentOptionalAttributesName
real(kind=wp), intent(in), dimension(6):: rv

[position,velocity] vector

Return Value real(kind=wp)


Called by

proc~~rdot~~CalledByGraph proc~rdot rdot proc~event_func event_func proc~event_func->proc~rdot

Contents

Source Code


Source Code

    pure function rdot(rv) result(rd)

    implicit none

    real(wp),dimension(6),intent(in) :: rv !! [position,velocity] vector
    real(wp) :: rd !! \( \dot{r} \)

    real(wp),dimension(3) :: r !! position vector
    real(wp),dimension(3) :: v !! velocity vector
    real(wp) :: rmag !! position vector magnitude

    r = rv(1:3)
    v = rv(4:6)
    rmag = norm2(r)

    rd = dot_product(r,v) / rmag

    end function rdot