From 52815742bbff9d839f8f78b106683368418ca40b Mon Sep 17 00:00:00 2001
From: Kristian Bendiksen <kristian.bendiksen@gmail.com>
Date: Wed, 24 Jun 2020 11:00:51 +0200
Subject: [PATCH] #6031 Initial porosity in pore compressibility calc must use
 geostatic time step data

---
 .../RigFemPartResultCalculatorPoreCompressibility.cpp         | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/ApplicationCode/GeoMech/GeoMechDataModel/RigFemPartResultCalculatorPoreCompressibility.cpp b/ApplicationCode/GeoMech/GeoMechDataModel/RigFemPartResultCalculatorPoreCompressibility.cpp
index 3b00f116c8..a3a998412d 100644
--- a/ApplicationCode/GeoMech/GeoMechDataModel/RigFemPartResultCalculatorPoreCompressibility.cpp
+++ b/ApplicationCode/GeoMech/GeoMechDataModel/RigFemPartResultCalculatorPoreCompressibility.cpp
@@ -152,7 +152,7 @@ RigFemScalarResultFrames*
         const std::vector<float>& referenceVerticalStrainData = verticalStrainDataFrames->frameData( referenceFrameIdx );
         const std::vector<float>& youngsModuliData            = youngsModuliFrames->frameData( fIdx );
         const std::vector<float>& poissonRatioData            = poissonRatioFrames->frameData( fIdx );
-        const std::vector<float>& voidRatioData               = voidRatioFrames->frameData( referenceFrameIdx );
+        const std::vector<float>& voidRatioData               = voidRatioFrames->frameData( 0 );
         const std::vector<float>& referencePorFrameData       = srcPORDataFrames->frameData( referenceFrameIdx );
         const std::vector<float>& porFrameData                = srcPORDataFrames->frameData( fIdx );
 
@@ -222,7 +222,7 @@ RigFemScalarResultFrames*
                             double bulkModulusFrame = youngsModuli / ( 3.0 * ( 1.0 - 2.0 * poissonRatio ) );
                             double bulkModulus      = bulkModulusFrame / ( 1.0 - biotCoefficient );
 
-                            // Calculate porosity
+                            // Calculate initial porosity (always from geostatic timestep)
                             double voidr    = voidRatioData[elmNodResIdx];
                             double porosity = voidr / ( 1.0 + voidr );