Exemple #1
0
def test_ecl_run():
    """Test import an eclrun with dates and export to roff after a diff"""

    dates = [19991201, 20030101]
    rprops = ["PRESSURE", "SWAT"]

    gg = Grid(REEKROOT,
              fformat="eclipserun",
              restartdates=dates,
              restartprops=rprops)

    # get the property object:
    pres1 = gg.get_prop_by_name("PRESSURE_20030101")
    tsetup.assert_almostequal(pres1.values.mean(), 308.45, 0.001)

    pres1.to_file(os.path.join(TMPDIR, "pres1.roff"))

    pres2 = gg.get_prop_by_name("PRESSURE_19991201")

    if isinstance(pres2, GridProperty):
        pass

    logger.debug(pres1.values)
    logger.debug(pres2.values)

    pres1.values = pres1.values - pres2.values
    # logger.debug(pres1.values)
    # logger.debug(pres1)
    avg = pres1.values.mean()
    # ok checked in RMS:
    tsetup.assert_almostequal(avg, -26.073, 0.001)

    pres1.to_file(os.path.join(TMPDIR, "pressurediff.roff"),
                  name="PRESSUREDIFF")