def test_hybridgrid2():
"""Making a hybridgrid for Emerald case in region"""
logger.info('Read grid...')
grd = Grid(EMEGFILE)
logger.info('Read grid... done, NLAY is {}'.format(grd.nlay))
reg = GridProperty()
reg.from_file(EMERFILE, name='REGION')
nhdiv = 40
grd.convert_to_hybrid(nhdiv=nhdiv, toplevel=1650, bottomlevel=1690,
region=reg, region_number=1)
grd.to_file('TMP/test_hybridgrid2.roff')
def test_hybridgrid2(tmpdir):
"""Making a hybridgrid for Emerald case in region"""
logger.info("Read grid...")
grd = Grid(EMEGFILE)
logger.info("Read grid... done, NLAY is {}".format(grd.nlay))
reg = GridProperty()
reg.from_file(EMERFILE, name="REGION")
nhdiv = 40
grd.convert_to_hybrid(
nhdiv=nhdiv, toplevel=1650, bottomlevel=1690, region=reg, region_number=1
)
grd.to_file(join(tmpdir, "test_hybridgrid2.roff"))
def test_hybridgrid1(tmpdir):
"""Making a hybridgrid for Emerald case"""
logger.info("Read grid...")
grd = Grid(EMEGFILE)
assert grd.subgrids is not None # initially, prior to subgrids
logger.info("Read grid... done, NZ is %s", grd.nlay)
grd.to_file(join(tmpdir, "test_hybridgrid1_asis.bgrdecl"), fformat="bgrdecl")
logger.info("Convert...")
nhdiv = 40
newnlay = grd.nlay * 2 + nhdiv
grd.convert_to_hybrid(nhdiv=nhdiv, toplevel=1700, bottomlevel=1740)
logger.info("Hybrid grid... done, NLAY is now %s", grd.nlay)
assert grd.nlay == newnlay, "New NLAY number"
assert grd.subgrids is None
dzv = grd.get_dz()
assert dzv.values3d.mean() == pytest.approx(1.395, abs=0.01)
logger.info("Export...")
grd.to_file(join(tmpdir, "test_hybridgrid1.roff"))
logger.info("Read grid 2...")
grd2 = Grid(join(tmpdir, "test_hybridgrid1_asis.bgrdecl"))
logger.info("Read grid... done, NLAY is %s", grd2.nlay)
logger.info("Convert...")
nhdiv = 40
newnz = grd2.nlay * 2 + nhdiv
grd2.convert_to_hybrid(nhdiv=nhdiv, toplevel=1700, bottomlevel=1740)
logger.info("Hybrid grid... done, NZ is now %s", grd2.nlay)
assert grd2.nlay == newnz, "New NLAY number"
dzv = grd2.get_dz()
assert dzv.values3d.mean() == pytest.approx(1.395, abs=0.01)
def test_inactivate_thin_cells():
"""Make hybridgrid for Emerald case in region, and inactive thin cells"""
logger.info('Read grid...')
grd = Grid(EMEGFILE)
logger.info('Read grid... done, NLAY is {}'.format(grd.nlay))
reg = GridProperty()
reg.from_file(EMERFILE, name='REGION')
nhdiv = 40
grd.convert_to_hybrid(nhdiv=nhdiv, toplevel=1650, bottomlevel=1690,
region=reg, region_number=1)
grd.inactivate_by_dz(0.001)
grd.to_file('TMP/test_hybridgrid2_inact_thin.roff')
def test_hybridgrid1():
"""Making a hybridgrid for Emerald case (ROFF and GRDECL"""
logger.info("Read grid...")
grd = Grid(EMEGFILE)
logger.info("Read grid... done, NZ is %s", grd.nlay)
grd.to_file(join(TDP, "test_hybridgrid1_asis.grdecl"), fformat="grdecl")
logger.info("Convert...")
nhdiv = 40
newnlay = grd.nlay * 2 + nhdiv
grd.convert_to_hybrid(nhdiv=nhdiv, toplevel=1700, bottomlevel=1740)
logger.info("Hybrid grid... done, NLAY is now %s", grd.nlay)
assert grd.nlay == newnlay, "New NLAY number"
dzv = grd.get_dz()
assert dzv.values3d.mean() == pytest.approx(1.395, abs=0.01)
logger.info("Export...")
grd.to_file("TMP/test_hybridgrid1.roff")
logger.info("Read grid 2...")
grd2 = Grid("TMP/test_hybridgrid1_asis.grdecl")
logger.info("Read grid... done, NLAY is %s", grd2.nlay)
logger.info("Convert...")
nhdiv = 40
newnz = grd2.nlay * 2 + nhdiv
grd2.convert_to_hybrid(nhdiv=nhdiv, toplevel=1700, bottomlevel=1740)
logger.info("Hybrid grid... done, NZ is now %s", grd2.nlay)
assert grd2.nlay == newnz, "New NLAY number"
dzv = grd2.get_dz()
assert dzv.values3d.mean() == pytest.approx(1.395, abs=0.01)
