This repository contains the output from a global suite of subduction zone models used in: "An introductory review of the thermal structure of subduction zones: II—numerical approach and validation", Wilson & van Keken, PEPS, 2023 and "An introductory review of the thermal structure of subduction zones: III—Comparison between models and observations", van Keken & Wilson, PEPS, 2023. Directories & Files The subdirectories of `vankeken_wilson_peps_2023.zip` are: `D80`: a reproduction of the "D80" Syracuse et al., PEPI, 2010 suite of subduction zones with a "dynamic" slab (Stokes equations solved in the slab portion of the domain as in van Keken & Wilson, PEPS, 2023 and `D80new` below) and performed using Sepran `D80new`: a new global suite based on Syracuse et al., PEPI, 2010 but using a lower mantle potential temperature and an error function boundary condition and performed using both Sepran and TerraFERMA. Input files provided for TerraFERMA. `benchmark`: a new subduction benchmark performed using both Sepran and TerraFERMA. Input files provided for TerraFERMA. `Figure_III_4_models`: the results used in Figure 4, Part III. Files include: `Table2_Syracuse2010_updated.pdf` and `Table2_Syracuse2010_updated.xlsx`: updated parameter values for the Syracuse et al., PEPI, 2010 global suite. `Table3_part_Syracuse2010_updated.pdf` and `Table3_part_Syracuse2010_updated.xlsx`: updated parameter values for the Syracuse et al., PEPI, 2010 global suite. `Figure_comp_D80_WW09_all.pdf`: a comparison between the Sepran models presented in `D80` and `D80new` and those in Wada & Wang, G-cubed, 2009. `Figure_comp_D80new_TF-Sepran_all.pdf`: comparison between the TerraFERMA and Sepran models presented in `D80new`. `Figure_III_6_all_D80.pdf`: complete set of model results for Figure 6, Part III. `Subduction_parameters.pdf` and `Subduction_parameters.json`: a description of the parameters used in the global suite. `PEPS_part_II_preprint.pdf`: preprint of part II. `PEPS_part_III_preprint.pdf`: preprint of part III. `vankeken_wilson_peps_2023_TF_lowres_minimal.zip`: contains a set of low resolution vtu files from TerraFERMA with the final potential temperature output from steady state and time dependent simulations. Intended for fast downloading (compared to `vankeken_wilson_peps_2023.zip`) and comparison. Running a simulation Input files are provided for the TerraFERMA (the Transparent Finite Element Rapid Model Assembler) simulations in this repository (`TF` subdirectories). Running them requires a working installation of TerraFERMA. If one is not available then consider using the docker image provided for this paper at: https://github.com/users/cianwilson/packages/container/package/vankeken_wilson_peps_2023 This docker image contains a complete installation of TerraFERMA and its dependencies, PETSc, FEniCS and SPuD, within an Ubuntu 20.04LTS OS. For a full description of TerraFERMA please refer to the webpage: http://terraferma.github.io