From 92bc55385d80bb063f80f6683409132ccf38f597 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Atgeirr=20Fl=C3=B8=20Rasmussen?= <atgeirr@sintef.no>
Date: Thu, 14 Jun 2012 21:53:56 +0200
Subject: [PATCH] Update porosities before calling transport solver.

---
 examples/polymer_reorder.cpp     | 5 +++++
 opm/polymer/SimulatorPolymer.cpp | 5 +++++
 2 files changed, 10 insertions(+)

diff --git a/examples/polymer_reorder.cpp b/examples/polymer_reorder.cpp
index 105542b66..110fddedc 100644
--- a/examples/polymer_reorder.cpp
+++ b/examples/polymer_reorder.cpp
@@ -767,6 +767,11 @@ main(int argc, char** argv)
             }
         } while (!well_control_passed);
 
+        // Update pore volumes if rock is compressible.
+        if (rock_comp->isActive()) {
+            computePorevolume(*grid->c_grid(), props->porosity(), *rock_comp, state.pressure(), porevol);
+        }
+
         // Process transport sources (to include bdy terms and well flows).
         Opm::computeTransportSource(*grid->c_grid(), src, state.faceflux(), 1.0,
                                     wells->c_wells(), well_state.perfRates(), reorder_src);
diff --git a/opm/polymer/SimulatorPolymer.cpp b/opm/polymer/SimulatorPolymer.cpp
index e7289ce39..d48067b52 100644
--- a/opm/polymer/SimulatorPolymer.cpp
+++ b/opm/polymer/SimulatorPolymer.cpp
@@ -296,6 +296,11 @@ namespace Opm
                 ptime += pt;
             } while (false);
 
+            // Update pore volumes if rock is compressible.
+            if (rock_comp_props_ && rock_comp_props_->isActive()) {
+                computePorevolume(grid_, props_.porosity(), *rock_comp_props_, state.pressure(), porevol);
+            }
+
             // Process transport sources (to include bdy terms and well flows).
             Opm::computeTransportSource(grid_, src_, state.faceflux(), 1.0,
                                         wells_, well_state.perfRates(), transport_src);