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

◆ set_surface_neumann_bc()

void pism::energy::enthSystemCtx::set_surface_neumann_bc ( double  dE)

Set enthalpy flux at the surface.

This method should probably be used for debugging only. Its purpose is to allow setting the enthalpy flux even if K == 0, i.e. in a "pure advection" setup.

Definition at line 222 of file enthSystem.cc.

References m_B_ks, m_D_ks, pism::columnSystemCtx::m_dt, pism::columnSystemCtx::m_dx, pism::columnSystemCtx::m_dy, pism::columnSystemCtx::m_dz, m_E_e, m_E_n, m_E_s, m_E_w, m_Enth, m_ice_density, pism::columnSystemCtx::m_ks, m_L_ks, m_lambda, m_margin_exclude_horizontal_advection, m_margin_exclude_strain_heat, m_marginal, m_nu, m_R, m_strain_heating, pism::columnSystemCtx::m_u, m_U_ks, pism::columnSystemCtx::m_v, pism::columnSystemCtx::m_w, and pism::energy::upwind().

Referenced by set_surface_heat_flux().