| Type | Visibility | Attributes | Name | Initial | |||
|---|---|---|---|---|---|---|---|
| integer, | public, | parameter | :: | n_eoms | = | 6 |
size of EOM derivative vector [x,y,z,vx,vy,vz] |
| real(kind=wp), | public | :: | mu_earth | = | 3.9860043543609593E+05_wp |
|
|
| real(kind=wp), | public | :: | mu_moon | = | 4.9028000661637961E+03_wp |
|
|
| real(kind=wp), | public | :: | mu_sun | = | 1.3271244004193938E+11_wp |
|
|
| real(kind=wp), | public | :: | mu_jupiter | = | 1.266865349218008E+08_wp |
-- from https://naif.jpl.nasa.gov/pub/naif/generic_kernels/pck/gm_de431.tpc |
|
| real(kind=wp), | public | :: | rad_moon | = | 1737.4_wp |
radius of the Moon |
|
| real(kind=wp), | public | :: | rad_sun | = | 695700.0_wp |
sun radius (km) |
|
| real(kind=wp), | public | :: | rad_earth | = | 6378.0_wp |
earth radius (km) |
|
| integer, | public | :: | maxnum | = | 10000 |
integrator max steps |
|
| integer, | public | :: | grav_n | = | 8 |
max grav degree |
|
| integer, | public | :: | grav_m | = | 8 |
max grav order |
|
| integer, | public | :: | grav_frame | = | 1 |
1: iau_moon, 2: moon_pa (splined) |
|
| real(kind=wp), | public, | dimension(:), allocatable | :: | xscale_x0 | |||
| real(kind=wp), | public, | dimension(:), allocatable | :: | fscale_xf | |||
| real(kind=wp), | public | :: | fscale_rdot | = | 1.0_wp |
scale value for the rdot=0 constraint |
|
| logical, | public | :: | use_battin_gravity | = | .false. |
use Battin gravity formulation |
|
| character(len=:), | public, | allocatable | :: | mkspk_path |
for mkspk: path to mkspk executable. e.g. |
||
| integer, | public | :: | object_id | = | -50000 |
for mkspk: OBJECT_ID |
|
| character(len=:), | public, | allocatable | :: | object_name |
for mkspk: OBJECT_NAME |
||
| character(len=:), | public, | allocatable | :: | leapseconds_file |
for mkspk: LEAPSECONDS_FILE |
||
| character(len=:), | public, | allocatable | :: | segment_id |
for mkspk: SEGMENT_ID |
||
| integer, | public | :: | polynom_degree | = | 9 |
for mkspk: POLYNOM_DEGREE |
|
| integer, | public | :: | output_spk_type | = | 9 |
for mkspk: OUTPUT_SPK_TYPE |
|
| real(kind=wp), | public | :: | min_export_time_step | = | 0.0_wp |
minimum time (sec) step for exporting to mkspk file [to try and avoid spice interpolation issues for very small time steps] |
a CR3BP Halo state, to be used as an initial guess for the full force model.
| Type | Visibility | Attributes | Name | Initial | |||
|---|---|---|---|---|---|---|---|
| real(kind=wp), | public | :: | t | = | zero |
Time from Periapse (days from periapsis) |
|
| real(kind=wp), | public, | dimension(6) | :: | rv | = | zero |
Rx (km) Ry (km) Rz (km) Vx (km/s) Vy (km/s) Vz (km/s) [moon-centered earth-moon rotating frame] |