version = 10000

type = '2+1'

data = 'setup'

name = 'mercury-moon-sun'

comment = 'setup for mercury+moon system'

time_yr = 0

ref = 'mercury orbit from Wikipedia'

angles = 'absolute'

orbits = ['moon orbit', 'mercury orbit']
objects = ['mercury', 'moon', 'sun']

[scales]
mercurymass_kg = 3.30104e23
moondensity = 3.344
moonmass_mercurymass = 0.01
moonradius_ = '( 3 * moonmass / (4 * pi * moondensity))**(1/3)'

# moon orbit
['moon orbit']
angles = 'relative'
name = 'moon orbit'
# an_AU = 0.00256955 # for earth moon; mercury hill radius is 0.0012 AU
an_mAU = 0.5
en = 0.0
anomaly_deg = 0
ascending_deg = 0.0
inclination_deg = 0.0
periapsis_deg = 0.0

# mercury orbit
['mercury orbit']
name = 'mercury orbit'
an_km = 57909050.0
# P_day = 87.9691
en = 0.20563593
# needed for relative alignment to spin of the sun
# mean_deg = 174.796
# inclination_deg = 7.005
# ascending_deg = 48.331
# periapsis_deg = 29.124
mean_deg = 0.0
inclination_deg = 0.0
ascending_deg = 0.0
periapsis_deg = 0.0

# mercury
[mercury]
angles = 'relative'
name = 'mercury'
# euler_deg = [281.01,  0.034, 0.0]
euler_deg = [0.0,  0.0, 0.0]
M_ = 'mercurymass'
I = 0.346
S_km = 2439.7
k = 0.1
Q = 100
P_d = 58.646

# moon
[moon]
name = 'moon'
# euler_deg = [90.0, 6.687, 0.0]
euler_deg = [0.0, 0.0, 0.0]
# M_kg = 7.342e22 # earh moon too large
M_ = 'moonmass'
I = 0.3929
S_ = 'moonradius'
k = 0.030
Q = 27
P_d = 27.321661

# sun
[sun]
name = 'sun'
# euler_deg = [ 0.0,  7.25, 0.0]
euler_deg = [ 0.0,  0.0, 0.0]
M_kg = 1.9885e30
I = 0.070
S_km = 696342.0
k = 0.013
Q = 1e5
P_d = 25.05