tramp.doc

利用fortran写的实现2D射线追踪的源代码！运行于linux下！

   path (default: 0)

itx -  plot the calculated travel time-distance  values   as lines
   (1) or symbols of height symht (2); no calculated travel times 
   are plotted if itx=0 (default: 1)

idata - plot the  observed travel times  from the file  tx.in with 
   the calculated travel times; the observed data are represented 
   by vertical  bars with  a  height of  twice  the corresponding  
   uncertainty included in the file tx.in (default: 0)
 
isep - plot the model and rays and travel times in the same plot
   frame (0), separate plot frames (1), or plot the model and rays 
   in separate plot frames for each shot followed by the travel times 
   in a single plot frame all shots (3) (default: 0)
   Note: the value of isep may be changed during a program run with
         entry from the keyboard of the values 0, 1, or 3 between
         plot frames when a <CR> is needed to continue; entering "s"
         will terminate the program run

iszero - plot the calculated and observed  travel times as if all 
   shot points were at 0 km and  the rays traveled left to right  
   (default: 0)

itxout - write the  calculated travel time-distance  pairs to the  
   file tx.out;  the  integer  code  identifying  each  pick  is  
   determined by the array ibreak; itx=1 or 2 must be used
   (default: 0)

iamout -  write the  calculated amplitude-distance  pairs  to the  
   file amp.out;  the  integer  code  identifying  each  pick is  
   determined by the array ibreak; iamp=2 must be used if iamout=1 
   is used (default: 0)

ivz -  plot velocity-depth  profiles, ivz=1  for  P-wave velocity  
   profiles, ivz=-1 for S-wave velocity profiles, ivz=2 for both  
   P- and S-wave  velocity profiles;  the velocity  profiles are  
   output to the file m.out if idump=1 (default: 0)

ivzp - plot  all velocity-depth profiles  in a single  plot frame  
   (default: 1)

ivrms -  plot  the RMS  velocity  of the  model  versus distance,  
   ivrms=1 for P-wave velocity, ivrms=-1 fro S-wave velocity; the 
   RMS velocity profile is  output to the  file m.out if  idump=1 
   (default:: 0)

igrid - output the  velocity model sampled  on a uniform  grid to 
   the file m.out; igrid=1 for P-wave velocity, igrid=2 for P-wave 
   velocity gradient, igrid=-1 for S-wave velocity, igrid=-2 for  
   S-wave  velocity  gradient,  igrid=-3   for  Poisson's  ratio   
   (default: 0)

idump  -   write  detailed   velocity  model   and   ray  tracing   
   information to  the file  r2.out  and the  namelist parameter  
   values to the file n.out; mainly for debugging (default: 0)
   Note: using idump=1 will greatly slow down the run time

isum -  write  a short  summary  about the  rays  traced  and the  
   velocity model to the screen and to the file r1.out (default:  
   1)

istep - pause after each ray is plotted writing its take-off angle
        to the screen and waiting until a <CR> is entered before
        plotting the next ray; this is intended to be used in
        conjunction with irays=1 to determine the cause of the
        failure of the search mode for a particular ray group;
        entering "s" will terminate the program run (default: 0)

b) Other plotting parameters:

symht - height of symbols (mm)  used when plotting each point 
   defining a  ray path,  single  point travel  time  curves (if  
   symht=0, single point travel time curves will not be plotted), 
   and the calculated travel times if itx=2 (default: 0.5)

velht - the maximum  height of the velocity  values (mm) when  
   ivel=1 or -1 is equal to velht*albht (default: 0.7)

nskip - the first nskip points  (step lengths) of all ray paths  are 
   not plotted (default: 0)

nrskip - plot every nrskip ray traced (default: 1)

vred -  reducing  velocity  (km/s)  of  travel  time-distance  plot;  
   observed data in the file  tx.in is assumed to be unreduced (a
   value  of zero is  permitted for no  reduction; vred must be less 
   than 10) (default: 8)

ttunc - uncertainty  of the  calculated travel times  (s) written  
   to the file tx.out if itxout=1 (default: 0.01)

ampunc - relative  uncertainty of  the calculated  amplitudes (%)  
   written to the file am.out if iamout=1 (default: 10)

vcntr - an  array containing  the isovelocity contours  (km/s) to  
   be plotted  if  icntr=1  or  -1; if  a  negative  velocity  is 
   specified, the contour corresponding to the absolute value will 
   be plotted but not labelled

xcinc -  the  spacing (km)  of  points at  which  the isovelocity  
   contours are calculated (default: (xmax-xmin)/100)

xclab - an  array containing  the model  distances (km)  at which  
   the   isovelocity    contours    are    labelled    (default:    
   (xmin+xmax)/2, for the first element only)

ncsmth - number of applications of  a 3-point smoothing filter to 
   the isovelocity contours (default: 0)

xvz -  an  array  containing  the  x-coordinates  (km)  at  which  
   velocity-depth  profiles  are  plotted  if  ivz=1,  -1  or  2   
   (default: (xmin+xmax)/2, for the first element only)

vrmstl, vrmsbl  - the  top and  bottom layer,  respectively, over  
   which the RMS velocity of the model is calculated if ivrms=1 or 
   -1 (default: 1, number of model layers)

xrinc - the spacing (km)  of points at which the  RMS velocity is 
   calculated (default: (xmaxt-xmint)/100)

nrsmth - number of applications of  a 3-point smoothing filter to 
   the RMS velocity versus distance curve (default: 0)

xgrid, zgrid - the distance  and depth spacing (km)  of points at 
   which the velocity model is sampled  if igrid=1, 2, -1, -2, 3  
   (default:: (xmax-xmin)/100, (zmax-zmin)/50)

2) Axes parameters (AXEPAR namelist):

a) Switches  (usually 0 = off, 1 = on):

iaxlab - plot axes labels (default: 1)

itrev -  draw the  time-distance  plot with  the  time increasing
   downward (default: 0)

b) Other axes parameters:

albht - height of axes labels (mm) (default: 2.5)
 
xmin, xmax - minimum and maximum distance (km) of velocity model
xmm - length of distance axis (mm)
xtmin, xtmax - minimum and maximum values (km) of tick marks
ntickx - number of intervals separated by tick marks along distance axis
ndecix - number of digits after decimal in distance axis annotation;
   -1 for integer annotation (defaults: 0, none, 250; values for
   xtmin, xtmax, ntickx and ndecix depend on values of xmin and xmax)

zmin, zmax - minimum and maximum depth (km) of velocity model
zmm - length of depth axis (mm)
ztmin, ztmax - minimum and maximum values (km) of tick marks
ntickz - number of intervals separated by tick marks along depth axis
ndeciz - number of digits after decimal in depth axis annotation;
    -1 for integer annotation (defaults: 0, 50, 75; values for
   ztmin, ztmax, ntickz and ndeciz depend on values of zmin and zmax)
 
xmint, xmaxt - minimum and maximum distance (km) of time-distance,
               amplitude-distance and RMS velocity plots
xmmt - length of distance axis (mm)
xtmint, xtmaxt - minimum and maximum values (km) of tick marks
ntckxt - number of intervals separated by tick marks along distance axis
ndecxt - number of digits after decimal in distance axis annotation;
   -1 for integer annotation (defaults: xmin, xmax, xmm; values for
   xtmint, xtmaxt, ntckxt and ndecxt depend on values of xmint and xmaxt)
 
tmin, tmax - minimum and maximum time (s) of time-distance plot
tmm - length of time axis (mm)
ttmin, ttmax - minimum and maximum values (s) of tick marks
ntickt - number of intervals separated by tick marks along time axis
ndecit - number of digits after decimal in time axis annotation;
   -1 for integer annotation (defaults: 0, 10, 75; values for
   ttmin, ttmax, ntickt and ndecit depend on values of tmin and tmax)
 
ampmin, ampmax - minimum and maximum amplitude (log base 10) of amplitude-
   distance plot
amm - length of amplitude axis (mm)
atmin, atmax - minimum and maximum values of tick marks
nticka - number of intervals separated by tick marks along amplitude axis
ndecia - number of digits after decimal in amplitude axis annotation; 
   -1 for integer annotation (defaults: 0, 8, 75, values for 
   atmin, atmax, nticka and ndecia depend on values of ampmin and ampmax) 
 
vmin, vmax - minimum and maximum velocity (km/s) of 1-D velocity
   profiles
vmm - length of velocity axis (mm)
vtmin, vtmax - minimum and maximum values (km/s) of tick marks (not
               implemented yet)
ntickv - number of intervals separated by tick marks along velocity axis
ndeciv - number of digits after decimal in velocity axis annotation; 
   -1 for integer annotation (defaults: 0, 8, 75, values for 
   vtmin, vtmax, ntickv and ndeciv depend on values of vmin and vmax) 
 
vrmin, vrmax - minimum and maximum velocity (km/s) of RMS 
   velocity-distance plot
vrmm - length of RMS velocity axis (mm)
vrtmin, vrtmax - minimum and maximum values (km/s) of tick marks (not
                 implemented yet)
ntckvr - number of intervals separated by tick marks along RMS velocity axis
ndecir - number of digits after decimal in RMS velocity  axis annotation; 
   -1 for integer annotation (defaults: 2, 7, 75, values for 
   vrtmin, vrtmax, ntckvr and ndecir depend on values of vrmin and vrmax)

orig - lower left origin of plot is (orig,orig) mm (default: 12.5)

sep - separation (mm) between the model and travel time digrams
   if isep=0 (default: 7.5)

3) Ray tracing parameters (TRAPAR namelist):

a) Switches  (usually 0 = off, 1 = on):

ishot - an array specifying the directions rays are to be
traced
 from the shot points listed in the arrays xshot and zshot; the
  following code is used:
   (1) 0 - no rays are traced
   (2) -1 - rays are traced to the left only
   (3) 1 - rays are traced to the right only
   (4) 2 - rays are traced to the left and right
   (default: 0)

iraysl - use the  array irayt to  select which ray  groups listed 
   in the array ray are to be traced for a particular shot point  
   (default: 0)

irayt - an array selecting  those ray groups listed  in the array 
   ray which are active for a particular shot point listed in the 
   arrays xshot and zshot.  If there are n  ray groups listed  in 
   the array ray,  then the  jth ray group  for the  ith shot is  
   referred to in  the ((n-1)*i+j)th element  of irayt (default:  
   1)

ifast - use a Runge-Kutta routine without error control to solve
   the ray tracing equations and a look-up/interpolation routine
   to evaluate certain trigonometric functions; with ifast=0, a
   Runge-Kutta method with error control and the intrinsic
   trigonometric functions are used. Using ifast=1 is about 30-40%
   faster and provides essentially the same accuracy as ifast=0
   (default: 1)  
   Note: ifast=1 is intended to be used for all routine modelling
         and ifast=0 only to test a final model or if ifast=1 
         obviously fails

iturn - an array  corresponding to the  ray groups listed  in the 
   array ray to trace  rays only to  their turning or  reflection 
   point in the search mode; if a refracted or head wave ray group 
   is to contain a reflection(s) before it turns (specified using 
   nrbnd and rbnd), then iturn=0 must be used (default: 1)

isrch - search again  for the take-off  angles of a  ray group if 
   the same ray code is listed more than once for the same shot in 
   the array ray; isrch=1 must be used if the take-off angles of 
   two or more ray  groups with the same  ray code are  different 
   because of multiple reflections and/or conversions specified by 
   nrbnd, rbnd, ncbnd,and cbnd (default: 0)

istop - stop  tracing any  ray which reflects  off a  boundary not
   specified in the arrays ray or rbnd; if istop=2, rays are also
   stopped if they enter a layer deeper than that specified by the
   ray code in the array ray, i.e., the layer number is greater
   than L and the ray code is L.n, where n=0, 1, 2, or 3
   (default: 1)

ibsmth - apply  a simulation  of smooth  layer boundaries  (1) or  
   apply the  simulation and  plot  the smoothed  boundaries (2);
    ibsmth=2 cannot be used if isep>1 (default: 0)

imodf - use the velocity  model in the file v.in  instead of the 
   model in part (3) of the file r.in (default: 0)

b) Other ray tracing parameters:

xshot - an  array containing  the x-coordinates  (km) of the  shot 
   points (default: xmin, for the first element only)

zshot - an  array containing  the z-coordinates  (km) of the  shot 
   points (default: a very small distance below the model surface)

ray -  an  array  containing the  ray  groups  to be  traced;  the  
   following code is used:
   (1) L.1 - rays which refract (turn) in the Lth layer
   (2) L.2 - rays which reflect off the bottom of the Lth layer
   (3) L.3 -  rays which  travel as  head waves along  the bottom