Unique composition and evolution histories of low velocity mantle domains: Data and analysis

Dataset and Jupyter notebook accompanying 'Unique composition and evolution histories of low velocity mantle domains'. 

Dataset includes:

• Present day properties of simulated mantle for simulations RCY, B=0.22, B=0.44, visc2, visc3, CMB2600, CMB2800, COMP, PRIM.
• Present day predicted seismic properties for simulations RCY, B=0.22, B=0.44, visc2, visc3, CMB2600, CMB2800, COMP, PRIM.
• Present day delta Vs for simulations RCY, B=0.22, B=0.44, visc2, visc3, CMB2600, CMB2800, COMP, PRIM, filtered using the resolution of seismic tomography model S40RTS.
• Properties of simulated mantle at 100 Myr intervals from 900 Ma - 100 Ma inclusive for simulation RCY. 
• P-T tables with predicted abundance of post-perovskite for different mantle lithologies (harzburgite, lherzolite and basalt - as defined in the paper).

Jupyter notebook `s-llvps.ipynb` contains code for identifying simulated large low-velocity provinces (S-LLVPs), extracting assoicated model properties and plotting results. Python module files terra_utils.py and ppv.py are also included and required by the code in the notebook. 

There are a number of pre-requisite packages that will need to be installed in order to run the Jupyter notebook, including terratools, a software package written specifically for reading and postprocessing outputs from TERRA simulations. Installation instructions can be found on the GitHub repository. 

Due to the TERRA code pre-dating open source licensing, we do not currently have permission to publicly share all aspects of the code. In code_pieces.F90 we include code snippets which were implemented for this study. 

Simulations were conducted using ARCHER2, the UK's national super-computing service.