diff --git a/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C b/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C
index 4473b8391dee292b5237eff8ea4fb1c3f7347eeb..3ec45ddfc91707e7cf16d0d7c9d3900215b06fc8 100644
--- a/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C
+++ b/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C
@@ -104,7 +104,7 @@ Foam::tmp<Foam::volScalarField> calcK
 
     // Do not correct BC
     // - operation may use inconsistent fields wrt these local manipulations
-    // k.correctBoundaryConditions();
+    //k.correctBoundaryConditions();
     correctProcessorPatches(k);
 
     Info<< "Writing k\n" << endl;
@@ -132,6 +132,7 @@ Foam::tmp<Foam::volScalarField> calcEpsilon
     scalar ce0 = ::pow(Cmu, 0.75)/kappa;
     epsilon = (1 - mask)*epsilon + mask*ce0*k*sqrt(k)/min(y, ybl);
     epsilon.max(SMALL);
+    epsilon.rename("epsilon");
 
     // Do not correct BC
     // - operation may use inconsistent fields wrt these local manipulations
@@ -297,6 +298,7 @@ void calcIncompressible
     // Calculate nut - reference nut is calculated by the turbulence model
     // on its construction
     volScalarField& nut = tnut();
+
     volScalarField S(mag(dev(symm(fvc::grad(U)))));
     nut = (1 - mask)*nut + mask*sqr(kappa*min(y, ybl))*::sqrt(2)*S;