Compute radial velocity magnitude
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=wp), | intent(in), | dimension(6) | :: | rv | [position,velocity] vector |
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