HR: 17:15h
AN: NG34A-06    [Abstracts]
TI: The Programming Language Python In Earth System Simulations
AU: * Gross, L
EM: gross@esscc.uq.edu.au
AF: Earth Systems Science Computational Center, he University of Queensland, Brisbane, QLD 4071 Australia
AU: Imranullah, A
EM: imran@esscc.uq.edu.au
AF: Earth Systems Science Computational Center, he University of Queensland, Brisbane, QLD 4071 Australia
AU: Mora, P
EM: morap@esscc.uq.edu.au
AF: Earth Systems Science Computational Center, he University of Queensland, Brisbane, QLD 4071 Australia
AU: Saez, E
EM: saez@esscc.uq.edu.au
AF: Earth Systems Science Computational Center, he University of Queensland, Brisbane, QLD 4071 Australia
AU: Smillie, J
EM: jgs@esscc.uq.edu.au
AF: Earth Systems Science Computational Center, he University of Queensland, Brisbane, QLD 4071 Australia
AU: Wang, C
EM: cynthaw@esscc.uq.edu.au
AF: Earth Systems Science Computational Center, he University of Queensland, Brisbane, QLD 4071 Australia
AB: Mathematical models in earth sciences base on the solution of systems of coupled, non-linear, time-dependent partial differential equations (PDEs). The spatial and time-scale vary from a planetary scale and million years for convection problems to 100km and 10 years for fault systems simulations. Various techniques are in use to deal with the time dependency (e.g. Crank-Nicholson), with the non-linearity (e.g. Newton-Raphson) and weakly coupled equations (e.g. non-linear Gauss-Seidel). Besides these high-level solution algorithms discretization methods (e.g. finite element method (FEM), boundary element method (BEM)) are used to deal with spatial derivatives. Typically, large-scale, three dimensional meshes are required to resolve geometrical complexity (e.g. in the case of fault systems) or features in the solution (e.g. in mantel convection simulations). The modelling environment escript allows the rapid implementation of new physics as required for the development of simulation codes in earth sciences. Its main object is to provide a programming language, where the user can define new models and rapidly develop high-level solution algorithms. The current implementation is linked with the finite element package finley as a PDE solver. However, the design is open and other discretization technologies such as finite differences and boundary element methods could be included. escript is implemented as an extension of the interactive programming environment python (see www.python.org). Key concepts introduced are Data objects, which are holding values on nodes or elements of the finite element mesh, and linearPDE objects, which are defining linear partial differential equations to be solved by the underlying discretization technology. In this paper we will show the basic concepts of escript and will show how escript is used to implement a simulation code for interacting fault systems. We will show some results of large-scale, parallel simulations on an SGI Altix system. Acknowledgements: Project work is supported by Australian Commonwealth Government through the Australian Computational Earth Systems Simulator Major National Research Facility, Queensland State Government Smart State Research Facility Fund, The University of Queensland and SGI.
UR: http://www.esscc.uq.edu.au/Research/EscriptFinley
DE: 7260 Theory and modeling
DE: 3200 MATHEMATICAL GEOPHYSICS (New field)
SC: Nonlinear Geophysics [NG]
MN: 2004 AGU Fall Meeting