diff --git a/applications/utilities/postProcessing/turbulence/createTurbulenceFields/createTurbulenceFields.C b/applications/utilities/postProcessing/turbulence/createTurbulenceFields/createTurbulenceFields.C
index ce074ccc54..be63b2a5f7 100644
--- a/applications/utilities/postProcessing/turbulence/createTurbulenceFields/createTurbulenceFields.C
+++ b/applications/utilities/postProcessing/turbulence/createTurbulenceFields/createTurbulenceFields.C
@@ -50,7 +50,6 @@ int main(int argc, char *argv[])
     instantList timeDirs = timeSelector::select0(runTime, args);
 
     #include "createMesh.H"
-    #include "createFields.H"
 
     forAll(timeDirs, timeI)
     {
@@ -58,7 +57,7 @@ int main(int argc, char *argv[])
 
         Info<< "Time = " << runTime.timeName() << endl;
 
-        // Cache the turbulence fields
+        #include "createFields.H"
 
         Info<< "\nRetrieving field k from turbulence model" << endl;
         const volScalarField k(RASModel->k());