Stephen,<br><br>Your suspicion was correct.  For some reason, this simulation has a number of particles with ids that are -1.  I went all the way back to the very first output of that simulation and found 126 out of 256^3 particles to have ids of -1.  I have absolutely no idea why this is.  This simulation data is more than a year old, so it is possible there was a bug in enzo at the time.  I just tested ParallelHOP out on a much newer dataset and it worked without a problem.<br>
<br>Thanks a lot for your help and thanks again for writing this code.<br><br>Britton<br><br><div class="gmail_quote">On Tue, Dec 1, 2009 at 12:22 PM, Stephen Skory <span dir="ltr"><<a href="mailto:stephenskory@yahoo.com">stephenskory@yahoo.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Britton,<br>
<div class="im"><br>
>>Ok, I got one of these:<br>
><br>
>Real index -1!!!<br>
><br>
>I didn't get a "Padded index -1!!!".  You're suspicion seems to have been right.<br>
<br>
</div>So it seems to me that you probably have a particle with index == -1. There are a few things we can try. If you have just one with -1, we can change the default 'empty' value in those arrays to something else. I suspect you won't have just one, however, at which point you're kind of up a creek, because that's absolutely required for parallel HOP.<br>

<br>
It may be worth your time to fix the -1 IDs to something unique in the enzo dataset.<br>
<br>
I should add a check that screens for bad particle index values, like -1. I've learned something today.<br>
<br>
Let me know if there's anything else I can help you with!<br>
<div><div></div><div class="h5"><br>
 _______________________________________________________<br>
<a href="mailto:sskory@physics.ucsd.edu">sskory@physics.ucsd.edu</a>           o__  Stephen Skory<br>
<a href="http://physics.ucsd.edu/%7Esskory/" target="_blank">http://physics.ucsd.edu/~sskory/</a> _.>/ _Graduate Student<br>
________________________________(_)_\(_)_______________<br>
<br>
_______________________________________________<br>
Yt-dev mailing list<br>
<a href="mailto:Yt-dev@lists.spacepope.org">Yt-dev@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org</a><br>
</div></div></blockquote></div><br>