PISM, A Parallel Ice Sheet Model  stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800

◆ 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 C, pism::ParallelSection::check(), pism::ParallelSection::failed(), pism::array::Array3D::get_column(), pism::array::Array::grid(), pism::k, and pism::array::max().

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