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 and refracted arrivals angle degrees. If the angle is , 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
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:
|+ =||emergent compression|
|- =||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'