I'm working with a dataset from a particle-based simulation with 64^3 particles and am trying to create a covering grid so that I can create a 3D visualization (using volume rendering). I have a problem with the dimensions of that covering grid. A ProjectionPlot from the covering grid dataset only produces a 2x2x2 pixel plot even though the original dataset ProjectionPlot is a 64x64x64 plot. I tried defining the dimensions of the covering grid to be [64,64,64], which results in a 64 pixel wide plot however data is only present in the first 2x2x2 pixels and the rest of the plot is blank. The way I define my covering grid is shown below. Could someone please suggest a solution to this problem? cov_grid = ds.covering_grid(level=0, left_edge=[-1,-1,-1], dims=ds.domain_dimensions, fields=('deposit','PartType1_density')) #alternate version cov_grid = ds.covering_grid(level=32, left_edge=[-1,-1,-1], dims=[64,64,64], fields=('deposit','PartType1_density')) ds_new = cov_grid.save_as_dataset() grid_ds = yt.load(ds_new)