if (runTime.outputTime())
{
  //K.write();

  volScalarField epsilonEq
  (
      IOobject
      (
          "epsilonEq",
          runTime.timeName(),
          mesh,
          IOobject::NO_READ,
          IOobject::AUTO_WRITE
      ),
      sqrt((2.0/3.0)*magSqr(dev(epsilon)))
  );

  Info<< "Max epsilonEq = " << max(epsilonEq).value() << endl;

  volScalarField sigmaEq
  (
      IOobject
      (
          "sigmaEq",
          runTime.timeName(),
          mesh,
          IOobject::NO_READ,
          IOobject::AUTO_WRITE
      ),
      sqrt((3.0/2.0)*magSqr(dev(sigma)))
  );

  Info<< "Max sigmaEq = " << max(sigmaEq).value() << endl;

//  runTime.write();
}
