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_refine_vertically(tmpdir): """Do a grid refinement vertically.""" logger.info("Read grid...") grd = Grid(EMEGFILE) logger.info("Read grid... done, NLAY is {}".format(grd.nlay)) logger.info("Subgrids before: %s", grd.get_subgrids()) avg_dz1 = grd.get_dz().values3d.mean() # idea; either a scalar (all cells), or a dictionary for zone wise grd.refine_vertically(3) avg_dz2 = grd.get_dz().values3d.mean() assert avg_dz1 == pytest.approx(3 * avg_dz2, abs=0.0001) logger.info("Subgrids after: %s", grd.get_subgrids()) grd.inactivate_by_dz(0.001) grd.to_file(join(tmpdir, "test_refined_by_3.roff"))