PBPointwiseIsostasy Class Reference

Pointwide isostasy bed deformation model. More...

#include <PISMBedDef.hh>

Inheritance diagram for PBPointwiseIsostasy:

[legend]

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 48 of file PISMBedDef.hh.

---

Constructor & Destructor Documentation

PBPointwiseIsostasy	(	IceGrid &	g,
		const NCConfigVariable &	conf
	)

Definition at line 21 of file PBPointwiseIsostasy.cc.

References allocate(), IceGrid.com, and PISMComponent.grid.

virtual ~PBPointwiseIsostasy

(

)

[virtual]

Definition at line 51 of file PISMBedDef.hh.

---

Member Function Documentation

PetscErrorCode allocate

(

)

[protected]

Definition at line 33 of file PBPointwiseIsostasy.cc.

References IceModelVec2S.create(), 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 IceGrid.com, IceModelVec.copy_to(), PISMComponent.grid, PISMBedDef.thk, thk_last, and verbPrintf().

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 56 of file PBPointwiseIsostasy.cc.

References IceModelVec.add(), PISMBedDef.compute_uplift(), PISMComponent.config, IceModelVec.copy_to(), PISMComponent.dt, e, NCConfigVariable.get(), PISMComponent.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 56 of file PISMBedDef.hh.

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

---

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

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