This Section presents the various solvers developed to solve glaciological applications using the finite element code Elmer. A solver is a Fortran subroutine that solves a differential equation on a given domain for all the nodal variables belonging in this domain. Very often, this results in solving a linear system. In certain cases, when the equation is only local, this can be done more simply by updating and iterating over all the nodes of the domain.
Solvers in the Elmer/Ice package can be called using the object file name ElmerIceSolvers:
Procedure = File "ElmerIceSolvers" "NameSolver"
NameSolver is the name of the solver you want to execute.
The source code of the solvers of the Elmer/Ice package can be obtained from the Elmer svn in
Non-distributed solvers must be linked to
libelmersolver.so, for which a wrapper script
elmerf90 MySolver.f90 -o MySolver
In the SIF file (Solver section), the solver has to point on this object file
Procedure = File "MySolver" "NameProcedureInMySolver"
NameProcedureInMySolver is the name of the fortran procedure in the file