New issue 854: Different dimensions for coordinates and fields in `smoothed_covering_grid` https://bitbucket.org/yt_analysis/yt/issue/854/different-dimensions-for-coor... John ZuHone: This script: ``` #!python import yt import numpy as np pf = yt.load("GasSloshing/sloshing_nomag2_hdf5_plt_cnt_0150") L = 617. # width in kpc dx = 2.*pf.index.get_smallest_dx().in_units("kpc").v nx = int(L/dx) le = pf.arr([-0.5*nx*dx]*3, "kpc").in_units("code_length").v dims = np.array([nx]*3) cube = pf.smoothed_covering_grid(pf.index.max_level-1, le, dims) print cube["x"].shape print cube["density"].shape ``` Gives this output: ``` #!python (131, 131, 131) (132, 132, 132) ``` So the shapes of the coordinates and the fields have gotten out of sync somehow. This only happens for certain values of `L` and `nx`, setting `L = 600.` results in both dimensions being the same.