PISM, A Parallel Ice Sheet Model
stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800
|
#include <BedrockColumn.hh>
Public Member Functions | |
BedrockColumn (const std::string &prefix, const Config &config, double dz, unsigned int M) | |
~BedrockColumn ()=default | |
void | solve (double dt, double Q_bottom, double T_top, const double *T_old, double *result) |
void | solve (double dt, double Q_bottom, double T_top, const std::vector< double > &T_old, std::vector< double > &result) |
Private Attributes | |
double | m_D |
double | m_k |
double | m_dz |
unsigned int | m_M |
TridiagonalSystem | m_system |
Tridiagonal linear system for conservation of energy in vertical columns of the bedrock thermal layer.
The top surface uses the Dirichlet boundary condition, the bottom surface a heat flux (Neumann) BC.
The implementation uses a second-order discretization in space and the backward-Euler (first-order, fully implicit) time-discretization.
Definition at line 39 of file BedrockColumn.hh.