diffusion

This module contains the code for the diffusion scheme models.

code to handle diffusion (heat, enthalpy) calls solver Draws from Numerical Heat Transfer and Heat Flow (Patankar, 1980) Enthalpy from Voller, Swaminathan, and Thomas (1990) Isotope diffusion now has its own class.

diffusion.enthalpyDiff(self, iii)

enthalpy diffusion function, new 1/30/19 - method from Voller and Swaminathan LWC is in volume (m^3) thermal diffusivity: alpha = K_firn / (rho*c_firn)

diffusion.firnConductivity(self, iii, K_ice)

Function to set the firn’s thermal conductivity based on one of a number of parameterizations.

Choose your favorite! Default is Calonne et al. 2019. References are provided at the end of this script.

diffusion.heatDiff(self, iii)

Heat diffusion function

Parameters

• z –

• dz –

• Ts –

• rho –

Returns self.Tz

Returns self.T10m

thermal diffusivity: alpha = K_firn / (rho*c_firn)