diff --git a/src/functionObjects/field/turbulenceFields/turbulenceFieldsTemplates.C b/src/functionObjects/field/turbulenceFields/turbulenceFieldsTemplates.C
index 5f17e49434300a3e9c806561b52072fc66884759..c94b118b8439e0ac12f5343e80477a949d3308fa 100644
--- a/src/functionObjects/field/turbulenceFields/turbulenceFieldsTemplates.C
+++ b/src/functionObjects/field/turbulenceFields/turbulenceFieldsTemplates.C
@@ -118,14 +118,15 @@ Foam::functionObjects::turbulenceFields::I
     const Model& model
 ) const
 {
-    // Assume k is available
-    const volScalarField uPrime(sqrt((2.0/3.0)*model.k()));
-    const dimensionedScalar U0("U0", dimVelocity, SMALL);
+    // (P:p. 183)
+    // root-mean-square of velocity fluctuations - isotropic turbulence
+    tmp<volScalarField> uPrime = sqrt((2.0/3.0)*model.k());
+    const dimensionedScalar U0(dimVelocity, SMALL);
 
     return tmp<volScalarField>::New
     (
         "I.tmp",
-        uPrime/max(max(uPrime, mag(model.U())), U0)
+        uPrime/max(mag(model.U()), U0)
     );
 }