Class Mountains

The Mountains class is defined by the following equations:

eps*w_i,aa+w_i,33=0 where 0<=eps<<1 and a=1,2 and w_i is the extension of the surface velocity where w_i(x_3=1)=v_i.
Integration of topography PDE using Taylor-Galerkin upwinding to stabilize the advection terms H^(t+dt)=H^t+dt*w_3+w_hat*dt*[(div(w_hat*H^t)+w_3)+(dt/2)+H^t], where w_hat=w*[1,1,0], dt<0.5*d/max(w_i), d is a characteristic element size; H(x_3=1)=lambda (?)

Instance Methods

__init__(self, domain, eps=0.01)
Sets up the level set method.

getDomain(self)
Returns the domain.

float

getSafeTimeStepSize(self)
Returns the time step value.

getSolverOptionsForSmooting(self)
returns the solver options for the smoothing/extrapolation

getSolverOptionsForUpdate(self)
returns the solver options for the topograthy update

getTopography(self)
returns the current topography.

getVelocity(self)
returns the smoothed/extrapolated velocity :rtype: vector Data

setTopography(self, H=None)
set the topography to H where H defines the vertical displacement.

setVelocity(self, v=None)
set a new velocity.

update(self, dt=None, allow_substeps=True)
Sets a new W and updates the H function.

Method Details

Sets up the level set method.

Parameters:

returns the current topography. :rtype: scalar Data

set the topography to H where H defines the vertical displacement. H is defined for the entire domain.

Parameters:

set a new velocity. v is define on the entire domain but only the surface values are used.

Parameters:

Sets a new W and updates the H function.

Parameters:

dt (positve float which is less or equal than the safe time step size.) - time step forward. If None the save time step size is used.