If you're going to use any of my software, I'd appreciate an e-mail, in order to keep track of how far it's spread -- thanks.
Ray-theoretical modelling of teleseismic waves in dipping, anisotropic structures
Download it here (zipfile). An example with necessary files to run it is here.
The current version is 1.3, updated Dec. 15, 2022. I've fixed a bug whereby the very top layer counldn't be isotropic; also, I've fixed another issue -- not a bug, but a wrong assumption -- whereby the multiples generated excluded some double-converted phases with significant amplitude. This had the side effect of making Moho multiples appear to be overly high amplitude, because a phase that cancelled some of that amplitude was excluded. Nonetheless, for complex uses of Raysum you should probably consider creating your own phase list. Thanks to Wasja Bloch for useful discussions during this update -- his Python wrapper for Raysum is available on Github.
See here for what this does. This is only the forward-modelling part of the project; the inversion code is based on Malcolm Sambridge's original implementation of his neighbourhood algorithm and cannot be distributed without his permission (which isn't hard to get -- e-mail me if you're interested). The program is written in Fortran 77 (with the END DO and IMPLICIT NONE extensions) and should be compilable just about anywhere; instructions are included.
Neighbourhood-algorithm inversion using this code -- it's a bit rough and not very well documented.
Other codes not yet on the web -- email me for info: 1) Scattering tomography for receiver functions, 2) directionally balanced shear-wave splitting, 3) transfer-function modelling of teleseismic waves without deconvolutio.