From e33600e6b74c2e8b87cd82e58ab21039be625f4e Mon Sep 17 00:00:00 2001
From: Kai Bao <kai.bao@sintef.no>
Date: Sat, 11 Dec 2021 00:13:18 +0100
Subject: [PATCH] adding compressFactor to FluidStateCompositionModules

---
 .../fluidstates/FluidStateCompositionModules.hpp     | 12 +++++++++++-
 1 file changed, 11 insertions(+), 1 deletion(-)

diff --git a/opm/material/fluidstates/FluidStateCompositionModules.hpp b/opm/material/fluidstates/FluidStateCompositionModules.hpp
index d8447548c..c4ca00d4a 100644
--- a/opm/material/fluidstates/FluidStateCompositionModules.hpp
+++ b/opm/material/fluidstates/FluidStateCompositionModules.hpp
@@ -125,6 +125,15 @@ public:
         }
     }
 
+    void setCompressFactor(unsigned phaseIdx, const Scalar& value) {
+        Valgrind::CheckDefined(value);
+        Z_[phaseIdx] = value;
+    }
+
+    Scalar compressFactor(unsigned phaseIdx) const {
+        return Z_[phaseIdx];
+    }
+
     /*!
      * \brief Retrieve all parameters from an arbitrary fluid
      *        state.
@@ -209,7 +218,7 @@ public:
         const auto& p_crit = FluidSystem::criticalPressure(compIdx);
         const auto& p = asImp_().pressure(0); //for now assume no capillary pressure
 
-        const auto& tmp = exp(5.3727 * (1+acf) * (1-T_crit/T)) * (p_crit/p);
+        const auto& tmp = exp(5.37 * (1+acf) * (1-T_crit/T)) * (p_crit/p);
         return tmp;
     }
 
@@ -220,6 +229,7 @@ protected:
     std::array<std::array<Scalar,numComponents>,numPhases> moleFraction_;
     std::array<Scalar,numPhases> averageMolarMass_;
     std::array<Scalar,numPhases> sumMoleFractions_;
+    std::array<Scalar,numPhases> Z_;
     std::array<Scalar,numComponents> K_;
     Scalar L_;
     bool twophaseflag_;