Understanding Sub-Lithospheric Small-Scale Convection By Linking Models Of Grain Size Evolution, Mantle Convection and Seismic Tomography: Data

Understanding Sub-Lithospheric Small-Scale Convection By Linking Models Of Grain Size Evolution, Mantle Convection and Seismic Tomography: Data This repository accompanies the paper Dannberg, J., Eilon, Z., Russell, J.B., Gassmöller, R.: Understanding Sub-Lithospheric Small-Scale Convection By Linking Models Of Grain Size Evolution, Mantle Convection and Seismic Tomography. This repository contains instructions on how to compile, execute, and postprocess the geodynamic models. Additionally, it contains scripts to convert the model output to seismic observables and analyze the resulting seismic models. To reproduce the workflow of the publication follow these steps: - Compile the version of the ASPECT geodynamic modeling software found in the `aspect` folder. You can find instructions on how to compile ASPECT in `aspect/doc/sphinx/user/install`. To make sure you use the same (or similar) underlying libraries, you should see the following (or similar) output when executing ASPECT: ```------------------------------------------------------------------------------- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion.-- . version 2.5.0-pre (grain_size_ridge_06_2023, 58cdddba8)-- . using deal.II 9.4.2-- . with 64 bit indices and vectorization level 3 (512 bits)-- . using Trilinos 13.2.0-- . using p4est 2.3.2-- . running in OPTIMIZED mode-- . running with 1536 MPI processes-----------------------------------------------------------------------------``` - Execute the geodynamic models found in the `geodynamic_models` folder. This step will require at least about 500 MPI processes (dependent on the amount of memory on your system) and several days of time per model. - The folder `geodynamic_postprocessing` contains a number of Python scripts and the necessary parts of the original model output to recreate the analysis of the geodynamic model output. - The folder `seismic_postprocessing` contains Matlab and Python scripts that convert the geodynamic model output and perform seismic analyses. The folder contains a separate file `README.txt` that explains the process.