PISM, A Parallel Ice Sheet Model  stable v2.1.1 committed by Constantine Khrulev on 2024-12-04 13:36:58 -0900

◆ flow_step()

void pism::GeometryEvolution::flow_step ( const Geometry geometry,
double  dt,
const array::Vector advective_velocity,
const array::Staggered diffusive_flux,
const array::Scalar thickness_bc_mask 
)
Parameters
[in]geometryice geometry
[in]dttime step, seconds
[in]advective_velocityadvective (SSA) velocity
[in]diffusive_fluxdiffusive (SIA) flux
[in]velocity_bc_valuesadvective velocity Dirichlet B.C. values
[in]thickness_bc_maskice thickness Dirichlet B.C. mask

Results are stored in internal fields accessible using getters.

Definition at line 250 of file GeometryEvolution.cc.

References pism::array::Array2D< T >::add(), pism::GeometryEvolution::Impl::area_specific_volume, pism::Geometry::bed_elevation, pism::GeometryEvolution::Impl::bed_elevation, pism::Profiling::begin(), pism::GeometryEvolution::Impl::cell_type, pism::GeometryCalculator::compute(), compute_flux_divergence(), compute_interface_fluxes(), pism::GeometryEvolution::Impl::conservation_error, pism::array::Array2D< T >::copy_from(), pism::array::Staggered::copy_from(), pism::Profiling::end(), ensure_nonnegativity(), pism::GeometryEvolution::Impl::flux_divergence, pism::GeometryEvolution::Impl::flux_staggered, pism::GeometryEvolution::Impl::gc, pism::Geometry::ice_area_specific_volume, pism::GeometryEvolution::Impl::ice_area_specific_volume_change, pism::Geometry::ice_thickness, pism::GeometryEvolution::Impl::ice_thickness, pism::GeometryEvolution::Impl::input_velocity, pism::Component::m_grid, m_impl, pism::make_nonnegative_preserving(), pism::Component::profiling(), pism::GeometryEvolution::Impl::sea_level, pism::Geometry::sea_level_elevation, pism::GeometryEvolution::Impl::surface_elevation, pism::GeometryEvolution::Impl::thickness_change, pism::array::Array::update_ghosts(), and update_in_place().