next up previous contents index
Next: 23.1.6 Making synthetic amplitudes Up: 23.1 FOCMEC Previous: 23.1.4 Local earthquakes   Contents   Index

23.1.5 Program operation

The program makes a grid-search and finds how many polarities and amplitude ratios fit each possible solution. All solutions with less than a given number of wrong polarities and/or amplitude ratios within given error limits, are then written out and can be plotted. With a cursor, the user can then select the preferred solution, which can be stored in the input file or the database. The program is intended to work from within EEV (option F), however it can also work independently (see below). The program uses an input file called focmec.inp (automatically generated). This is a Nordic format file. Direct waves have angle \bgroup\color{black}$ >90$\egroup and refracted arrivals angle \bgroup\color{black}$ <90$\egroup degrees. If the angle is \bgroup\color{black}$ >90$\egroup , the polarity is plotted at an azimuth+180. If the user wants to use FOCMEC as a freestanding program, the angle of incidence information may have to be put in manually in a standard CAT-file, which is then renamed focmec.inp. This can be done automatically by FOCMEC if a hyp.out and corresponding print.out file is available. FOCMEC can also be used to convert angles, like dip, strike and rake to T and P-axis, simply say 'focmec a', where argument a stands for angles and you will be prompted for input.

When the program runs, all amplitude information and corresponding corrections are listed. First there is a question of which type of amplitudes to be used. In the example below, manual amplitudes are chosen.

============ FOCMEC ============

Number of polarities:            10
Amplitude types:   Manual:       8   Automatic:   8   Spectral:    8

Amplitude to use:  Manual(1), Automatic(2), Spectral(3) ?
1

No FOCMEC.DEF file, use defaults

Q: Local: Qp= 100.0**1.00  Qs= 100.0** 1.0   Global: t*(P)=1.00  t*(S)=4.00

STAT  C PH       AMP    PER TRTIME   QCOR ANGINC ANGEMG Fcor   AZ  DIST
SNART Z PG      1582   0.16   12.6    1.2    100     79  0.6  301    77
SNART Z SG      9397   0.19   21.8    1.3    100     79 -0.3  301    77
SNART T SG     10577   0.09   21.8    1.3    100     79  2.0  301    77
MUD   Z PG        53   0.10   26.3    1.4     94     85  0.3  163   179
MUD   Z SG       197   0.15   45.5    1.7     94     85 -0.2  163   179
MUD   T SG       209   0.22   45.5    1.6     94     85  2.0  163   179
BLS5  Z PG       749   0.28   28.0    1.3     94     85  0.3  326   192
BLS5  T SG      1102   0.10   49.8    1.9     94     85  2.0  326   192
BLS5  Z SG       662   0.10   49.8    1.9     94     85 -0.2  326   192

 STAT  Ratio type  T     Amp 1    Amp 2  Fcor LogRat
 SNART SV(Z)/P(Z)  V      9397     1582   1.0   0.80
 SNART SH(T)/P(Z)  H     10577     1582   0.3   0.34
 SNART SV(Z)/SH(T) S      9397    10577   3.5   0.47
 MUD   SV(Z)/P(Z)  V       197       53   1.3   0.75
 MUD   SH(T)/P(Z)  H       209       53   0.2  -0.12
 MUD   SV(Z)/SH(T) S       197      209   7.5   0.87
 BLS5  SH(T)/P(Z)  H      1102      749   0.2  -0.45
 BLS5  SV(Z)/P(Z)  V       662      749   1.3   0.23

The abbreviations are STAT: Station code, C: Component, PH: Phase, AMP: Amplitude in count, PER: Period in sec, TRTIME: Travel time in sec, QCOR: Log Q-correction, ANGINC: Angle of incidence at the source, ANGEMG: Angle of emergence at the station, Fcorr: Free surface correction for this amplitude, Az: Azimuth from the event to the station, DIST: Epicentral distance in km., Ratio type (see text), T: indicator of ratio type, Amp1 and Amp2: The two amplitudes (count) in the ratio, Fcor is the free surface correction in the amplitude ratio (to be multiplied with ratio) and LogRat is the logarithm of the corrected amplitude ratio used.

Note that for station SNART, amplitudes were also read on the radial component so more then 3 amplitude ratios were used.

Following, the user get the choices:

Stop                        (0)
Plot saved solution(s)      (1)
Plot new solutions          (2)
Plot selected solution      (3)
Find new solutions          (4)
-1, -2, -3 also plot station

  1. This is the solution(s) already stored in the data base (S-file). See secetion "Storing and selecting fault plane solutions" below.
  2. Plotting new solution after having used option 4
  3. Plotting the selected solution after using option 4
    Using e.g. -1 instead of 1, also plots the stations to help identify them on the plot, see Figure 23.1
  4. Starting a search for new solutions

    Option 4 gives the following information and questions:

    There are   10 polarity readings
    Maximum number of allowed polarity errors or -1 to show best solutions only
    

    Depending on number of data values, 0-5 is a good answer. To let the program find the minimum number of polarity errors, type '-1', which is particular useful if there is a significant minimum number of polarity errors.

    There are  8 amp ratio readings
    Maximum number of allowed amplitude ratio errors
    

    Equivalent for ratios to 'Maximum number of polarity errors', however, error is defined by amplitude ratio error. Number of errors depends on number of observations. For 9 observations 1-2 errors is reasonable.

    Maximum amplitude ratio error,  return for default of .2
    

    Give maximum allowed difference between observed and computed log amplitude ratio, default is 0.2, which often is a good value.

    Degree increment in search, enter for default 2
    

The program will now start the searching and write out on the screen (and in a file) the solutions which fit the requirement of number of misfits. The maximum number of solutions is limited to 100 as a default, or to the value defined by `FOCMEC MAXSOL' in SEISAN.DEF. At the end, the number of acceptable solutions is written out as well as the minimum number of bad fits. This can then be used for the next search. Now option 0 to 4 can be used again.

When plotting the solution with option 2, the cursor comes up. Also, the solutions will be printed in text form to the screen, see Figure 23.1.

The abbreviations are Pol: Number of polarity errors for P, SV(not used) and SH(not used), Rat Err: Number of ratio errors, RMS RErr: The RMS error for the ratios used, RErr (All): The RMS error for all ratios.

The polarities and amplitude ratios can be plotted on the focal sphere using the same convention as the original FOCMEC program, which is:

o = compression
+ = emergent compression
\bgroup\color{black}$ \Delta$\egroup = dilatation
- = emergent dilatation
V = amplitude ratio SV/P
S = amplitude ratio SV/SH
H = amplitude ratio SH/P

The user can select a preferred solution by moving the cursor near one of the letters T or P (T and P axis). By pressing T, the program will find the nearest T axis (same for P and nearest P-axis) and corresponding fault plane solution, which can be stored in the database and/or plotted with option 3. If no solution is to be selected, press q for quit. If a solution has been selected, the user will be asked if it is to be saved or not after selecting option 0. The saved solution goes into the focmec.out and from there into the S-file (type F-line) in the database if FOCMEC is operated from EEV and the solutiosn will also be written to fps.out.

When working from EEV, the event will always be located before the FOCMEC program starts up. In the Nordic format the solution is stored simply as strike, dip, rake and number of bad polarities (3f10.1,I5). Aki and Richards convention is used. In addition, the name FOCMEC will be written near the end of the line to indicate that the fault plane solution was made by FOCMEC. The line type is F.

The following files are created:

focmec.dat: Input parameters to FOCMEC_EXE.
focmec.log: Log of the FOCMEC_EXE run.
focmec.lst: More details on solutions
focmec.out: Gives input parameters and solutions
focmec.eps: A Postscript plot file of LATEST plot
focmec.run: Run parameters for FOCMEC_EXE, you can re-run FOCMEC by `focmec_exe < focmec.run'


next up previous contents index
Next: 23.1.6 Making synthetic amplitudes Up: 23.1 FOCMEC Previous: 23.1.4 Local earthquakes   Contents   Index
Peter Voss : Mon Feb 27 10:16:12 UTC 2017