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

◆ solve() [1/2]

void pism::energy::BedrockColumn::solve ( double  dt,
double  Q_bottom,
double  T_top,
const double *  T_old,
double *  T_new 
)

Advance the heat equation in time.

Parameters
[in]dttime step length
[in]Q_bottomheat flux into the column through the bottom surface
[in]T_toptemperature at the top surface
[in]T_oldcurrent temperature in the column
[out]T_newoutput

Note: T_old and T_new may point to the same location.

Definition at line 53 of file BedrockColumn.cc.

References pism::TridiagonalSystem::D(), G, pism::k, pism::TridiagonalSystem::L(), m_D, m_dz, m_k, m_M, m_system, pism::TridiagonalSystem::RHS(), pism::TridiagonalSystem::solve(), and pism::TridiagonalSystem::U().

Referenced by solve().