A comparison between parameterized and two dimesional finite element reservoir cooling models is described. The models feature strongly temperature dependent Newtonian rheology and thermal boundary conditions simulating sills cooled from above. Cooling histories predicted by the parameterized scheme are found to be in excellent agreement with the more exact finite calculations for viscosity variations of up to 105 occurring at a given time in the reservoir.