diff --git a/applications/solvers/incompressible/simpleFoam/Allwmake b/applications/solvers/incompressible/simpleFoam/Allwmake
new file mode 100755
index 0000000000000000000000000000000000000000..577d11018fd0905d54cae6bba8e43bcb5f4f474c
--- /dev/null
+++ b/applications/solvers/incompressible/simpleFoam/Allwmake
@@ -0,0 +1,8 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # run from this directory
+set -x
+
+wmake
+wmake porousSimpleFoam
+
+# ----------------------------------------------------------------- end-of-file
diff --git a/applications/solvers/incompressible/simpleFoam/porousSimpleFoam/createPorousZones.H b/applications/solvers/incompressible/simpleFoam/porousSimpleFoam/createPorousZones.H
new file mode 100644
index 0000000000000000000000000000000000000000..d14adbb493da24b6c82d5a18cef2cd330163fa96
--- /dev/null
+++ b/applications/solvers/incompressible/simpleFoam/porousSimpleFoam/createPorousZones.H
@@ -0,0 +1,26 @@
+    porousZones pZones(mesh);
+    Switch pressureImplicitPorosity(false);
+
+    // nUCorrectors used for pressureImplicitPorosity
+    int nUCorr = 0;
+    if (pZones.size())
+    {
+        // nUCorrectors for pressureImplicitPorosity
+        if (mesh.solutionDict().subDict("SIMPLE").found("nUCorrectors"))
+        {
+            nUCorr = readInt
+            (
+                mesh.solutionDict().subDict("SIMPLE").lookup("nUCorrectors")
+            );
+        }
+
+        if (nUCorr > 0)
+        {
+            pressureImplicitPorosity = true;
+            Info<< "Using pressure implicit porosity" << endl;
+        }
+        else
+        {
+            Info<< "Using pressure explicit porosity" << endl;
+        }
+    }
diff --git a/src/engine/include/createEngineMesh.H b/src/engine/include/createEngineMesh.H
index 88a9e06686e181d6ed10c0569e55018a6a909651..263ebed5ff64edf6050930a3cd537e5b833ddf07 100644
--- a/src/engine/include/createEngineMesh.H
+++ b/src/engine/include/createEngineMesh.H
@@ -9,7 +9,8 @@
             (
                 engineMesh::defaultRegion,
                 runTime.timeName(),
-                runTime
+                runTime,
+                Foam::IOobject::MUST_READ
             )
         )
     );
diff --git a/tutorials/combustion/engineFoam/kivaTest/system/fvSchemes b/tutorials/combustion/engineFoam/kivaTest/system/fvSchemes
index 7d9e15bae0202334efd1ef07b4655299108cbe59..0a28cdbe5810ab4d6084efb2531f725eb340a106 100644
--- a/tutorials/combustion/engineFoam/kivaTest/system/fvSchemes
+++ b/tutorials/combustion/engineFoam/kivaTest/system/fvSchemes
@@ -60,7 +60,7 @@ laplacianSchemes
     laplacian(DepsilonEff,epsilon) Gauss linear limited 0.5;
     laplacian(DREff,R) Gauss linear limited 0.5;
     laplacian((rho*(1|A(U))),p) Gauss linear limited 0.5;
-    laplacian(muEff,b) Gauss linear limited 0.5;
+    laplacian(alphaEff,b) Gauss linear limited 0.5;
     laplacian(muEff,ft) Gauss linear limited 0.5;
     laplacian(alphaEff,h) Gauss linear limited 0.5;
     laplacian(alphaEff,hu) Gauss linear limited 0.5;