const dimensionedScalar& rho1f(rho1);
const dimensionedScalar& rho2f(rho2);