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<double> threshold_pressures = thresholdPressures(parseMode, eclipseState, grid);
+    std::map<std::pair<int, int>, double> maxDp;
+    computeMaxDp(maxDp, deck, eclipseState, grid, state, props, gravity[2]);
+    std::vector<double> 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<double> threshold_pressures = thresholdPressures(parseMode, eclipseState, grid);
+    std::map<std::pair<int, int>, double> maxDp;
+    computeMaxDp(maxDp, deck, eclipseState, grid, state, props, gravity[2]);
+    std::vector<double> threshold_pressures = thresholdPressures(deck, eclipseState, grid, maxDp);
 
     SimulatorFullyImplicitBlackoilSolvent< Grid >  simulator(param,
                                                       grid,