<div dir="ltr">Sorry for all the questions, but this is very puzzling. Can you share the exact script you're using? If the dataset is small enough and it's ok to share, it would also help to see the dataset as well. We've used dropbox or google drive to share these sorts of files in the past with a lot of success.</div>
<div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Jul 10, 2014 at 3:48 AM, 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">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="HOEnZb"><div class="h5"><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><font color="#888888">
<div><br></div><div>-d</div></font></span></div><div><div><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" 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>