halfsp - Computes motion stress vectors for Rayleigh waves in a half space
halfsp
Compatible with Basic Seismic Utilities (BSU), this program computes Rayleigh wave motion stress vectors in an elastic, half-space earth. The program is interactive, and prompts the user for the input below. Fortran 77 Version.
Input
ENTER RHO,VP,VS
These are the half-space properties which are entered on a line separated by commas RHO=mass density (kg/m^3) VP=P-wave velocity (m/s) VS=S-wave velocity (m/s) Example: 1600, 800, 200
ENTER
FREQ,NZ,Z0,ZEND
These are the frequency, number of depths, top depth, bottom
depth for motion-stress vector computation. FREQ=frequency
in Hz NZ=number of depth computations Z0=shallow-most depth
to calculate, usually 0 ZEND=maximum depth to calculate
Example: 10, 50, 0, 20
Output
Output includes Rayleighwave phase velocity (stdout) and
file halfsp.tmp which is a text listing of the computed
motion-stress vectors.
no input
standard output
phase velocity of Rayleigh Wave
halfsp.tmp
Text listing of motion-stress vectors
gendis(1), genwave(1), waves(1), showmdl(1), halfsp(1)
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>