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