bed_smoother compute_3d_horizontal_velocity compute_diffusive_flux compute_diffusivity compute_I compute_surface_gradient diagnostics_impl diffusivity init interglacial m_bed_smoother m_D m_delta_0 m_delta_1 m_e_factor m_e_factor_interglacial m_eemian_end m_eemian_start m_event_sia m_h_x m_h_y m_holocene_start m_seconds_per_year m_stencil_width m_work_2d_0 m_work_2d_1 m_work_3d_0 m_work_3d_1 SIAFD surface_gradient_eta surface_gradient_haseloff surface_gradient_mahaffy surface_gradient_x surface_gradient_y update ~SIAFD

compute_3d_horizontal_velocity() void pism::stressbalance::SIAFD::compute_3d_horizontal_velocity ( const Geometry &  geometry, const array::Staggered &  h_x, const array::Staggered &  h_y, const array::Vector &  sliding_velocity, array::Array3D &  u_out, array::Array3D &  v_out  ) protectedvirtual 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.$$ Note This is one of the places where "hybridization" is done. Parameters [in] h_x the X-component of the surface gradient, on the staggered grid [in] h_y the Y-component of the surface gradient, on the staggered grid [in] sliding_velocity the thickness-advective velocity from the underlying stress balance module [out] u_out the X-component of the resulting horizontal velocity field [out] v_out the Y-component of the resulting horizontal velocity field 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().