Is it possible to get grid size from an enzo AMR data?
Dear all, Is it possible to take grid size from an enzo AMR data? Thanks -- Reju Sam John
Hi Reju, If what you want is to get the size of each AMR grid, then for a given dataset, ds, you can access the list of all AMR grids as ds.index.grids. Then, for example, you can get the shape of the first grid with ds.index.grids[0].shape. Britton On Wed, Nov 26, 2014 at 10:45 AM, Reju Sam John <rejusamjohn@gmail.com> wrote:
Dear all, Is it possible to take grid size from an enzo AMR data? Thanks
-- Reju Sam John
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Thanks Britton. If I want to define a field like * (grid size)*x *, where 'x' is another user defined field, is it possible with the above method? Thank you On Wed, Nov 26, 2014 at 4:42 PM, Britton Smith <brittonsmith@gmail.com> wrote:
Hi Reju,
If what you want is to get the size of each AMR grid, then for a given dataset, ds, you can access the list of all AMR grids as ds.index.grids. Then, for example, you can get the shape of the first grid with ds.index.grids[0].shape.
Britton
On Wed, Nov 26, 2014 at 10:45 AM, Reju Sam John <rejusamjohn@gmail.com> wrote:
Dear all, Is it possible to take grid size from an enzo AMR data? Thanks
-- Reju Sam John
_______________________________________________ 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
-- Reju Sam John
Hi Reju, When you say (grid size), do you mean the number of cells in a single grid? If so, that may be a bit tricky. When you create a data container, a field isn't necessarily evaluated on a grid object, but on the cells contained within the container. Can you explain what it is you are trying to accomplish? There may be another way to do it. Britton On Wed, Nov 26, 2014 at 11:20 AM, Reju Sam John <rejusamjohn@gmail.com> wrote:
Thanks Britton. If I want to define a field like * (grid size)*x *, where 'x' is another user defined field, is it possible with the above method? Thank you
On Wed, Nov 26, 2014 at 4:42 PM, Britton Smith <brittonsmith@gmail.com> wrote:
Hi Reju,
If what you want is to get the size of each AMR grid, then for a given dataset, ds, you can access the list of all AMR grids as ds.index.grids. Then, for example, you can get the shape of the first grid with ds.index.grids[0].shape.
Britton
On Wed, Nov 26, 2014 at 10:45 AM, Reju Sam John <rejusamjohn@gmail.com> wrote:
Dear all, Is it possible to take grid size from an enzo AMR data? Thanks
-- Reju Sam John
_______________________________________________ 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
-- Reju Sam John
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Hi Britton, Actually I was confused with what a gird is. Now I understood that, want I actually need is the cell size. Thanks Britton for your kind information. What I am trying to accomplish? I am trying to define a new field which has the dimension of flux (say CRF_flux) . I think the following definition would work. @derived_field(name = "CRF_flux",take_log=False) def my_new_field(field, data): return data['dx']*data['dx']*data['CRFraction'] # along with some unit conversion constants. Where 'CRFraction' is one of my user defined field. I think data['dx'] is equivalent to the cell size, isn't it ? One more doubt, what will be the unit of the field 'dx'? Is it in code unit? Once again thanks for your suggestions and information. Thank You. On Fri, Nov 28, 2014 at 8:14 PM, Britton Smith <brittonsmith@gmail.com> wrote:
Hi Reju,
When you say (grid size), do you mean the number of cells in a single grid? If so, that may be a bit tricky. When you create a data container, a field isn't necessarily evaluated on a grid object, but on the cells contained within the container. Can you explain what it is you are trying to accomplish? There may be another way to do it.
Britton
On Wed, Nov 26, 2014 at 11:20 AM, Reju Sam John <rejusamjohn@gmail.com> wrote:
Thanks Britton. If I want to define a field like * (grid size)*x *, where 'x' is another user defined field, is it possible with the above method? Thank you
On Wed, Nov 26, 2014 at 4:42 PM, Britton Smith <brittonsmith@gmail.com> wrote:
Hi Reju,
If what you want is to get the size of each AMR grid, then for a given dataset, ds, you can access the list of all AMR grids as ds.index.grids. Then, for example, you can get the shape of the first grid with ds.index.grids[0].shape.
Britton
On Wed, Nov 26, 2014 at 10:45 AM, Reju Sam John <rejusamjohn@gmail.com> wrote:
Dear all, Is it possible to take grid size from an enzo AMR data? Thanks
-- Reju Sam John
_______________________________________________ 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
-- Reju Sam John
_______________________________________________ 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
-- Reju Sam John
Hi Reju, Yes, the "dx" field is the cell width, so that will do what you want. It will be in units of length and will likely default to code_length, but you can convert it to anything you like. Additionally, I recommend the following procedure for creating your field: def my_new_field(field, data): return data['dx']*data['dx']*data['CRFraction'] ds = yt.load(...) ds.add_field("CRF_flux", function=my_new_field, take_log=False, units=SOME_UNITS) Adding the field directly to your loaded dataset is now the preferred method. Notice, you can also specify the units that you want the field returned in. Britton On Fri, Nov 28, 2014 at 9:49 PM, Reju Sam John <rejusamjohn@gmail.com> wrote:
Hi Britton, Actually I was confused with what a gird is. Now I understood that, want I actually need is the cell size. Thanks Britton for your kind information. What I am trying to accomplish? I am trying to define a new field which has the dimension of flux (say CRF_flux) . I think the following definition would work.
@derived_field(name = "CRF_flux",take_log=False) def my_new_field(field, data): return data['dx']*data['dx']*data['CRFraction'] # along with some unit conversion constants.
Where 'CRFraction' is one of my user defined field. I think data['dx'] is equivalent to the cell size, isn't it ? One more doubt, what will be the unit of the field 'dx'? Is it in code unit? Once again thanks for your suggestions and information. Thank You.
On Fri, Nov 28, 2014 at 8:14 PM, Britton Smith <brittonsmith@gmail.com> wrote:
Hi Reju,
When you say (grid size), do you mean the number of cells in a single grid? If so, that may be a bit tricky. When you create a data container, a field isn't necessarily evaluated on a grid object, but on the cells contained within the container. Can you explain what it is you are trying to accomplish? There may be another way to do it.
Britton
On Wed, Nov 26, 2014 at 11:20 AM, Reju Sam John <rejusamjohn@gmail.com> wrote:
Thanks Britton. If I want to define a field like * (grid size)*x *, where 'x' is another user defined field, is it possible with the above method? Thank you
On Wed, Nov 26, 2014 at 4:42 PM, Britton Smith <brittonsmith@gmail.com> wrote:
Hi Reju,
If what you want is to get the size of each AMR grid, then for a given dataset, ds, you can access the list of all AMR grids as ds.index.grids. Then, for example, you can get the shape of the first grid with ds.index.grids[0].shape.
Britton
On Wed, Nov 26, 2014 at 10:45 AM, Reju Sam John <rejusamjohn@gmail.com> wrote:
Dear all, Is it possible to take grid size from an enzo AMR data? Thanks
-- Reju Sam John
_______________________________________________ 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
-- Reju Sam John
_______________________________________________ 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
-- Reju Sam John
Hello Britton, Thank you very much..... On Sat, Nov 29, 2014 at 2:51 PM, Britton Smith <brittonsmith@gmail.com> wrote:
Hi Reju,
Yes, the "dx" field is the cell width, so that will do what you want. It will be in units of length and will likely default to code_length, but you can convert it to anything you like.
Additionally, I recommend the following procedure for creating your field:
def my_new_field(field, data): return data['dx']*data['dx']*data['CRFraction']
ds = yt.load(...) ds.add_field("CRF_flux", function=my_new_field, take_log=False, units=SOME_UNITS)
Adding the field directly to your loaded dataset is now the preferred method. Notice, you can also specify the units that you want the field returned in.
Britton
On Fri, Nov 28, 2014 at 9:49 PM, Reju Sam John <rejusamjohn@gmail.com> wrote:
Hi Britton, Actually I was confused with what a gird is. Now I understood that, want I actually need is the cell size. Thanks Britton for your kind information. What I am trying to accomplish? I am trying to define a new field which has the dimension of flux (say CRF_flux) . I think the following definition would work.
@derived_field(name = "CRF_flux",take_log=False) def my_new_field(field, data): return data['dx']*data['dx']*data['CRFraction'] # along with some unit conversion constants.
Where 'CRFraction' is one of my user defined field. I think data['dx'] is equivalent to the cell size, isn't it ? One more doubt, what will be the unit of the field 'dx'? Is it in code unit? Once again thanks for your suggestions and information. Thank You.
On Fri, Nov 28, 2014 at 8:14 PM, Britton Smith <brittonsmith@gmail.com> wrote:
Hi Reju,
When you say (grid size), do you mean the number of cells in a single grid? If so, that may be a bit tricky. When you create a data container, a field isn't necessarily evaluated on a grid object, but on the cells contained within the container. Can you explain what it is you are trying to accomplish? There may be another way to do it.
Britton
On Wed, Nov 26, 2014 at 11:20 AM, Reju Sam John <rejusamjohn@gmail.com> wrote:
Thanks Britton. If I want to define a field like * (grid size)*x *, where 'x' is another user defined field, is it possible with the above method? Thank you
On Wed, Nov 26, 2014 at 4:42 PM, Britton Smith <brittonsmith@gmail.com> wrote:
Hi Reju,
If what you want is to get the size of each AMR grid, then for a given dataset, ds, you can access the list of all AMR grids as ds.index.grids. Then, for example, you can get the shape of the first grid with ds.index.grids[0].shape.
Britton
On Wed, Nov 26, 2014 at 10:45 AM, Reju Sam John <rejusamjohn@gmail.com
wrote:
Dear all, Is it possible to take grid size from an enzo AMR data? Thanks
-- Reju Sam John
_______________________________________________ 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
-- Reju Sam John
_______________________________________________ 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
-- Reju Sam John
-- Reju Sam John
participants (2)
-
Britton Smith -
Reju Sam John