<div dir="ltr">Hi Desika,<div><br></div><div>If I read the log right the latest change is:</div><div><br></div><div><div>changeset: 12732:c9e843865b2a</div><div>branch: yt</div><div>tag: tip</div><div>parent: 12499:1b93a1840f38</div>
<div>user: dcollins4096</div><div>date: Sun Jul 06 20:08:39 2014 -0400</div></div><div><br></div><div>Thanks,</div><div>Daniel</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On 10 July 2014 10:08, 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">Hi Daniel,<div><br></div><div>Which changeset are you on? (Matt can correct me on this) but I think there was a time when the default was that over_refine_factor = 1. Though that wouldn't explain equal cell and particle numbers...</div>
<span class="HOEnZb"><font color="#888888">
<div><br></div><div>-d</div></font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Jul 9, 2014 at 3:29 PM, Daniel <span dir="ltr"><<a href="mailto:cunnama@gmail.com" target="_blank">cunnama@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">I was comparing density and temperature to the cell_volume.</div><div><div><div class="gmail_extra">
<br><br><div class="gmail_quote">On 9 July 2014 20:04, Nathan Goldbaum <span dir="ltr"><<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Which fields were you comparing with?<div><div><span></span><br><br>On Wednesday, July 9, 2014, Daniel <<a href="mailto:cunnama@gmail.com" target="_blank">cunnama@gmail.com</a>> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">Hi Matt,<div><br></div><div>No, just loading a standard GadgetDataset:</div><div><br></div><div><div>bbox = [[0.0, 1000000.0],</div><div> [0.0, 1000000.0],</div><div> [0.0, 1000000.0]]</div><div>
<br></div><div>unit_base = {'UnitLength_in_cm' : 3.08568e+21,</div><div> 'UnitMass_in_g' : 1.989e+43,</div><div> 'UnitVelocity_in_cm_per_s' : 100000}</div>
<div><br></div><div>pf = GadgetDataset("snapshot_050", unit_base = unit_base, bounding_box = bbox)<br></div></div><div><br></div><div>Cheers,</div><div>Daniel</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On 9 July 2014 17:58, Matthew Turk <span dir="ltr"><<a>matthewturk@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 Daniel,<br>
<div><br>
On Wed, Jul 9, 2014 at 10:57 AM, Daniel <<a>cunnama@gmail.com</a>> wrote:<br>
> Hi Nathan,<br>
><br>
> Thanks for getting back to me.<br>
><br>
> That is odd, because then if my region includes all the particles I have<br>
> equal elements in cell_volume and all the other gas quantities?<br>
><br>
<br>
</div>Yup, that is a bit odd. How are you loading the data? Are you<br>
specifying either n_ref or over_refine_factor when doing so?<br>
<br>
-Matt<br>
<div><div><br>
> Kind Regards,<br>
> Daniel<br>
><br>
><br>
> On 9 July 2014 15:06, Nathan Goldbaum <<a>nathan12343@gmail.com</a>> wrote:<br>
>><br>
>> Hi Daniel,<br>
>><br>
>> The cell_volume field is an alias for ('gas, 'cell_volume'), which for<br>
>> Gadget data corresponds to the volume of octree cells. In general there will<br>
>> be fewer octree cells than particles. I'm a little surprised the numbers<br>
>> are so close since the default option is to refine octree cells when more<br>
>> than 64 particles are in a cell, although the precise structure of the<br>
>> octree will depend on the particle distribution.<br>
>><br>
>> Nathan<br>
>><br>
>><br>
>> On Wednesday, July 9, 2014, Daniel <<a>cunnama@gmail.com</a>> wrote:<br>
>>><br>
>>> Hi there,<br>
>>><br>
>>> I am extracting a region from a Gadget dataset using the sphere function<br>
>>> and then attempting to access the "cell_volume" field.<br>
>>><br>
>>> For some reason the number of elements in the cell volume field is not<br>
>>> equal to the number of particles in the region.<br>
>>><br>
>>> For example I have 4341760 particles but only 4228930 cell volume<br>
>>> elements.<br>
>>><br>
>>> I have extracted a sphere which encompassed all the particles and the<br>
>>> number of cell volume elements is then equal.<br>
>>><br>
>>> Any ideas or advice would be greatly appreciated.<br>
>>><br>
>>> Kind Regards,<br>
>>> Daniel<br>
>><br>
>><br>
>> _______________________________________________<br>
>> yt-users mailing list<br>
>> <a>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>
><br>
><br>
> _______________________________________________<br>
> yt-users mailing list<br>
> <a>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>
_______________________________________________<br>
yt-users mailing list<br>
<a>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>
</blockquote>
</div></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><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><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></div>