Gpsh

Synopsis

gpsh [OPTIONS] [FILE]

Description

Computes the SH transfer function for layered models given through stdin or FILE.

The input is the soil column defined as listed below. The output is the Fourier transfer function for VERTICALLY INCIDENT plane S waves

Format for layered models:

 Line 1  : <number of layers n including half-space for first model>
 Line 2  : <thickness (m)> <Vp (m/s)> <Vs (m/s)> <Density (kg/m3)>[ <Qp> <Qs>]
 ....
 Line n+1 : 0 <Vp (m/s)> <Vs (m/s)> <Density (kg/m3)>[ <Qp> <Qs>]
 Line n+2 : <number of layers m including half-space for a second model - if wished>
 Line n+3 : <thickness (m)> <Vp (m/s)> <Vs (m/s)> <Density (kg/m3)>[ <Qp> <Qs>]
 ....
 Line n+m+2 : 0 <Vp (m/s)> <Vs (m/s)> <Density (kg/m3)>[ <Qp> <Qs>]
 Line n+m+3 : <number of layers l including half-space for a third model - if wished>
 ....

Quality factors are mandatory. Any number of models can be given as input.

Gpsh options

• -n <count> Number of frequency samples (default=100)
• -max <max> Maximum frequency up to which the transfer function s to be computed (default=20 Hz)
• -angle <angle> Incidence angle (default=90, vertical incidence)

Generic options:

• -help <SECTION> Show help about options. SECTION may be empty or:

                           all, examples, generic, gpsh

• -version Show version information
• -app-version Show short version information
• -reportbug Start bug report dialog, information about bug is

                           passed through stdin. This option is used
                           internally to report bugs if option -nobugreport is
                           not specified.

• -reportint Start bug report dialog, information about

                           interruption is passed through stdin. This option
                           is used internally to report interruptions if
                           option -nobugreport is not specified.

• -nobugreport Do not generate bug report in case of error

Examples

        gpsh < test.model

 Calculate SH transfer function from 0.2 Hz to 20 Hz for model 'test.model'.

        gpsh < test.model | figue -c

 Calculate the same function and plot it.

        gpsh < test.model | figue -c -m sh.mkup

 Show the same function on a log-log plot. 'sh.mkup' is a tar.gz file
 containing an xml description of the graphic format, it can be generated from
 figue's interface.

Copyright

Copyright (c) 2002-2008 Marc Wathelet. All rights reserved. See http://www.geopsy.org for more information.