diff --git a/src/scheme/HybridHLLcRusanovEulerianCompositeSolver_v2.cpp b/src/scheme/HybridHLLcRusanovEulerianCompositeSolver_v2.cpp
index 373dc1c25bbfd8acf252d8bd2a5c9fd54e773b07..1a30528b9e6a38b63fe39841ed244cd053977f86 100644
--- a/src/scheme/HybridHLLcRusanovEulerianCompositeSolver_v2.cpp
+++ b/src/scheme/HybridHLLcRusanovEulerianCompositeSolver_v2.cpp
@@ -1230,7 +1230,7 @@ class HybridHLLcRusanovEulerianCompositeSolver_v2
                 for (size_t d = 2; d < (Dimension + 1); ++d){
                   DL[d] = Uj[d-1];
                 }
-                DL[Dimension + 1] = Ej + (SC - uL) * (SC + (Pressionj / (rhoj * diffVelocityj)));
+                DL[Dimension + 1] = Ej + rhoj * (SC - uL) * (SC + (Pressionj / (rhoj * diffVelocityj)));
                   
                 const Rp UCL = rhoj * (diffVelocityj / (MinVpNormjk - SC)) * DL;
                 const Rp diffStates = MinVpNormjk * l2Norm(Cjf_loc) * (UCL - StateAtFace[j][l]);
@@ -1250,7 +1250,7 @@ class HybridHLLcRusanovEulerianCompositeSolver_v2
                 for (size_t d = 2; d < (Dimension + 1); ++d){
                   DR[d] = UK[d-1];            
                 }
-                DR[Dimension + 1] = EK + (SC - uR) * (SC + (PressionK / (rhoK * diffVelocityK)));
+                DR[Dimension + 1] = EK + rhoK * (SC - uR) * (SC + (PressionK / (rhoK * diffVelocityK)));
                     
                 const Rp UCR = rhoK * (diffVelocityK / (MaxVpNormjk - SC)) * DR;
                 const Rp diffStates = MaxVpNormjk * l2Norm(Cjf_loc) * (UCR - StateAtFace[K][R]);