From 4d800cad3d4d2c6c08ef2b47e631c0c9ff4b2433 Mon Sep 17 00:00:00 2001
From: Xavier Raynaud <xavier.raynaud@sintef.no>
Date: Wed, 29 Aug 2012 10:54:21 +0200
Subject: [PATCH] Added functionality to compute porosity for compressible
 fluid.

---
 opm/core/utility/miscUtilities.cpp | 19 +++++++++++++++++++
 opm/core/utility/miscUtilities.hpp | 11 +++++++++++
 2 files changed, 30 insertions(+)

diff --git a/opm/core/utility/miscUtilities.cpp b/opm/core/utility/miscUtilities.cpp
index 16c59b17..ac9552dc 100644
--- a/opm/core/utility/miscUtilities.cpp
+++ b/opm/core/utility/miscUtilities.cpp
@@ -70,6 +70,25 @@ namespace Opm
         }
     }
 
+    /// @brief Computes porosity of all cells in a grid, with rock compressibility effects.
+    /// @param[in]  grid               a grid
+    /// @param[in]  porosity_standard  array of grid.number_of_cells porosity values (at standard conditions)
+    /// @param[in]  rock_comp          rock compressibility properties
+    /// @param[in]  pressure           pressure by cell
+    /// @param[out] porosity           porosity (at reservoir condition)
+    void computePorosity(const UnstructuredGrid& grid,
+                         const double* porosity_standard,
+                         const RockCompressibility& rock_comp,
+                         const std::vector<double>& pressure,
+                         std::vector<double>& porosity)
+    {
+        int num_cells = grid.number_of_cells;
+        porosity.resize(num_cells);
+        for (int i = 0; i < num_cells; ++i) {
+            porosity[i] = porosity_standard[i]*rock_comp.poroMult(pressure[i]);
+        }
+    }
+
 
     /// @brief Computes total saturated volumes over all grid cells.
     /// @param[in]  pv        the pore volume by cell.
diff --git a/opm/core/utility/miscUtilities.hpp b/opm/core/utility/miscUtilities.hpp
index 37b252bc..de452fba 100644
--- a/opm/core/utility/miscUtilities.hpp
+++ b/opm/core/utility/miscUtilities.hpp
@@ -54,6 +54,17 @@ namespace Opm
                            const std::vector<double>& pressure,
                            std::vector<double>& porevol);
 
+    /// @brief Computes porosity of all cells in a grid, with rock compressibility effects.
+    /// @param[in]  grid               a grid
+    /// @param[in]  porosity_standard  array of grid.number_of_cells porosity values (at standard conditions)
+    /// @param[in]  rock_comp          rock compressibility properties
+    /// @param[in]  pressure           pressure by cell
+    /// @param[out] porosity           porosity (at reservoir condition)
+    void computePorosity(const UnstructuredGrid& grid,
+                         const double* porosity_standard,
+                         const RockCompressibility& rock_comp,
+                         const std::vector<double>& pressure,
+                         std::vector<double>& porosity);
 
     /// @brief Computes total saturated volumes over all grid cells.
     /// @param[in]  pv        the pore volume by cell.