crevprop.temperature_field¶

A two-dimensional thermal model used to solve for temperature within the ice block defining model geometry.

Solvers used include a semi-implicit finite-difference scheme staggered leapfrog method for vertical advection upward motion scaled linearly from -b at the ice surface to 0 at the bed. the entire domain is advected horizontally at 200 m/a Plug flow with Lagrangian reference frame

For the plug flow simulation the model’s domain expands for each timestep at the upstream end of the horizontal model domain. The amount of ice added at the upstream end is determined by the annual ice velocity u, which when applied, pushes the model domain downstream. Once the model’s horzontal domain exceeds a length of 500 m, we remove the uppermost 200 m of domain. This model configuration allows the model’s domain to track the crevasse field as it advects downglacier and evolvs thermo-mechanically

Thermal model components include: - horizontal diffusion - vertical diffusion - latent heat transfer from refreezing

Functions

flatten(nested)

return flatened list

Classes

ThermalModel(iceblock_geometry, dt_T, ...[, ...])

Thermal model used to create IceBlock's temperature field.