Skip to content

Home

Jump to bottom
Stephanie Durand edited this page Dec 19, 2017 · 2 revisions

SeisTomoPy

What is SeisTomoPy ?

SeisTomoPy is a new Python tool that facilitates the use of a suite of tomographic models available to the public, with a single programme. SeisTomoPy provides six tools that allow to visualize tomographic models, compare them and extract information for further scientific purposes. The tool comes with a graphical interface with intuitive buttons and simple parameters but the same information can also be gained by using the Python classes that can be run routinely in Python scripts. SeisTomoPy is suited for global and spherical tomographic models and is provided with a default list of recent tomographic models. For detailed information on each model the reader is referred to the original publication. However, the user can also upload his own model if desired.

Associated to this first library of models, we provide six tools that will aid extracting information about each of these models:

  • Cross section for generating cross sections and extracting values of Vs, Vp and density as input for AXISEM to generate synthetic seismograms,
  • Map for producing global tomographic maps at a given depth,
  • Spectrum for computing the amplitude spectrum at any depth,
  • Correlation for computing the correlation at any depth between the different models,
  • Path for plotting wave paths on the top of tomographic models,
  • Travel Time for computing travel times of body waves for the different tomographic models.

SeisTomoPy has been developed for saving the generated results from each of the functions in output files or figures that can be used for further scientific studies.
It is also possible to upload a model if desired. A fortran routine is provided in order to create the model files. Our choice was to develop models in spherical harmonics which:

  • enables to store the models quite efficently, and avoids to deal with huge model files,
  • is an efficient way to compute of any section, map, spectrum, etc, in a reasonable amount of time.

This means that certain global models won’t be suitable to be included in the tool, in particular those with irregular parametrization and regional models.

Clone this wiki locally