add THELCOEF field property

this is the thermal expansion coefficient for each cell in a mechanics model
write this to the .INIT file for relevant models

opm/input/eclipse/EclipseState/Grid/FieldProps.hpp

@@ -152,6 +152,7 @@ static const std::unordered_map<std::string, keyword_info<double>> double_keywor
                                                                                       {"BIOTCOEF", keyword_info<double>{}.unit_string("1")},
                                                                                       {"POELCOEF", keyword_info<double>{}.unit_string("1")},
                                                                                       {"THERMEXR", keyword_info<double>{}.unit_string("1/AbsoluteTemperature")},
+                                                                                      {"THELCOEF", keyword_info<double>{}.unit_string("Pressure/AbsoluteTemperature")},
                                                                                       {"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)},

src/opm/input/eclipse/share/keywords/900_OPM/T/THELCOEF

@@ -0,0 +1,11 @@
+{
+  "name": "THELCOEF",
+  "sections": [
+    "GRID"
+  ],
+  "description": "The THELCOEF item is used to set the thermal expansion coefficient for a cell in mechanics models.",
+  "data": {
+    "value_type": "DOUBLE",
+    "dimension": "1"
+  }
+}

src/opm/input/eclipse/share/keywords/keyword_list.cmake

@@ -1154,6 +1154,7 @@ set( keywords
      900_OPM/S/STRESSEQUILNUM
      900_OPM/S/SUREA
      900_OPM/S/SWOFLET
+     900_OPM/T/THELCOEF
      900_OPM/T/THERMEXR
      900_OPM/T/TLPMIXPA
      900_OPM/V/VAPWAT

src/opm/output/eclipse/WriteInit.cpp

@@ -448,6 +448,7 @@ namespace {
             {"POELCOEF" , ::Opm::UnitSystem::measure::identity},
             {"PRATIO"   , ::Opm::UnitSystem::measure::identity},
             {"THERMEXR" , ::Opm::UnitSystem::measure::identity}, // 1/(temperature difference)
+            {"THELCOEF" , ::Opm::UnitSystem::measure::identity}, // Pressure/Temperature
             {"YMODULE"  , ::Opm::UnitSystem::measure::ymodule},
         };