PISM, A Parallel Ice Sheet Model 2.2.1-cd005eec8 committed by Constantine Khrulev on 2025-03-07
|
◆ update_impl()
Remove "icebergs" using the finite element notion of connectivity: two elements are connected if they share a boundary.
Once "iceberg" elements are labeled we need to remove nodes that belong to icebergs but do not belong to any elements connected to grounded ice.
Reimplemented from pism::calving::IcebergRemover. Definition at line 67 of file IcebergRemoverFEM.cc. References pism::fem::Element2::add_contribution(), pism::array::Array2D< T >::array(), pism::array::Scalar::as_int(), pism::mask::grounded(), pism::mask::icy(), pism::connected_components::label_isolated(), pism::Component::m_grid, pism::calving::IcebergRemover::m_iceberg_mask, m_mask, pism::MASK_ICE_FREE_OCEAN, n, pism::fem::q1::n_chi, pism::fem::Element2::nodal_values(), pism::fem::Element2::reset(), pism::array::Array::set(), pism::array::Array::stencil_width(), and pism::array::Array::update_ghosts(). |