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_I() void pism::stressbalance::SIAFD::compute_I ( const Geometry &  geometry ) protectedvirtual Compute I. This computes \[ I(z) = \int_b^z\delta(s)ds.\] Uses the trapezoidal rule to approximate the integral. See compute_diffusive_flux() for the definition of \(\delta\). The result is stored in work_3d[0,1] and is used to compute the SIA component of the 3D-distributed horizontal ice velocity. Definition at line 781 of file SIAFD.cc. References pism::Geometry::cell_type, pism::ParallelSection::check(), pism::ParallelSection::failed(), pism::array::Array3D::get_column(), I, pism::Geometry::ice_surface_elevation, pism::Geometry::ice_thickness, pism::k, m_bed_smoother, m_delta_0, m_delta_1, pism::Component::m_grid, m_work_2d_0, m_work_3d_0, m_work_3d_1, pism::stressbalance::BedSmoother::smoothed_thk(), and pism::array::Array::stencil_width(). Referenced by compute_3d_horizontal_velocity().