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

◆ compute_cts()

void pism::energy::compute_cts ( const array::Array3D ice_enthalpy,
const array::Scalar ice_thickness,
array::Array3D result 
)

Compute the CTS field, CTS = E/E_s(p), from ice_enthalpy and ice_thickness, and put in result.

The actual cold-temperate transition surface (CTS) is the level set CTS = 1.

Does not communicate ghosts for array::Array3D result.

Definition at line 179 of file utilities.cc.

References pism::array::Array3D::get_column(), pism::array::Array::grid(), pism::array::Array::inc_state_counter(), pism::k, pism::VariableMetadata::long_name(), pism::array::Array::metadata(), pism::array::Array::set_name(), pism::VariableMetadata::set_name(), and pism::VariableMetadata::units().

Referenced by pism::diagnostics::CTS::compute_impl().