[FLASH-USERS] Continuing to evolve particles that leave the domain

[Klaus Weide]

On Fri, 27 Mar 2015, James Guillochon wrote:

> Hi all, I'm using active particles in my calculation, and I would like to
> continue to evolve the particle trajectories even once the particles have
> left the domain. I see that particles that leave the domain can be saved
> after they leave by enabling a flag, but it appears that their evolution at
> this point is frozen. It looks like a few relatively simple changes in
> Particles_advance may be able to continue their evolution, but I wasn't
> sure if evolving these particles would cause issues.

James,

Have you tried to set the keepLostParticles runtime parameter?

I have just noticed an apparent bug in Grid/GridParticles/gr_ptInit.F90.
That file contains the line

   call RuntimeParameters_get("keepLostParticles",gr_ptKeepLostParticles)

but later on another line

   gr_ptKeepLostParticles=.false.

that makes the runtime parameter ineffective.

Does fixing this help with your problem?

Klaus