Finite difference time domain code using operator splitting to solve coupled flow and heat equations.