Planet plysical and Hill's radius
This is a #523 (closed) spin off.
Also, we will need to inject Hill's spheres and sizes (rpl in the code)
I would propose something like
PlanetarySystem {
....
Integrator Symba7 {
MaxDeltaTime ...
...
}
Planets {
pl1 {
Name pl1
Init
{
Mass ...
SemiMajorAxis ...
Eccentricity ...
Inclination ...
...
Radius ... ; optional if the integrator is not symba7, necessary if it is symba7
HillSphere ... ; optional; if it's not given by the user, we calculate it as a*(1-e)*(mpl/(3M*))^(1/3)
}
...
}
pl2 {
...
}
...
}
}
Another option is
PlanetarySystem {
....
Integrator Symba7 {
MaxDeltaTime ...
...
}
Planets {
pl1 {
Name pl1
Init
{
Mass ...
SemiMajorAxis ...
Eccentricity ...
Inclination ...
...
}
Symba7Init ; necessary if the integrator is symba7
{
Radius ...
HillSphere ... ; optional; if it's not given by the user, we calculate it as a*(1-e)*(mpl/(3M*))^(1/3)
}
...
}
pl2 {
...
}
...
}
}