move volScalarField& T = thermo.T(); const volScalarField& psi = thermo.psi(); to createFieldRefs.H (and make T const)