From 217f8f7b6014de58a31f6e5e323b54a3e9f1e5f8 Mon Sep 17 00:00:00 2001
From: hnil <Halvor.M.Nilsen@sintef.no>
Date: Thu, 17 Aug 2023 08:23:59 +0200
Subject: [PATCH] add YMODULE field property

this is the Young's modulus for each cell in a mechanics model
write this to the .INIT file for relevant models
---
 opm/input/eclipse/EclipseState/Grid/FieldProps.hpp    |  1 +
 .../input/eclipse/share/keywords/900_OPM/Y/YMODULE    | 11 +++++++++++
 .../input/eclipse/share/keywords/keyword_list.cmake   |  1 +
 src/opm/output/eclipse/WriteInit.cpp                  |  1 +
 4 files changed, 14 insertions(+)
 create mode 100644 src/opm/input/eclipse/share/keywords/900_OPM/Y/YMODULE

diff --git a/opm/input/eclipse/EclipseState/Grid/FieldProps.hpp b/opm/input/eclipse/EclipseState/Grid/FieldProps.hpp
index ed6f7cb07..207813233 100644
--- a/opm/input/eclipse/EclipseState/Grid/FieldProps.hpp
+++ b/opm/input/eclipse/EclipseState/Grid/FieldProps.hpp
@@ -147,6 +147,7 @@ static const std::unordered_map<std::string, keyword_info<double>> double_keywor
                                                                                       {"THCOIL",  keyword_info<double>{}.unit_string("Energy/AbsoluteTemperature*Length*Time")},
                                                                                       {"THCGAS",  keyword_info<double>{}.unit_string("Energy/AbsoluteTemperature*Length*Time")},
                                                                                       {"THCWATER",keyword_info<double>{}.unit_string("Energy/AbsoluteTemperature*Length*Time")},
+                                                                                      {"YMODULE", keyword_info<double>{}.unit_string("Giga*Pascal")},
                                                                                       {"MULTX",   keyword_info<double>{}.init(1.0).mult(true)},
                                                                                       {"MULTX-",  keyword_info<double>{}.init(1.0).mult(true)},
                                                                                       {"MULTY",   keyword_info<double>{}.init(1.0).mult(true)},
diff --git a/src/opm/input/eclipse/share/keywords/900_OPM/Y/YMODULE b/src/opm/input/eclipse/share/keywords/900_OPM/Y/YMODULE
new file mode 100644
index 000000000..e8c05bb65
--- /dev/null
+++ b/src/opm/input/eclipse/share/keywords/900_OPM/Y/YMODULE
@@ -0,0 +1,11 @@
+{
+  "name": "YMODULE",
+  "sections": [
+    "GRID"
+  ],
+  "description": "The YMODULE item is used to set the Young's modulus for a cell in mechanics models.",
+  "data": {
+    "value_type": "DOUBLE",
+    "dimension": "Ymodule"
+  }
+}
diff --git a/src/opm/input/eclipse/share/keywords/keyword_list.cmake b/src/opm/input/eclipse/share/keywords/keyword_list.cmake
index cd2794845..87644a6a5 100644
--- a/src/opm/input/eclipse/share/keywords/keyword_list.cmake
+++ b/src/opm/input/eclipse/share/keywords/keyword_list.cmake
@@ -1152,6 +1152,7 @@ set( keywords
      900_OPM/S/SWOFLET
      900_OPM/T/TLPMIXPA
      900_OPM/V/VAPWAT
+     900_OPM/Y/YMODULE
      900_OPM/W/WATJT
      900_OPM/W/WELTRAJ
      900_OPM/W/WINJCLN
diff --git a/src/opm/output/eclipse/WriteInit.cpp b/src/opm/output/eclipse/WriteInit.cpp
index 270e52264..278fb9c3e 100644
--- a/src/opm/output/eclipse/WriteInit.cpp
+++ b/src/opm/output/eclipse/WriteInit.cpp
@@ -444,6 +444,7 @@ namespace {
             {"PERMY" , ::Opm::UnitSystem::measure::permeability },
             {"PERMZ" , ::Opm::UnitSystem::measure::permeability },
             {"NTG"   , ::Opm::UnitSystem::measure::identity },
+            {"YMODULE"   , ::Opm::UnitSystem::measure::ymodule },
         };
 
         // The INIT file should always contain the NTG property, we