Surface Wave Processing

The BSU package includes some software for Rayleigh wave analysis. These programs are as follows:

• bvax.f A fortran program which extracts a dispersion curve from a bsegy data set.
• FwdR1.m A Octave program which does the forward problem and displays a measured dispersion curve. Where invR1.m is automated, this program is interactive, requiring the user to make changes to the model with each iteration.
• invR1.m An automated Octave inversion program which reads the dispersion curve from a bvax run. The input file to the Octave program is bvax.his. If multiple runs of bvax are to be made, the file bvax.his should be deleted before each new run, since it is appended to during a bvax execution.

Example programs:

• rayleigh.m This is an example program which shows how to link the dispersion computing module disper.oct (shared object) with the octave session. The fortran subroutine in file rwv.f is based on the fortran program disper.f included in BSU. When running rayleigh.m, the subroutine file, rwv.f must be in the directory where one is launching Octave. On execution of rayleigh.m, a check is made to see if a compiled version exists, and if not a compilation step occurs. The object module is dynamically linked to the Octave procedure to produce high speed computations. In addition to the Fortran file, rwv.f, two other files should also be in the path or directory of computation. These are wrapper.cpp and build_disper_oct.
• moho.m Similar to rayleigh.m, but computes dispersion as a function of period rather than frequency. Also an example program.