[yt-users] cannot SlicePlot() using cut_region() filter on data from yt.load_uniform_grid()
Dale Braden
genghis17 at comcast.net
Sun Dec 10 21:43:06 PST 2017
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
More information about the yt-users
mailing list