PISM, A Parallel Ice Sheet Model 2.2.1-cd005eec8 committed by Constantine Khrulev on 2025-03-07
|
◆ compute_3d_horizontal_velocity()
Compute horizontal components of the SIA velocity (in 3D). Recall that \mathbf{U}(z) = -2 \nabla h \int_b^z F(s)P(s)ds + \mathbf{U}_b, which can be written in terms of I(z) defined in compute_I(): \mathbf{U}(z) = -I(z) \nabla h + \mathbf{U}_b.
Definition at line 862 of file SIAFD.cc. References compute_I(), pism::array::Array3D::get_column(), pism::k, pism::Component::m_grid, m_work_3d_0, m_work_3d_1, and pism::array::Array::update_ghosts(). Referenced by update(). |