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

◆ set_ocean_input_fields()

void pism::ocean::Pico::set_ocean_input_fields ( const PicoPhysics physics,
const array::Scalar ice_thickness,
const array::CellType1 mask,
const array::Scalar basin_mask,
const array::Scalar shelf_mask,
const std::vector< double > &  basin_temperature,
const std::vector< double > &  basin_salinity,
array::Scalar Toc_box0,
array::Scalar Soc_box0 
) const
private

Set ocean ocean input from box 0 as boundary condition for box 1.

Set ocean temperature and salinity (Toc_box0, Soc_box0) from box 0 (in front of the ice shelf) as inputs for box 1, which is the ocean box adjacent to the grounding line.

We enforce that Toc_box0 is always at least the local pressure melting point.

Definition at line 461 of file Pico.cc.

References pism::array::Scalar::as_int(), pism::GlobalSum(), pism::Component::m_grid, pism::Component::m_log, m_n_basins, m_n_shelves, pism::MASK_FLOATING, pism::MASK_ICE_FREE_OCEAN, pism::array::max(), pism::ocean::PicoPhysics::pressure(), pism::array::Array2D< T >::star(), and pism::ocean::PicoPhysics::theta_pm().

Referenced by update_impl().