Class implementing the bed deformation model described in [BLKfastearth]. More...

#include <deformation.hh>

Public Member Functions
	BedDeformLC ()
	~BedDeformLC ()
PetscErrorCode	settings (const NCConfigVariable &config, PetscTruth myinclude_elastic, PetscInt myMx, PetscInt myMy, PetscScalar mydx, PetscScalar mydy, PetscInt myZ, PetscScalar myicerho, Vec myHstart, Vec mybedstart, Vec myuplift, Vec myH, Vec *mybed)
PetscErrorCode	alloc ()
PetscErrorCode	uplift_init ()
PetscErrorCode	step (const PetscScalar dtyear, const PetscScalar yearFromStart)
Protected Attributes
PetscTruth	include_elastic
PetscInt	Mx
PetscInt	My
PetscScalar	dx
PetscScalar	dy
PetscInt	Z
PetscScalar	icerho
PetscScalar	rho
PetscScalar	eta
PetscScalar	D
Private Attributes
PetscScalar	standard_gravity
PetscTruth	settingsDone
PetscTruth	allocDone
PetscInt	Nx
PetscInt	Ny
PetscInt	Nxge
PetscInt	Nyge
PetscInt	i0_plate
PetscInt	j0_plate
PetscScalar	Lx
PetscScalar	Ly
PetscScalar	fatLx
PetscScalar	fatLy
PetscScalar *	cx
PetscScalar *	cy
Vec *	H
Vec *	bed
Vec *	Hstart
Vec *	bedstart
Vec *	uplift
Vec	Hdiff
Vec	dbedElastic
Vec	platefat
Vec	plateoffset
Vec	vleft
Vec	vright
Vec	lrmE
fftw_complex *	bdin
fftw_complex *	bdout
fftw_plan	bdplanfor
fftw_plan	bdplanback

Detailed Description

Class implementing the bed deformation model described in [BLKfastearth].

This class implements the [LingleClark] bed deformation model by a Fourier spectral collocation method, as described in [BLKfastearth]. (The former reference is where the continuum model arose, and a flow-line application is given. The latter reference describes a new, fast method and gives verification results. See also [BLK2006earth] if more technical detail and/or Matlab programs are desired.)

Both a viscous half-space model (with elastic lithosphere) and a spherical elastic model are computed. They are superposed because the underlying earth model is linear.

The class assumes that the supplied Petsc Vecs are *sequential*. It is expected to be run only on processor zero (or possibly by each processor once each processor owns the entire 2D gridded ice thicknesses and bed elevations.)

This class SHOULD! include the scatter structures necessary to make this work in parallel.

A test program for this class is pism/src/verif/tryLCbd.cc.

Called by IceModel in pism/src/iMbeddef.cc.

Definition at line 51 of file deformation.hh.

Constructor & Destructor Documentation

BedDeformLC ( )

Definition at line 32 of file deformation.cc.

References allocDone, and settingsDone.

~BedDeformLC ( )

Definition at line 38 of file deformation.cc.

References allocDone, bdin, bdout, bdplanback, bdplanfor, cx, cy, dbedElastic, Hdiff, lrmE, platefat, plateoffset, vleft, and vright.

Member Function Documentation

PetscErrorCode alloc ( )

Definition at line 104 of file deformation.cc.

References allocDone, bdin, bdout, bdplanback, bdplanfor, cx, cy, dbedElastic, dblquad_cubature(), dx, ge_params::dx, dy, ge_params::dy, e, fatLx, fatLy, ge_integrand(), H, Hdiff, include_elastic, lrmE, Nx, Nxge, Ny, Nyge, ge_params::p, pi, platefat, plateoffset, ge_params::q, settingsDone, vleft, and vright.

Referenced by PBLingleClark::allocate(), and main().

PetscErrorCode settings	(	const NCConfigVariable &	config,
		PetscTruth	myinclude_elastic,
		PetscInt	myMx,
		PetscInt	myMy,
		PetscScalar	mydx,
		PetscScalar	mydy,
		PetscInt	myZ,
		PetscScalar	myicerho,
		Vec *	myHstart,
		Vec *	mybedstart,
		Vec *	myuplift,
		Vec *	myH,
		Vec *	mybed
	)