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

◆ construct_surface_temp()

void pism::icebin::IBIceModel::construct_surface_temp ( pism::array::Scalar deltah,
double  default_val,
double  timestep_s,
pism::array::Scalar surface_temp 
)

Merges surface temperature derived from m_ice_enthalpy into any NaN values in the vector provided.

Parameters
deltahIN: Input from Icebin (change in enthalpy of each grid cell over the timestep) [W m-2].
default_valThe value that deltah(i,j) will have if no value is listed for that grid cell
timestep_sLength of the current coupling timestep [s]
surface_tempOUT: Resulting surface temperature to use as the Dirichlet B.C.

Merges surface temperature derived from m_ice_enthalpy into any NaN values in the vector provided.

Parameters
deltahIN: Input from Icebin (change in enthalpy of each m_grid cell over the timestep) [W m-2].
default_valThe value that deltah(i,j) will have if no value is listed for that m_grid cell
timestep_sLength of the current coupling timestep [s]
surface_tempOUT: Resulting surface temperature to use as the Dirichlet B.C.

Definition at line 420 of file IBIceModel.cc.

References pism::IceModel::ctx(), pism::Geometry::ice_thickness, pism::IceModel::m_config, pism::IceModel::m_energy_model, pism::IceModel::m_geometry, pism::IceModel::m_grid, and pism::printf().