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)