*** UEqns.H	2010-09-05 21:22:27.175381992 -0500
--- UEqnsPatched.H	2010-09-06 15:53:47.426940167 -0500
***************
*** 4,14 ****
  {
      {
          volTensorField gradUaT = fvc::grad(Ua)().T();
-         volTensorField Rca
-         (
-             "Rca",
-             ((2.0/3.0)*I)*(sqr(Ct)*k + nuEffa*tr(gradUaT)) - nuEffa*gradUaT
-         );
  
          if (kineticTheory.on())
          {
--- 4,9 ----
***************
*** 19,24 ****
--- 14,25 ----
          {
              nuEffa = sqr(Ct)*nutb + nua;
          }
+         
+         volTensorField Rca
+         (
+             "Rca",
+             ((2.0/3.0)*I)*(sqr(Ct)*k + nuEffa*tr(gradUaT)) - nuEffa*gradUaT
+         );
  
          if (kineticTheory.on())
          {