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

◆ compute_grounded_cell_fraction()

void pism::compute_grounded_cell_fraction ( double  ice_density,
double  ocean_density,
const array::Scalar1 sea_level,
const array::Scalar1 ice_thickness,
const array::Scalar1 bed_topography,
array::Scalar result 
)
Parameters
[in]ice_densityice density, kg/m3
[in]ocean_densityocean_density, kg/m3
[in]sea_levelsea level (flotation) elevation, m
[in]ice_thicknessice thickness, m
[in]bed_topographybed elevation, m
[out]resultgrounded cell fraction, between 0 (floating) and 1 (grounded)

Compute grounded cell fractions by splitting control volumes into triangles and treating the flotation criterion as a linear function on each triangle.

Definition at line 257 of file grounded_cell_fraction.cc.

References pism::array::Array2D< T >::box(), pism::stencils::Box< T >::c, pism::ParallelSection::check(), clip(), pism::stencils::Box< T >::e, F(), pism::ParallelSection::failed(), pism::array::Array::grid(), grounded_area_fraction(), pism::stencils::Box< T >::n, pism::stencils::Box< T >::ne, pism::stencils::Box< T >::nw, S(), pism::stencils::Box< T >::s, pism::stencils::Box< T >::se, pism::stencils::Box< T >::sw, and pism::stencils::Box< T >::w.

Referenced by pism::diagnostics::IceAreaFractionGrounded::compute_impl(), and pism::Geometry::ensure_consistency().