diff --git a/src/cpu/VirtualFluidsCore/BoundaryConditions/NoSlipBCAlgorithmMultiphase.cpp b/src/cpu/VirtualFluidsCore/BoundaryConditions/NoSlipBCAlgorithmMultiphase.cpp
index 7ea78b5269565bb2e6030f9736b7dd3e4e163f9d..d32a42d9a362e20d6985474397268a8679be4717 100644
--- a/src/cpu/VirtualFluidsCore/BoundaryConditions/NoSlipBCAlgorithmMultiphase.cpp
+++ b/src/cpu/VirtualFluidsCore/BoundaryConditions/NoSlipBCAlgorithmMultiphase.cpp
@@ -70,7 +70,7 @@ void NoSlipBCAlgorithmMultiphase::applyBC()
    LBMReal heq[D3Q27System::ENDF+1];
    distributions ->getDistributionInv(f, x1, x2, x3);
    distributionsH->getDistributionInv(h, x1, x2, x3);
-   LBMReal phi, rho, vx1, vx2, vx3, p1;
+   LBMReal phi, vx1, vx2, vx3, p1;
    
    D3Q27System::calcDensity(h, phi);
    
@@ -80,7 +80,7 @@ void NoSlipBCAlgorithmMultiphase::applyBC()
    //rho = phi + (1.0 - phi)*1.0/densityRatio;
    LBMReal rhoH = 1.0;
    LBMReal rhoL = 1.0/densityRatio;
-   rho = rhoH + (rhoH - rhoL)*(phi - phiH)/(phiH - phiL);
+   //rho = rhoH + (rhoH - rhoL)*(phi - phiH)/(phiH - phiL);
    
    calcMacrosFct(f, p1, vx1, vx2, vx3);
    /*vx1/=(rho*c1o3);