diff --git a/dumux/implicit/mpnc/mpncpropertydefaults.hh b/dumux/implicit/mpnc/mpncpropertydefaults.hh
index 3067f972bc595ac3520856eeeeea9ca7590d0276..0ca7b873491fcdd7fdc99715fbc89f317ebb0581 100644
--- a/dumux/implicit/mpnc/mpncpropertydefaults.hh
+++ b/dumux/implicit/mpnc/mpncpropertydefaults.hh
@@ -234,6 +234,11 @@ SET_PROP(MPNC, FluidState){
         typedef Dumux::CompositionalFluidState<Scalar, FluidSystem> type;
 };
 
+//! Set the default pressure formulation to the pressure of the (most) wetting phase
+SET_INT_PROP(MPNC,
+             PressureFormulation,
+             MpNcPressureFormulation::mostWettingFirst);
+
 //! default value for the forchheimer coefficient
 // Source: Ward, J.C. 1964 Turbulent flow in porous media. ASCE J. Hydraul. Div 90.
 //        Actually the Forchheimer coefficient is also a function of the dimensions of the