nchoksi4456

Thanks Dan. What exactly is REBOUND doing then, when one applies the line `sim.particles[1].e += delta_e`? Because REBOUND does not throw an error. My impression based on the outputs was...

Thanks. Could you point me to where in the REBOUND source code these changes are applied?

Hi Hanno and Dan. I've uploaded a full working example "test.py" [here](http://ugastro.berkeley.edu/~nchoksi/test.py) as requested. One issue is that I am using my own version of modify_orbits_forces -- I've also included...

Hi Hanno and Dan. FYI the "test.py" script I uploaded ~40 minutes ago used the wrong random seed. It should be np.random.seed(0) to get the same initial conditions as I've...

Thanks Hanno and Dan, I'll try out these suggestions. Could I ask one more question right now, which might or might not be related? To record output, I was using...

Thanks Hanno, that's probably it... a(t) is wild.

Will do! And yes, I'm using the latest version of REBOUND, but haven't made any changes to it.

Here's another diagnostic test I tried which continues to add to my confusion. If I just replace ``` Noutputs = int(( max_time/dt)) # I've been using dt = 1 yr...

Hanno, you were right that the output timestep was less than the simulation timestep -- your suggested code caught that, thanks. But (somewhat surprisingly) this doesn't seem to fix the...

I modified your if statement (see the code snippet in my last post). The new if statement avoids integrating backwards by continuing to the next iteration of the loop (until...