[Yt-dev] Clarification of Enzo Parameters, TopGridRank and TopGridDimensions

[Chris Malone]

Hi,

In attempting to develop yt support for an additional computer code,
MAESTRO, I'm finding a little confusion in relating the Enzo parameters to
some parameters used in MAESTRO.  I have checked
http://lca.ucsd.edu/software/enzo/v1.0.1/amr_guide/parameters-sorted.htmland
compared, for example, TopGridRank and TopGridDimensions to what is
used
in the yt code.  My understanding is as follows:

(1) TopGridRank = < int >
     essentially the spatial dimensionality of the problem;  As far as I can
tell, yt supports only 3-d data sets, so this /should/ be set to 3 for
current use?
(2) TopGridDimensions = TopGridRank-vector of < int >'s
     essentially the number of zones (not counting ghosts/guards) in each
direction of the problem;  Again, as far as I can tell this should be
something like:
          TopGridDimensions = number_zones_x number_zones_y number_zones_z


Is this correct?

If it is, then I'm confused about how the Orion frontend loads these
parameters.  For that frontend, there is a dictionary which translates
between the Orion (i.e. BoxLib) parameters and those of Enzo:
orion2enzoDict.  It is defined in frontends/orion/definitions.py as:

orion2enzoDict = {"amr.n_cell": "TopGridRank",
                  "materials.gamma": "Gamma",
                  "amr.ref_ratio": "RefineBy"
                  }

I don't have a sample Orion parameter file, but I do have one for Castro,
which also uses the BoxLib data format.  For Castro, the amr.n_cell line of
a parameter file looks like:

 amr.n_cell           = 32   32   32

If I make the assumption (and maybe this is where I'm going wrong...) that
the grid parameters used in a Castro parameter file are the same as those in
an Orion parameter file, then the OrionStaticOutput._parse_parameter_file
method is mapping TopGridRank to a list of < int >'s, which is what I
thought TopGridDimensions was given (2) above.  In other words, if (1) and
(2) are correct, shouldn't

TopGridRank = len(amr.n_cell.split())
TopGridDimensions = map(int, amr.n_cell.split())

?

Any clarification would be helpful!

Chris
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-dev-spacepope.org/attachments/20101219/0f8150bb/attachment.htm>