need to work on enthalpy-dependent diffusivity