From a1de02ef61560e4b66459e5960f5ba7b032edf2e Mon Sep 17 00:00:00 2001
From: Halvor M Nilsen <Halvor.M.Nilsen@sintef.no>
Date: Tue, 3 Dec 2024 10:10:05 +0100
Subject: [PATCH] -add blockoutput of stress

---
 opm/simulators/flow/OutputBlackoilModule.hpp | 38 +++++++++++++++++++-
 1 file changed, 37 insertions(+), 1 deletion(-)

diff --git a/opm/simulators/flow/OutputBlackoilModule.hpp b/opm/simulators/flow/OutputBlackoilModule.hpp
index 213900de794..ccc0203340d 100644
--- a/opm/simulators/flow/OutputBlackoilModule.hpp
+++ b/opm/simulators/flow/OutputBlackoilModule.hpp
@@ -152,7 +152,7 @@ class OutputBlackOilModule : public GenericOutputBlackoilModule<GetPropType<Type
         this->forceDisableFipresvOutput_ =
             Parameters::Get<Parameters::ForceDisableResvFluidInPlaceOutput>();
 
-        if (! Parameters::Get<Parameters::OwnerCellsFirst>()) {
+        if (!Parameters::Get<Parameters::OwnerCellsFirst>()){// && !Parameters::Get<Parameters::EnableEclOutput>()) {
             const std::string msg = "The output code does not support --owner-cells-first=false.";
             if (collectToIORank.isIORank()) {
                 OpmLog::error(msg);
@@ -882,6 +882,42 @@ class OutputBlackOilModule : public GenericOutputBlackoilModule<GetPropType<Type
                         else if (FluidSystem::phaseIsActive(waterPhaseIdx))
                             val.second = getValue(fs.temperature(waterPhaseIdx));
                     }
+                    else if (key.first == "BSTRSSXX")
+                        {
+                            const auto& model = problem.geoMechModel();
+                            auto stress = model.stress(globalDofIdx,/*include_fracture*/true);
+                            val.second = stress[0];
+                        }
+                    else if (key.first == "BSTRSSYY")
+                        {
+                            const auto& model = problem.geoMechModel();
+                            auto stress = model.stress(globalDofIdx,/*include_fracture*/true);
+                            val.second = stress[1];
+                        }
+                    else if (key.first == "BSTRSSZZ")
+                        {
+                            const auto& model = problem.geoMechModel();
+                            auto stress = model.stress(globalDofIdx,/*include_fracture*/true);
+                            val.second = stress[2];
+                        }
+                    else if (key.first == "BSTRSSXY")
+                        {
+                            const auto& model = problem.geoMechModel();
+                            auto stress = model.stress(globalDofIdx,/*include_fracture*/true);
+                            val.second = stress[5];
+                        }
+                    else if (key.first == "BSTRSSXZ")
+                        {
+                            const auto& model = problem.geoMechModel();
+                            auto stress = model.stress(globalDofIdx,/*include_fracture*/true);
+                            val.second = stress[4];
+                        }
+                    else if (key.first == "BSTRSSYZ")
+                        {
+                            const auto& model = problem.geoMechModel();
+                            auto stress = model.stress(globalDofIdx,/*include_fracture*/true);
+                            val.second = stress[3];
+                        }
                     else if (key.first == "BWKR" || key.first == "BKRW")
                         val.second = getValue(intQuants.relativePermeability(waterPhaseIdx));
                     else if (key.first == "BGKR" || key.first == "BKRG")