https://doi.org/10.5194/gmd-2020-372 Preprint. Discussion started: 16 December 2020 c Author(s) 2020. CC BY 4.0 License. LARGE 0.2.0: 2D numerical modelling of geodynamic problems Nicola Creati1,* and Roberto Vidmar1 1National Institute of Oceanography and Applied Geophysics - OGS, Borgo Grotta Gigante 42/c, Sgonico, Trieste, Italy Correspondence: Nicola Creati (
[email protected]) Abstract. We present here LARGE 0.2.0 (Lithosphere AsthenospheRe Geodynamic Evolution) a geodynamic modelling Python package that implements a flexible and user friendly tool for the geodynamic/modelling community. It simulates 2D large scale geodynamic processes by solving the conservation equations of mass, momentum, and energy by a finite difference method with the moving tracers technique. LARGE uses advanced modern numerical libraries and algorithms but unlike com- 5 mon simulation code written in Fortran or C this code is written entirely in Python. Simulations are driven by configuration files that define thoroughly the lithologies and the parameters that distinguish the model. Documentation for them and for all the modules is included in the package together with a complete set of examples and utilities. The package can be used to reproduce results of published studies and models or to experiment new simulations. LARGE can run in serial mode on desktop computers but can take advantage of MPI to run in parallel on multi node HPC systems. 10 1 Introduction Numerical geodynamics aims to propose and to verify hypothesis and transform observations to prediction. Numerical sim- ulations help to constrain the spatial and temporal evolution of the surface and the internal structure of the Earth since some quantities cannot be measured, In the late 1960s the theory of plate tectonic pushed geodynamics to move from a descrip- tive to a quantitative view (McKenzie and Parker, 1967; Minear and Toksöz, 1970).