Loading [MathJax]/extensions/tex2jax.js
PISM, A Parallel Ice Sheet Model 2.2.2-d6b3a29ca committed by Constantine Khrulev on 2025-03-28
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages

◆ cryo_hydrologic_warming_flux()

void pism::energy::cryo_hydrologic_warming_flux ( double  k,
double  R,
const array::Scalar ice_thickness,
const array::Array3D ice_enthalpy,
const array::Array3D ch_enthalpy,
array::Array3D result 
)

Compute the heat flux corresponding to the cryo-hydrologic warming.

Q = (k / R**2) * (T_ch - T_ice),

where k is the thermal conductivity of ice and R us the average spacing of channels in the cryo-hydrologic system.

Definition at line 251 of file CHSystem.cc.

References pism::ParallelSection::check(), pism::ParallelSection::failed(), pism::array::Array3D::get_column(), pism::array::Array::grid(), and pism::k.

Referenced by pism::CHHeatFlux::compute_impl(), and pism::IceRegionalModel::energy_step().