bootstrap bootstrap_impl compute_basal_melt_rate compute_overburden_pressure compute_surface_input_rate define_model_state_impl diagnostics_impl enforce_bounds flux Hydrology init init_impl initialization_message m_basal_melt_rate m_conservation_error_change m_flow_change m_flow_change_incremental m_grounded_margin_change m_grounding_line_change m_input_change m_no_model_mask_change m_Pover m_Q m_surface_input_rate m_total_change m_W m_Wtill mass_change mass_change_at_domain_boundary mass_change_at_grounded_margin mass_change_at_grounding_line mass_change_due_to_conservation_error mass_change_due_to_input mass_change_due_to_lateral_flow overburden_pressure restart restart_impl subglacial_water_thickness surface_input_rate till_water_thickness update update_impl write_model_state_impl ~Hydrology

enforce_bounds() void pism::hydrology::Hydrology::enforce_bounds ( const array::CellType &  cell_type, const array::Scalar *  no_model_mask, double  max_thickness, double  ocean_water_thickness, array::Scalar &  water_thickness, array::Scalar &  grounded_margin_change, array::Scalar &  grounding_line_change, array::Scalar &  conservation_error_change, array::Scalar &  no_model_mask_change  ) protected Correct the new water thickness based on boundary requirements. At ice free locations and ocean locations we require that water thicknesses (i.e. both the transportable water thickness \(W\) and the till water thickness \(W_{till}\)) be zero at the end of each time step. Also we require that any negative water thicknesses be set to zero (i.e. we do projection to enforce lower bound). This method should be called once for each thickness field which needs to be processed. This method alters the field water_thickness in-place. Parameters [in] cell_type cell type mask [in] no_model_mask (optional) mask of zeros and ones, zero within the modeling domain, one outside [in] max_thickness maximum allowed water thickness (use a zero or a negative value to disable) [in,out] water_thickness adjusted water thickness (till storage or the transport system) [in,out] grounded_margin_change change in water thickness at the grounded margin [in,out] grounding_line_change change in water thickness at the grounding line [in,out] conservation_error_change change in water thickness due to mass conservation errors [in,out] no_model_mask_change change in water thickness outside the modeling domain (regional models) Definition at line 669 of file Hydrology.cc. References pism::array::Array2D< T >::add(), pism::array::CellType::ice_free_land(), pism::array::CellType::ice_free_ocean(), pism::Component::m_config, pism::Component::m_grid, and pism::array::CellType::ocean(). Referenced by pism::hydrology::Distributed::update_impl(), pism::hydrology::NullTransport::update_impl(), and pism::hydrology::Routing::update_impl().