diff --git a/LDDsimulation/domainPatch.py b/LDDsimulation/domainPatch.py
index 92d2aef7de3fe16f2dcb528aba2f9486f871fc6a..098089a829a23133d5023f7946a6cdc42744712f 100644
--- a/LDDsimulation/domainPatch.py
+++ b/LDDsimulation/domainPatch.py
@@ -171,7 +171,7 @@ class DomainPatch(df.SubDomain):
         if self.isRichards:
             self.has_phases = ['wetting']
         else:
-            self.has_phases = ['nonwetting', 'wetting'] #['wetting', 'nonwetting']
+            self.has_phases = ['wetting', 'nonwetting']  #['nonwetting', 'wetting']
         # dictionary holding the dof indices corresponding to an interface of
         # given interface. see self._calc_dof_indices_of_interfaces()
         self._dof_indices_of_interface = dict()