[yt-users] 2D Axisymmetric BoxLib Dataset Issue

[Kacper Kowalik]

On 01/18/2016 11:59 AM, Max Katz wrote:
> Hi all,
>
> I'm trying to use yt for analysis of a 2D axisymmetric (R-z) dataset using
> the BoxLib frontend. After I load the dataset I am trying to make a
> covering grid with the following commands:
>
>      ds = yt.load(pltfile)
>
>      problo = ds.domain_left_edge.v
>      probhi = ds.domain_right_edge.v
>      dim    = ds.domain_dimensions
>
>      grid = ds.covering_grid(level=0, left_edge = problo, dims = dim)
>
> That leads to the following output at the terminal:
>
> yt : [INFO     ] 2016-01-18 01:35:47,516 Parameters:
> current_time              = 0.0
> yt : [INFO     ] 2016-01-18 01:35:47,516 Parameters:
> domain_dimensions         = [1024 2048]
> yt : [INFO     ] 2016-01-18 01:35:47,517 Parameters:
> domain_left_edge          = [  0.00000000e+00  -5.12000000e+09
> 0.00000000e+00]
> yt : [INFO     ] 2016-01-18 01:35:47,517 Parameters:
> domain_right_edge         = [  5.12000000e+09   5.12000000e+09
> 6.28318531e+00]
> Traceback (most recent call last):
>    File "analyze.py", line 78, in <module>
>      wdmerger.rho_T_scatterplot(eps_filename, output_dir + pltfile)
>    File "/home/mkatz/codes/wdmerger/analysis/wdmerger.py", line 615, in
> rho_T_scatterplot
>      grid = ds.covering_grid(level=0, left_edge = problo, dims = dim)
>    File
> "/gpfs/home/mkatz/codes/yt-x86_64/src/yt-hg/yt/data_objects/construction_data_containers.py",
> line 503, in __init__
>      self.ActiveDimensions = self._sanitize_dims(dims)
>    File
> "/gpfs/home/mkatz/codes/yt-x86_64/src/yt-hg/yt/data_objects/construction_data_containers.py",
> line 552, in _sanitize_dims
>      "Length of dims must match the dimensionality of the dataset")
> RuntimeError: Length of dims must match the dimensionality of the dataset
>
> Thoughts? It seems like the code has some inconsistency in what the
> dimensionality of the dataset is. The Castro plot file in question doesn't
> provide that third angular dimension of size 2*pi.

Could you try dims=[1024, 2048, 1] ?

Cheers,
Kacper

> Max Katz
> Ph.D. Candidate, Department of Physics and Astronomy
> Stony Brook University
> http://astro.sunysb.edu/mkatz/
>
>
>
> _______________________________________________
> 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