added temperature-dependent diffusivity function