<meta http-equiv="Content-Type" content="text/html; charset=utf-8"><div dir="ltr">I'm not able to reproduce this with any of the 2D datasets I have on hand.<div><br></div><div>Would it be ok to share the dataset you're having trouble with or a test dataset that triggers the same error? You can use the curldrop Kacper set up to do this, which allows you to upload a file using:</div><div><br></div><div>$ curl -T <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__file.name&d=BQMFaQ&c=8hUWFZcy2Z-Za5rBPlktOQ&r=hgcBC3x6dKFoTrmFmMYYbKNfiHZlGLKliIidd1LwmHI&m=7et9dS598tPYZ-vRP5VZixdqUw7OqEIHoYfL3-rIqlo&s=aph5FFRIk3SDz9B7tMMK47Er4Ssm0JO5pMc6UJnyu9c&e=">file.name</a> <span class=""><a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__use.yt_upload_&d=BQMFaQ&c=8hUWFZcy2Z-Za5rBPlktOQ&r=hgcBC3x6dKFoTrmFmMYYbKNfiHZlGLKliIidd1LwmHI&m=7et9dS598tPYZ-vRP5VZixdqUw7OqEIHoYfL3-rIqlo&s=OoPRYvsMNHTTsztYkXUIDU-cVRV2mRQDlk5F24lvirs&e=">http://use.yt/upload/</a></span><br></div><div><p class="">You could also use e.g. dropbox or google drive.</p><p class="">-Nathan</p></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jan 18, 2016 at 3:40 PM, Max Katz <span dir="ltr"><<a href="mailto:maxpkatz@gmail.com" target="_blank">maxpkatz@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 dir="ltr">That results in the following:<br><br><div style="margin-left:40px"><span class="">Traceback (most recent call last):<br> File "analyze.py", line 78, in <module><br> wdmerger.rho_T_scatterplot(eps_filename, output_dir + pltfile)<br></span> 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"><span class=""><br clear="all"><div><div><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=7et9dS598tPYZ-vRP5VZixdqUw7OqEIHoYfL3-rIqlo&s=dlbHbQtrSA_hyr5jwVSZ7QHABfueggr9a4D9GiAroYE&e=" target="_blank">http://astro.sunysb.edu/mkatz/</a></div></div></div></div></div></div></div>
<br></span><div><div class="h5"><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><div>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>
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=7et9dS598tPYZ-vRP5VZixdqUw7OqEIHoYfL3-rIqlo&s=dlbHbQtrSA_hyr5jwVSZ7QHABfueggr9a4D9GiAroYE&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=7et9dS598tPYZ-vRP5VZixdqUw7OqEIHoYfL3-rIqlo&s=PDhEpLgBFcIg3IuDyaHJAKl8c6_wJFT_Fk8JuX3EJyk&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=7et9dS598tPYZ-vRP5VZixdqUw7OqEIHoYfL3-rIqlo&s=PDhEpLgBFcIg3IuDyaHJAKl8c6_wJFT_Fk8JuX3EJyk&e=" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
</blockquote></div><br></div></div></div>
<br>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">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=7et9dS598tPYZ-vRP5VZixdqUw7OqEIHoYfL3-rIqlo&s=PDhEpLgBFcIg3IuDyaHJAKl8c6_wJFT_Fk8JuX3EJyk&e=" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>