volScalarField& p = thermo.p();
const volScalarField& T = thermo.T();
const volScalarField& psi = thermo.psi();
const volScalarField& mu = thermo.mu();
const label inertIndex(composition.species()[inertSpecie]);
bool inviscid(true);
if (max(mu.primitiveField()) > 0.0)
{
    inviscid = false;
}
