def test_geometrics_reek():
"""Import Reek and test geometrics."""
grd = Grid(REEKFILE, fformat="egrid")
geom = grd.get_geometrics(return_dict=True, cellcenter=False)
for key, val in geom.items():
logger.info("%s is %s", key, val)
# compared with RMS info:
tsetup.assert_almostequal(geom["xmin"], 456510.6, 0.1, "Xmin")
tsetup.assert_almostequal(geom["ymax"], 5938935.5, 0.1, "Ymax")
# cellcenter True:
geom = grd.get_geometrics(return_dict=True, cellcenter=True)
tsetup.assert_almostequal(geom["xmin"], 456620, 1, "Xmin cell center")
def test_reduce_to_one_layer():
"""Crop a grid."""
logger.info("Read grid...")
grd1 = Grid(EMEGFILE2)
geom1v1 = grd1.get_geometrics(allcells=True, cellcenter=True)
geom1v2 = grd1.get_geometrics(allcells=True, cellcenter=True, _ver=2)
grd2 = Grid(EMEGFILE2)
grd2.reduce_to_one_layer()
geom2v1 = grd2.get_geometrics(allcells=True, cellcenter=True)
geom2v2 = grd2.get_geometrics(allcells=True, cellcenter=True, _ver=2)
print("XXXX1V1", geom1v1)
print("XXXX1V2", geom1v2)
print("XXXX2V1", geom2v1)
print("XXXX2V2", geom2v2)
def test_eclgrid_import3(tmp_path):
"""Eclipse GRDECL import and translate."""
grd = Grid(BRILGRDECL, fformat="grdecl")
mylist = grd.get_geometrics()
xori1 = mylist[0]
# translate the coordinates
grd.translate_coordinates(translate=(100, 100, 10), flip=(1, 1, 1))
mylist = grd.get_geometrics()
xori2 = mylist[0]
# check if origin is translated 100m in X
tsetup.assert_equal(xori1 + 100, xori2, txt="Translate X distance")
grd.to_file(tmp_path / "g1_translate.roff", fformat="roff_binary")
grd.to_file(tmp_path / "g1_translate.bgrdecl", fformat="bgrdecl")
def test_reverse_row_axis_eme(tmpdir):
"""Reverse axis for emerald grid"""
grd1 = Grid(EMEGFILE)
assert grd1.ijk_handedness == "left"
grd1.to_file(join(tmpdir, "eme_left.roff"), fformat="roff")
grd2 = grd1.copy()
geom1 = grd1.get_geometrics(return_dict=True)
grd2.reverse_row_axis()
assert grd2.ijk_handedness == "right"
grd2.to_file(join(tmpdir, "eme_right.roff"), fformat="roff")
geom2 = grd2.get_geometrics(return_dict=True)
tsetup.assert_almostequal(geom1["avg_rotation"], geom2["avg_rotation"], 0.01)
