diff --git a/applications/solvers/multiphase/compressibleMultiphaseInterFoam/pEqn.H b/applications/solvers/multiphase/compressibleMultiphaseInterFoam/pEqn.H
index 649063da7..f0d7e05ed 100644
--- a/applications/solvers/multiphase/compressibleMultiphaseInterFoam/pEqn.H
+++ b/applications/solvers/multiphase/compressibleMultiphaseInterFoam/pEqn.H
@@ -106,6 +106,8 @@
                 phase().dgdt() =
                     pos0(phase())
                   *(p_rghEqnComps[phasei] & p_rgh)/phase().thermo().rho();
+
+                phasei++;
             }
 
             phi = phiHbyA + p_rghEqnIncomp.flux();