FEniCS-SZ

Authors: Cian Wilson, Cameron Seebeck, Kidus Teshome, Nathan Sime, Peter van Keken

Welcome to FEniCS Subduction Zones (FEniCS-SZ), a python-based package for the thermal modeling of subduction zones!

This package was developed by undergraduate interns Kidus Teshome and Cameron Seebeck at the Carnegie Science Earth & Planets Laboratory. It is based on Wilson & van Keken, PEPS, 2023 (II), which is part II of a three part introductory review of the thermal structure of subduction zones by Peter van Keken and Cian Wilson.

Our goal is both to develop software tools for kinematic-slab thermal models of subduction zones in as open-source a way as possible and to provide an easily accessible and modifiable source of the global suite of subduction zone models as described in Wilson & van Keken, PEPS, 2023 (II) and van Keken & Wilson, PEPS, 2023 (III). We do this by basing all of our development in jupyter notebooks and making them accessible through the FEniCS-SZ website.

This website is published as a jupyter book. Each page is a jupyter notebook that can be run interactively in the browser. To start such an interactive session using binder click the -symbol in the top right corner of the relevant page. For more details on running an interactive session see the usage section.

Comments and corrections to the package should be submitted to the issue tracker by going to the relevant page in the jupyter book, then clicking the -symbol in the top right corner and “open issue”.

Acknowledgments

We make heavy use of the finite element library FEniCSx and linear algebra package PETSc. Simulations are run in parallel using ipyparallel and we rely on python’s extensive scientific ecosystem, such as numpy, scipy, matplotlib, and pyvista, combined with project jupyter’s web-based, interactive, annotated notebooks.

This jupyter book is based on the FEniCSx Tutorial by Jørgen S. Dokken, which is an excellent resource for learning how to use FEniCS in a similar interactive jupyter book.