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

◆ residual_redistribution_iteration()

void pism::GeometryEvolution::residual_redistribution_iteration ( const array::Scalar bed_topography,
const array::Scalar sea_level,
array::Scalar1 ice_surface_elevation,
array::Scalar ice_thickness,
array::CellType1 cell_type,
array::Scalar area_specific_volume,
array::Scalar residual,
bool &  done 
)
protected

Perform one iteration of the residual mass redistribution.

Parameters
[in]bed_topographybed elevation
[in]sea_levelsea level elevation
[in,out]ice_surface_elevationsurface elevation; used as temp. storage
[in,out]ice_thicknessice thickness; updated
[in,out]cell_typecell type mask; used as temp. storage
[in,out]area_specific_volumearea specific volume; updated
[in,out]residualice volume that still needs to be distributed; updated
[in,out]doneresult flag: true if this iteration should be the last one

Definition at line 825 of file GeometryEvolution.cc.

References pism::GeometryCalculator::compute(), pism::GeometryCalculator::compute_mask(), pism::array::Array2D< T >::copy_from(), pism::East, pism::GeometryEvolution::Impl::gc, pism::GlobalSum(), pism::array::CellType::ice_free_ocean(), pism::mask::ice_free_ocean(), pism::Component::m_grid, m_impl, pism::North, pism::part_grid_threshold_thickness(), pism::array::Array::set(), pism::South, pism::array::Array2D< T >::star(), pism::array::CellType1::star_int(), pism::GeometryEvolution::Impl::thickness, pism::array::Array::update_ghosts(), and pism::West.

Referenced by update_in_place().