Package esys :: Package escript :: Module mountains :: Class Mountains
[hide private]
[frames] | no frames]

Class Mountains


The Mountains class is defined by the following equations:

(1) w_i,aa+(1/eps)*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 (?) 
    

Instance Methods [hide private]
 
__init__(self, domain, eps=0.01, reduced=False)
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
scalar Data
getTopography(self)
returns the current topography.
vector Data
getVelocity(self)
returns the smoothed/extrapolated velocity
 
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)
Sets a new W and updates the H function.
Method Details [hide private]

__init__(self, domain, eps=0.01, reduced=False)
(Constructor)

 

Sets up the level set method.

Parameters:
  • domain - the domain where the mountains is used
  • eps - the smoothing parameter for (1)
  • reduced - if True topography reduced order is used for reconstruction of internal velocity and topography

setTopography(self, 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(self, 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(self, dt=None)

 

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.