diff --git a/applications/solvers/modules/multiphaseEuler/multiphaseEuler/thermophysicalPredictor.C b/applications/solvers/modules/multiphaseEuler/multiphaseEuler/thermophysicalPredictor.C
index 280e6d8575..0389c451de 100644
--- a/applications/solvers/modules/multiphaseEuler/multiphaseEuler/thermophysicalPredictor.C
+++ b/applications/solvers/modules/multiphaseEuler/multiphaseEuler/thermophysicalPredictor.C
@@ -121,9 +121,10 @@ void Foam::solvers::multiphaseEuler::thermophysicalPredictor()
             compositionPredictor();
             energyPredictor();
 
-            forAll(phases, phasei)
+            forAll(fluid.anisothermalPhases(), anisothermalPhasei)
             {
-                phaseModel& phase = phases[phasei];
+                phaseModel& phase =
+                    fluid.anisothermalPhases()[anisothermalPhasei];
 
                 Info<< phase.name() << " min/max T "
                     << min(phase.thermo().T()).value()