Tridiagonal linear system for conservation of energy in vertical column of ice enthalpy. More...

#include <enthSystem.hh>

Inheritance diagram for pism::energy::enthSystemCtx:

Public Member Functions
	enthSystemCtx (const std::vector< double > &storage_grid, const std::string &prefix, double dx, double dy, double dt, const Config &config, const array::Array3D &Enth3, const array::Array3D &u3, const array::Array3D &v3, const array::Array3D &w3, const array::Array3D &strain_heating3, EnthalpyConverter::Ptr EC)

	~enthSystemCtx ()=default

void	init (int i, int j, bool ismarginal, double ice_thickness)

double	k_from_T (double T) const

void	set_surface_heat_flux (double heat_flux)
	Set the top surface heat flux into the ice.

void	set_surface_neumann_bc (double dE)
	Set enthalpy flux at the surface.

void	set_surface_dirichlet_bc (double E_surface)

void	set_basal_dirichlet_bc (double E_basal)
	Set coefficients in discrete equation for \(E = Y\) at base of ice.

void	set_basal_heat_flux (double heat_flux)
	Set coefficients in discrete equation for Neumann condition at base of ice.

void	set_basal_neumann_bc (double dE)
	Set enthalpy flux at the base.

virtual void	save_system (std::ostream &output, unsigned int M) const

void	solve (std::vector< double > &result)
	Solve the tridiagonal system, in a single column, which determines the new values of the ice enthalpy.

double	lambda () const

double	Enth (size_t i) const

double	Enth_s (size_t i) const

Public Member Functions inherited from pism::columnSystemCtx
	columnSystemCtx (const std::vector< double > &storage_grid, const std::string &prefix, double dx, double dy, double dt, const array::Array3D &u3, const array::Array3D &v3, const array::Array3D &w3)
	A column system is a kind of a tridiagonal system.

	~columnSystemCtx ()

void	save_to_file (const std::vector< double > &x)
	Write system matrix, right-hand-side, and (provided) solution into Python script. Constructs file name from m_prefix.

void	save_to_file (const std::string &filename, const std::vector< double > &x)

unsigned int	ks () const

double	dz () const

const std::vector< double > &	z () const

void	fine_to_coarse (const std::vector< double > &input, int i, int j, array::Array3D &output) const

Protected Member Functions
void	compute_enthalpy_CTS ()
	Compute the CTS value of enthalpy in an ice column.

double	compute_lambda ()
	Compute the lambda for BOMBPROOF.

void	assemble_R ()
	Assemble the R array. The R value switches at the CTS.

void	checkReadyToSolve () const

Protected Member Functions inherited from pism::columnSystemCtx
void	init_column (int i, int j, double ice_thickness)

void	reportColumnZeroPivotErrorMFile (unsigned int M)
	Write system matrix and right-hand-side into an Python script. The file name contains ZERO_PIVOT_ERROR.

void	init_fine_grid (const std::vector< double > &storage_grid)

void	coarse_to_fine (const array::Array3D &input, int i, int j, double *output) const

Protected Attributes
std::vector< double >	m_Enth

std::vector< double >	m_Enth_s

std::vector< double >	m_E_ij

std::vector< double >	m_E_n

std::vector< double >	m_E_e

std::vector< double >	m_E_s

std::vector< double >	m_E_w

std::vector< double >	m_strain_heating
	strain heating in the ice column

std::vector< double >	m_R
	values of \( k \Delta t / (\rho c \Delta x^2) \)

double	m_ice_density

double	m_ice_c

double	m_ice_k

double	m_p_air

double	m_nu

double	m_R_cold

double	m_R_temp

double	m_R_factor

double	m_ice_thickness

double	m_lambda

double	m_D0
	implicit FD method parameter

double	m_U0

double	m_B0

double	m_L_ks

double	m_D_ks

double	m_U_ks

double	m_B_ks

bool	m_marginal

bool	m_k_depends_on_T

bool	m_margin_exclude_horizontal_advection

bool	m_margin_exclude_vertical_advection

bool	m_margin_exclude_strain_heat

const array::Array3D &	m_Enth3

const array::Array3D &	m_strain_heating3

EnthalpyConverter::Ptr	m_EC

Protected Attributes inherited from pism::columnSystemCtx
TridiagonalSystem *	m_solver

ColumnInterpolation *	m_interp

unsigned int	m_ks
	current system size; corresponds to the highest vertical level within the ice

int	m_i
	current column indexes

int	m_j

double	m_dx

double	m_dy

double	m_dz

double	m_dt

std::vector< double >	m_u
	u-component of the ice velocity

std::vector< double >	m_v
	v-component of the ice velocity

std::vector< double >	m_w
	w-component of the ice velocity

std::vector< double >	m_z
	levels of the fine vertical grid

const array::Array3D &	m_u3
	pointers to 3D velocity components

const array::Array3D &	m_v3

const array::Array3D &	m_w3

Detailed Description

Tridiagonal linear system for conservation of energy in vertical column of ice enthalpy.

See the page documenting bombproofenth. The top of the ice has a Dirichlet condition.

Definition at line 38 of file enthSystem.hh.

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

src/energy/enthSystem.hh
src/energy/enthSystem.cc

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description