rayinvr.doc

实现地震射线追踪和反演的frotran程序

Output file: m.out  contains the  two-way travel  times of  model 
   boundaries and RMS velocity profile corresponding to the plots 
   generated by VMODEL

Output file: p.out contains all plot commands for the run used as 
   input by the program RAYPLOT

                         Input parameters

1) Plotting parameters (PLTPAR namelist):

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

iplot - generate the plot  during the run (1), or  write all plot 
   commands to the file p.out (0), or do both (2) (default: 1)

imod - plot  the 2-D  model as  distance versus depth  or two-way  
   travel time (default: 0)

idash - plot model boundaries as dashed lines (default: 1)

ivel -  plot  the velocity  values  (km/s)  at their  appropriate  
   locations within model (default: 0)

ivp - plot 1-D  velocity profiles versus depth  or two-way travel  
   time; the velocity  profiles are  written to  the file vm.out  
   (default: 0)

izort -  plot  the 1-D  velocity  profiles and  2-D  model versus  
   two-way travel time; if izort=0, they are plotted versus depth 
   (default: 0)

ivrms - plot the RMS  velocity of the model  versus distance; the 
   RMS velocity profile is written to  the file m.out if idump=1  
   (default: 0)

idump - write  the two-way travel  times of model  boundaries and  
   the RMS  velocity  profile  as  plotted  to  the  file  m.out  
   (default: 0)

ifrefl - plot floating reflectors from the file f.in in the 2-D
   model (default: 0)

b) Other plotting parameters:

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

xtinc - the spacing  (km) of points  at which the  two-way travel 
   time of model boundaries is  calculated if imod=1 and izort=1  
   (default: (xmax-xmin)/100)

ntsmth - number of applications of  a 3-point averaging filter to 
   the model boundaries when plotted  versus two-way travel time  
   (default: 0)

xvp - an  array containing  the x-coordinates  (km) at  which 1-D  
   velocity   profiles   are   plotted    if   ivp=1   (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  
   (default: 1, number of model layers)

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

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

ibcol, ifcol - background and foreground colours (defaults: 0, 1)

2) Axes parameters (AXEPAR namelist):

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

iaxlab - plot axes labels (default: 1)

b) Other axes parameters:

albht - height of axes labels (mm) (default: 0.2)

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)

tmin, tmax - minimum and maximum two-way travel time (s) of 1-D
             velocity profiles or 2-D model boundary plots
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)

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
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 profile
vrmm - length of RMS velocity axis (mm)
vrtmin, vrtmax - minimum and maximum values (km/s) of tick marks
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) 

3) Model parameters (MODPAR namelist):

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

iorder - re-order all boundary  and velocity nodes from  left to 
   right; if iorder=0, a  message is written to  the screen and  
   the file vm.out to indicate the locations where boundaries or 
   velocities are not listed left to right (default: 0)

izint  -  an  array  containing  the  layer   boundaries  to  be  
   re-sampled to increase  or decrease  the number  of boundary  
   nodes

ivuint, ivlint -  an array  containing the  layers in which  the 
   upper and  lower  velocity  values  are  to  be  re-sampled,  
   respectively, to increase or decrease the number of velocity 
   points

iswit - when interpolating  the boundaries and  velocities using 
   izint, ivuint and ivlint, the corresponding switches equal to 
   0, 1 or -1 (see lines  c, f and i in  the description of the 
   file r.in for the program RAYINVR) will be set equal to 1 if 
   either of  the switches  across which  the  interpolation is  
   performed is equal to 1 and it will be set equal to -1 if the 
   two switches equal 0  and -1; if  iswit=0, the corresponding  
   switch will be set equal  to 0 if the  switches across which 
   the interpolation is performed are different (default: 1)

izsmt - an array containing the layer  boundaries to be smoothed 
   using a three-point averaging filter

ivusmt, ivlsmt -  an array  containing the  layers in which  the 
   upper  and  lower  velocity  values   are  to  be  smoothed,  
   respectively,  using   a   three-point   averaging   filter;   
   smoothing is  performed  laterally  on  the  upper and  lower  
   velocities separately.

b) Other model parameters:

nzint - an  array containing  the number of  points at  which to 
   uniformly re-sample  the layer  boundaries specified  in the  
   array izint;  if  nzint=1,  the  average  boundary depth  is  
   calculated; if nzint  equals the current  number of boundary  
   nodes, no re-sampling is performed  (default: none; however, 
   the default  is the  first element  of nzint  for  all layer  
   boundaries if only one value is specified)

nvuint, nvlint -  an array  containing the  number of points  at 
   which to  uniformly  re-sample  the  upper  and lower  layer  
   velocities  specified  in  the  arrays  ivuint  and  ivlint,  
   respectively; if nvuint=1 or nvlint=1,  the average upper or 
   lower layer  velocity  is calculated;  if  nvuint  or nvlint  
   equals the current number of upper or lower velocity points, 
   no re-sampling  is  performed (default:  none;  however, the  
   default is the  first element of  nvuint and nvlint  for all  
   layers if only one value is specified)

nzsmt -  an array  containing the  number of  applications of  a  
   three-point  averaging  filter   to  the   layer  boundaries   
   specified in  the  array  izsmt  (default:  1; however,  the  
   default  is  the  first  element  of  nzsmt  for  all  layer  
   boundaries if only one value is specified)

nvusmt, nvlsmt - an array containing the  number of applications 
   of a  three-point averaging  filter to  the upper  and lower  
   layer velocities specified in the  arrays ivusmt and ivlsmt, 
   respectively (default: 1; however, the  default is the first 
   element of nvusmt and nvlsmt for all layers if only one value 
   is specified)

velmin, velmax -  write a  message to  the screen  and the  file 
   vm.out indicating the location and  value of velocities less 
   than  or  greater  than  velmin   and  velmax,  respectively  
   (defaults: 0.1, 7)

dsmax -  write  a message  to  the screen  and  the file  vm.out  
   indicating the location  and value of  slope changes between  
   three successive boundary nodes greater than dsmax (degrees), 
   (default: 10)

dvvmax, dlvmax -  write a  message to  the screen  and the  file 
   vm.out indicating the location and value of vertical velocity 
   gradients greater  than  dvvmax (s-1)  and  lateral velocity  
   gradients greater than dlvmax (s-1), (defaults: 1.0, 0.1)

                        Additional Notes

Order of operations:  the following list  is the order  in which 
   the main operations are performed by VMODEL:
      1) check  that boundary  nodes and  velocities are  listed 
        left to right, or re-order them
      2) re-sample boundary nodes and velocities to  increase or 
        decrease the number of model parameters
      3) smooth layer boundaries and lateral velocity variations
      4)  check  for  small  and  large  velocities,  boundaries  
        crossing, low-velocity zones, negative vertical velocity 
        gradients, large changes  in boundary  slope, and  large 
        vertical and lateral velocity gradients
      5) plot  the model,  velocity  profiles and  RMS  velocity 
        profile
      6) write velocity model to  the file v.out if it  has been 
       changed

Low-velocity zones and pinchouts: low-velocity  zones that occur 
   across one or more layers that have been pinched out will not 
   be detected.

External plot routines: VMODEL is self-contained except  for the
   graphics calls inside the following CALCOMP or CALCOMP-like plot 
   routines: PLOTS, PLOT, NUMBER, SYMBOL, LINE, EMPTY, ERASE, PLOTND. 
   The file pltlib.f contains ALL these routines (and a few other 
   less important ones used by the graphics system where the 
   code was developed) and it will be necessary to replace the "call" 
   statements within each of these with their equivalents in the 
   local graphics package.

                           Array Sizes

The parameters specifying the size of all arrays within VMODEL are 
contained in the  file vmodel.par and  may be altered  to suit  a 
particular problem and memory requirements.

model layers - player

points defining a single layer boundary (must be a multiple of 10) 
   - ppcntr

points at which the upper or lower layer velocity is defined (must 
   be a multiple of 10) - ppvel

1-D velocity profiles - pnvp

points defining layer boundary versus two-way travel time or RMS 
   velocity profile - ppnpts

                    Warning and Error Messages

The following  is a  list and  explanation of  warning and  error 
messages generated by VMODEL:

***  xmax  not  specified  ***  -  the  value  of  xmax  must  be  
   specified in the AXEPAR namelist.

*** xmin not  specified ***  - the  value of xmin  is needed  and 
   must be specified in the  AXEPAR namelist since it  cannot be 
   determined from the file v.in.

*** a value  of nzint  > ppcntr  *** -  all values  in the  array 
   nzint must  be less  than or  equal  to the  value of  ppcntr 
   specified in the file vmodel.par.

*** a value  of nvuint  > ppvel  *** -  all values  in the  array 
   nvuint must  be less  than or  equal  to the  value of  ppvel 
   specified in the file vmodel.par.

*** a value  of nvlint  > ppvel  *** -  all values  in the  array 
   nvlint must  be le