Modeling crustal structure by waveform-matching P receiver functions and autocorrelations

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. 

See the readme for instructions.

This is a modified version of the code used in my previous paper: Modeling Crustal Structure in the Permian Basin by Waveform-Matching P Receiver Functions and Autocorrelograms With Particle Swarm Optimization  (https://agupubs.onlinelibrary.wiley.com/doi/full/10.1029/2023JB026730). The previous version of the modeling code is also available on Zenodo. The biggest differences between this and the previous version of the code is that this version is driven by a basinhopping optimizer instead of a particle swarm optimizer, and the code has been updated for running variable numbers of crustal layers.