Software for finite-frequency tomography - Download from this page

 

 

General remarks

The software that is made available on this page has been developed by many authors. The most important contributions came from (in alphabetical order): Jean Charlety, Shu-Huei Hung, Raffaella Montelli, Guust Nolet, Karin Sigloch, Yue Tian and Jean Virieux.

We are still in the process of writing more extensive manuals, for now you'll have to make do with what is offered. At times you may have to consult the program itself (e.g. for output formats).

For an extensive background see Nolet, G., A Breviary for Seismic Tomography, Cambridge Univ. Press (2008). Click here for a running list of errata. You can also download the Solutions to the Excercises, courtesy of Yue Tian.

In the following, equation numbers refer to the equations in this book, sometimes referred to as Breviary.

Dynamic raytracing in a spherically symmetric Earth

There are two packages of software. The first one, program Raydyntrace.f is a stand-alone program that can be used to program your own banana-dougnut kernels in a spherical Earth. The latest version (Jan 15, 2014) has improved crustal corrections in oceanic areas (including for OBS's) and computes also reflection- and transmission coeficients.

Those who have already developed their own tomography software, and wish to move from single-frequency, ray-theoretical inversions to multiple-frequency inversions using banana-doughnut kernels, shall wish to use this program. It provides all variables needed to program eq. (7.28) or its simplified version (7.20) for cross-correlation delay times, eq. (7.30) for delays obtained by cross-correlation waveform fitting, or the expression for waveform maxima picks (see page 136 of Breviary).

It can be used with the paraxial approximation (computing detour times for scatterers off the geometrical ray path assuming the wavefront has a parabolic shape), but it can also be used to provide a table of all necessary parameters to be used in a table look-up. The tar file contains extensive documentation as well as examples.

The package of programs that we use ourselves to compute the matrix and solve the system is still in continuous evolution, most recently because we have abandoned irregular tetrahedrons as parameters, and replaced this with a cubic Earth representation. We therefore do not put the software up on the web, but if you are interested in acting as a guinea pig and help testing whatever is available, please contact Guust Nolet (nolet@geoazur.unice.fr).

Dynamic raytracing and matrix computation in 3D shallow models

The latest version has the programs grafbdyn.f, matrix3D.f and assemblematrix3D.f, as well as a number of utilities to plot or analyse the output. These result in a fully assembled matrix, including correction factors.

The final MPI program mpisolvetomo3D.f that does the inversion is still subject to changes and improvements. It is available in its beta-version to whoever wishes to collaborate with us. If you are interested in this, please contact Guust Nolet (nolet@geoazur.unice.fr)

The package BD3D.19Apr2013.tar.gz has all you need to compute banana-dougnut kernels and the matrix in a 3D Cartesian model.

The BD3D manual is currently being written, the latest version can be downloaded separately.

References

1. Nolet, G., A Breviary for Seismic Tomography, Cambridge Univ. Press, 2008.

2. Tian, Y., G. Nolet, F.A. Dahlen, S.-H. Hung and R. Montelli, Dynamic ray tracing and traveltime corrections in finite-frequency tomography, J. Comp. Phys., 226, 672-687, 2007.

3. Tian, Y.,R. Montelli, G. Nolet, and F.A. Dahlen, Computing traveltime and amplitude sensitivity kernels in finite-frequency tomography, J. Comp. Phys., 226, 2271-2288, 2007a.

 

Return to the GLOBALSEIS home page