OCTAVE FwdR1.m - Compute forward model Rayleigh wave dispersion (model.txt, bvax.his)
FwdR1 prompts for:
IF exists
disper.oct
Enter name of model file (ex. model.txt)
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
Enter name of model file (ex. model.txt)
ENDIF
Basic Seismic Utilities (BSU) program run in Octave. Within an octave session, type FwdR1. The program runs and generates a plot for the soil profile in model.txt or a file name of your choosing. 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). If file bvax.his is included from a bvax run, then the observed data will be plotted with the calculated data.
Plot of
frequency vs. phase velocities for soil profile in model.txt
Plot of velocity vs. wavelength
Plot of the model (Vs, Vp, Density) Title has VSave=time
averaged velocity top to bottom
NOTE: if bottom of soil profile definition is 30 meters,
VSave=Vs30
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)
Optional: bvax.his if available (code tests for presence of
file)
disper(1)
waves(1) gendis(1) genwav(1)
OCTAVE: FwdR1.m(7) invR1.m.(7) moho.m(7) mastercurve.m(7)
rayleigh.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]>