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_interface_fluxes() void pism::GeometryEvolution::compute_interface_fluxes ( const array::CellType1 &  cell_type, const array::Scalar &  ice_thickness, const array::Vector &  velocity, const array::Staggered &  diffusive_flux, array::Staggered &  output  ) protectedvirtual Combine advective velocity and the diffusive flux on the staggered grid with the ice thickness to compute the total flux through cell interfaces. Uses first-order upwinding to compute the advective flux. Limits the diffusive flux to prevent SIA-driven flow in the ocean and ice-free areas. Reimplemented in pism::RegionalGeometryEvolution. Definition at line 575 of file GeometryEvolution.cc. References pism::ParallelSection::check(), pism::ParallelSection::failed(), pism::mask::icy(), pism::limit_advective_flux(), pism::limit_diffusive_flux(), pism::Component::m_grid, pism::array::max(), and n. Referenced by pism::RegionalGeometryEvolution::compute_interface_fluxes(), and flow_step().