diff --git a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/pEqn.H b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/pEqn.H
index bbb8765e1457861ced59a080e25accc437da8d1f..ea3e825448c1f644eda5055c20fd57c510988e24 100644
--- a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/pEqn.H
+++ b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/pEqn.H
@@ -84,23 +84,12 @@
               + (alpha2/rho2)*pEqnComp2()
             ) +
             pEqnIncomp,
-            mesh.solver
-            (
-                p.select(pimple.finalInnerIter(corr, nonOrth))
-                /*
-                (
-                    oCorr == pimple.nOuterCorr()-1
-                 && corr == nCorr-1
-                 && nonOrth == nNonOrthCorr
-                )
-                */
-            )
+            mesh.solver(p.select(pimple.finalInnerIter(corr, nonOrth)))
         );
 
         if (nonOrth == pimple.nNonOrthCorr())
         {
             surfaceScalarField mSfGradp = pEqnIncomp.flux()/Dp;
-
             phi1 += rAlphaAU1f*mSfGradp/fvc::interpolate(rho1);
             phi2 += rAlphaAU2f*mSfGradp/fvc::interpolate(rho2);
             phi = alpha1f*phi1 + alpha2f*phi2;
diff --git a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/readTwoPhaseEulerFoamControls.H b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/readTwoPhaseEulerFoamControls.H
index eec13762cc7eb982aca4e5a355ed34e152cc82d9..9767d7c845f7b55f05995fae8a8bf32d44ebf69f 100644
--- a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/readTwoPhaseEulerFoamControls.H
+++ b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/readTwoPhaseEulerFoamControls.H
@@ -1,4 +1,3 @@
     #include "readTimeControls.H"
 
     int nAlphaCorr(readInt(pimple.dict().lookup("nAlphaCorr")));
-    Switch correctAlpha(pimple.dict().lookup("correctAlpha"));