From 4f23ef92e3c8a476b19d3cee65670f7aa264e2f8 Mon Sep 17 00:00:00 2001 From: Andreas Lauser Date: Tue, 20 Oct 2015 13:34:53 +0200 Subject: [PATCH] adapt to the threshold pressure related API changes to be able to determine the threshold pressure from the initial condition it needs to be able to access the initial condition, i.e., the initial simulator state, material properties object and gravity constant need to be passed to the thresholdPressure() function. --- examples/flow.cpp | 4 +++- examples/flow_solvent.cpp | 4 +++- 2 files changed, 6 insertions(+), 2 deletions(-) diff --git a/examples/flow.cpp b/examples/flow.cpp index 091685aa0..5127ec1a9 100644 --- a/examples/flow.cpp +++ b/examples/flow.cpp @@ -367,7 +367,9 @@ try // initialize variables simtimer.init(timeMap); - std::vector threshold_pressures = thresholdPressures(parseMode, eclipseState, grid); + std::map, double> maxDp; + computeMaxDp(maxDp, deck, eclipseState, grid, state, props, gravity[2]); + std::vector threshold_pressures = thresholdPressures(deck, eclipseState, grid, maxDp); SimulatorFullyImplicitBlackoil< Grid > simulator(param, grid, diff --git a/examples/flow_solvent.cpp b/examples/flow_solvent.cpp index 528d7b8b6..7ee180b76 100644 --- a/examples/flow_solvent.cpp +++ b/examples/flow_solvent.cpp @@ -377,7 +377,9 @@ try // initialize variables simtimer.init(timeMap); - std::vector threshold_pressures = thresholdPressures(parseMode, eclipseState, grid); + std::map, double> maxDp; + computeMaxDp(maxDp, deck, eclipseState, grid, state, props, gravity[2]); + std::vector threshold_pressures = thresholdPressures(deck, eclipseState, grid, maxDp); SimulatorFullyImplicitBlackoilSolvent< Grid > simulator(param, grid,