get_orbit Subroutine

  subroutine get_orbit( this, r, phi, ur, En, Lz )
  !! Routine that gets all the orbit variables except for \(\chi\).
    class(g_orbit), intent(inout) :: this
    !! The routine is called on this [[g_orbit]] object.
    real(wp), intent(out) :: r
    !! The radial coordinate, \(r\) (in Schwarzschild coordinates).
    real(wp), intent(out) :: phi
    !! The azimuthal angle, \(\phi\).
    real(wp), intent(out) :: ur
    !! The radial component of the 4-velocity, \(u^r\).
    real(wp), intent(out) :: En
    !! The energy per unit mass of the orbit, \(E\).
    real(wp), intent(out) :: Lz
    !! The angular momentum per unit mass of the orbit, \(L_z\).

    r = this%r
    phi = this%phi
    ur = this%ur 
    En = this%En
    Lz = this%Lz
  end subroutine get_orbit