UWGeodynamics: A teaching and research tool for numerical geodynamic modelling

International audienceThe UWGeodynamics module facilitates development of 2D and 3D thermo-mechanicalgeodynamic models (Subduction, Rift, Passive Margins, Orogenic systems etc.). It isdesigned to be used for research and teaching, and combined the flexibility of the Under-world Application Programming Interface, (Moresi, Dufour, & Mühlhaus, 2002, Moresi,Dufour, & Mühlhaus (2003), Moresi et al. (2007)) with a structured workflow.Designing geodynamic numerical models can be a daunting task which often requiresgood understanding of the numerical code. UWGeodynamics provides a simple interfacewith examples to get you started with development of numerical models. Users can startdesigning their models without any pre-existing knowledge of programming. Expert userscan easily modify the framework and adapt it to more specific needs. The code can be runin parallel on multiple CPUs on personal computers and/or High Performance Computingsystems.Although UWGeodynamics has been primarily designed to address geodynamic problems,it can also be used to teach fluid dynamics and material mechanics.UWGeodynamics uses the flexibility of the Python language and the Jupyter Notebookenvironment, which allows leveraging the wide range of scientific libraries available fromthe Python community. It also facilitates the coupling with existing scientific Pythonmodules such as Badlands (Salles, Ding, & Brocard, 2018).The functionalities include:•Dimensional input values, using user’s choice of physical units.•Automated and transparent scaling of dimensional values.•Sets of predefined geometries that can be combined to define the initial geometryof a model.•Handles Newtonian and non-Newtonian rheologies (Viscous, Visco-plastic andVisco-elasto-plastic).•Database of common rheologies used in geodynamics, which can be personalised /extended by users.•Simple definition of kinematic, stress, and thermal boundary conditions.•Lithostatic pressure calculation•Thermal equilibrium (steady-state) calculation.•Pseudo Isostasy using a range of kinematic or stress boundary conditions.•Partial melt calculation and associated change in viscosity / heat production.•Simple definition of passive tracers and grid of tracers.•Simple Phase changes•2-way coupling with the surface processes model pyBadlands (Salles et al., 2018).UWGeo comes with a series of examples, benchmarks and tutorial setups that can be usedas cookbook recipes. They provide a wide range of teaching materials useful to introducenumerical geodynamic modeling to students.New functionalities are constantly being added to the code and contributions are morethan welcomed. You can access the full documentation online athttps://uwgeodynamics.readthedocs.i