def test_import_restart(): """Import Restart""" g = Grid() g.from_file(GFILE1, fformat="egrid") x = GridProperties() names = ["PRESSURE", "SWAT"] dates = [19991201, 20010101] x.from_file(RFILE1, fformat="unrst", names=names, dates=dates, grid=g) # get the object pr = x.get_prop_by_name("PRESSURE_19991201") swat = x.get_prop_by_name("SWAT_19991201") logger.info(x.names) logger.info(swat.values3d.mean()) logger.info(pr.values3d.mean()) txt = "Average PRESSURE_19991201" assert pr.values.mean() == pytest.approx(334.52327, abs=0.0001), txt txt = "Average SWAT_19991201" assert swat.values.mean() == pytest.approx(0.87, abs=0.01), txt pr = x.get_prop_by_name("PRESSURE_20010101") logger.info(pr.values3d.mean()) txt = "Average PRESSURE_20010101" assert pr.values.mean() == pytest.approx(304.897, abs=0.01), txt
def test_import_init(): """Import INIT Reek""" g = Grid() g.from_file(GFILE1, fformat="egrid") x = GridProperties() names = ["PORO", "PORV"] x.from_file(IFILE1, fformat="init", names=names, grid=g) # get the object poro = x.get_prop_by_name("PORO") logger.info("PORO avg {}".format(poro.values.mean())) porv = x.get_prop_by_name("PORV") logger.info("PORV avg {}".format(porv.values.mean())) assert poro.values.mean() == pytest.approx(0.1677402, abs=0.00001)
def test_import_restart_gull(): """Import Restart Reek""" g = Grid() g.from_file(GFILE1, fformat="egrid") x = GridProperties() names = ["PRESSURE", "SWAT"] dates = [19991201] x.from_file(RFILE1, fformat="unrst", names=names, dates=dates, grid=g) # get the object pr = x.get_prop_by_name("PRESSURE_19991201") swat = x.get_prop_by_name("SWAT_19991201") logger.info(x.names) logger.info(swat.values3d.mean()) logger.info(pr.values3d.mean())
def test_import_soil(): """SOIL need to be computed in code from SWAT and SGAS""" g = Grid() g.from_file(GFILE1, fformat="egrid") x = GridProperties() names = ["SOIL", "SWAT", "PRESSURE"] dates = [19991201] x.from_file(RFILE1, fformat="unrst", names=names, dates=dates, grid=g) logger.info(x.names) # get the object instance soil = x.get_prop_by_name("SOIL_19991201") logger.info(soil.values3d.mean()) logger.debug(x.names) txt = "Average SOIL_19850101" assert soil.values.mean() == pytest.approx(0.121977, abs=0.001), txt