OpenFOAM-4.x-lab/applications/solvers/lagrangian/ThermoMPPICFoam/createFieldRefs.H

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