output

这是一个TOMO TO GMT的源码

   ----  FATOMO Ver.08/11/99  ---
  start at Fri Aug 20 09:49:02 1999
  
  
  FATOMO performs an inversion of Vp-velocities
   using fat rays. Up to now earthquake relocation
   has to be done external with SIMULPS. Therefore,
   the new velocity model after inversion is printed
   in the file "f88" for direct use in SIMULPS.
  
 FATOMO has been set up to use: 
  max. no. of events :   800
  max. no. of obs per event :   50
  max. no. of stations :   60
  
 control parameters: 

     origin (lower left corner) of numerical grid
      relative to short distance origin:
          ox =  80.00 km  oy =   0.00 km
     elevation of top of the model:
          baselev =   4.00 km

     gridspacing for FD modelling
          h =  0.50  km

     number of earthquakes
          neqs =   50

     doing Vp velocity inversion
          nitmax =  1

     damping of velocity parameters:
          vthet =   10.0

     fat ray width (s):
          delta_t =  0.07 s

     full resolution matrix will be computed:
          ires =  2

     output parameters for ray points and fat rays:
          rout = 0  fout = 0

     short distance origin:   latitude  longitude
           38 30           -21   0

     short distance conversion factors
          one min lat   1.8501 km
          one min lon   1.4539 km

 station   latitude   longitude   elev   x    y   z                pdl s-pdl  nfixst
    1   R001  38N35.39    20E53.11    100   69.99   9.97   3.90 0.00 0.00  1
    2   R002  38N46.17    20E53.11    100   70.00  29.92   3.90 0.00 0.00  1
    3   R003  38N51.56    20E53.11    100   70.01  39.89   3.90 0.00 0.00  1
    4   R004  38N56.95    20E53.11    100   70.01  49.86   3.90 0.00 0.00  1
    5   R005  39N 7.73    20E53.11    100   70.03  69.80   3.90 0.00 0.00  1
    6   R006  38N35.39    20E39.34    100   49.98   9.97   3.90 0.00 0.00  1
    7   R007  38N46.17    20E39.34    100   50.02  29.92   3.90 0.00 0.00  1
    8   R008  38N51.56    20E39.34    100   50.04  39.89   3.90 0.00 0.00  1
    9   R009  38N56.95    20E39.34    100   50.06  49.86   3.90 0.00 0.00  1
   10   R010  39N 7.73    20E39.34    100   50.09  69.80   3.90 0.00 0.00  1
   11   R011  38N35.39    20E32.45    100   39.97   9.97   3.90 0.00 0.00  1
   12   R012  38N46.17    20E32.45    100   40.02  29.92   3.90 0.00 0.00  1
   13   R013  38N51.56    20E32.45    100   40.04  39.89   3.90 0.00 0.00  1
   14   R014  38N56.95    20E32.45    100   40.07  49.86   3.90 0.00 0.00  1
   15   R015  39N 7.73    20E32.45    100   40.12  69.80   3.90 0.00 0.00  1
   16   R016  38N35.39    20E25.56    100   29.96   9.97   3.90 0.00 0.00  1
   17   R017  38N46.17    20E25.56    100   30.02  29.92   3.90 0.00 0.00  1
   18   R018  38N51.56    20E25.56    100   30.05  39.89   3.90 0.00 0.00  1
   19   R019  38N56.95    20E25.56    100   30.08  49.86   3.90 0.00 0.00  1
   20   R020  39N 7.73    20E25.56    100   30.15  69.80   3.90 0.00 0.00  1
   21   R021  38N35.39    20E11.79    100    9.95   9.97   3.90 0.00 0.00  1
   22   R022  38N46.17    20E11.79    100   10.04  29.92   3.90 0.00 0.00  1
   23   R023  38N51.56    20E11.79    100   10.08  39.89   3.90 0.00 0.00  1
   24   R024  38N56.95    20E11.79    100   10.13  49.86   3.90 0.00 0.00  1
   25   R025  39N 7.73    20E11.79    100   10.21  69.80   3.90 0.00 0.00  1
  
 INVERSION GRID

   inversion grid size (no. of edges):
   inv_nx =  15 inv_ny =  15 inv_nz =  12

   no. of inversion blocks :
   in x-dir. = 14
   in y-dir. = 14
   in z-dir. = 11


   no. of parameter to be inverted :  2156

 inv. block boundaries :

  xgrid
  80.0  72.5  62.5  57.5  52.5  47.5  42.5  37.5  32.5  27.5  22.5  17.5  12.5   5.0   0.0

  ygrid
   0.0   5.0  12.5  17.5  22.5  27.5  32.5  37.5  42.5  47.5  52.5  57.5  62.5  72.5  80.0

  zgrid
  -4.0   1.0   3.0   4.5   7.5  12.5  17.5  22.5  27.5  32.5  37.5  40.0
  
  SEISMIC GRID (gridnodes) 

   no. of gridnodes:
   sei_nx =  14 sei_ny =  14 sei_nz =  11

 position of grid nodes:

  xgrid
  80.0  65.0  60.0  55.0  50.0  45.0  40.0  35.0  30.0  25.0  20.0  15.0  10.0   0.0

  ygrid
   0.0  10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0  55.0  60.0  65.0  80.0

  zgrid
  -4.0   0.0   2.0   4.0   5.0  10.0  15.0  20.0  25.0  30.0  40.0
  
 First blocks in inv. and seis. layer :
  
 seis-layer   1 no. first block   1 in inv. layer   1
 inv-layer   1 no. of first block   1
 seis-layer   2 no. first block   197 in inv. layer   1
 inv-layer   1 no. of first block   1
 seis-layer   3 no. first block   393 in inv. layer   2
 inv-layer   2 no. of first block   197
 seis-layer   4 no. first block   589 in inv. layer   3
 inv-layer   3 no. of first block   393
 seis-layer   5 no. first block   785 in inv. layer   4
 inv-layer   4 no. of first block   589
 seis-layer   6 no. first block   981 in inv. layer   5
 inv-layer   5 no. of first block   785
 seis-layer   7 no. first block   1177 in inv. layer   6
 inv-layer   6 no. of first block   981
 seis-layer   8 no. first block   1373 in inv. layer   7
 inv-layer   7 no. of first block   1177
 seis-layer   9 no. first block   1569 in inv. layer   8
 inv-layer   8 no. of first block   1373
 seis-layer   10 no. first block   1765 in inv. layer   9
 inv-layer   9 no. of first block   1569
 seis-layer   11 no. first block   1961 in inv. layer   10
 inv-layer   10 no. of first block   1765
  
 initial velocity values on three dimensional grid

layer   1    P-velocity    z =  -4.0
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00

layer   2    P-velocity    z =   0.0
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00
 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00

layer   3    P-velocity    z =   2.0
 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40
 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40
 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40
 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40 3.40