def test_crop_grid_after_copy():
"""Copy a grid, then crop and check number of active cells."""
logger.info("Read grid...")
grd = Grid(EMEGFILE2)
grd.describe()
zprop = GridProperty(EMEZFILE2, name="Zone", grid=grd)
grd.describe(details=True)
logger.info(grd.dimensions)
grd2 = grd.copy()
grd2.describe(details=True)
logger.info("GRD2 props: %s", grd2.props)
assert grd.propnames == grd2.propnames
logger.info("GRD2 number of active cells: %s", grd2.nactive)
act = grd.get_actnum()
logger.info(act.values.shape)
logger.info("ZPROP: %s", zprop.values.shape)
grd2.crop((1, 30), (40, 80), (23, 46))
grd2.describe(details=True)
def test_import_grdecl_and_bgrdecl():
"""Eclipse import of GRDECL and binary GRDECL"""
grd1 = Grid(REEKFIL2, fformat="grdecl")
grd1.describe()
assert grd1.dimensions == (40, 64, 14)
assert grd1.nactive == 35812
# get dZ...
dzv1 = grd1.get_dz()
grd2 = Grid(REEKFIL3, fformat="bgrdecl")
grd2.describe()
assert grd2.dimensions == (40, 64, 14)
assert grd2.nactive == 35812
# get dZ...
dzv2 = grd2.get_dz()
tsetup.assert_almostequal(dzv1.values.mean(), dzv2.values.mean(), 0.001)
