OCTAVE moho.m - Sample program to compute Rayleigh wave dispersion (period x-axis) Crustal Model
moho prompts for:
IF exists
disper.oct
Nothing. Just run it and see the plot.
ELSE
ELF LSB shared object disper.oct will need to be
compiled
CLICK OK
from within octave, the object will be built, provided that
the
following files exist in the path
rwv.f Fortran Source, wrapper.cpp, build_disper_oct script
ENDIF
Basic Seismic Utilities (BSU) program run in Octave. Within an octave session, type moho. The program runs and generates a plot for the built in soil profile. This is a demonstration example that links an ELF LSB shared object, disper.oct with a running Octave session. One must first run the build script, build_disper_oct which compiles rwv.f a Fortran subroutine. Calling on wrapper.cpp, the result is the shared object, disper.oct. You need only build disper.oct once, and then can link it to any program you may write, or to the other BSU *.m files that depend on it (see ALSO below).
Plot of moho phase velocities for built in model
One first
compile the Fortran subroutine, rwv.f using the files:
rwv.f
wrapper.cpp
build_disper_oct
disper.oct (created by running build_disper_oct)
disper(1)
waves(1) gendis(1) genwav(1)
OCTAVE: FwdR1.m(7) invR1.m.(7) rayleigh.m(7)
mastercurve.m(7)
saswv.m(7) SASW.m(7)
No known bugs.
Copyright © 2024 by Paul Michaels
This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
P. Michaels, PE. <[email protected]>