PISM, A Parallel Ice Sheet Model
stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800
|
◆ compute_enthalpy_cold()
Compute ice enthalpy from temperature temperature by assuming the ice has zero liquid fraction. First this method makes sure the temperatures is at most the pressure-melting value, before computing the enthalpy for that temperature, using zero liquid fraction. Because of how EnthalpyConverter::pressure() works, the energy content in the air is set to the value that ice would have if it a chunk of it occupied the air; the atmosphere actually has much lower energy content. It is done this way for regularity (i.e. dEnth/dz computations). Definition at line 48 of file utilities.cc. References pism::array::Array3D::get_column(), pism::array::Array::grid(), pism::array::Array::inc_state_counter(), pism::k, and pism::array::Array::update_ghosts(). Referenced by bootstrap_ice_enthalpy(), pism::energy::TemperatureModel::bootstrap_impl(), pism::energy::EnergyModel::init_enthalpy(), pism::energy::TemperatureModel::initialize_impl(), pism::energy::TemperatureModel_Verification::initialize_impl(), pism::energy::TemperatureModel::restart_impl(), and pism::energy::TemperatureModel::update_impl(). |