Trigger WELL_THP_UPDATE also for WTMULT with THP control

WTMULT with THP control multiplies the well THP limit, so it must
cancel a THP limit imposed dynamically by the simulator like the other
THP re-specifications do.

Author: GitPaean

opm/input/eclipse/Schedule/Action/SimulatorUpdate.hpp

@@ -73,11 +73,12 @@ struct SimulatorUpdate
     std::unordered_set<std::string> welpi_wells{};
 
     /// Wells whose THP limit and/or VFP table has been (re)specified by
-    /// WCONPROD or WCONHIST for the named wells, or by WELTARG with THP or
-    /// VFP control. Unlike affected_wells, this set also includes wells for
-    /// which the entered values are unchanged, since re-specifying either
-    /// input cancels a THP limit imposed dynamically by the simulator (e.g.
-    /// by network balancing). See also the WELL_THP_UPDATE schedule event.
+    /// WCONPROD or WCONHIST for the named wells, by WELTARG with THP or VFP
+    /// control, or by WTMULT with THP control. Unlike affected_wells, this
+    /// set also includes wells for which the entered values are unchanged,
+    /// since re-specifying either input cancels a THP limit imposed
+    /// dynamically by the simulator (e.g. by network balancing). See also
+    /// the WELL_THP_UPDATE schedule event.
     std::unordered_set<std::string> thp_respec_wells{};
 
     /// New well connections created as a result of a geomechanical

opm/input/eclipse/Schedule/Events.hpp

@@ -115,11 +115,12 @@ namespace Opm
             REQUEST_OPEN_COMPLETION = (UINT64_C(1) << 24),
 
             /// The well THP limit and/or VFP table has been (re)specified
-            /// by WCONPROD or WCONHIST for the named wells, or by WELTARG
-            /// with THP or VFP control.  Unlike PRODUCTION_UPDATE, the event
-            /// is also triggered when the entered values are unchanged,
-            /// since re-specifying either input cancels a THP limit imposed
-            /// dynamically by the simulator (e.g. by network balancing).
+            /// by WCONPROD or WCONHIST for the named wells, by WELTARG with
+            /// THP or VFP control, or by WTMULT with THP control.  Unlike
+            /// PRODUCTION_UPDATE, the event is also triggered when the
+            /// entered values are unchanged, since re-specifying either
+            /// input cancels a THP limit imposed dynamically by the
+            /// simulator (e.g. by network balancing).
             WELL_THP_UPDATE = (UINT64_C(1) << 25),
         };
     } // namespace ScheduleEvents

opm/input/eclipse/Schedule/Well/WellPropertiesKeywordHandlers.cpp

@@ -633,6 +633,18 @@ void handleWTMULT(HandlerContext& handlerContext)
                 properties->handleWTMULT(cmode, factor.get<double>());
 
                 well.updateProduction(properties);
+
+                // The THP limit is re-specified (multiplied).
+                if (cmode == Well::WELTARGCMode::THP) {
+                    handlerContext.state().events()
+                        .addEvent(ScheduleEvents::WELL_THP_UPDATE);
+
+                    handlerContext.state().wellgroup_events()
+                        .addEvent(well_name, ScheduleEvents::WELL_THP_UPDATE);
+
+                    handlerContext.thp_respec_well(well_name);
+                }
+
                 if (update_well) {
                     handlerContext.state().events()
                         .addEvent(ScheduleEvents::PRODUCTION_UPDATE);

tests/parser/ScheduleTests.cpp

@@ -7288,6 +7288,14 @@ DATES        -- 5
  10 OKT 1998 /
 /
 
+WTMULT
+    'P1' 'THP' 2.0 /
+/
+
+DATES        -- 6
+ 10 NOV 1998 /
+/
+
 WCONHIST
     'P1' 'OPEN' 'ORAT' 100 0 0 /
 /
@@ -7307,6 +7315,8 @@ WCONHIST
     // WELTARG VFP re-specifies the VFP table, also for an unchanged value
     BOOST_CHECK( schedule[4].wellgroup_events().hasEvent("P1", ScheduleEvents::WELL_THP_UPDATE));
     BOOST_CHECK(!schedule[4].wellgroup_events().hasEvent("P1", ScheduleEvents::PRODUCTION_UPDATE));
-    // WCONHIST re-specifies the VFP table
+    // WTMULT with THP control re-specifies (multiplies) the THP limit
     BOOST_CHECK( schedule[5].wellgroup_events().hasEvent("P1", ScheduleEvents::WELL_THP_UPDATE));
+    // WCONHIST re-specifies the VFP table
+    BOOST_CHECK( schedule[6].wellgroup_events().hasEvent("P1", ScheduleEvents::WELL_THP_UPDATE));
 }