pism::stressbalance Namespace Reference

Stress balance models and related diagnostics. More...

Classes
class	BedSmoother
	PISM bed smoother, plus bed roughness parameterization, based on Schoof (2003). More...
class	Blatter
class	BlatterISMIPHOM
class	BlatterMod
class	BlatterTestCFBC
class	BlatterTestHalfar
class	BlatterTestvanderVeen
class	BlatterTestXY
class	BlatterTestXZ
class	ConstantInColumn
	The trivial Shallow Stress Balance modifier. More...
struct	DeviatoricStresses
class	Inputs
class	PrescribedSliding
struct	PrincipalStrainRates
class	PSB_bfrict
	Computes basal frictional heating. More...
class	PSB_deviatoric_stresses
	Reports the vertically-integrated (2D) deviatoric stresses. More...
class	PSB_flux
	Computes uflux and vflux, components of vertically-integrated horizontal flux of ice. More...
class	PSB_flux_mag
	Computes flux_mag, the magnitude of vertically-integrated horizontal flux of ice. More...
class	PSB_pressure
	Reports the pressure within the ice (3D). More...
class	PSB_strain_rates
	Reports the vertically-integrated (2D) principal strain rates. More...
class	PSB_strainheat
	Reports the volumetric strain heating (3D). More...
class	PSB_tauxz
	Reports the xz component of the shear stress within the ice (3D), according to the SIA formula. More...
class	PSB_tauyz
	Reports the yz component of the shear stress within the ice (3D), according to the SIA formula. More...
class	PSB_uvel
	Computes the x-component of the horizontal ice velocity. More...
class	PSB_velbar
	Computes the vertically-averaged ice velocity. More...
class	PSB_velbar_mag
	Computes velbar_mag, the magnitude of vertically-integrated horizontal velocity of ice and masks out ice-free areas. More...
class	PSB_velbase
	Computes horizontal ice velocity at the base of ice. More...
class	PSB_velbase_mag
	Computes velbase_mag, the magnitude of horizontal velocity of ice at base of ice and masks out ice-free areas. More...
class	PSB_velsurf
	Computes velsurf, the horizontal velocity of ice at ice surface. More...
class	PSB_velsurf_mag
	Computes velsurf_mag, the magnitude of horizontal ice velocity at the surface. More...
class	PSB_vonmises_stress
	tensile von Mises stress More...
class	PSB_vvel
	Computes the y-component of the horizontal ice velocity. More...
class	PSB_wvel
	Computes vertical ice velocity (relative to the geoid). More...
class	PSB_wvel_rel
	Computes vertical velocity of ice, relative to the bed directly below. More...
class	PSB_wvelbase
	Computes wvelbase, the vertical velocity of ice at the base of ice. More...
class	PSB_wvelsurf
	Computes wvelsurf, the vertical velocity of ice at ice surface. More...
class	ShallowStressBalance
	Shallow stress balance (such as the SSA). More...
class	SIAFD
class	SIAFD_diffusivity
	Compute diffusivity of the SIA flow. More...
class	SIAFD_diffusivity_staggered
	Compute diffusivity of the SIA flow (on the staggered grid). More...
class	SIAFD_h_x
	Reports the x-component of the ice surface gradient on the staggered grid as computed by SIAFD. More...
class	SIAFD_h_y
	Reports the y-component of the ice surface gradient on the staggered grid as computed by SIAFD. More...
class	SIAFD_Regional
	A version of the SIA stress balance with tweaks for outlet glacier simulations. More...
class	SIAFD_schoofs_theta
	Computes the multiplier \(\theta\) in Schoof's (2003) theory of the effect of bed roughness on the diffusivity of the SIA. More...
class	SIAFD_thksmooth
	Computes the thickness relative to the smoothed bed elevation in Schoof's (2003) theory of the effect of bed roughness on the SIA. More...
class	SIAFD_topgsmooth
	Computes the smoothed bed elevation from Schoof's (2003) theory of the effect of bed roughness on the SIA. More...
class	SSA
	PISM's SSA solver. More...
class	SSAFD
	PISM's SSA solver: the finite difference implementation. More...
class	SSAFD_nuH
	Reports the nuH (viscosity times thickness) product on the staggered grid. More...
class	SSAFD_residual_mag
	Computes the magnitude of the driving shear stress at the base of ice (diagnostically). More...
class	SSAFD_SNES
class	SSAFD_taud
	Computes the driving shear stress at the base of ice (diagnostically). More...
class	SSAFD_taud_mag
	Computes the magnitude of the driving shear stress at the base of ice (diagnostically). More...
class	SSAFDBase
class	SSAFEM
	PISM's SSA solver: the finite element method implementation written by Jed and David. More...
class	SSAStrengthExtension
	Gives an extension coefficient to maintain ellipticity of SSA where ice is thin. More...
class	SSATestCase
class	SSATestCaseCFBC
class	SSATestCaseConst
class	SSATestCaseExp
class	SSATestCaseI
class	SSATestCaseJ
class	SSATestCasePlug
class	SSB_beta
class	SSB_Modifier
	Shallow stress balance modifier (such as the non-sliding SIA). More...
class	SSB_taub
	Computes the basal shear stress \( \tau_b \). More...
class	SSB_taub_mag
	Computes the magnitude of the basal shear stress (diagnostically). More...
class	SSB_taud
	Computes the gravitational driving stress (diagnostically). More...
class	SSB_taud_mag
	Computes the magnitude of the gravitational driving stress (diagnostically). More...
class	StressBalance
	The class defining PISM's interface to the shallow stress balance code. More...
class	WeertmanSliding
class	ZeroSliding
	Returns zero velocity field, zero friction heating, and zero for D^2. More...

Typedefs
typedef SSA *(*	SSAFactory) (std::shared_ptr< const Grid >)

Enumerations
enum	ISMIPHOMTest { HOM_A , HOM_B , HOM_C , HOM_D }

Functions
static PetscErrorCode	blatter_restriction_hook (DM fine, Mat mrestrict, Vec rscale, Mat inject, DM coarse, void *ctx)
static PetscErrorCode	blatter_coarsening_hook (DM dm_fine, DM dm_coarse, void *ctx)
static void	disable_ew (::SNES snes, double rtol)
static void	enable_ew (::SNES snes)
static double	A_surface (double x, double y, double L)
static double	A_bed (double x, double y, double L)
static double	B_bed (double x, double y, double L)
static double	C_surface (double x, double y, double L)
static double	C_bed (double x, double y, double L)
std::shared_ptr< StressBalance >	create (const std::string &model, std::shared_ptr< const Grid > grid, bool regional)
static void	copy_staggered_vec (const array::Staggered &input, array::Staggered &output)
PetscErrorCode	SSAFDSNESConvergenceTest (SNES snes, PetscInt it, PetscReal xnorm, PetscReal gnorm, PetscReal f, SNESConvergedReason *reason, void *ctx)
static int	weight (bool margin_bc, int M_ij, int M_n, double h_ij, double h_n, int N_ij, int N_n)
static double	diff_uphill (double L, double C, double R)
static double	diff_centered (double L, double, double R)
static bool	is_marginal (int i, int j, const array::CellType1 &cell_type, bool ssa_dirichlet_bc)
	Checks if a cell is near or at the ice front.
static double	H_exact (double V0, double H0, double C, double x)
static double	u_exact (double V0, double H0, double C, double x)
std::shared_ptr< Grid >	ssa_test_cfbc_grid (std::shared_ptr< Context > ctx, int Mx, int My)
std::shared_ptr< Grid >	ssa_test_const_grid (std::shared_ptr< Context > ctx, int Mx, int My)
std::shared_ptr< Grid >	ssa_test_i_grid (std::shared_ptr< Context > ctx, int Mx, int My)
std::shared_ptr< Grid >	ssa_test_j_grid (std::shared_ptr< Context > ctx, int Mx, int My)
std::shared_ptr< Grid >	ssa_test_linear_grid (std::shared_ptr< Context > ctx, int Mx, int My)
std::shared_ptr< Grid >	ssa_test_plug_grid (std::shared_ptr< Context > ctx, int Mx, int My)
static double	D2 (double u_x, double u_y, double u_z, double v_x, double v_y, double v_z)
void	compute_2D_principal_strain_rates (const array::Vector1 &V, const array::CellType1 &mask, array::Array2D< PrincipalStrainRates > &result)
	Compute eigenvalues of the horizontal, vertically-integrated strain rate tensor.
void	compute_2D_stresses (const rheology::FlowLaw &flow_law, const array::Vector1 &velocity, const array::Scalar &hardness, const array::CellType1 &cell_type, array::Array2D< DeviatoricStresses > &result)
	Compute 2D deviatoric stresses.
static void	zero_above_ice (const array::Array3D &F, const array::Scalar &H, array::Array3D &result)

Variables
static const Vector2d	u_exterior = {0.0, 0.0}
const double	m_schoof = 10
const double	L_schoof = 40e3
const double	aspect_schoof = 0.05
const double	H0_schoof = aspect_schoof * L_schoof
const double	B_schoof = 3.7e8

Detailed Description

Stress balance models and related diagnostics.