Supplementary material for "Dynamic restrengthening and stress heterogeneity explain megathrust earthquake complexity"

This repository contains the supplementary files of the "Coseismic restrengthening and stress heterogeneity explain megathrust earthquake complexity" by Jeremy Wing Ching Wong (wcwong@ucsd.edu), Alice-Agnes Gabriel, and Wenyuan Fan. 

[Code availability] 
- The open-source SeisSol is available at [https://github.com/SeisSol/SeisSol]
- We use the commit code '8cddb43' on the master branch to perform all dynamic rupture simulations presented in this study. 
- Instructions for downloading, installing, and running the code are available in the SeisSol documentation.
- Downloading and compiling instructions are here.
- Instructions for setting up and running simulations are here.
- Quickstart containerized installations and introductory materials are provided in the docker container and Jupyter Notebooks at the SeisSol Training Repository.
- Example problems and model configuration files are provided at the SeisSol Example repository,  many of which reproduce the SCEC 3D Dynamic Rupture benchmark problems. 

[File structure]
├── Readme.md
├── Tohoku_FL103_plasticity.par :  SeisSol parameter file for fast-velocity weakening rate-state friction
├── Tohoku_fault.yaml  :  yaml file defining frictional properties
├── Tohoku_initial_stress.yaml  :  yaml file defining initial stress conditions
├── Tohoku_lithostress_gamma.yaml  :  yaml file defining lithostatic stress and pore-pressure level
├── Tohoku_material_1D_C0_1MPa.yaml  :  yaml file for 1D velocity structure and plasticity  
├── Tohoku_nucleation_stress_R0.8.yaml :  yaml file for nucleation stress
├── FL33_stress_change_kinematic_model.yaml  :  yaml file for stress heterogeneity based on kinematic slip distribution
├── M_stress_change_1km_1Hz_1000.nc  :  NC grid file for 3D stress change tensor based on kinematic slip 
├── mesh : mesh files
  │   ├── JIVSM_v7_15d_1000_1D_1Hz.puml.h5
  │   ├── JIVSM_v7_15d_1000_1D_1Hz.xdmf
  │   └── JIVSM_v7_15d_1000_1D_1Hz.xml
├── receivers :  on-fault and off-fault receivers
  │   ├── Fault_10km_node.txt
  │   ├── Tohoku_receivers.csv
  │   └── Tohoku_receivers.dat
└── run_Seissol_general_o4.sh : SuperMUC job submission file 

[Computational Resources]
-  SeisSol simulations were conducted on the supercomputing cluster SuperMUC-NG, hosted at the Leibniz Supercomputing Center (https://www.lrz.de/), and Frontera, hosted at the Texas Advanced Computing Center (https://tacc.utexas.edu/).