Use functions to return the subgrid T and rho rather than directly read it from the particle properties.