I am simulating electron precipitation into a planetary atmosphere with a magnetic field. When I turn off the magnetic field everything works as expected, the electrons slow down in the atmosphere and eventually stop. When I turn the magnetic field on, however, this does not happen. Instead my electrons stop seemingly arbitrarily. They aren’t cut off by a fixed number of steps or some low energy threshold or (as far as I can tell) by some distance traveled.
I am using the DormandPrince745 stepper and have set the largest allowable step size to be several times the radius of my planet (as recommended in the docs). Has anyone encountered this problem before and knows how to solve it?
I am using Geant v11.0.0 on Mac OS X 11.4