| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer(kind=I4B), | intent(in) | :: | nbody |
number of bodies in the system) |
||
| logical, | intent(in), | dimension(:) | :: | lmask |
logical mask for active bodies in the system |
|
| real(kind=DP), | intent(in), | dimension(:) | :: | mu |
semi-major axis, mass, radius, and mu of the particle |
|
| real(kind=DP), | intent(in), | dimension(:) | :: | mass |
semi-major axis, mass, radius, and mu of the particle |
|
| real(kind=DP), | intent(in), | dimension(:) | :: | radius |
semi-major axis, mass, radius, and mu of the particle |
|
| real(kind=DP), | intent(in), | dimension(:, :) | :: | r_vec |
position and velocity vectors of the particle |
|
| real(kind=DP), | intent(in), | dimension(:, :) | :: | v_vec |
position and velocity vectors of the particle |
|
| real(kind=DP), | intent(inout), | dimension(:, :) | :: | acc |
Acceleration vector for all bodies |
|
| real(kind=DP), | intent(in), | dimension(:, :) | :: | rot |
rotation vector of the particle |
|
| real(kind=DP), | intent(in), | dimension(:) | :: | a |
semi-major axis, mass, radius, and mu of the particle |
|
| real(kind=DP), | intent(in), | dimension(:) | :: | emissivity |
particle characteristics for Yarkovsky calculations |
|
| real(kind=DP), | intent(in), | dimension(:) | :: | gamma |
particle characteristics for Yarkovsky calculations |
|
| real(kind=DP), | intent(in), | dimension(:) | :: | albedo |
particle characteristics for Yarkovsky calculations |
|
| real(kind=DP), | intent(in), | dimension(:) | :: | rot_k |
particle characteristics for Yarkovsky calculations |
|
| real(kind=DP), | intent(in) | :: | L_SUN_sys |
constants and parameters needed for Yarkovsky calculations |
||
| real(kind=DP), | intent(in) | :: | inv_c2 |
constants and parameters needed for Yarkovsky calculations |
||
| real(kind=DP), | intent(in) | :: | sigma_sys |
constants and parameters needed for Yarkovsky calculations |
||
| real(kind=DP), | intent(in) | :: | yark_radius_threshold_sys |
constants and parameters needed for Yarkovsky calculations |
Loop over all bodies to calculate the Yarkovsky effect.