invR1.m

NAME
SYNOPSIS
DESCRIPTION
OUTPUT
REQUIRES:
SEE ALSO
BUGS
COPYRIGHT
AUTHOR

NAME

OCTAVE invR1.m - Invert Rayleigh wave dispersion data (model.txt, bvax.his) (depth,Vs) points

SYNOPSIS

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

DESCRIPTION

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.

OUTPUT

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

REQUIRES:

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

SEE ALSO

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)

BUGS

No known bugs.

COPYRIGHT

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.

AUTHOR

P. Michaels, PE. <paulmichaels@.boisestate.edu>