def test_banal7(xtgshow):
"""Create a simple property in a small grid box"""
grd = Grid(BANAL7)
assert grd.dimensions == (4, 2, 3)
disc = GridProperty(BANAL7, name="DISC")
assert disc.dimensions == (4, 2, 3)
assert disc.values.mean() == pytest.approx(0.59091, abs=0.001)
gprops = grd.get_gridquality_properties()
mix = gprops.get_prop_by_name("minangle_sides")
assert mix.values.mean() == pytest.approx(81.31036, abs=0.001)
if xtgshow:
lay = 2
layslice = xtgeo.plot.Grid3DSlice()
layslice.canvas(title=f"Layer {lay}")
layslice.plot_gridslice(
grd,
prop=mix,
mode="layer",
index=lay,
window=None,
linecolor="black",
)
layslice.show()
def test_gridquality_properties(show_plot):
"""Get grid quality props."""
grd1 = Grid(GRIDQC1)
props1 = grd1.get_gridquality_properties()
minang = props1.get_prop_by_name("minangle_topbase")
assert minang.values[5, 2, 1] == pytest.approx(71.05561, abs=0.001)
if show_plot:
lay = 1
layslice = xtgeo.plot.Grid3DSlice()
layslice.canvas(title=f"Layer {lay}")
layslice.plot_gridslice(
grd1,
prop=minang,
mode="layer",
index=lay + 1,
window=None,
linecolor="black",
)
layslice.show()
grd2 = Grid(GRIDQC2)
props2 = grd2.get_gridquality_properties()
neg = props2.get_prop_by_name("negative_thickness")
assert neg.values[0, 0, 0] == 0
assert neg.values[2, 1, 0] == 1
grd3 = Grid(EMEGFILE)
props3 = grd3.get_gridquality_properties()
concp = props3.get_prop_by_name("concave_proj")
if show_plot:
lay = 23
layslice = xtgeo.plot.Grid3DSlice()
layslice.canvas(title=f"Layer {lay}")
layslice.plot_gridslice(
grd3,
prop=concp,
mode="layer",
index=lay + 1,
window=None,
linecolor="black",
)
layslice.show()
def test_gridquality_properties():
"""Get grid quality props."""
grd1 = Grid(GRIDQC1)
props1 = grd1.get_gridquality_properties()
minang = props1.get_prop_by_name("minangle_topbase")
assert minang.values[5, 2, 1] == pytest.approx(71.05561, abs=0.001)
grd2 = Grid(GRIDQC2)
props2 = grd2.get_gridquality_properties()
neg = props2.get_prop_by_name("negative_thickness")
assert neg.values[0, 0, 0] == 0
assert neg.values[2, 1, 0] == 1
grd3 = Grid(EMEGFILE)
props3 = grd3.get_gridquality_properties()
concp = props3.get_prop_by_name("concave_proj")
assert concp.values.sum() == 7949
