next up previous contents index
Next: 29. Making synthetic seismograms Up: 28. Merge events in Previous: 28.3 Catalogue merging and   Contents   Index


28.4 COMPARE_HYP, compare hypocenters of two cat-files

This program compares two Nordic event files with the same events in both. The intention is to be able to see the effect of locating the same data set using a different model and or different magnitude scales. The input requires two files: The first file is the reference and the second file is the modified reference file. The compared content is origin time, rms, hypocenter and magnitudes. For each parameter, the average difference with standard deviations is calculated. The difference is calculated as the value in the reference minus the value in the modified file. In order to compare the files, they must have the same number of events and it is assumed that the events are the same and come in the same order, however the only thing checked is that the files contain the same number of events. If an event has no data (no hypocenter), the event is skipped. The the list of skipped events is given in a file compare_hyp.skip. The number of magnitudes for each event does not have to be the same so not all events need to have magnitudes. Only magnitudes available in both data sets are compared and listed. The magnitudes are always listed in order Ml, Mc, Mb, MB, Ms, MS, Mw. For the magnitudes there is no check if the agency is the same. If a magnitude appears 2 times, the last one will be used. Note that magnitude in position 3 on header line will not be overwritten when relocating unless it has the agency identical to the agency in the STATION0.HYP file, this can create problems when comparing so make sure the agency in that position corresponds to the agency used when locating. There is no check if epicenter is fixed since normally this is what will be compared. Since the depth very often is fixed, fixed depths are not used for comparison. The output is a file called compare_hyp.out. It lists each event with the data from the reference event and the difference with the relocated event. Events with a difference of more than user defined values for the parameters, have their differences marked with a star, see example below. Events with less than a given number of stations are skipped. The skipped events are listed in file compare_hyp.skip. An example run is shown below. The default minimum differences used indicating with a star are 5 km in lat-lon and depth, 1 s in origin time, 0.5 s in rms and 0.2 magnitude units.

Example run

c:\seismo\wor>compare\_hyp
 Give input file 1, the reference
eev.out
 2016  721 1517 46.7 D  41.340  44.040  5.0F SC3 20 1.5 4.2LSC3 5.7BSC3 4.3bSC3
 2016  721 1628 47.7 D -23.860-175.560 10.0F SC3 43 1.9 6.0LSC3 5.8BSC3 5.3bSC3
 2016  721 23 9  3.5 D  40.780-124.070 10.0F SC3 33 2.0 4.9LSC3 4.9BSC3 4.7bSC3
 2016  722  040 26.3 D  39.360  25.970  7.0  SC3 34 1.4 3.1LSC3 3.0SSC3
 2016  722 1131 49.5 D  19.900 120.900 10.0FSSC3 17 1.8 5.1BSC3 4.7bSC3
 2016  722 2140 10.2 D  42.400  35.210 23.0  SC3 49 1.0 4.4LSC3 3.9SSC3 4.6BSC3
 2016  722 2344 57.8 D  36.940  26.680  8.0  SC3 15 1.0 2.7LSC3
 2016  723  1 0 20.0 D  47.680 147.050395.0  SC3214 1.0 6.9LSC3 5.7BSC3 5.7bSC3
 2016  723  226 20.3 D  -7.410 128.920132.0  SC3 19 1.1 5.4LSC3 5.3BSC3 5.0bSC3
 2016  723  252 12.1 D  30.140 142.110 61.0  SC3 14 0.7 4.9LSC3 4.9BSC3 4.8bSC3
 2016  723  741 45.1 D  -6.710 101.270 10.0F SC3 83 1.5 6.0LSC3 5.8BSC3 5.3bSC3
  End of s-file
 Give input file 2, the file to compare
hyp.out
 2016  721 1517 46.5 D  41.342  44.053  5.0F SC3 20 1.6 4.2LSC3 4.6bSC3
 2016  721 1628 47.6 D -23.873-175.548 10.0F SC3 43 1.8 6.1LSC3 5.2bSC3 5.9BSC3
 2016  721 2309  3.2 D  40.770-124.070 10.0F SC3 33 2.0 4.9LSC3 5.0bSC3 5.1BSC3
 2016  722 0040 27.1 D  39.365  25.988  7.3  SC3 25 1.4 3.0LSC3
 2016  722 1131 49.4 D  19.911 120.895 10.0FSSC3 17 1.7 4.7bSC3 5.3BSC3
 2016  722 2140 11.9 D  42.343  35.206 22.2  SC3 42 1.3 4.4LSC3-4.1SSC3
 2016  722 2344 58.5 D  36.951  26.683 10.5  SC3 13 0.9 2.7LSC3
 2016  723 0100 19.3 D  47.772 147.040385.8  SC3214 0.8 7.3LSC3 5.4bSC3 5.6BSC3
 2016  723 0226 20.3 D  -7.409 128.924135.3  SC3 19 1.1 5.6LSC3 4.8bSC3 5.0BSC3
 2016  723 0252 12.1 D  30.140 142.097 61.9  SC3 14 0.7 5.1LSC3 4.4bSC3 4.9BSC3
 2016  723 0741 45.1 D  -6.701 101.262 10.0F SC3 83 1.5 6.1LSC3 5.2bSC3 5.9BSC3
  End of s-file

 Min difference to mark, lat-lon, depth, origin, rms, mag
 Defaults (enter) are 5km, 5km, 1s, 0.5s, 0.2s

 Min number of stations to compare, def 1

           Origin time RMS     Lat      Lon Depth          Ml   Mb   MB  
   Average diff    0.4 0.0  -0.005   -0.001   0.5        -0.1  0.1 -0.0
   Standard dev    1.0 0.1   0.033    0.010   4.5         0.1  0.3  0.2
   Number of values                      11     6          10    8    7


 Output file name is compare_hyp.out

All the input headers are listed for the two files and at the end a summary of the differences is given. In the example above, note that there were 11 events. Five events had fixed depth so only 6 were used for depth comparison. Also note that not all events had all the magnitude types. This example is comparing locations done with SeisComp3 and SEISAN. Below is shown the output file.

    #      Origin time RMS     Lat      Lon Depth    N     Ml   Mb   MB  
    1 20160721151746.7 1.5  41.340   44.040   5.0   20    4.2  4.3              
                   0.2-0.1  -0.002   -0.013               0.0 -0.3*             
    2 20160721162847.7 1.9 -23.860 -175.560  10.0   43    6.0  5.3  5.8         
                   0.1 0.1   0.013   -0.012              -0.1  0.1 -0.1         
    3 20160721230903.5 2.0  40.780 -124.070  10.0   33    4.9  4.7  4.9         
                   0.3 0.0   0.010    0.000               0.0 -0.3*-0.2         
    4 20160722004026.3 1.4  39.360   25.970   7.0   34    3.1                   
                  -0.8 0.0  -0.005   -0.018  -0.3         0.1                   
    5 20160722113149.5 1.8  19.900  120.900  10.0   17         4.7  5.1         
                   0.1 0.1  -0.011    0.005                    0.0 -0.2*        
    6 20160722214010.2 1.0  42.400   35.210  23.0   49    4.4            
                  -1.7*0.3   0.057*   0.004   0.8         0.0            
    7 20160722234457.8 1.0  36.940   26.680   8.0   15    2.7                   
                  -0.7 0.1  -0.011   -0.003  -2.5         0.0                   
    8 20160723010020.0 1.0  47.680  147.050 395.0  214    6.9  5.7  5.7         
                   0.7 0.2  -0.092*   0.010   9.2*       -0.4* 0.3* 0.1         
    9 20160723022620.3 1.1  -7.410  128.920 132.0   19    5.4  5.0  5.3         
                   0.0 0.0  -0.001   -0.004  -3.3        -0.2  0.2  0.3*        
   10 20160723025212.1 0.7  30.140  142.110  61.0   14    4.9  4.8  4.9         
                   0.0 0.0   0.000    0.013  -0.9        -0.2  0.4* 0.0         
   11 20160723074145.1 1.5  -6.710  101.270  10.0   83    6.0  5.3  5.8         
                   0.0 0.0  -0.009    0.008              -0.1  0.1 -0.1         

           Origin time RMS     Lat      Lon Depth          Ml   Mb   MB   MS    
   Average diff    0.4 0.0  -0.005   -0.001   0.5        -0.1  0.1 -0.0         
   Standard dev    1.0 0.1   0.033    0.010   4.5         0.1  0.3  0.2         
   Number of values                      11     6          10    8    7

For each event, the difference in value is shown. Note that the difference value has a star behind it if it exceeds the user specified values. The event numbers are also listed.


next up previous contents index
Next: 29. Making synthetic seismograms Up: 28. Merge events in Previous: 28.3 Catalogue merging and   Contents   Index
Peter Voss : Wed Aug 9 08:16:56 UTC 2017