[yt-users] getting the average density inside a data container
Britton Smith
brittonsmith at gmail.com
Sun Dec 1 04:51:07 PST 2013
Hi Kearn,
Have a look at this section in the documentation:
http://yt-project.org/docs/2.6/analyzing/objects.html#processing-objects-derived-quantities
That should tell you what you need to know.
Britton
On Sat, Nov 30, 2013 at 11:49 AM, <k.grisdale at surrey.ac.uk> wrote:
> Hi Britton,
>
> Thanks for the reply. Your advice with the region command helped and I
> now have that working.
>
> I have tried to implement the commands you suggest to find the average
> density in side my object, but I don’t really understand how to apply them.
> Could you elaborate on them?
>
> Thanks
>
> Kearn
>
> On 30 Nov 2013, at 09:23, Britton Smith <brittonsmith at gmail.com> wrote:
>
> Hi Kearn,
>
> I believe the issue with the region is that the arguments are center,
> left corner, right corner. Each of those needs to be a triple.
>
> On the other issue, is AveragedDensity a field you created? If not, I
> don't believe it exists. If you want to take the average inside some
> object, you should do something like:
> object.quantities["WeightedAverageQuantity"](field, weight_field)
>
> Britton
>
>
> On Sat, Nov 30, 2013 at 7:12 AM, <k.grisdale at surrey.ac.uk> wrote:
>
>> Hi yt Users
>>
>> I am trying to create a data container that is a cube from which I can
>> extract the average density of the gas inside the cube. I am using Ramses
>> data.
>>
>> When I create a data container using:
>>
>> p6 = load("output_00006/info_00006.txt", fields =
>> ["Density","x-velocity", "y-velocity", "z-velocity", "Pressure"])
>> cen = [0, 0, 0]
>> box = p6.h.region(cen, 0.2/p6['cm'], 0.3/p6['cm’])
>> print box["Density”]
>>
>> I get this error:
>>
>> TypeError: 'float' object has no attribute ‘__getitem__'
>>
>> I have tried creating a sphere as a data container and this works for
>> the Density, but then when I try to take the AverageDensity in the sphere:
>>
>> sphere = p6.h.sphere(cen, 0.7/p6['cm’])
>> print sphere['AveragedDensity’]
>>
>> I get this error:
>>
>> enerationInProgress Traceback (most recent call
>> last)
>> <ipython-input-28-574d78c4fb5c> in <module>()
>> ----> 1 print sphere['AveragedDensity']
>>
>> /Applications/Code/yt/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc
>> in __getitem__(self, key)
>> 187 return self.field_data[f]
>> 188 else:
>> --> 189 self.get_data(f)
>> 190 # Note that this is less succinct so that we can account
>> for the case
>> 191 # when there are, for example, no elements in the object.
>>
>> /Applications/Code/yt/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc
>> in get_data(self, fields)
>> 514 return
>> 515 elif self._locked == True:
>> --> 516 raise GenerationInProgress(fields)
>> 517 # At this point, we want to figure out *all* our
>> dependencies.
>> 518 fields_to_get = self._identify_dependencies(fields_to_get)
>>
>> Any idea what I am doing wrong or if there is another way to go about
>> this?
>>
>> Thanks
>>
>> Kearn
>>
>> _______________________________________________
>> yt-users mailing list
>> yt-users at lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>
>>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20131201/00f9d9d1/attachment.htm>
More information about the yt-users
mailing list