Software available to the community

I maintain several software packages which provide functionality which may be useful to the geophysicist. Bug reports and pull requests to add functionality are welcome, and I will try to see to these in a reasonable time frame.

The following codes—and other useful packages—are hosted at my GitHub pages:

• Seis.jl – a fast, flexible package for dealing with seismic data

  Seis.jl is a Julia package for dealing with seismic data. It is fast (many times faster than Obspy for most operations), flexible, and open source.

  Seis.jl provides the basics building blocks to process data (e.g., filtering), and other packages build on top of Seis.jl's basic types to provide other functionality:

  • SeisSplit.jl: Shear wave splitting analysis.
  • Beamforming.jl: Beamforming with seismic arrays.
  • SeisRequests.jl: Download seismic data from IRIS and other data centres.

• CIJ.jl – Anisotropy in Julia

  A package to deal with elastic constants for geophysics problems, including computing velocities from arbitrary elastic constants.

• SeisModels.jl – Global seismic models

  Obtain seismic and other physical properties at any point in the Earth for various Global seismic models, including computing gravity, moment of intertia, and so on.

• Fortran geophysics utility libraries

  General utility modules for dealing with various geophysical problems, including:

  • anisotropy_ajn.f90, for dealing with anisotropic sets of elastic constants
  • global_1d_models.f90, for getting 1D Earth properties from PREM and AK135
  • plate_motion.f90, to compute absolute and relative plate motion from various models
  • spherical_splines, allowing fitting of arbitrary functions onto spherical B-splines
  • splitwave, modelling shear wave splitting in anisotropic media along precomputed ray paths

• Elasticity – Anisotropic modelling on the command line

  Making use of seismo-fortran, this repo gives contains many small command-line programs for manipulating and plotting first-order general elastic constants.

  

  3D view of the shear wave splitting produced by a kimberlite as given by Ben Ismaïl and Mainprice, Tectonophysics, 1998, produced using CIJ_plot from the repo.

• f90rad3 – Command-line GPR

  Fortran codes to deal with RAD3-format data as produced by, for example, Malå ground-penetrating radar units. Many proprietary or graphical programs exist for processing such data, but few offer any automated way to create plots and filter the data. These codes aim to make that process easier.

  Output options include making sections and writing to NetCDF for visualisation in your 3D program of choice (such as Paraview).

  

  3D GPR volume, viewed in Paraview, and created by f90rad3.