<div dir="ltr">Hi Anshu, evolution of the particles is done via the Sinks module currently, which seems to have its own routines for evolving the particles. At the moment I am not using the sink functionality (removing gas from the grid), but in principle it's something I'd like to be able to do in the future, which is why I'm using that module rather than just the active module.<div><br></div><div>Anyway, if we ignore the sink functionality for now, I guess I'm not understanding why my current solution isn't working. Right now when the LOST flag would be applied, I instead just assign it to the first block on a processor within gr_ptLocalMatch:</div><div><br></div><div><font face="monospace, monospace">if(outside)blockID=gr_ptBlkList(1)</font><br></div><div><br></div><div>I've also tried assigning it to its original block, which doesn't work either. And I make sure that the LOST flag is not applied in gr_ptOneFaceBC (I just comment out the "particle(gr_ptBlk)=LOST" line).</div><div><br></div><div>Otherwise, there are no other locations in the code where particles can be marked as "LOST", as far as I can tell. Particles can be marked as "NONEXISTENT," but if this is happening I would think the particles wouldn't even be saved in the lost list, which would be a bug.</div><div><br></div><div>My guess would be is that the Sinks module is doing something else that is removing particles that leave the domain?</div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Mar 28, 2015 at 12:38 PM, Anshu <span dir="ltr"><<a href="mailto:adubey@lbl.gov" target="_blank">adubey@lbl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="auto"><div>Try adding another category such as OUTSIDE, with that particle should not get deleted.</div><div>Give it a number > maxblocks so that sort puts those particles at the end. When you pass the data structure to particle routines make sure the local count is such that these particles are not</div><div>processed. And then add a routine to treat them separately. It has been a couple of years since</div><div>I worked with this code, so I might have a few things wrong and a few details missing, but I think </div><div>in general this should work. </div><div><br>Best,<div>Anshu</div></div><div><div><div><br>On Mar 28, 2015, at 9:16 AM, James Guillochon <<a href="mailto:jfg@ucolick.org" target="_blank">jfg@ucolick.org</a>> wrote:<br><br></div><blockquote type="cite"><div><div dir="ltr"><div style="font-size:12.8000001907349px">Hi all,</div><div style="font-size:12.8000001907349px"><br></div><div style="font-size:12.8000001907349px">I'm having a bit of trouble getting this to work. I need to leave a particle assigned to a block at all times, and not marked as NONEXISTENT or LOST, resulting in its deletion. Unfortunately it seems like there's a number of places where particles can be mark as NONEXISTENT in the code, and changing them all seems difficult (there's a only a few places where the particle is marked as "LOST", which I've disabled).<br></div><div style="font-size:12.8000001907349px"><br></div><div style="font-size:12.8000001907349px">Anyone have any additional insight here? At the moment, particles are deleted as soon as they leave the domain, even with the changes I've made.</div><div style="font-size:12.8000001907349px"><br></div><div style="font-size:12.8000001907349px">Thanks,</div><div style="font-size:12.8000001907349px">- James</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Mar 28, 2015 at 12:15 PM, James Guillochon <span dir="ltr"><<a href="mailto:guillochon@gmail.com" target="_blank">guillochon@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi all,<div><br></div><div>I'm having a bit of trouble getting this to work. I need to leave a particle assigned to a block at all times, and not marked as NONEXISTENT or LOST, resulting in its deletion. Unfortunately it seems like there's a number of places where particles can be mark as NONEXISTENT in the code, and changing them all seems difficult (there's a only a few places where the particle is marked as "LOST", which I've disabled).<br></div><div><br></div><div>Anyone have any additional insight here? At the moment, particles are deleted as soon as they leave the domain, even with the changes I've made.</div><div><br></div><div>Thanks,</div><div>- James</div></div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Mar 27, 2015 at 4:58 PM, Anshu Dubey <span dir="ltr"><<a href="mailto:adubey@lbl.gov" target="_blank">adubey@lbl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">If you don't need them to interact with the mesh in anyway I don't see<div>any conceptual problem with advancing them. You will still have to make </div><div>sure that they are not treated the same as particles within domain</div><div>in any function. </div><div><br></div><div>Anshu</div></div><div class="gmail_extra"><div><div><br><div class="gmail_quote">On Fri, Mar 27, 2015 at 1:48 PM, James Guillochon <span dir="ltr"><<a href="mailto:jfg@ucolick.org" target="_blank">jfg@ucolick.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><span style="font-size:12.8000001907349px">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.</span><div><span style="font-size:12.8000001907349px"><br></span></div><div><span style="font-size:12.8000001907349px">Has anyone tried to do something like this before?</span><div style="font-size:12.8000001907349px"><br></div><div style="font-size:12.8000001907349px">Cheers,</div><div style="font-size:12.8000001907349px">- James</div><span><font color="#888888"><div><br></div>-- <br><div><div dir="ltr">James Guillochon<br><span style="font-family:arial;font-size:small">Einstein Fellow at the Harvard-Smithsonian CfA</span><br><a href="mailto:jguillochon@cfa.harvard.edu" target="_blank">jguillochon@cfa.harvard.edu</a></div></div>
</font></span></div></div>
</blockquote></div><br><br clear="all"><div><br></div></div></div><span><font color="#888888">-- <br><div>Anshu Dubey <<a href="mailto:adubey@lbl.gov" target="_blank">adubey@lbl.gov</a>><br>Mailstop: 50A1148<br>Building-Room 050A-2154 <br><a href="tel:510-486-5242" value="+15104865242" target="_blank">510-486-5242</a> </div>
</font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div><div dir="ltr">James Guillochon<br><span style="font-family:arial;font-size:small">Einstein Fellow at the Harvard-Smithsonian CfA</span><br><a href="mailto:jguillochon@cfa.harvard.edu" target="_blank">jguillochon@cfa.harvard.edu</a><br></div></div>
</div>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div><div dir="ltr">James Guillochon<br><span style="font-family:arial;font-size:small">Einstein Fellow at the Harvard-Smithsonian CfA</span><br><a href="mailto:jguillochon@cfa.harvard.edu" target="_blank">jguillochon@cfa.harvard.edu</a></div></div>
</div>
</div></blockquote></div></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr">James Guillochon<br><span style="font-family:arial;font-size:small">Einstein Fellow at the Harvard-Smithsonian CfA</span><br><a href="mailto:jguillochon@cfa.harvard.edu" target="_blank">jguillochon@cfa.harvard.edu</a></div></div>

</div></div>