resol_main.out

这是一个TOMO TO GMT的源码

 



  PROGRAM RESOL:,
  ==============

 Version 6.0 of 20. February 1996 by E.Kissling
 
  file name of print output:
resolantofo.out                                                                 
 
  file name of binary model input:
antofo.bin                                                                      
 
  file name of input earthquake data:
antofo3d.cnv                                                                    
 
  file name of steplength weight output:
antofo.slw                                                                      
 
  file name of binary model output:
antofotomo.cnt                                                                  
 
  file name of selected events output:
antofotomo.eve                                                                  
 
 
 >>>>> Start of program RESOL <<<<<
 
Antofogasta Preliminary Tomo SHEK97                                             
 
 



  INPUT - MODEL AS READ BY SUBR. INPUT10:



  CHILE ANTOFOGASTA PRELIMINARY 3D MODEL (SHEK0897)                              

  NUMBER OF VELO-MODELS(NMOD)=    1
     NUMBER OF LAYERS(NL)= 16


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





 
 
 SHORT DISTANCE CONVERSION on ELLIPSOIDE of WORLD GEODETIC SYSTEM 1972 (WGS72)
 =============================================================================
 
 Radius at equator (REARTH)= 6378.13500  km
   1. / (ellipticity)      =    298.260
 
 Origin of cartesian coordinates [degrees]:
  -23.5000000 N       70.0000000 W
 
  without rotation of grid, 
              pos. x-axis toward WEST
              pos. y-axis toward NORTH
 
  Rotation of y-axis from North anticlock-
  wise with pos. angle given in degrees
 
  Rotation of grid anticlockwise by
     0. degrees
 
 Radius of sphere at OLAT =   6374.773 km
 
 Conversion of GEOGRAPHICAL LATITUDE to GEOCENTRICAL LATITUDE:
 RLATC = TAN(GEOCENTR.LAT) / TAN(GEOGRAPH.LAT)
 RLATC =   0.99330157
 
          Short distance conversions:
          one min lat ~  1.8458 km 
          one min lon ~  1.7024 km 

 
 

  CNTRL DATA FOR EVENT CHECK:
     RMS1= 5.00  RMS2= 5.00     AMAG1= 0.00  AMAG2= 0.00     IGAP1= 180  IGAP2= 180    IOBSM= 8   MINW= 3

     DM1=15.00  DM2=15.00     ERH1= 1.00  ERH2= 1.00     ERZ1= 2.00  ERZ2= 2.00


  KFIRST= 0   NBL= 11408      NEQS=   603      NSHOT=   0    NMIN=   1    NSTA= 57

  DELTAT=   0.10SEC  AA=    1.8458  BB=    1.8543  RPHI=   0.00     VTHET=  10.000



  VELOCITY MODEL:
  ===============


  NUMBER OF VELO-MODELS(NMOD)=    1
     NUMBER OF LAYERS(NL)= 16


  IGEO=  1          IPLOT=  1

  NR STATION      P-DELAY  ISTA LSTA  X-COOR.    Y-COOR.   Z-COOR.

    1   MJN       0.000     0    1     51.545     44.112     -0.500
    2   MJC       0.000     0    1     54.403     23.689     -0.150
    3   MJS       0.000     0    1     54.478      3.356     -0.250
    4   MJP       0.000     0    1     32.695     36.027     -0.150
    5   CAC       0.000     0    1     31.211     55.464     -0.100
    6   MEX       0.000     0    1      6.420     23.845     -1.200
    7   DES       0.000     0    1     33.695     16.801     -0.450
    8   PRT       0.000     0    1     23.016      4.773     -0.851
    9   ERC       0.000     0    1    -12.866     39.903     -1.368
   10   AME       0.000     0    1    -24.669    -46.150     -0.850
   11   NAV       0.000     0    1      4.510    -30.091     -0.675
   12   PA1       0.000     0    1     29.571    -49.096     -0.646
   13   COL       0.000     0    1     41.720    -31.630     -0.311
   14   QC1       0.000     0    1     34.906    -13.722     -0.535
   15   PEL       0.000     0    1    -13.210    -16.424     -0.990
   16   TRO       0.000     0    1    -14.304      8.617     -0.940
   17   MIC       0.000     0    1     16.647     88.445     -0.900
   18   NEU       0.000     0    1     -4.016     90.328     -1.300
   19   VAL       0.000     0    1    -29.089     94.415     -1.500
   20   ELE       0.000     0    1    -35.715    134.374     -1.200
   21   NEG       0.000     0    1    -64.364    135.714     -1.450
   22   SOM       0.000     0    1    -62.694    111.659     -1.350
   23   PIM       0.000     0    1    -65.357     68.544     -1.660
   24   GAT       0.000     0    1     22.044    113.674     -0.370
   25   TOC       0.000     0    1     11.917    149.358     -1.250
   26   CER       0.000     0    1    -18.756    136.743     -1.450
   27   CRP       0.000     0    1    -16.468    117.240     -1.830
   28   GRA       0.000     0    1      8.943     52.695     -1.100
   29   GUA       0.000     0    1    -38.156     52.374     -1.360
   30   PUN       0.000     0    1    -59.316     42.417     -1.550
   31   VEN       0.000     0    1    -39.558     14.497     -1.400
   32   AMA       0.000     0    1    -62.016     -6.711     -1.724
   33   PSR       0.000     0    1    -50.778    -44.965     -1.300
   34   PA2       0.000     0    1     25.776    -51.898     -0.877
   35   SUR       0.000     0    1     41.475    -82.542     -1.605
   36   QC2       0.000     0    1     36.216    -10.887     -1.713
   37   TET       0.000     0    1      5.572    -81.490     -1.250
   38   MUL       0.000     0    1     45.394    -57.712     -0.900
   39   PRE       0.000     0    1     33.077   -104.748     -1.950
   40   PAS       0.000     0    1     37.248   -137.364     -1.900
   41   PRI       0.000     0    1    -15.350    -61.312     -1.036
   42   LCA       0.000     0    1    -19.047   -113.121     -1.550
   43   MTO       0.000     0    1     28.985     80.815     -0.250
   44   LDO       0.000     0    1    -12.142   -152.481     -2.190
   45   BSJ       0.000     0    1     42.469   -167.765     -0.800
   46   ALG       0.000     0    1    -19.822     70.947     -1.380
   47   OFE       0.000     0    1    -10.476    -67.115     -1.100
   48   OSG       0.000     0    1    -21.575    -79.187     -1.046
   49   LTO       0.000     0    1     -1.240   -148.893     -2.130
   50   ZAA       0.000     0    1    120.521    -82.941      5.055
   51   ZBA       0.000     0    1     95.200    -64.590      2.238
   52   ZCA       0.000     0    1     71.377    -41.587      1.501
   53   ZDA       0.000     0    1     91.900    -18.473      3.124
   54   ZEA       0.000     0    1    114.222      7.464      5.218
   55   ZFA       0.000     0    1     88.701     27.688      3.260
   56   ZHA       0.000     0    1     58.071     64.490      0.696
   57   ZIA       0.000     0    1     86.523     73.550      3.238

  LSMAX=    1


  NUMBER OF BLOCKS(STARTING NETWORK) IN:  X-DIR.(NX)=   21     Y-DIR.(NY)=   34

  X-DIR. BLOCK-BOUNDARIES (KM)
   130.0   120.0   110.0   100.0    90.0    80.0    70.0    60.0    50.0    40.0
    30.0    20.0    10.0     0.0   -10.0   -20.0   -30.0   -40.0   -50.0   -60.0
   -70.0   -80.0
  Y-DIR. BLOCK-BOUNDARIES (KM)
   160.0   150.0   140.0   130.0   120.0   110.0   100.0    90.0    80.0    70.0
    60.0    50.0    40.0    30.0    20.0    10.0     0.0   -10.0   -20.0   -30.0
   -40.0   -50.0   -60.0   -70.0   -80.0   -90.0  -100.0  -110.0  -120.0  -130.0
  -140.0  -150.0  -160.0  -170.0  -180.0


  DXYMIN=     10.000



  X-DIR. CENTERS OF BLOCKS  (KM)

   125.0   115.0   105.0    95.0    85.0    75.0    65.0    55.0    45.0    35.0
    25.0    15.0     5.0    -5.0   -15.0   -25.0   -35.0   -45.0   -55.0   -65.0
   -75.0

  Y-DIR. CENTERS OF BLOCKS  (KM)

   155.0   145.0   135.0   125.0   115.0   105.0    95.0    85.0    75.0    65.0
    55.0    45.0    35.0    25.0    15.0     5.0    -5.0   -15.0   -25.0   -35.0
   -45.0   -55.0   -65.0   -75.0   -85.0   -95.0  -105.0  -115.0  -125.0  -135.0
  -145.0  -155.0  -165.0  -175.0

  BLOCKNUMBER OF FIRST BLOCK IN EACH LAYER,NNB(L):
 NUMBERING OF BLOCKS IN FINAL MODEL, VP-ARRAY!

     1   714  1427  2140  2853  3566  4279  4992  5705  6418
  7131  7844  8557  9270  9983 10696 11409

  NUMBERS OF BLOCKS USED FOR INVERSION (NBINV(K)):

  NUMBER OF BLOCK IN UPPER LEFT CORNER FOR EACH INVERSION-LAYER:
  BLOCKNR.   INV.LAYERNR.

         1         1
       714         2
      1427         3
      2140         4
      2853         5
      3566         6
      4279         7
      4992         8
      5705         9


  NBINV(NLINV+1)=   6418          TOTAL NR. OF INVERSION BLOCKS; NBMAX=   6417
     TOT. NR. OF INV. LAYER; NLINV=    9


 MATRIX VELOCITY-MODEL  AREA=FULL AREA UNDER STUDY TO BEGIN WITH

  SEIS-LAYER     INV-LAYER     VP(KM/S)  DEPTH(KM)  THICKNESS(KM)