PISM, A Parallel Ice Sheet Model  stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800

◆ update()

void pism::Isochrones::update ( double  t,
double  dt,
const array::Array3D u,
const array::Array3D v,
const array::Scalar ice_thickness,
const array::Scalar top_surface_mass_balance,
const array::Scalar bottom_surface_mass_balance 
)

Update layer thicknesses.

Parameters
[in]ttime step start, seconds
[in]dttime step length, seconds
[in]ux-component of the ice velocity, m/s
[in]vy-component of the ice velocity, m/s
[in]icethickness, m
[in]top_surface_mass_balancetotal top surface mass balance over the time step, meters
[in]bottom_surface_mass_balancetotal bottom surface mass balance over the time step, meters

For mass balance inputs, positive corresponds to mass gain.

Definition at line 557 of file Isochrones.cc.

References pism::stencils::Star< T >::c, pism::Time::date(), pism::details::deposition_times(), pism::stencils::Star< T >::e, pism::array::Array3D::interpolate(), pism::k, pism::Component::m_config, pism::Component::m_grid, m_layer_thickness, pism::Component::m_log, m_tmp, m_top_layer_index, pism::array::max(), pism::stencils::Star< T >::n, S(), pism::stencils::Star< T >::s, pism::Component::time(), top_layer_deposition_time(), and pism::stencils::Star< T >::w.