Inclusion of temperature variations in a numerical scheme for partial equilibrium flow