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

compute_surface_and_basal_mass_balance() void pism::GeometryEvolution::compute_surface_and_basal_mass_balance ( double  dt, const array::Scalar &  thickness_bc_mask, const array::Scalar &  ice_thickness, const array::CellType &  cell_type, const array::Scalar &  smb_flux, const array::Scalar &  basal_melt_rate, array::Scalar &  effective_SMB, array::Scalar &  effective_BMB  ) protectedvirtual Compute effective surface and basal mass balance. Parameters [in] dt time step, seconds [in] thickness_bc_mask mask specifying ice thickness Dirichlet B.C. locations [in] ice_thickness ice thickness, m [in] thickness_change thickness change due to flow, m [in] cell_type cell type mask [in] smb_rate top surface mass balance rate, kg m-2 s-1 [in] basal_melt_rate basal melt rate, m s-1 [out] effective_smb effective top surface mass balance, m [out] effective_bmb effective basal mass balance, m This computation is purely local. Positive outputs correspond to mass gain. Reimplemented in pism::RegionalGeometryEvolution. Definition at line 1018 of file GeometryEvolution.cc. References pism::array::Scalar::as_int(), pism::ParallelSection::check(), pism::effective_change(), pism::ParallelSection::failed(), pism::GeometryEvolution::Impl::ice_density, pism::array::CellType::ice_free_ocean(), pism::Component::m_grid, m_impl, and pism::GeometryEvolution::Impl::use_bmr. Referenced by pism::RegionalGeometryEvolution::compute_surface_and_basal_mass_balance(), and source_term_step().