PISM, A Parallel Ice Sheet Model  stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800
Public Member Functions | Private Attributes | List of all members
pism::energy::BedrockColumn Class Reference

#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
 

Detailed Description

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.


The documentation for this class was generated from the following files: