set ws.thp = thp_limit for wells under THP control

opm/simulators/wells/StandardWellEval.cpp

@@ -572,6 +572,12 @@ updateThp(WellState& well_state,
         ws.thp = 0;
         return;
     }
+    bool thp_controled = baseif_.isInjector() ? ws.injection_cmode == Well::InjectorCMode::THP:
+                                              ws.production_cmode == Well::ProducerCMode::THP;
+    if (thp_controled) {
+        return;
+    }
+
 
     // the well is under other control types, we calculate the thp based on bhp and rates
     std::vector<double> rates(3, 0.0);