apply_flux_divergence apply_mass_fluxes area_specific_volume_change_due_to_flow bottom_surface_mass_balance compute_flux_divergence compute_interface_fluxes compute_surface_and_basal_mass_balance conservation_error diagnostics_impl ensure_nonnegativity flow_step flux_divergence flux_staggered GeometryEvolution init init_impl m_impl reset residual_redistribution_iteration set_no_model_mask set_no_model_mask_impl source_term_step thickness_change_due_to_flow top_surface_mass_balance update_in_place ~GeometryEvolution

ensure_nonnegativity() void pism::GeometryEvolution::ensure_nonnegativity ( const array::Scalar &  ice_thickness, const array::Scalar &  area_specific_volume, array::Scalar &  thickness_change, array::Scalar &  area_specific_volume_change, array::Scalar &  conservation_error  ) protectedvirtual Correct thickness_change and area_specific_volume_change so that applying them will not result in negative ice_thickness and area_specific_volume. Compute the the amount of ice that is added to preserve non-negativity of ice thickness. Parameters [in] ice_thickness ice thickness (m) [in] area_specific_volume area-specific volume (m3/m2) [in,out] thickness_change "proposed" thickness change (m) [in,out] area_specific_volume_change "proposed" area-specific volume change (m3/m2) [in,out] conservation_error computed conservation error (m) This computation is purely local. Definition at line 955 of file GeometryEvolution.cc. References pism::ParallelSection::check(), conservation_error(), pism::ParallelSection::failed(), and pism::Component::m_grid. Referenced by flow_step().