<div dir="ltr">Hmm.  This may have something to do with a lagrangian code too.  Perhaps someone who has more experience with SPH in yt should weigh in on this one.  Sorry I don't have the answer!<div><br></div><div>Cameron</div>

</div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Jun 25, 2014 at 5:52 AM, Desika Narayanan <span dir="ltr"><<a href="mailto:dnarayan@haverford.edu" target="_blank">dnarayan@haverford.edu</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">Hey Cameron,<div><br></div><div>Thanks for checking in!  </div><div><br></div><div>I had originally wondered about this too, and you're right that the snaps considered in the notebook are slightly different z (z=1.3 vs 1.5).   You're also right that reducing the size of the filtered bbox for the second snap (the one that throws currently) fixes the problem.</div>




<div><br></div><div>The part of this train of thought that stumped me is that all of the bounds of the filtered bbox are within the edges of the first domain that I construct -   i.e. when I first load the snapshot with some huge boxsize (1e5 kpc), then even with the re-centering the new smaller bbox that I'm using for filtering lies within the ds1.domain_left_edge, right_edge etc.   It could be that I'm understanding how the bbox filtering works incorrectly though.</div>




<div><br></div><div>thanks</div><span class="HOEnZb"><font color="#888888"><div>-d</div></font></span><div><div class="h5"><div><br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Jun 24, 2014 at 8:11 PM, Cameron Hummels <span dir="ltr"><<a href="mailto:chummels@gmail.com" target="_blank">chummels@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">Desika,<div><br></div><div>Is it possible that the bounds that you're providing to yt are in absolute units (not in comoving), and that in the different dataset, those bounds are no longer valid (e.g. if it were an early dataset at higher z so the volume is smaller in absolute coordinates)?  Just an idea.</div>






<div><br></div><div>Also, when you get this working, it might be a good notebook to include in the cookbook of the documentation for how to explore an SPH dataset.</div><div><br></div><div>Cameron</div></div><div class="gmail_extra">






<br><br><div class="gmail_quote"><div><div>On Mon, Jun 23, 2014 at 7:57 AM, Desika Narayanan <span dir="ltr"><<a href="mailto:dnarayan@haverford.edu" target="_blank">dnarayan@haverford.edu</a>></span> wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div>

<div dir="ltr">Hi All,<div><br></div><div>I'm playing around with filtering out particles from an SPH data set using the bounding box.  I have a (somewhat lengthy but not complicated) notebook that highlights some of the capabilities (and an issue that I've found) here:</div>








<div><br></div><div><a href="http://nbviewer.ipython.org/gist/dnarayanan/517b736c9bb0cd72a6fa" target="_blank">http://nbviewer.ipython.org/gist/dnarayanan/517b736c9bb0cd72a6fa</a><br></div><div><br></div><div>The notebook does 3 things:</div>








<div><br></div><div>1. It first loads a snapshot from a cosmological simulation, and does a projection plot.</div><div><br></div><div>2. It then sets a smaller bounding box, and attempts to zoom in on that region (centered on the max density).   This works great.</div>








<div><br></div><div>3. Perplexingly, though, when I try basically the exact same code (the section after the second image) but with a different snapshot in the same simulation, I get a particle out of bounds error.</div>







<div>
<br></div><div>One relevant clue is that this error only comes when trying to access smoothed quantities (i.e. the snippets where I try to access [as an example] the SFR gives no problems).    Probably obvious, but this is the bleeding edge yt (most recent update).</div>








<div><br></div><div>Any thoughts would definitely be welcome.</div><div><br></div><div>thanks,</div><div>desika</div></div>
<br></div></div>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><span><font color="#888888"><br><br clear="all"><div><br></div>-- <br>Cameron Hummels<div>Postdoctoral Researcher</div><div>Steward Observatory</div><div>University of Arizona</div><div>
<a href="http://chummels.org" target="_blank">http://chummels.org</a></div>


</font></span></div>
<br>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div></div></div></div>
<br>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>Cameron Hummels<div>Postdoctoral Researcher</div><div>Steward Observatory</div><div>University of Arizona</div><div><a href="http://chummels.org" target="_blank">http://chummels.org</a></div>


</div>