[yt-users] cannot SlicePlot() using cut_region() filter on data from yt.load_uniform_grid()

[Dale Braden]

Hi Nathan,

Here is a script whose output will show the same problem.  I'm using 
version 3.4.0

--------------------
import yt
import numpy as np

arr = np.random.random(size=(41,56,47))
data = dict(density = (arr,"g/cm**3"))
bbox = np.array([[-4.0,4.2],[-5.5,5.7],[-4.175872,5.22413]])
ds = yt.load_uniform_grid(data, arr.shape, length_unit="cm", bbox=bbox)

ad = ds.all_data()
print ad.quantities.extrema([("gas", "density")])

den = ad.cut_region(['obj["density"] > 0.0'])
print den.quantities.extrema([("gas", "density")])

p = yt.SlicePlot(ds, 'x','density', center="m", data_source=den)
----------------------

thanks for investigating,
Dale