diff --git a/src/equil/MultiPhaseEquil.cpp b/src/equil/MultiPhaseEquil.cpp
index 4d3ba51f5..be846cb65 100644
--- a/src/equil/MultiPhaseEquil.cpp
+++ b/src/equil/MultiPhaseEquil.cpp
@@ -66,6 +66,7 @@ MultiPhaseEquil::MultiPhaseEquil(MultiPhase* mix, bool start, int loglevel) : m_
     m_temp = mix->temperature();
 
     index_t m, k;
+    m_force = true;
     m_nel = 0;
     m_nsp = 0;
     m_eloc = 1000;