Description of the program  B P L O T

    Program BPLOT  is designed for  plotting synthetic seismogram
sections from any  file which has the same  structure as the file
LU8 generated for this purpose by programs SYNFAN.

Short description of the program BPLOT

    Program  BPLOT is  a modified  version of  the equally  named
program from the program package BEAM87 by V.Cerveny.
    Synthetic seismogram section is plotted  in a frame, in which
horizontal axis corresponds either to the range (if receivers are
distributed along the surface or an internal interface) or to the
depth (if  receivers are distributed  along a vertical  profile).
The vertical axis  corresponds to the travel time  or the reduced
travel  time  axis.  It  is  possible  to  use  various  forms of
amplitude scaling  along the profile, i.e.  along the  horizontal
axis. The scaling of amplitudes along time axis is not  used, but 
may  be simply  introduced. A selection of receiver positions, at 
which  the  synthetic  seismograms  should  be  computed  can  be 
performed.  Positive  polarity  of  impulses  is  to  the left of 
vertical time axis.

                                                     
Description of input and output data

    Main input data are read from the standard input by list-directed
input (free format) and consist of a single line containing following
data:
    'LIN' 'LOU' 'LU8' 'LU4'/
Here:
    'LIN' is the name of the input data file LIN.
    'LOU' is the name of the output log file LOU.
    'LU8' is the name of the intput data file LU8,
          generated by program SYNFAN.
    'LU4' is the name of the output data file LU4.
    / is a slash recomended in batch and script files to enable future
        extensions.
Defaults:
    'LIN'='bplot.dat'
    'LOU'='bplot.out'
    'LU8'='lu8.dat'
    'LU4'='lu4.dat'
EXAMPLE OF THE MAIN INPUT DATA:
    'bplot.ani' /
    'bplot.qi' /

    Input  data consist  partially of  the data  obtained in  the
program SYNFAN (or any other program which generates  the data in
the required form) and are stored  in the file LU8, and partially
of the  additional input  data  read by list-directed input (free 
format), data specified by the user in the file  LIN. Output data  
describing the computations  are stored in  the file  LOU. Output 
data for plotting particle motion diagrams are stored in the file 
LU4. The program  generates a  postscript file  with the desired 
plot.

                                                      
The data stored in LU8

    The data are stored in a formatted form in the file LU8. For 
details  see  the  description of the content of the file LU8 in
program SYNFAN.

1) MPRINT                                FORMAT(A)
2) IPRINT                                FORMAT(A)
3) TITLE                                 FORMAT(A)
4) XSOUR,YSOUR,ZSOUR,TSOUR,RSTEP,DT,DF   FORMAT(5F10.5,2E15.7)
5) NDST,NT,MCOMP,ILOC                    FORMAT(16I5)
6) DIST,TO,AREDUC,NT                     FORMAT(2F10.3,1E12.5,I5)
7) IS(I),I=1,NT                          FORMAT(20I4)


                                                      
The additional input data in the file LIN

     These data are specified by  the user. They control plotting
of synthetic seismograms. The places, where the data from LU8 are
read in are denoted by **LU8/1, **LU8/2, etc.

1) TEXT

   TEXT...    arbitrary  alphanumeric  text describing  the  data
              set. Default 'BPLOT'

2) Switch, which controls the output into the file LOU and optional
   name of the plot file.

   IPR,PSTEXT

    IPR...    controls the printout into the file LOU and the type
              of the name of the postscript file with a plot.
              IPR.GE.0... automatic generation of postscript file
              names: plot00.ps, plot01.ps, etc.
              IPR.LT.0... postscript file name specified by the user,
              see PSTEXT.
              IPR=0... standard print (default value).
              IABS(IPR)=1... also a table of coordinates of receivers,
              with maximum amplitudes.
              IABS(IPR)=2...  also  synthetic  seismograms,  in reduced
              form.
              For more details see the section on output tables.
    PSTEXT... the name of the postscript file specified by the user;
              it is used only if IPR.LT.0.

**LU8/1
**LU8/2
**LU8/3
**LU8/4
**LU8/5

3) Various switches.

   MRED,MEPIC,NTICX,NTICY,NDX,NDY

    MRED...   controls the reduction of travel times.
              MRED=0... no reduction (default value).
              MRED=1... reduction applied, see input data No.6.
    MEPIC...  controls the selection of receiver positions.
              MEPIC=0... no selection (default value).
              MEPIC=1... selection applied, see input data No.4
    NTICX...  the  number  of  marked  intervals  along the range
              (depth)  axis  between  two  neighbouring tics with
              corresponding coordinate values (default value 1).
    NTICY...  the same as NTICY, but for the travel-time axis 
              (default value 1).
    NDX,NDY...control  the  precision  of  numbers describing the
              coordinate  axes in  the plots.  NDX corresponds to
              range  axis.  NDY  corresponds  to  the travel-time
              axis.
              ND.GT.0: The  number of digits to  the right of the
              decimal point.
              ND=0:  Only integer  portions of  the numbers  with
              decimal points (default value).
              ND.LT.0:  Integers.

4) Selection of receivers positions.
   Included only when MEPIC.NE.0.

   NEPIC,(IEP(I),I=1,NEPIC)

    NEPIC...  number  of  selected  receiver  positions  at which
              synthetic seismograms are to be constructed.
    IEP(1),IEP(2)...  sequential  numbers  of  selected  receiver
              positions.

    Note: Similar  selection can be  done already in  the program
    SYNFAN.  The selection here  is made  among  receivers, which
    remained after the selection in program SYNFAN.

5) Description of axes of the plot.

   XMIN,XMAX,XLEN,DTICX,YMIN,YMAX,YLEN,DTICY

    XMIN,XMAX... minimum and maximum  values along the range axis
              (in the users length units).
    XLEN...   length of the range/depth axis (in cm).
    DTICX...  the distance between two  neighbouring tics on  the
              range/depth   axis   which   are   denoted  by  the
              corresponding coordinate values (in the user length
              units).
              DTICX.GT.0.0:  Tic  marks  starting  from  XMIN and
              appearing  at  the  subsequent  points  XMIN+DTICX,
              XMIN+2.0*DTICX,...
              DTICX.LT.0.0:  Tic marks  start and  continue to be
              plotted   from  the   first  integer   multiple  of
              ABS(DTICX) greater than XMIN.
    YMIN,YMAX,YLEN,DTICY...  the  same  for  the  time  axis. For
              MRED=1  (reduced  time  along  the  time axis), the
              quantities correspond to the reduced travel times.

6) Parameters of the time reduction and of the amplitude scaling
   along the range/depth axis.

   VRED,AMP,B1,EPICS,EPS,SC

    VRED...   reduction velocity (in km/s) default VRED=6.0.
    AMP,B1,EPICS,EPS...  quantities   controlling  the  amplitude
              scaling. Default values 0., 1., 10. and 0., respectively.
              The  scaling  factor  may  generally  depend on the
              receiver position.  Denote the x-coordinate  of the
              receiver  by XX,  and the  average distance between
              neighbouring receiver  positions by DDX  (i.e., the
              average distance between individual traces, in plot
              units). Then the following  five systems of scaling
              can be used:

              1) FACTOR=B1*DDX/SMAXI,     when AMP=0.
              Normalization of maximum amplitudes in all traces:
              For  B1=1,  the  maximum  amplitudes  in all traces
              equal  to the  average distance  between individual
              traces.

              2) FACTOR=B1*DDX/SMAXIM,    when AMP=-1 and EPS=0.
              Real amplitudes: For B1=1, the maximum amplitude in
              the  whole  plot  equals  to  the  distance between
              individual traces.

              3) FACTOR=B1,               when AMP=1 and EPS=0.
              Real   amplitudes:  Manual   determination  of  the
              scaling factor by specifying B1.

              4) FACTOR=B1*DDX*((ABS(XX-XSOUR)/EPICS)**EPS)/SMAXIM
              when AMP=-1. and EPS.GT.0.
              Power scaling of amplitudes: When EPICS corresponds
              to the  epicentral distance of  the trace with  the
              maximum  amplitude  and   B1=1,  then  the  maximum
              amplitude  of  the  trace  at  EPICS  equals to the
              average distance between individual traces.

              5) FACTOR=B1*(ABS(XX-XSOUR)/EPICS)**EPS,
              when AMP=1. and EPS.GT.0.
              Power scaling of amplitudes:  Fully manual, all the
              parameters must be specified.

              Note: The scaling systems 4 and 5 give for XX=XSOUR
              traces with zero amplitudes.

    SC...     scaling factor, controls scaling of tics and
              alphanumeric texts. For SC=1.0, the tics are 0.15 
              cm long and coordinates and text describing the 
              plot are 0.4 and 0.45 cm high, respectively.

**LU8/6
**LU8/7

Example of data LIN for anisotropic model
Example of data LIN for isotropic model for QI computations

Termination of computations.

The data LU8/6 and LU8/7 are  read in NDST times (see **LU8/5 for
NDST). Afterwards, the computations terminate.

                                                      
Output to the file LOU

    All  the  input  data  are  reproduced  in  the file LOU. The
storage  of other  data is  controlled by  the parameter IPR, see
input data line No.2.
    For IPR=0, also the data LU8/1-LU8/5 are stored.
    For IPR=1, the same data as for IPR=0 are stored. In addition, 
for each receiver position the following data are stored:
    XX,SMAXI,FACTOR,SFMAX,
where XX is the coordinate of  the receiver, SMAXI is the maximum
amplitude  in   the  corresponding  seismogram,   FACTOR  is  the
amplitude  scaling  factor  corresponding  to  this receiver, and
SFMAX is the maximum amplitude  in the synthetic seismogram after
scaling.
    For IPR=2,  also the synthetic  seismograms are stored.  Each
synthetic seismogram has the heading:
XX,TSTART,SMAXI,FACTOR,SFMAX.
The quantities  XX,SMAXI,FACTOR,SFMAX have the  meaning as above.
TSTART is  the minimum time,  corresponding to the  first sample.
Synthetic  seismograms are  printed  in  a reduced  form, maximum
equals to unity.

                                                      
Output to the file LU4

The file LU4 contains the computed ray synthetic seismograms
in the form required for plotting particle motion diagrams in
program POLAR. The zero samples in the initial and final parts of 
each seismogram are not stored. The data are stored in the formatted 
form to allow for an inspection of computations. To perform the 
plotting of particle motion diagrams from the file LU4, the program 
POLAR, included in this package, can be used. The data on LU4 are 
stored in the following order:

1) MTEXT                                       FORMAT(A)
  Arbitrary alphanumeric text describing the computations. this
  text will appear under the plots. The text was specified in
  the program ANRAY.
2) MMD,MRED,MCOMP,ITPR,VRED,RSTEP,XSOUR,YSOUR,DT  FORMAT(4I5,5F10.5)
     MMD...    the number of selected receiver positions.
     MRED...   MRED=0: non reduced travel times,
               MRED=0: reduced travel times.
     MCOMP...  MCOMP=0: vertical component,
               MCOMP=1: component along the x-axis,
               MCOMP=2: component along the y-axis.
     ITPR...   ITPR=0: receivers along the earth's surface,
               ITPR=1: receivers along a vertical profile,
               ITPR.GT.100: receivers along the (itpr-100)th 
                            INTERFACE.
     VRED...   reduction velocity.
     RSTEP...  average difference between coordinates of
               neighbouring receiver positions.
     XSOUR, YSOUR...  x- and y-coordinates of the source.
     DT...     time step in synthetic seismogram.
3) XMX,SMAXIM                             FORMAT(22X,F10.5,9X,F15.9)
     XMX...    coordinate of the receiver at which maximum
               amplitude smaxim was recorded.
     SMAXIM... the maximum amplitude from all stored synthetic
               seismograms.
The following lines 4,5 are successively repeated for all receiver 
positions. They contain information about individual seismograms.
Only parts of seismograms are stored, starting with the first 
non-zero sample and ending with the last non-zero sample:
4) XX,SMAX,TM,NPS                    FORMAT(F10.5,E15.8,F10.5,I5)
     XX...     coordinate of the receiver on the profile
               (radius from the vertical line through the source
               for receivers along a surface or an interface
               profile, 'model' z coordinate for receivers along
               a vertical profile)
     SMAX...   maximum amplitude in the synthetic seismogram
     TM...     TM is the time corresponding to the first non-zero 
               sample in the seismogram. if MRED.NE.0, TM
               is a reduced time.
     NPS...    number of points in synthetic seismograms.
5) Synthetic seismograms in normalized form:
   JS(I), I=1,NPS                              FORMAT(20I4)
     JS(I)..   I-th sample, JS(I)=IFIX(999.1*SEIS(I)/SMAX),
               where SEIS(I) is computed value of i-th non-zero
               sample.

Note: The maximum length of the synthetic seismograms is 3001
points.