<div>Stephen,</div><div><br></div>I'm printing the warning message when along with the <a href="http://halo.id">halo.id</a> as the exception is raised, so it appears that halo 5 in this case has a radius smaller than dx, but what it really should be checking is the maximum instead of the virial radius instead if I'm not mistaken.<div>
<br></div><div>Matt,</div><div><br></div><div>I just want the field attributes of where the halo is located, and I'm comparing the attributes to the region surrounding the halo (up to 3x the virial or maximum radius) to see if there's anything that is substantially different.  I want to see what's different in the field quantity that suppresses star formation in lower mass haloes.  In this analysis the detailed shape and size of the sphere doesn't really matter, I'm just using it as a data container.</div>
<div><br></div><div>Right now I'm trying to output Density averaged values of each halo but whenever I encounter a halo with radius smaller than a cell (triggering YTSphereTooSmall), I output 0 for all attributes.  So I expect lots of zeros to occur on the small end.  However I'm seeing some zeros scattered even in the high mass end with large maximum radius, and I think the high mass halos with zeros is the result of it having a small virial radius, at least that's my current hypothesis.  </div>
<div><br></div><div>I'm going to try to go down the list from the high mass end to see if the second halo giving me zeros on the high end also have a virial radius smaller than dx.</div><div><br></div><div>I guess a work around for this problem temporarily is to use pf.h.sphere(halo.center_of_mass(), halo.maximum_radius()) instead of calling halo.get_sphere().</div>
<div><br></div><div>From</div><div>G.S.<br><br><div class="gmail_quote">On Tue, May 22, 2012 at 6:16 PM, Stephen Skory <span dir="ltr"><<a href="mailto:s@skory.us" target="_blank">s@skory.us</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Geoffrey and Matt,<br>
<div class="im"><br>
> While I think Stephen might know why the difference between<br>
> maximum_radius / virial_radius is confusing the system here,<br>
<br>
</div>I don't quite know what's going on. halo.get_sphere() should give you<br>
a YT sphere centered on the center of mass of the halo, with radius of<br>
.maximum_radius(). Perhaps it's not halo 5 that's giving you the<br>
error? Or there's an error in what you're printing?<br>
<div class="im"><br>
> I want to<br>
> suggest you back up a bit.  Do you really want to analyze halos that<br>
> are O(dx) across?<br>
<br>
</div>I agree with Matt that you should take a close look at what you are<br>
asking from these halos, and if one or just a few cells are sufficient<br>
to give you that.<br>
<span class="HOEnZb"><font color="#888888"><br>
--<br>
Stephen Skory<br>
<a href="mailto:s@skory.us">s@skory.us</a><br>
<a href="http://stephenskory.com/" target="_blank">http://stephenskory.com/</a><br>
<a href="tel:510.621.3687" value="+15106213687">510.621.3687</a> (google voice)<br>
</font></span><div class="HOEnZb"><div class="h5">_______________________________________________<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>
</div></div></blockquote></div><br></div>