From 30dea503923879ed89999b9474c44ac91b588ef5 Mon Sep 17 00:00:00 2001
From: Kristian Bendiksen <kristian.bendiksen@gmail.com>
Date: Wed, 8 Jul 2020 09:38:31 +0200
Subject: [PATCH] #6090 Fracture Model: export pressure at top of zone to frk
 file.

---
 .../ProjectDataModel/RimFractureModelPlot.cpp | 39 ++++++++++++++++++-
 .../ProjectDataModel/RimFractureModelPlot.h   |  4 ++
 2 files changed, 42 insertions(+), 1 deletion(-)

diff --git a/ApplicationCode/ProjectDataModel/RimFractureModelPlot.cpp b/ApplicationCode/ProjectDataModel/RimFractureModelPlot.cpp
index f6d70eb541..959372cd2f 100644
--- a/ApplicationCode/ProjectDataModel/RimFractureModelPlot.cpp
+++ b/ApplicationCode/ProjectDataModel/RimFractureModelPlot.cpp
@@ -191,6 +191,19 @@ void RimFractureModelPlot::computeAverageByLayer( const std::vector<std::pair<si
     }
 }
 
+//--------------------------------------------------------------------------------------------------
+///
+//--------------------------------------------------------------------------------------------------
+void RimFractureModelPlot::extractTopOfLayerValues( const std::vector<std::pair<size_t, size_t>>& layerBoundaryIndexes,
+                                                    const std::vector<double>&                    inputVector,
+                                                    std::vector<double>&                          result )
+{
+    for ( size_t i = 0; i < layerBoundaryIndexes.size(); i++ )
+    {
+        result.push_back( findValueAtTopOfLayer( inputVector, layerBoundaryIndexes, i ) );
+    }
+}
+
 //--------------------------------------------------------------------------------------------------
 ///
 //--------------------------------------------------------------------------------------------------
@@ -254,6 +267,30 @@ std::vector<double> RimFractureModelPlot::findCurveAndComputeLayeredAverage( Ria
     return result;
 }
 
+//--------------------------------------------------------------------------------------------------
+///
+//--------------------------------------------------------------------------------------------------
+std::vector<double> RimFractureModelPlot::findCurveAndComputeTopOfLayer( RiaDefines::CurveProperty curveProperty ) const
+{
+    RimWellLogExtractionCurve* curve = findCurveByProperty( curveProperty );
+    if ( !curve )
+    {
+        QString curveName = caf::AppEnum<RiaDefines::CurveProperty>::uiText( curveProperty );
+        RiaLogging::error( QString( "No curve for '%1' found" ).arg( curveName ) );
+        return std::vector<double>();
+    }
+
+    std::vector<std::pair<double, double>> layerBoundaryDepths;
+    std::vector<std::pair<size_t, size_t>> layerBoundaryIndexes;
+    calculateLayers( layerBoundaryDepths, layerBoundaryIndexes );
+
+    const RigWellLogCurveData* curveData = curve->curveData();
+    std::vector<double>        values    = curveData->xValues();
+    std::vector<double>        result;
+    extractTopOfLayerValues( layerBoundaryIndexes, values, result );
+    return result;
+}
+
 //--------------------------------------------------------------------------------------------------
 ///
 //--------------------------------------------------------------------------------------------------
@@ -267,7 +304,7 @@ std::vector<double> RimFractureModelPlot::calculatePorosity() const
 //--------------------------------------------------------------------------------------------------
 std::vector<double> RimFractureModelPlot::calculateReservoirPressure() const
 {
-    return findCurveAndComputeLayeredAverage( RiaDefines::CurveProperty::PRESSURE );
+    return findCurveAndComputeTopOfLayer( RiaDefines::CurveProperty::PRESSURE );
 }
 
 //--------------------------------------------------------------------------------------------------
diff --git a/ApplicationCode/ProjectDataModel/RimFractureModelPlot.h b/ApplicationCode/ProjectDataModel/RimFractureModelPlot.h
index 42a7127511..7c60dd9058 100644
--- a/ApplicationCode/ProjectDataModel/RimFractureModelPlot.h
+++ b/ApplicationCode/ProjectDataModel/RimFractureModelPlot.h
@@ -78,6 +78,7 @@ public:
 protected:
     std::vector<double> findCurveAndComputeLayeredAverage( RiaDefines::CurveProperty curveProperty ) const;
     std::vector<double> findCurveXValuesByProperty( RiaDefines::CurveProperty curveProperty ) const;
+    std::vector<double> findCurveAndComputeTopOfLayer( RiaDefines::CurveProperty curveProperty ) const;
 
     void                       calculateLayers( std::vector<std::pair<double, double>>& layerBoundaryDepths,
                                                 std::vector<std::pair<size_t, size_t>>& layerBoundaryIndexes ) const;
@@ -93,6 +94,9 @@ protected:
     static void   computeAverageByLayer( const std::vector<std::pair<size_t, size_t>>& layerBoundaryIndexes,
                                          const std::vector<double>&                    inputVector,
                                          std::vector<double>&                          result );
+    static void   extractTopOfLayerValues( const std::vector<std::pair<size_t, size_t>>& layerBoundaryIndexes,
+                                           const std::vector<double>&                    inputVector,
+                                           std::vector<double>&                          result );
 
     void defineUiOrdering( QString uiConfigName, caf::PdmUiOrdering& uiOrdering ) override;