resol_main.out

这是一个TOMO TO GMT的源码

95 924 2211  8.24 23.5717S  70.4727W  35.36 24     48.27   -7.94   35.36
     event nr=   521   obs nr= 10 residuum =    1.13
95 923 2216 40.85 23.5642S  70.2336W  45.17 21     23.85   -7.11   45.17
95 923 2234 51.75 23.9540S  70.4391W  35.04 15     44.77  -50.28   35.04
95 925 2352 38.63 23.3749S  70.3961W  37.83 23     40.48   13.85   37.83
95 925 2354 45.65 23.1886S  70.4014W  34.91 17     41.05   34.49   34.91
95 926  028 49.92 23.9383S  70.3655W  38.43 19     37.27  -48.54   38.43
95 926  052 38.67 23.3746S  70.3948W  37.42 22     40.34   13.89   37.42
95 926  054 45.67 23.1868S  70.4006W  34.75 18     40.97   34.69   34.75
95 926  128 49.90 23.9362S  70.3642W  38.58 18     37.14  -48.31   38.58
95 926  4 9 57.08 23.1992S  70.1753W  44.78 26     17.93   33.31   44.78
95 925  6 4  1.76 23.2085S  70.1793W  42.68 14     18.33   32.28   42.68
95 925  714 21.33 23.4199S  70.5333W  32.34 32     54.49    8.87   32.34
95 926  724 59.95 23.7658S  70.2886W  41.73 25     29.45  -29.44   41.73
95 926  8 4 25.65 23.2437S  70.4540W  35.22 28     46.42   28.39   35.22
95 926  831 10.00 24.0993S  70.4241W  38.80 19     43.22  -66.37   38.80
95 925  910 39.12 23.9733S  70.1139W  47.84 17     11.61  -52.42   47.84
95 925  947 47.37 23.2375S  70.4651W  34.28 12     47.55   29.07   34.28
95 925 10 6 21.80 23.9070S  70.6023W  29.81 31     61.43  -45.08   29.81
95 925 1050 44.19 23.2607S  70.5249W  31.17 18     53.66   26.50   31.17
95 925 1153 28.06 23.7095S  70.2953W  40.89 24     30.14  -23.20   40.89
95 926 1231  4.90 24.4951S  70.6095W  36.74 13     62.02 -110.21   36.74
95 926 13 0 50.86 23.7212S  70.2869W  42.30 12     29.28  -24.50   42.30
95 925 13 1 24.01 23.0928S  69.8551W  39.61 14    -14.82   45.10   39.61
95 926 1331 34.66 24.1494S  70.4142W  37.90 11     42.20  -71.92   37.90
95 926 1421  5.32 23.4100S  70.4456W  35.43 10     45.53    9.97   35.43
95 926 1434 42.79 24.0081S  70.3874W  40.46 10     39.49  -56.27   40.46
95 925 1435 18.46 23.2337S  70.5101W  32.58 14     52.15   29.49   32.58
95 925 1536 52.01 23.2466S  70.4313W  35.79 13     44.10   28.06   35.79
95 925 1542  2.17 24.2055S  70.7786W  20.99 15     79.32  -78.13   20.99
95 925 16 2 35.38 23.2399S  70.4527W  35.34 29     46.28   28.81   35.34
95 926 1638 40.91 23.6141S  70.2972W  44.90 10     30.34  -12.64   44.90
95 925 1716 11.52 24.1733S  70.5972W  34.68 22     60.84  -74.57   34.68
95 925 1821 52.50 23.5359S  70.5181W  34.72 20     52.91   -3.98   34.72
95 926 1854 14.75 23.3865S  70.5068W  32.29 11     51.79   12.57   32.29
95 925 19 5  5.76 23.4258S  70.4514W  34.37 21     46.12    8.22   34.37
95 925 1916  5.77 24.0588S  70.3439W  40.57 21     35.05  -61.89   40.57
95 926 1950 24.51 23.6040S  70.4446W  38.97 28     45.39  -11.52   38.97
95 926 2038 48.73 23.3882S  70.3655W  38.99 14     37.35   12.38   38.99
95 926 21 9 13.84 23.1760S  70.1013W  47.71 14     10.36   35.88   47.71
95 926 2112 46.81 23.1789S  70.1028W  45.49 10     10.51   35.56   45.49
95 926 2144 55.02 23.6588S  70.4231W  48.22 20     43.19  -17.59   48.22
95 925 22 1 49.89 23.3756S  70.3974W  36.92 19     40.61   13.78   36.92
95 925 2232 10.17 24.0390S  70.6150W  30.53 15     62.69  -59.69   30.53
95 926 2333  3.10 24.0086S  70.2080W  48.14 11     21.20  -56.33   48.14
95 925 2346 10.91 23.5970S  70.3496W  41.00 19     35.70  -10.74   41.00
95 927 2331 35.73 24.0807S  70.4077W  41.13 22     41.55  -64.31   41.13
95 928  0 1 12.72 23.1866S  70.4048W  36.96 15     41.40   34.71   36.96
95 928  036 41.49 23.4829S  70.3605W  41.07 22     36.82    1.89   41.07
95 927  1 4 47.46 23.9697S  70.3795W  41.36 10     38.69  -52.02   41.36
95 928  128 49.72 23.6017S  70.4080W  40.06 22     41.66  -11.26   40.06
95 927  143 25.64 23.8530S  70.3308W  40.75 10     33.74  -39.10   40.75
95 928  219 50.22 23.3378S  70.5160W  31.43 16     52.74   17.96   31.43
95 927  247 29.79 23.9880S  70.2335W  46.35 28     23.81  -54.05   46.35
95 927  326 34.61 23.3673S  70.3526W  40.87 28     36.03   14.70   40.87
95 927  333  2.70 23.9022S  70.1710W  49.29 16     17.44  -44.54   49.29
95 927  335 29.13 23.3727S  70.3553W  39.22 10     36.31   14.10   39.22
95 928  350 22.63 24.0160S  70.4540W  37.70 24     46.28  -57.15   37.70
95 928  457 22.50 23.8471S  70.3337W  39.95 21     34.04  -38.44   39.95
95 928  523  7.54 23.6893S  70.4066W  38.52 22     41.50  -20.97   38.52
95 928  549 20.05 23.2641S  69.6300W  65.39 27    -37.83   26.13   65.39
95 927  655 46.77 23.4966S  70.3474W  41.42 20     35.48    0.38   41.42
95 928  7 3 22.51 24.1463S  69.9791W  53.59 10     -2.13  -71.58   53.59
95 927  742 50.98 23.4499S  70.4784W  36.73 29     48.87    5.55   36.73
95 927  812 58.98 23.5150S  70.4076W  39.36 13     41.63   -1.66   39.36
95 928  828 39.97 23.8673S  70.3271W  37.41 13     33.36  -40.68   37.41
95 928  9 5 59.69 23.3986S  70.3983W  36.38 18     40.70   11.23   36.38
95 928  929 14.46 23.2405S  70.4525W  36.90 32     46.26   28.74   36.90
95 928  933  2.90 23.2418S  70.4508W  35.30 29     46.09   28.60   35.30
95 928 11 9 17.67 23.3745S  70.3510W  38.87 26     35.87   13.90   38.87
95 928 1140 34.07 23.2671S  70.4353W  35.28 17     44.50   25.79   35.28
95 927 12 9 13.87 23.4753S  70.5270W  32.67 10     53.83    2.74   32.67
95 928 1218 50.56 23.9969S  70.3808W  37.77 20     38.82  -55.03   37.77
95 927 1224 45.40 23.4123S  70.4295W  34.21 10     43.88    9.71   34.21
95 928 1337 29.91 24.1237S  70.4142W  38.99 18     42.21  -69.08   38.99
95 927 1444  3.86 23.4074S  70.4483W  36.20 17     45.81   10.26   36.20
95 927 1612 57.02 23.1469S  70.0294W  52.94 11      3.01   39.11   52.94
95 927 1640 46.01 24.0199S  70.2597W  43.78 13     26.47  -57.58   43.78
95 928 1739  6.77 23.6980S  70.2441W  45.36 14     24.91  -21.93   45.36
95 927 1833 34.93 23.4794S  70.3585W  40.19 14     36.62    2.28   40.19
95 927 1934 27.52 22.7647S  69.7624W  57.34 12    -24.34   81.43   57.34
95 927 2053 55.83 23.8800S  70.1838W  46.08 16     18.75  -42.09   46.08
95 928 2146 18.06 23.8537S  70.4543W  33.12 27     46.34  -39.17   33.12
95 927 2237 38.00 23.2965S  70.4274W  36.06 23     43.69   22.54   36.06
  subr. wricnt: writing on *.cnt file:

  olat=    -23.5000  (degr.)  olon=     70.0000  (degr.)
     icoordsystem=  1
   rota=   0.00 degr.





  end of subr. wricnt 



  DATA - STATISTICS COMING OUT OF PROGRAM RESOL:
*************************************************

      EQUAL DATA STATISTICS GOING INTO GEOTOMOGRAPHY


  TOTAL NUMBER OF OBSERVATIONS (OUT OF KNOBSTO) NEGLECTED FOR VELOCITY INVERSION: NOBNE=       0
      REASON: SEE SWITCHES OF ISTA(K) FOR STATIONS



   ATTENTION: AVERAGE RESIDUAL (FOR EVENT) NOT REMOVED IN
  THE FOLLOWING VARIANCE CALCULATIONS, AS
  IT WILL BE DONE IN PROGRAM TOMOGR!




   also check if station corrections that were
   used for locating the earthquakes have been used in this
   run of RESOL. if NOT, the variance must be much larger
   than the variance for the same data calculated by RELOC.
   Since RESOL is the first program used for the 3D-inver-
   sion, where the station residuals will be back projected
   the normal case of usage of RESOL will be without
   station corrections, and consequently larger variances!






      VARIANCES CALCULATIONS WITHOUT ANY WEIGTHS
     --------------------------------------------------

  DATA-VARIANCES CALCULATED AS FOR 1D-MODEL,
   DAVAR=(DAVAR-TRES**2/KNOBSTO)/(KNOBSTO-NL-NSTA)

KNOBSTO=     12906     NL= 16     NSTA= 57

     knobsto = tot. nr. of observation = sum of all used either
     for hypo.-loc. only or for model and hypo. inversion.

  DATA VAR.=   0.07048          AVREAGE RES.=     0.010


  DATA VARIANCE AS FOR 1D-MODEL :   DAVAR= (DAVAR-TRES**2/KNOBSTO)/KNOBSTO
     DAVAR=   0.07008          STAND. DEV. (SQRT(DAVAR)) =  0.264729





      VARIANCES CALCULATIONS WITH READING WEIGTHS
     -----------------------------------------------

  DATA-VARIANCES CALCULATED AS FOR 1D-MODEL,
   DAVAR=(DAVAR-TRES1**2*(2*SW1TO-KNOBSTO)/(SW1TO-NL-NSTA)
  SUM OF ALL WEIGHTS (SW1TO)= 12906.000     NL= 16     NSTA= 57
  DATA VAR.=   0.06499          AVREAGE RES.=    -0.009


  DATA VARIANCE AS FOR 1D-MODEL :   DAVAR= (DAVAR-TRES1**2*(2*SW1TO-KNOBSTO)/SW1TO
     DAVAR=   0.06463          STAND. DEV. (SQRT(DAVAR)) =  0.254219





      VARIANCES CALCULATIONS WITH GEOMETRIC WEIGTHS
     -----------------------------------------------

  DATA-VARIANCES CALCULATED AS FOR 1D-MODEL,
   DAVAR=(DAVAR-TRES2**2*(2*SW2TO-KNOBSTO)/(SW2TO-NL-NSTA)
  SUM OF ALL WEIGHTS (SW2TO)=  8531.987     NL= 16     NSTA= 57
  DATA VAR.=      0.0621          AVREAGE RES.=    -0.002


  DATA VARIANCE AS FOR 1D-MODEL :   DAVAR= (DAVAR-TRES2**2*(2*SW2TO-KNOBSTO)/SW2TO
     DAVAR=      0.0615          STAND. DEV. (SQRT(DAVAR)) =  0.248042


 NUMBER OF EVENTS (NEV) =  603




  TOTAL NUMBER OF HITS PER INVERSION LAYER:
     LAYER     NR.HITS   NR.BLOCKS HIT   NR.OF BLOCKS (NMIN)

        1        24818      227              191
        2        25719      365              307
        3        45356      410              360
        4        42850      376              304
        5        24289      342              253
        6         7768      318              236
        7         1510      227              125
        8          181       21               10
        9            0        0                0


   TOTAL NUMBER OF INV-BLOCKS HIT BY AT LEAST NMIN RAYS: NMIB=  1786




 reminder:  NMIN was set in *.cnt-file as:     4

1 OBSERVATION MATRIX  (COMBINED BLOCKS ONLY)



      NHIT -MATRIX : (NMIN=   4)






 LAYER    1


 SEIS-LAYER NR:  1     INV-LAYER NR:  1
     NR OF FIRST SEIS-BLOCK:      1          NR. OF LAST SEIS-BLOCK:    713
 NR.OF FIRST INV-BLOCK:      1          NR.OF LAST INV-BLOCK:    713
          TOP OF LAYER=    -2.50 KM, LAYER THICKNESS (KM)=      7.80


         1     2     3     4     5     6     7     8     9    10    11    12    13    14    15    16    17    18    19    20

   1