const volScalarField& psi = thermo.psi();
//const volScalarField& T = thermo.T();