<meta http-equiv="Content-Type" content="text/html; charset=utf-8"><div dir="ltr">That results in the following:<br><br><div style="margin-left:40px">Traceback (most recent call last):<br> File "analyze.py", line 78, in <module><br> wdmerger.rho_T_scatterplot(eps_filename, output_dir + pltfile)<br> File "/home/mkatz/codes/wdmerger/analysis/wdmerger.py", line 617, in rho_T_scatterplot<br> grid = ds.covering_grid(level=0, left_edge = problo, dims = [1024, 2048, 1])<br> File "/gpfs/home/mkatz/codes/yt-x86_64/src/yt-hg/yt/data_objects/construction_data_containers.py", line 506, in __init__<br> rdx[np.where(np.array(dims) - 2 * num_ghost_zones <= 1)] = 1 # issue 602<br>IndexError: index 2 is out of bounds for axis 1 with size 2<br></div><br></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">Max Katz<br>Ph.D. Candidate, Department of Physics and Astronomy<br>Stony Brook University<br><a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__astro.sunysb.edu_mkatz_&d=BQMFaQ&c=8hUWFZcy2Z-Za5rBPlktOQ&r=hgcBC3x6dKFoTrmFmMYYbKNfiHZlGLKliIidd1LwmHI&m=cHTPBXK_eSljbf7el6Bgce3bRfuLZUq9zIwXN75WLms&s=yEEyu799_PkLvghZKW-tOQKE5rGgTrwKZQ3TFmrYnW0&e=" target="_blank">http://astro.sunysb.edu/mkatz/</a></div></div></div></div></div></div></div>
<br><div class="gmail_quote">On Mon, Jan 18, 2016 at 3:32 PM, Kacper Kowalik <span dir="ltr"><<a href="mailto:xarthisius.kk@gmail.com" target="_blank">xarthisius.kk@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="HOEnZb"><div class="h5">On 01/18/2016 11:59 AM, Max Katz wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi all,<br>
<br>
I'm trying to use yt for analysis of a 2D axisymmetric (R-z) dataset using<br>
the BoxLib frontend. After I load the dataset I am trying to make a<br>
covering grid with the following commands:<br>
<br>
ds = yt.load(pltfile)<br>
<br>
problo = ds.domain_left_edge.v<br>
probhi = ds.domain_right_edge.v<br>
dim = ds.domain_dimensions<br>
<br>
grid = ds.covering_grid(level=0, left_edge = problo, dims = dim)<br>
<br>
That leads to the following output at the terminal:<br>
<br>
yt : [INFO ] 2016-01-18 01:35:47,516 Parameters:<br>
current_time = 0.0<br>
yt : [INFO ] 2016-01-18 01:35:47,516 Parameters:<br>
domain_dimensions = [1024 2048]<br>
yt : [INFO ] 2016-01-18 01:35:47,517 Parameters:<br>
domain_left_edge = [ 0.00000000e+00 -5.12000000e+09<br>
0.00000000e+00]<br>
yt : [INFO ] 2016-01-18 01:35:47,517 Parameters:<br>
domain_right_edge = [ 5.12000000e+09 5.12000000e+09<br>
6.28318531e+00]<br>
Traceback (most recent call last):<br>
File "analyze.py", line 78, in <module><br>
wdmerger.rho_T_scatterplot(eps_filename, output_dir + pltfile)<br>
File "/home/mkatz/codes/wdmerger/analysis/wdmerger.py", line 615, in<br>
rho_T_scatterplot<br>
grid = ds.covering_grid(level=0, left_edge = problo, dims = dim)<br>
File<br>
"/gpfs/home/mkatz/codes/yt-x86_64/src/yt-hg/yt/data_objects/construction_data_containers.py",<br>
line 503, in __init__<br>
self.ActiveDimensions = self._sanitize_dims(dims)<br>
File<br>
"/gpfs/home/mkatz/codes/yt-x86_64/src/yt-hg/yt/data_objects/construction_data_containers.py",<br>
line 552, in _sanitize_dims<br>
"Length of dims must match the dimensionality of the dataset")<br>
RuntimeError: Length of dims must match the dimensionality of the dataset<br>
<br>
Thoughts? It seems like the code has some inconsistency in what the<br>
dimensionality of the dataset is. The Castro plot file in question doesn't<br>
provide that third angular dimension of size 2*pi.<br>
</blockquote>
<br></div></div>
Could you try dims=[1024, 2048, 1] ?<br>
<br>
Cheers,<br>
Kacper<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">
Max Katz<br>
Ph.D. Candidate, Department of Physics and Astronomy<br>
Stony Brook University<br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__astro.sunysb.edu_mkatz_&d=BQMFaQ&c=8hUWFZcy2Z-Za5rBPlktOQ&r=hgcBC3x6dKFoTrmFmMYYbKNfiHZlGLKliIidd1LwmHI&m=cHTPBXK_eSljbf7el6Bgce3bRfuLZUq9zIwXN75WLms&s=yEEyu799_PkLvghZKW-tOQKE5rGgTrwKZQ3TFmrYnW0&e=" rel="noreferrer" target="_blank">http://astro.sunysb.edu/mkatz/</a><br>
<br>
<br>
<br></span>
_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.spacepope.org_listinfo.cgi_yt-2Dusers-2Dspacepope.org&d=BQMFaQ&c=8hUWFZcy2Z-Za5rBPlktOQ&r=hgcBC3x6dKFoTrmFmMYYbKNfiHZlGLKliIidd1LwmHI&m=cHTPBXK_eSljbf7el6Bgce3bRfuLZUq9zIwXN75WLms&s=j__RAlG2YtU40I45jhVHQZHjWlqJJAFUDJ2CUTpY8NI&e=" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br>
</blockquote>
<br>
_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.spacepope.org_listinfo.cgi_yt-2Dusers-2Dspacepope.org&d=BQMFaQ&c=8hUWFZcy2Z-Za5rBPlktOQ&r=hgcBC3x6dKFoTrmFmMYYbKNfiHZlGLKliIidd1LwmHI&m=cHTPBXK_eSljbf7el6Bgce3bRfuLZUq9zIwXN75WLms&s=j__RAlG2YtU40I45jhVHQZHjWlqJJAFUDJ2CUTpY8NI&e=" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
</blockquote></div><br></div>