esys.escriptcore.mountains Package

Classes

class esys.escriptcore.mountains.Mountains(domain, eps=0.01)

Bases: object

The Mountains class is defined by the following equations:

  1. 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.
  2. 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 (?)
getDomain()

Returns the domain.

getSafeTimeStepSize()

Returns the time step value.

Return type:float
getSolverOptionsForSmooting()

returns the solver options for the smoothing/extrapolation

getSolverOptionsForUpdate()

returns the solver options for the topograthy update

getTopography()

returns the current topography. :rtype: scalar Data

getVelocity()

returns the smoothed/extrapolated velocity :rtype: vector Data

setTopography(H=None)

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

Parameters:H (scalar) – the topography. If None zero is used.
setVelocity(v=None)

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

Parameters:v (vector) – velocity field. If None zero is used.
update(dt=None, allow_substeps=True)

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.
class esys.escriptcore.mountains.SubSteppingException

Bases: exceptions.Exception

Thrown if the L{Mountains} class uses substepping.

args
message

Functions

Others

  • __builtins__
  • __copyright__
  • __doc__
  • __file__
  • __license__
  • __name__
  • __package__
  • __url__