PISM, A Parallel Ice Sheet Model
stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800
|
◆ solve()
First-order upwind scheme with implicit in the vertical: one column solve. The PDE being solved is \[ \frac{\partial \tau}{\partial t} + \frac{\partial}{\partial x}\left(u \tau\right) + \frac{\partial}{\partial y}\left(v \tau\right) + \frac{\partial}{\partial z}\left(w \tau\right) = 1. \] Definition at line 68 of file AgeColumnSystem.cc. References pism::RuntimeError::add_context(), pism::k, m_A, m_A_e, m_A_n, m_A_s, m_A_w, pism::columnSystemCtx::m_dt, pism::columnSystemCtx::m_dx, pism::columnSystemCtx::m_dy, pism::columnSystemCtx::m_i, pism::columnSystemCtx::m_j, pism::columnSystemCtx::m_ks, m_nu, pism::columnSystemCtx::m_solver, pism::columnSystemCtx::m_u, pism::columnSystemCtx::m_v, pism::columnSystemCtx::m_w, pism::columnSystemCtx::reportColumnZeroPivotErrorMFile(), and S(). Referenced by pism::AgeModel::update(). |