The scripts and files included here use a waveform matching algorithm to model crustal structure using radial receiver functions and vertical autocorrelations. The code is written in Python and was used to model crustal structure across the South Central US.

The data included here is radial receiver functions and teleseismic S wave arrivals for station US-NATX, a broadband station located in Nacogdoches, Texas. 

The script files are used to run the modeling codes.

script_computeAcorr_bin.py will compute autocorrelations, and then bin and stack the radial receiver functions and autocorrelations
script_modelRF.py will run the modeling code
script_plotSolutions.py will make scatter plots and corner plots
script_sWave_importPlot.py will import and plot S arrivals, along with synthetics based on the receiver function and autocorrelation modeling results

P_DATA contains the teleseismic P arrivals and receiver functions
S_DATA contains the teleseismic S arrivals
bins contains binned receiver functions and autocorrelations
models contains results from running modeling algorithm
sWave_Compare contains the comparison plots of the real and synthetic S arrivals

To run, this needs to be in a conda environment with telewavesim: https://paudetseis.github.io/Telewavesim/

for any questions, email me at sdlrbn@gmail.com