getting the average density inside a data container
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
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@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@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
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@gmail.com<mailto:brittonsmith@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@surrey.ac.uk<mailto:k.grisdale@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@lists.spacepope.org<mailto:yt-users@lists.spacepope.org> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org<mailto:yt-users@lists.spacepope.org> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Hi Kearn, Have a look at this section in the documentation: http://yt-project.org/docs/2.6/analyzing/objects.html#processing-objects-der... That should tell you what you need to know. Britton On Sat, Nov 30, 2013 at 11:49 AM, <k.grisdale@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@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@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@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Hi Britton, Just wanted to say thank you for the help. That link was exactly what I needed and have now got my code running. Thanks Kearn On 1 Dec 2013, at 12:51, Britton Smith <brittonsmith@gmail.com<mailto:brittonsmith@gmail.com>> wrote: Hi Kearn, Have a look at this section in the documentation: http://yt-project.org/docs/2.6/analyzing/objects.html#processing-objects-der... That should tell you what you need to know. Britton On Sat, Nov 30, 2013 at 11:49 AM, <k.grisdale@surrey.ac.uk<mailto:k.grisdale@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@gmail.com<mailto:brittonsmith@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@surrey.ac.uk<mailto:k.grisdale@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@lists.spacepope.org<mailto:yt-users@lists.spacepope.org> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org<mailto:yt-users@lists.spacepope.org> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org<mailto:yt-users@lists.spacepope.org> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org<mailto:yt-users@lists.spacepope.org> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Excellent! Glad to hear it! Britton On Tue, Dec 3, 2013 at 6:13 PM, <k.grisdale@surrey.ac.uk> wrote:
Hi Britton,
Just wanted to say thank you for the help. That link was exactly what I needed and have now got my code running.
Thanks
Kearn On 1 Dec 2013, at 12:51, Britton Smith <brittonsmith@gmail.com> wrote:
Hi Kearn,
Have a look at this section in the documentation:
http://yt-project.org/docs/2.6/analyzing/objects.html#processing-objects-der...
That should tell you what you need to know.
Britton
On Sat, Nov 30, 2013 at 11:49 AM, <k.grisdale@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@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@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@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
participants (2)
-
Britton Smith
-
k.grisdale@surrey.ac.uk