const volScalarField& rho = trho(); volScalarField& h = thermo.he(); const volScalarField& betav = *betavPtr;