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

◆ compute_node_type()

void pism::stressbalance::Blatter::compute_node_type ( double  min_thickness)
protected

Compute node type using domain thickness and the thickness threshold min_thickness.

An element contains ice if ice thickness at all its nodes equal or exceeds the min_thickness threshold.

A node is interior if all four elements it belongs to contain ice.

A node is exterior if it belongs to zero icy elements.

A node that is neither interior nor exterior is a boundary node.

Definition at line 56 of file Blatter.cc.

References pism::grid_transpose(), pism::k, pism::fem::Element2::local_to_global(), m_da, pism::Component::m_grid, m_parameters, pism::fem::q1::n_chi, pism::fem::Element2::nodal_values(), pism::NODE_BOUNDARY, pism::NODE_EXTERIOR, pism::NODE_INTERIOR, PISM_CHK, pism::fem::Element2::reset(), pism::array::Array::set(), and pism::array::Array::update_ghosts().

Referenced by init_2d_parameters().