Is it possible to make a slice or projection of a extracted, fixed resolution region? i.e. the resultant object of:
extractCube = pf.h.covering_grid(extract_level, left_edge=extractLE, right_edge=extractRE, # How many dimensions along each axis dims=extractDims, # And any fields to preload (this is optional!) fields=["PotentialField"], num_ghost_zones = 3)
Or of an array later made to the same size? e.g.
cloud = na.zeros(extractCube["x-velocity"].shape, dtype='float64')