diff --git a/applications/solvers/multiphase/icoReactingMultiphaseInterFoam/phasesSystem/phaseSystem/multiphaseSystem/multiphaseSystem.C b/applications/solvers/multiphase/icoReactingMultiphaseInterFoam/phasesSystem/phaseSystem/multiphaseSystem/multiphaseSystem.C
index f4a1414df33f6f1cde51e2489ec647f08eae8bc3..48b980e003bdbafd1d56f80d273f47854d9b6415 100644
--- a/applications/solvers/multiphase/icoReactingMultiphaseInterFoam/phasesSystem/phaseSystem/multiphaseSystem/multiphaseSystem.C
+++ b/applications/solvers/multiphase/icoReactingMultiphaseInterFoam/phasesSystem/phaseSystem/multiphaseSystem/multiphaseSystem.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2017 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 2017-2019 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -59,7 +59,7 @@ Foam::multiphaseSystem::multiphaseSystem
 )
 :
     phaseSystem(mesh),
-    cAlphas_(mesh.solverDict("alpha").lookup("cAlphas")),
+    cAlphas_(),
     ddtAlphaMax_(0.0),
     limitedPhiAlphas_(phaseModels_.size()),
     Su_(phaseModels_.size()),
@@ -73,6 +73,8 @@ Foam::multiphaseSystem::multiphaseSystem
         phases_.set(phasei++, &pm);
     }
 
+    mesh.solverDict("alpha").readEntry("cAlphas", cAlphas_);
+
     // Initiate Su and Sp
     forAllConstIters(phaseModels_, iter)
     {