PBPointwiseIsostasy Class Reference

Pointwide isostasy bed deformation model. More...

#include <PISMBedDef.hh>

Inheritance diagram for PBPointwiseIsostasy:

List of all members.

Public Member Functions
	PBPointwiseIsostasy (IceGrid &g, const NCConfigVariable &conf)
virtual	~PBPointwiseIsostasy ()
virtual PetscErrorCode	init (PISMVars &vars)
virtual PetscErrorCode	update (PetscReal t_years, PetscReal dt_years)
	Updates the pointwise isostasy model.
Protected Member Functions
PetscErrorCode	allocate ()
Protected Attributes
IceModelVec2S	thk_last
	last ice thickness

---

Detailed Description

Pointwide isostasy bed deformation model.

Definition at line 53 of file PISMBedDef.hh.

---

Constructor & Destructor Documentation

PBPointwiseIsostasy	(	IceGrid &	g,
		const NCConfigVariable &	conf
	)

Definition at line 21 of file PBPointwiseIsostasy.cc.

References allocate(), and PISMComponent::grid.

virtual ~PBPointwiseIsostasy

(

)

[inline, virtual]

Definition at line 56 of file PISMBedDef.hh.

---

Member Function Documentation

PetscErrorCode allocate

(

)

[protected]

Definition at line 33 of file PBPointwiseIsostasy.cc.

References PISMComponent::grid, and thk_last.

Referenced by PBPointwiseIsostasy().

PetscErrorCode init

(

PISMVars &

vars

)

[virtual]

Reimplemented from PISMBedDef.

Definition at line 42 of file PBPointwiseIsostasy.cc.

References PISMComponent::grid, PISMBedDef::thk, thk_last, PISMBedDef::topg, and PISMBedDef::topg_last.

PetscErrorCode update	(	PetscReal	t_years,
		PetscReal	dt_years
	)		[virtual]

Updates the pointwise isostasy model.

Our goal: topg = topg_last - f*(thk - thk_last)

Step 1: topg = topg_last - f*thk

Step 2: topg = topg + f*thk_last = (topg_last - f*thk) + f*thk_last = topg_last - f*(thk - thk_last)

This code is written this way to avoid allocating temp. storage for (thk - thk_last).

Finally, we need to update bed uplift, topg_last and thk_last.

Implements PISMBedDef.

Definition at line 57 of file PBPointwiseIsostasy.cc.

References PISMBedDef::compute_uplift(), PISMComponent::config, PISMComponent_TS::dt, e, PISMComponent_TS::t, PISMBedDef::t_beddef_last, PISMBedDef::thk, thk_last, PISMBedDef::topg, and PISMBedDef::topg_last.

---

Member Data Documentation

IceModelVec2S thk_last [protected]

last ice thickness

Definition at line 61 of file PISMBedDef.hh.

Referenced by allocate(), init(), and update().

---

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

• src/earth/PISMBedDef.hh
• src/earth/PBPointwiseIsostasy.cc