True. Actually smagorinsky has an intermediate k (and also epsilon) variable in its formulation, solved from a quadratic equation. I wonder if there's a switch in dictionary to write the generic turbulence->k() variable.. for now I simply had to recompile my solver with a variable to store k(), but would be much easier to just use a switch in turbulence model, since right now I am blocking memory unnecessarily, just to maintain compatibility with RAS and LES.