From c23c64893ec2a99df0979e24f38fba32fec0a216 Mon Sep 17 00:00:00 2001
From: Tor Harald Sandve <tosa@norceresearch.no>
Date: Fri, 21 Oct 2022 15:34:05 +0200
Subject: [PATCH] fix bug in boundary conditions for energy by adding the
 scaling factor

---
 opm/models/blackoil/blackoilboundaryratevector.hh | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/opm/models/blackoil/blackoilboundaryratevector.hh b/opm/models/blackoil/blackoilboundaryratevector.hh
index 090da47c8..6c9449e2f 100644
--- a/opm/models/blackoil/blackoilboundaryratevector.hh
+++ b/opm/models/blackoil/blackoilboundaryratevector.hh
@@ -161,7 +161,7 @@ public:
                 }
 
                 Evaluation enthalpyRate = density*extQuants.volumeFlux(phaseIdx)*specificEnthalpy;
-                EnergyModule::addToEnthalpyRate(*this, enthalpyRate);
+                EnergyModule::addToEnthalpyRate(*this, enthalpyRate*getPropValue<TypeTag, Properties::BlackOilEnergyScalingFactor>());
             }
         }
 
@@ -189,7 +189,7 @@ public:
 
         // heat conduction
         if constexpr (enableEnergy)
-            EnergyModule::addToEnthalpyRate(*this, extQuants.energyFlux());
+            EnergyModule::addToEnthalpyRate(*this, extQuants.energyFlux()*getPropValue<TypeTag, Properties::BlackOilEnergyScalingFactor>());
 
 #ifndef NDEBUG
         for (unsigned i = 0; i < numEq; ++i) {