advective_fluxes bootstrap_impl compute_conductivity compute_velocity define_model_state_impl diagnostics_impl init_impl initialization_message m_bottom_surface m_dx m_dy m_Kstag m_Qstag m_Qstag_average m_R m_rg m_Vstag m_Wnew m_Wstag m_Wtillnew max_timestep_W_cfl max_timestep_W_diff restart_impl Routing subglacial_water_pressure ts_diagnostics_impl update_impl update_W update_Wtill velocity_staggered W_change_due_to_flow water_thickness_staggered write_model_state_impl ~Routing

compute_conductivity() void pism::hydrology::Routing::compute_conductivity ( const array::Staggered &  W, const array::Scalar &  P, const array::Scalar &  bed_elevation, array::Staggered &  result, double &  KW_max  ) const protected Compute the nonlinear conductivity at the center of cell edges. Computes \[ K = K(W, \nabla P, \nabla b) = k W^{\alpha-1} |\nabla R|^{\beta-2} \] on the staggered grid, where \(R = P+\rho_w g b\). We denote \[ \Pi = |\nabla R|^2 \] internally; this is computed on a staggered grid by a Mahaffy-like ([Mahaffy]) scheme. This requires \(R\) to be defined on a box stencil of width 1. Also returns the maximum over all staggered points of \( K W \). Definition at line 436 of file Routing.cc. References pism::array::Array2D< T >::add(), pism::GlobalMax(), pism::k, pism::Component::m_config, m_dx, m_dy, pism::Component::m_grid, m_R, m_rg, pism::array::max(), and pism::array::Array::update_ghosts(). Referenced by pism::hydrology::Distributed::update_impl(), and update_impl().