The programming language python in earth system simulations

-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.

Anthropogenic and Natural Organohalogen Compounds in Blubber of Dolphins and Dugongs (Dugong dugon) from Northeastern Australia

A range of organohalogen compounds (10 polychlorinated biphenyl [PCB] congeners, DDT and metabolites, chlordane-related compounds, the potential natural organochlorine compound Q1, toxaphene, hexachlorobenzene, hexachlorocyclohexanes, dieldrin, and several yet unidentified brominated compounds) were detected in the blubber of four bottlenose dolphins (Tursiops truncatus), one common dolphin (Delphinus delphis), and seven dugongs (Dugong dugon), as well as in adipose tissue of a green turtle (Chelonia mydas) and a python (Morelia spilota) from northeast Queensland (Australia). The green turtle and dugongs accumulated lower organohalogen levels than the dolphins. Lower levels in dugongs were expected because this species is exclusively herbivorous. Highest PCB and DDT levels recorded in dugongs were 209 and 173 mug/kg lipids, respectively. Levels of the nonanthropogenic heptachlorinated compound Q1 (highest level in dugongs was 160 mug/kg lipids) were estimated using the ECD response factor of trans-nonachlor. Highest organohalogen levels were found in blubber of dolphins for sumDDT (575-52,500 mug/kg) and PCBs (600-25,500 mug/kg lipids). Furthermore, Q1 was a major organohalogen detected in all samples analyzed, ranging from 450 -9,100 mug/kg lipids. The highest concentration of Q1 determined in this study represents the highest concentration reported to date in an environmental sample. Levels of chlordane-related compounds were also high (280-7,700 mug/kg, mainly derived from trans-nonachlor), but concentrations of hexachlorobenzene, hexachlorocyclohexanes, dieldrin, and toxaphene were relatively low and contributed little to the overall organohalogen contamination. Furthermore, a series of three major (BC-1, BC-2, and BC-3) and six minor (BC-4 through BC-9) unknown brominated compounds were observable by extracting m/z 79 and m/z 81 from the GC/ECNI-MS full scan run. Structural proposals were made for the two major recalcitrant compounds (referred to as BC-1 and BC-2). BC-2 appears to be a tetrabromo-methoxy-diphenylether (512 u) and BC-1 has 14 u (corresponding with an additional CH2 group) more relative to BC-1. In general the organohalogen pattern observed in blubber of dolphins was different compared to similar samples from other locations in the world, which is apparent from the fact that the four major abundant signals in the GC/ECD chromatogram. of D. delphis originated from the four unknown compounds Q1, BC-1, BC-2, and BC-3.