A difference method for the two-dimensional quench front problem using temperature as an independent variable is described and used to compute velocities of traveling wave solutions. Results indicate that this method is more efficient than standard difference methods for large values of the Biot number.