Over the past 10 years I have generated numerous codes to analyze seismic data. These codes include the Multitaper Library (Fortran 90, Prieto et al., 2009), on-the-fly plotting in Fortran, using gnuplot and notes on computer issues I have faced in the past that may help you.

Over the past few years I have created the Multitaper Spectrum Estimation Library. It is written in Fortran 90, and there are some available adaptations to Python (http://pypi.python.org/pypi/mtspec) and Matlab (send email if you want to contact the developer). This Library is available for download. Some extra features are also available in this library, including Thomson’s Multitaper, Sine Multitaper, Quadratic Multitaper, and calculation of coherence, decolvolution, transfer functions, etc. Go to download Multitaper page. Download latest library here, mtspec_v6.0.tar.gz

If you decide to use these codes, please use the following reference. If you want to alter, modify and improve the codes, please feel free to do so but also send me an email to keep track of problems, issues or improvements that are always welcome.

Prieto, G. A., R. L. Parker, F. L. Vernon. (2009), A Fortran 90 library for multitaper spectrum analysis, Computers and Geosciences, 35, pp. 1701-1710. doi:10.1016/ j.cageo.2008.06.007. (PDF)

This is a plotting library that allows on the fly plots, using Gnuplot (http://www.gnuplot.info/) inside your Fortran programs. That means you can plot variables from inside any of your F90 codes. You do not need to save your files and then open Matlab to plot the data. In the past the codes used Bob Parker’s Plotxy (visit webpage) but is now using Gnuplot. It is possible to generate high quality plots of data in 2D, and 3D.

Visit gplot page

This notes are intended to help you (actually myself) remember how some simple Unix, and programming issues can be resolved. Here I discuss some simple examples of using Fortran 90 in various environments and in particular to use Antelope databased thru Fortran.

Visit Computer Notes page.