OCTAVE invR1.m - Invert Rayleigh wave dispersion data (model.txt, bvax.his) (depth,Vs) points
invR1 prompts for:
IF exists
disper.oct
Enter name of model file (ex. model.txt)
GUI Entry Box p=number singular values
deltz=layer thickness used between control points
locksw = fixes some parameters
locksw=0 free bottom control, velocity and depth
locksw=1 lock bottom depth, free bottom velocity
locksw=2 lock bottom velocity, free bottom depth
locksw=3 lock both bottom velocity and depth (default)
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)
GUI Entry Box p=number singular values
deltz=layer thickness used between control points
locksw = fixes some parameters
locksw=0 free bottom control, velocity and depth
locksw=1 lock bottom depth, free bottom velocity
locksw=2 lock bottom velocity, free bottom depth
locksw=3 lock both bottom velocity and depth (default)
ENDIF
Basic Seismic
Utilities (BSU) program run in Octave. Within an octave
session, type invR1. The program prompts as above.
Enter the text file with the starting model. There are 3
rows, number of control points, Vs of control, depths of
control.
EXAMPLE: 4 control points in model.txt
4
140 90 300 300
0 6.0 6.001 30
Enter any changes for the parameters, click OK. The program assumes the observed data are in a file, bvax.his. This file would have been created in an earlier bvax run.
Plot of
frequency vs. phase velocities for soil profile updated
during inversion.
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
Plot of singular values
FILES that can
be run with FwdR1.m
solution.txt
slow.txt (using uncertainties, slower limit on solution)
fast.txt (using uncertainties, faster limit on solution)
FILES that can
be run with Gnuplot
dispcrv.dat and dispcrv.gp
model.dat control.dat and model.gp
wavelen.dat and wavelen.gp
sval.dat and singlval.gp
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)
bvax.his
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 © 2017 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. <pm@cgiss.boisestate.edu>