next up previous contents
Next: ampcorr Up: Computer programs Previous: Computer programs   Contents

sti

The stress tensor inversion program sti implements the ideas that was presented in paper II. The program is written entirely in C and includes library functions for vector and matrix manipulation, angular conversion between different focal mechanism conventions and the stress estimation routines. Input data are taken either from the SIL data base or from a program native format, obtainable by the utility program conv2sti. The grid size can be chosen between increments of $ {2}^{\circ}$ to $ {20}^{\circ}$, with $ R$ in increments of 0.05 or 0.1. Output data are in the form of ASCII files with information on the estimated state of stress, confidence limits and chosen fault planes. I have a number of GMT (Generic Mapping Tool) [Wessel and Smith, 1998] scripts for the visualization of the results.
thor:~/sil/rst> sti -h

   STRESS TENSOR INVERSION OF EARTHQUAKE FOCAL MECHANISMS

sti [-h] data file [-c config file][-d num] [-f #] [-l #] [-n]
    [-o output file name] [-s]

    The program inverts earthquake focal mechanisms or geologic
    fault slip data for the stress tensor that caused the movement on
    the faults.
    The inversion process is described in Lund and Slunga, JGR, 1999.

    For optimal use of the program's error messages, run the program as:
    sti data_file [options] >& data_file.err &


    data file       The program accepts input data either in SIL format,
                    a *.grx or *.grs file, or in a program native format.
                    The native format can be obtained from strike,dip and
                    rake data with the program conv2sti.

    Options:
    -c      Configuration file with other name than rst.conf or
            in another directory. Default is rst.conf in the
            current directory.
    -d      Dump SIL data into other formats:
            1       A total dump in sti native format.
            2       Dump optimal planes in Gephart & Forsyth #7
            3       Dump final planes in Gephart & Forsyth #7
    -f      First event in inversion. Number in list of approved events,
            default is 1.
    -h      This text.
    -l      Last event in inversion. Number in list of approved events,
            default is last.
    -n      Write a file with the non-redundant events as defined by
            the gridding of Magee, 1997
    -o      Base name for output files. Default is the input file name
            extended with the inversion method
    -s      Read in the events and write out relevant info, then exit.

    Bjorn Lund, November 1996. Version 1.4, BL June 1999.



Bjorn Lund 2000-06-14