Hi Dave,<br><br>The best I have been able to do for this is to manually bin up a 3d data object with the cut_region function and then do the calculation myself bin by bin.<br>For example, if you have a sphere, you can do<br>
new_sphere = sphere.cut_region(['grid["Temperature"] > 1e5', 'grid["Temperature"] < 1e6'])<br>That'll give you access to all the cells within that bin for a variance calculation.  Note the quotes around those expressions.  The cuts are applied in yt using eval functions.<br>
<br>This will work, but it's slow and requires you to do a lot by hand.<br><br>Britton<br><br><div class="gmail_quote">On Wed, Oct 5, 2011 at 11:51 AM, david collins <span dir="ltr"><<a href="mailto:antpuncher@gmail.com">antpuncher@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Hi, Everybody!<br>
<br>
Does anyone out there have a technique for getting the variance out of<br>
a profile object?  A profile object is good at getting <X> vs. B, I'd<br>
then like to get < (X - <X>)^2 > vs B.  Matt and I had spittballed the<br>
possibility some time ago, but I was wondering if anyone out there had<br>
successfully done it.<br>
<br>
Thanks,<br>
d.<br>
<font color="#888888"><br>
--<br>
Sent from my computer.<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>
</font></blockquote></div><br>