Routine that sets all the orbit variables.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(g_orbit), | intent(inout) | :: | this | The routine is called on this g_orbit object. |
||
| real(kind=wp), | intent(in) | :: | r | The radial coordinate, (in Schwarzschild coordinates). |
||
| real(kind=wp), | intent(in) | :: | phi | The azimuthal angle, . |
||
| real(kind=wp), | intent(in) | :: | ur | The radial component of the 4-velocity, . |
||
| real(kind=wp), | intent(in) | :: | En | The energy per unit mass of the orbit, . |
||
| real(kind=wp), | intent(in) | :: | Lz | The angular momentum per unit mass of the orbit, . |
||
| real(kind=wp), | intent(in) | :: | chi | The osculaing orbits parameter, that changes by over a full radial cycle. |
subroutine set_orbit ( this, r, phi, ur, En, Lz, chi )
!! Routine that sets all the orbit variables.
class(g_orbit), intent(inout) :: this
!! The routine is called on this [[g_orbit]] object.
real(wp), intent(in) :: r
!! The radial coordinate, \(r\) (in Schwarzschild coordinates).
real(wp), intent(in) :: phi
!! The azimuthal angle, \(\phi\).
real(wp), intent(in) :: ur
!! The radial component of the 4-velocity, \(u^r\).
real(wp), intent(in) :: En
!! The energy per unit mass of the orbit, \(E\).
real(wp), intent(in) :: Lz
!! The angular momentum per unit mass of the orbit, \(L_z\).
real(wp), intent(in) :: chi
!! The osculaing orbits parameter, \(\chi\) that changes by \(2\pi\) over
!! a full radial cycle.
this%r = r
this%phi = phi
this%ur = ur
this%En = En
this%Lz = Lz
this%chi = chi
end subroutine set_orbit