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