diff --git a/src/cpu/VirtualFluidsCore/BoundaryConditions/NoSlipBCAlgorithmMultiphase.cpp b/src/cpu/VirtualFluidsCore/BoundaryConditions/NoSlipBCAlgorithmMultiphase.cpp
index d32a42d9a362e20d6985474397268a8679be4717..ac64d31f68808d091006fa977d9e2ac7c3824bac 100644
--- a/src/cpu/VirtualFluidsCore/BoundaryConditions/NoSlipBCAlgorithmMultiphase.cpp
+++ b/src/cpu/VirtualFluidsCore/BoundaryConditions/NoSlipBCAlgorithmMultiphase.cpp
@@ -78,8 +78,8 @@ void NoSlipBCAlgorithmMultiphase::applyBC()
    //LBMReal collFactorM = collFactorL + (collFactorL - collFactorG)*(phi - phiH)/(phiH - phiL);
    
    //rho = phi + (1.0 - phi)*1.0/densityRatio;
-   LBMReal rhoH = 1.0;
-   LBMReal rhoL = 1.0/densityRatio;
+   //LBMReal rhoH = 1.0;
+   //LBMReal rhoL = 1.0/densityRatio;
    //rho = rhoH + (rhoH - rhoL)*(phi - phiH)/(phiH - phiL);
    
    calcMacrosFct(f, p1, vx1, vx2, vx3);