hannorein / rebound

💫 An open-source multi-purpose N-body code.
https://rebound.readthedocs.io/
GNU General Public License v3.0
844 stars 219 forks source link

run solar system with and without planet to planet gravity forces #775

Closed tckosvic closed 3 months ago

tckosvic commented 3 months ago

Would like a tip on how to run full solar system model using WHFast both with body-to-body gravitation forces and also without body-to-body gravitational forces. Just out of curiosity, I'd like to investigate starting from now after alot of years how solar system would be different.

Is there a setting in WHFast to make it only consider gravity to act for each planet with the sun (major mass) and ignore interactions between planets?

Alternatively, might there be a setting to effectively turn off gravity between the planets and only have gravity act between the planet and the sun?

Alternatively again, could some planet property be used to generate the desired effect?

I think I could compare the differences by just looking at the orbit differences for sometime in the future in some manner.

hannorein commented 3 months ago

You can simply set sim.N_active = 1. See https://rebound.readthedocs.io/en/latest/ipython_examples/Testparticles/