diff --git a/src/coreComponents/physicsSolvers/fluidFlow/ThermalCompositionalMultiphaseBaseKernels.hpp b/src/coreComponents/physicsSolvers/fluidFlow/ThermalCompositionalMultiphaseBaseKernels.hpp
index a27d879c215..d39b7f51ddc 100644
--- a/src/coreComponents/physicsSolvers/fluidFlow/ThermalCompositionalMultiphaseBaseKernels.hpp
+++ b/src/coreComponents/physicsSolvers/fluidFlow/ThermalCompositionalMultiphaseBaseKernels.hpp
@@ -911,7 +911,7 @@ class ResidualNormKernel : public solverBaseKernels::ResidualNormKernelBase< 3 >
     // step 2: volume residual
 
     real64 const valVol = LvArray::math::abs( m_localResidual[stack.localRow + m_numComponents] ) / volumeNormalizer;
-    if( valVol > stack.localValue[0] )
+    if( valVol > stack.localValue[1] )
     {
       stack.localValue[1] = valVol;
     }
@@ -919,7 +919,7 @@ class ResidualNormKernel : public solverBaseKernels::ResidualNormKernelBase< 3 >
     // step 3: energy residual
 
     real64 const valEnergy = LvArray::math::abs( m_localResidual[stack.localRow + m_numComponents + 1] ) / energyNormalizer;
-    if( valEnergy > stack.localValue[1] )
+    if( valEnergy > stack.localValue[2] )
     {
       stack.localValue[2] = valEnergy;
     }