Custom operator on eccentricity and semimajor axis about reboundx

[dtamayo]

Hi,

Unfortunately it's not well defined to just set the eccentricity or semimajor axis (you have to choose what to do with the remaining elements), so REBOUND doesn't let you do that.

If you just want exponential damping of the eccentricity and semimajor axis, you can just use modify_orbits_forces that's already implemented, and just set the appropriate tau_a and tau_e using your tidal parameters (tau_e = 1/de/dt).

If that's not good enough, you're doing the right thing by implementing an operator, but you have to set xyzzy and vxvyvz for an operator. One option would be to follow appendix A of Lee and Peale 2002.

[huhell]

Hi,
Thank you very much for the clarification !