This repository contains the output from a global suite of subduction zone models used in: "An introductory review of thermal structure of subduction zones: II. Numerical approach, validation, and comparison", van Keken & Wilson, PEPS, 2023. Directories & Files The subdirectories are: `D80`: a reproduction of the "D80" Syracuse et al., PEPI, 2010 suite of subduction zones 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. 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. 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