output

这是一个射线追踪的源码

 Computation began at 26-Apr-99 15:04:36
 Program Simulps14 (20-Apr-99 FH)
Added paraxial raytracing.
Simulps13q (15-Mar-98 DMEP).  Solves for Vp and Vp/Vs; Input data is P travel-time and S-P time.
       Can vary relative weighting of S-P times in hypocenter location.

       OR can solve for Qp using tstar as input

       Allows fixed nodes (up to  7350 parameters, up to  5000 solution parameters);
       up to  600 events,   200 stations,  100 observations per event
       Psuedo-bending; Allows less curvature "moho" for initial arcuate paths.

 control parameters
 kout kout2 kout3
   3    0     0

 neqs nsht nbls wtsht nitloc wtsp   zmin   eigtol  rmscut hitct
   50   0   0     0.0   3    0.50   0.00   0.020    0.002    7.

 Only P velocities are included in inversion  (iusep= 1, iuses= 0)

    dvpmx  dvpvsmx  idmp     vpdamp    vpvsdmp    stadamp
     0.40    0.200   0       10.000      2.000      2.000
  ires nitmax snrmct dxmax rderr ercof
  1      1  0.03000  1.25  0.10  0.01
 # its. for 1-d vel. model =   0
 rms for term. = 0.005
 fix locations for iterations =   0
 distance weighting:  50.00 200.00; residual weighting: 0.30 1.00 2.00
 step length for integration: 0.500
 parameters for approximate ray tracer
      ndip iskip scale1 scale2
       9     0   1.00   0.50
  parameters for pseudo-bending:
  i3d xfac   tlim  nitpb(1) nitpb(2)
   4  1.90 0.0010    35        45
 i3d = 4: use ray shooting



  Control parameters for paraxial raytracing
    read from file RAYTRAC:


 iheter, epsob, epsca, ides, ampr, iterrai
   20    1.00000E-03    1.00000E-01  0    1.00000  1
 dxrt, dyrt, dzrt 
     1.00000    1.00000    1.00000

  computed machine constants eta and tol:    0.119209E-06    0.986076E-31

  origin :  latitude   longitude   rotation
             38 30.00   -21  0.00      0.00
     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   10.01   9.97  -0.10 0.00 0.00  0
    2   R002  38N46.17    20E53.11    100   10.00  29.92  -0.10 0.00 0.00  0
    3   R003  38N51.56    20E53.11    100    9.99  39.89  -0.10 0.00 0.00  0
    4   R004  38N56.95    20E53.11    100    9.99  49.86  -0.10 0.00 0.00  0
    5   R005  39N 7.73    20E53.11    100    9.97  69.80  -0.10 0.00 0.00  0
    6   R006  38N35.39    20E39.34    100   30.02   9.97  -0.10 0.00 0.00  0
    7   R007  38N46.17    20E39.34    100   29.98  29.92  -0.10 0.00 0.00  0
    8   R008  38N51.56    20E39.34    100   29.96  39.89  -0.10 0.00 0.00  0
    9   R009  38N56.95    20E39.34    100   29.94  49.86  -0.10 0.00 0.00  0
   10   R010  39N 7.73    20E39.34    100   29.91  69.80  -0.10 0.00 0.00  0
   11   R011  38N35.39    20E32.45    100   40.03   9.97  -0.10 0.00 0.00  0
   12   R012  38N46.17    20E32.45    100   39.98  29.92  -0.10 0.00 0.00  0
   13   R013  38N51.56    20E32.45    100   39.96  39.89  -0.10 0.00 0.00  0
   14   R014  38N56.95    20E32.45    100   39.93  49.86  -0.10 0.00 0.00  0
   15   R015  39N 7.73    20E32.45    100   39.88  69.80  -0.10 0.00 0.00  0
   16   R016  38N35.39    20E25.56    100   50.04   9.97  -0.10 0.00 0.00  0
   17   R017  38N46.17    20E25.56    100   49.98  29.92  -0.10 0.00 0.00  0
   18   R018  38N51.56    20E25.56    100   49.95  39.89  -0.10 0.00 0.00  0
   19   R019  38N56.95    20E25.56    100   49.92  49.86  -0.10 0.00 0.00  0
   20   R020  39N 7.73    20E25.56    100   49.85  69.80  -0.10 0.00 0.00  0
   21   R021  38N35.39    20E11.79    100   70.05   9.97  -0.10 0.00 0.00  0
   22   R022  38N46.17    20E11.79    100   69.96  29.92  -0.10 0.00 0.00  0
   23   R023  38N51.56    20E11.79    100   69.92  39.89  -0.10 0.00 0.00  0
   24   R024  38N56.95    20E11.79    100   69.87  49.86  -0.10 0.00 0.00  0
   25   R025  39N 7.73    20E11.79    100   69.79  69.80  -0.10 0.00 0.00  0


 velocity grid size:
bld = 1.0      nx = 16     ny = 16     nz = 12

 xgrid
     -30.0    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  110.0

 ygrid
     -30.0    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  110.0

 zgrid
      -4.0    0.0    2.0    4.0    5.0   10.0   15.0   20.0   25.0   30.0   40.0   80.0


 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  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00
  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00
  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00
  3.00  3.00  3.00  3.00  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  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00
  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00
  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00  3.00
  3.00  3.00  3.00  3.00  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  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
  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  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
  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

 layer  4     P velocity          z =    4.0
  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80
  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80
  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80
  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80
  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80
  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80
  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80
  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80
  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80
  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80
  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80
  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80
  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80
  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80
  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80
  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80  3.80

 layer  5     P velocity          z =    5.0
  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00
  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00
  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00
  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00
  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00
  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00  4.00