def test_reverse_row_axis_dual(tmpdir):
"""Reverse axis for distorted but small grid"""
grd = Grid(DUAL)
assert grd.ijk_handedness == "left"
grd.to_file(join(tmpdir, "dual_left.grdecl"), fformat="grdecl")
cellcorners1 = grd.get_xyz_cell_corners((5, 1, 1))
grd.reverse_row_axis()
assert grd.ijk_handedness == "right"
grd.to_file(join(tmpdir, "dual_right.grdecl"), fformat="grdecl")
cellcorners2 = grd.get_xyz_cell_corners((5, 3, 1))
assert cellcorners1[7] == cellcorners2[1]
def test_reverse_row_axis_box(tmpdir):
"""Crop a grid."""
grd = Grid()
grd.create_box(
origin=(1000, 4000, 300),
increment=(100, 100, 2),
dimension=(2, 3, 1),
rotation=0,
)
assert grd.ijk_handedness == "left"
grd.to_file(join(tmpdir, "reverse_left.grdecl"), fformat="grdecl")
grd.reverse_row_axis()
assert grd.ijk_handedness == "right"
grd.to_file(join(tmpdir, "reverse_right.grdecl"), fformat="grdecl")
def test_reverse_row_axis_dualprops():
"""Reverse axis for distorted but small grid with props"""
grd = Grid(DUALPROPS, fformat="eclipserun", initprops=["PORO", "PORV"])
poro = grd.gridprops.props[0]
logger.info(grd.gridprops.describe())
porowas = poro.copy()
tsetup.assert_almostequal(poro.values[1, 0, 0], 0.17777, 0.01)
assert grd.ijk_handedness == "left"
grd.reverse_row_axis()
tsetup.assert_almostequal(poro.values[1, 2, 0], 0.17777, 0.01)
assert poro.values[1, 2, 0] == porowas.values[1, 0, 0]
grd.reverse_row_axis()
assert poro.values[1, 0, 0] == porowas.values[1, 0, 0]
assert grd.ijk_handedness == "left"
grd.reverse_row_axis(ijk_handedness="left") # ie do nothing in this case
assert poro.values[1, 0, 0] == porowas.values[1, 0, 0]
assert grd.ijk_handedness == "left"
