********************************************
DESCRIPTION OF THE PROGRAM   S E I S P L O T
********************************************

  PROGRAM SEISPLOT IS DESIGNED FOR PLOTTING OF SYNTHETIC SEISMO-
GRAMS STORED IN THE FILE LU4, GENERATED IN THE PROGRAM SYNTPL
FROM PROGRAM PACKAGE SEIS83 OR IN THE PROGRAM SYNTAN FROM
PROGRAM PACKAGE ANRAY. IT CAN BE ALSO USED TO PLOT OTHER
SYNTHETIC SEISMOGRAMS, WHICH ARE STORED IN THE SAME FORM AS THOSE
GENERATED IN PROGRAMS SYNTPL AND SYNTAN


A SHORT DESCRIPTION OF THE PROGRAM
----------------------------------

  SYNTHETIC SEISMOGRAMS ARE PLOTTED IN A FRAME, IN WHICH THE
HORIZONTAL AXIS CORRESPONDS TO THE X-AXIS OF THE MODEL IN CASE OF
RECEIVERS SITUATED ALONG THE EARTH'S SURFACE OR ALONG AN INTERFACE,
AND TO THE Z-AXIS IN CASE OF RECEIVERS SITUATED ALONG A VERTICAL
PROFILE. THE VERTICAL AXIS OF THE FRAME CORRESPONDS TO THE TRAVEL
TIME OR THE REDUCED TRAVEL TIME. IT IS POSSIBLE TO USE VARIOUS FORMS
OF AMPLITUDE SCALING ALONG THE HORIZONTAL AXIS OF THE FRAME. THE
SCALING OF SYNTHETIC SEISMOGRAMS ALONG TIME AXIS IS NOT USED, BUT MAY
BE SIMPLY INTRODUCED. A SELECTION OF RECEIVERS AT WHICH THE SYNTHETIC
SEISMOGRAMS SHOULD BE COMPUTED CAN BE PERFORMED. POSITIVE POLARITY OF
SEISMOGRAMS IS TO THE LEFT OF VERTICAL TIME AXIS.
   PLOTTING ROUTINES PLOTS, PLOTN, PLOT, NUMBER AND SYMBOL ARE
REQUIRED.

                                                     
DESCRIPTION OF INPUT AND OUTPUT FILES
-------------------------------------
    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' 'LU4'/
HERE:
    'LIN' IS THE NAME OF THE INPUT DATA FILE LIN.
    'LOU' IS THE NAME OF THE OUTPUT LOG FILE LOU.
    'LU4' IS THE NAME OF THE INTPUT DATA FILE LU4,
          GENERATED BY PROGRAM SYNTAN.
    / IS A SLASH RECOMENDED IN BATCH AND SCRIPT FILES TO ENABLE FUTURE
        EXTENSIONS.
DEFAULTS:
    'LIN'='seispl.dat'
    'LOU'='seispl.out'
    'LU4'='lu4.dat'
EXAMPLE OF THE MAIN INPUT DATA:
    'seispl.chs' /

INPUT DATA CONSIST PARTIALLY OF THE DATA GENERATED BY THE PROGRAM
SYNTAN AND STORED IN A FORMATTED FORM IN THE FILE LU4, AND
PARTIALLY OF THE ADDITIONAL INPUT DATA READ BY LIST-DIRECTED
INPUT (FREE FORMAT). THE ADDITIONAL DATA CONTROL PLOTTING, ARE
PREPARED BY THE USER AND ARE STORED IN THE FILE LIN. OUTPUT DATA
DESCRIBING THE COMPUTATIONS ARE STORED IN THE FILE LOU. THE
PROGRAM GENERATES A POSTSCRIPT FILE WITH THE DESIRED PLOT, SEE ALSO
SEP PARAMETERS CONTROLLING THE FORM OF THE OUTPUT POSTSCRIPT FILE.

                                                      
THE DATA STORED IN LU4
----------------------

IN THE FILE LU4, THE COMPUTED RAY SYNTHETIC SEISMOGRAMS AT INDI-
VIDUAL RECEIVER POSITIONS ARE STORED, TOGETHER WITH SOME OTHER
RELEVANT IMPORTANT INFORMATION.
THE DATA IN LU4 ARE STORED IN THE FOLLOWING ORDER:
1) MTEXT                                  FORMAT(A)
   ARBITRARY ALPHANUMERIC TEXT DESCRIBING THE COMPUTATIONS. THIS
   TEXT WILL APPEAR UNDER THE PLOTS. IT WAS READ IN IN THE PRO-
   GRAM SEIS88 OR ANRAY.
2) MDIST,MRED,MCOMP,ITPR,VRED,RSTEP,XSOUR,YSOUR,DT  FORMAT(4I5,5F10.5)
      MDIST...  THE NUMBER OF RECEIVER POSITIONS.
      MRED...   MRED=0:NON REDUCED TRAVEL TIMES,
                MRED=1: 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 AN INTERFACE.
      VRED...   REDUCTION VELOCITY
      RSTEP...  AVERAGE DIFFERENCE BETWEEN X-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,E15.9)
      XMX...    COORDINATE OF OF THE RECEIVER AT WHICH THE
                MAXIMUM AMPLITUDE 'SMAXIM' WAS RECORDED.
      SMAXIM... THE MAXIMUM AMPLITUDE OF ALL TRACES.
THE FOLLOWING CARDS 4,5 ARE SUCCESSIVELY REPEATED FOR ALL
RECEIVER POSITIONS:
4) XX,SMAX,TMIN,NPTS            FORMAT(F10.5,E15.8,F10.5,I5)
      XX...     COORDINATE OF THE RECEIVER (X-COORDINATE WHEN
                RECEIVERS ARE SITUATED ALONG THE EARTH'S SURFACE
                OR ALONG AN INTERFACE AND Z-COORDINATE WHEN RE-
                CEIVERS ARE SITUATED ALONG A VERTICAL PROFILE.
      SMAX...   MAXIMUM AMPLITUDE IN THE SYNTHETIC SEISMOGRAM
                AT XX.
      TMIN...   THE TIME CORRESPONDING TO THE FIRST NON-ZERO
                POINT IN THE SYNTHETIC SEISMOGRAM.
      NPTS...   NUMBER OF POINTS IN SYNTHETIC SEISMOGRAM.
5) SYNTHETIC SEISMOGRAM                   FORMAT(20F4.0)
THE MAXIMUM LENGTH OF ONE SYNTHETIC SEISMOGRAM IS 3001 POINTS.
ONLY THE PARTS OF THE SEISMOGRAMS SITUATED BETWEEN THE FIRST
AND THE LAST NON-ZERO ARRIVAL ARE STORED. EACH SEISMOGRAM IS
STORED IN A NORMALIZED FORM, SEE THE FORMAT, WITH MAXIMUM
AMPLITUDE EQUAL 999.

                                                      
THE ADDITIONAL INPUT DATA IN THE FILE LIN
-----------------------------------------

THE FOLLOWING DATA CONTROL THE PLOTTING OF SYNTHETIC SEISMOGRAMS.
THE DATA FROM LU4 ARE READ AFTER ALL THE ADDITIONAL INPUT DATA
ARE READ IN.

1) ONE LINE
   IPRINT,XSHIFT,YSHIFTY,PSTEXT
        IPRINT...CONTROLS STORAGE OF SOME DATA IN THE FILE LOU
                 (SEE DETAILS IN THE SECTION ON OUTPUT TABLES)
                 AND THE TYPE OF THE NAME OF THE POSTSCRIPT FILE
                 WITH A PLOT.
                 IPRINT.GE.0... AUTOMATIC GENERATION OF POSTSCRIPT FILE
                 NAMES: PLOT00.PS, PLOT01.PS, ETC.
                 IPRINT.LT.0... POSTSCRIPT FILE NAME SPECIFIED BY THE USER,
                 SEE PSTEXT.
                 DEFAULT IPRINT=0.
        XSHIFT,YSHIFT... THE X- AND Y-SHIFTS WITH RESPECT TO
                 THE ORIGIN OF COORDINATES, IN CM. DEFAULT XSHIFT=3,
                 YSHIFT=6.

2) ONE LINE, VARIOUS SWITCHES
  MCONT,MEPIC,NTICX,NTICY,NDX,NDY
         MCONT... MCONT=0: THE COMPUTATION TERMINATES. DEFAULT
                  VALUE.
                  MCONT=1: A NEW PLOT OF SYNTHETIC SEISMOGRAMS
                  WILL BE MADE FROM THE SAME FILE AS THE PRECE-
                  DING PLOT.
         MEPIC... CONTROLS THE SELECTION OF RECEIVER POSITIONS
                  MEPIC=0: NO SELECTION. DEFAULT VALUE.
                  MEPIC=1: SELECTION APPLIED, SEE LINE NO.3
         NTICX... THE NUMBER OF MARKED INTERVALS ALONG THE X-AXIS
                  BETWEEN TWO ADJOINING TICS DENOTED BY COORDI-
                  NATE VALUES. DEFAULT NTICX=1.
         NTICY... THE SAME AS NTICY, BUT FOR THE TRAVEL-TIME AXIS.
                  DEFAULT NTICY=1.
         NDX,NDY... CONTROL THE PRECISION OF NUMBERS DESCRIBING
                  THE COORDINATE AXES IN THE PLOTS.
                  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.

3) ONE LINE, SELECTION OF RECEIVER 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),...,IEP(NEPIC)... SEQUENTIAL NUMBERS OF
                  SELECTED RECEIVER POSITIONS (FROM THE LEFT TO
                  THE RIGHT).
NOTE: SIMILAR SELECTION CAN BE DONE ALREADY IN THE PROGRAM
SYNTPL. THE SELECTION HERE IS MADE AMONG THOSE RECEIVERS WHICH
REMAINED AFTER THE SELECTION IN SYNTPL.

4) ONE LINE, DESCRIPTION OF AXES IN THE PLOT
  XMIN,XMAX,XLEN,DTICX,YMIN,YMAX,YLEN,DTICY
        XMIN,XMAX...THE MINIMUM AND MAXIMUM VALUES AT THE X-AXIS
                  (IN THE USERS LENGTH UNITS)
        XLEN...   LENGTH OF THE X-AXIS (IN CM)
        DTICX...  THE DISTANCE BETWEEN TWO NEIGHBOURING MARKS ON
                  THE X-AXIS WHICH ARE DENOTED BY CORRESPONDING
                  COORDINATE VALUES (IN THE USERS LENGTH UNITS).
                   DTICX.GT.0.: TIC MARKS STARTING FROM XMIN AND
                  APPEARING AT THE SUBSEQUENT POINTS XMIN+DTICX,
                  XMIN+2.*DTICX,...
                   DTICX.LT.0.: TIC MARKS START AND CONTINUED 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
                   REDUCED TIME AXIS (MRED=1), THE QUANTITIES
                   CORRESPOND TO THE REDUCED TRAVEL TIMES

5) ONE LINE. PARAMETERS OF THE AMPLITUDE SCALING ALONG THE X-AXIS
    AMP,B1,EPICS,EPS,SC
    DEFAULT AMP=0, B1=1, EPICS=10, EPS=0, SC=1.
    THE PARAMETERS AMP,B1,EPICS,EPS CONTROL THE SCALING OF SEIS-
    MOGRAMS. THE SCALING FACTOR OF SYNTHETIC SEISMOGRAMS MAY GENE-
    RALLY DEPEND ON THE RECEIVER POSITION. DENOTE THE COORDINATE
    OF THE RECEIVER BY XX, AND THE AVERAGE DISTANCE BETWEEN NEIGH-
    BOURING RECEIVERS 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. (THE WHOLE LINE 5 MAY BE BLANK IN THIS CASE).
       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.
       THE WHOLE LINE 5 WITH THE EXCEPTION OF AMP MAY BE BLANK IN
       THIS CASE.
       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 PARA-
       METERS MUST BE SPECIFIED.
NOTE THAT THE SCALING SYSTEMS 4 AND 5 GIVE FOR XX=XSOUR TRACES
WITH ZERO AMPLITUDES.
       SC... CONTROLS SCALING OF TICS AND ALPHANUMERICAL TEXTS.
       FOR SC=1, THE TICS ARE 0.15 CM LONG AND COORDINATES AND
       TEXT DESCRIBING THE PLOTS ARE 0.4 AND 0.45 CM HIGH, RES-
       PECTIVELY.

Example of data LIN for anisotropic model (isosurface approximation)


TERMINATION OF COMPUTATIONS
---------------------------

IF MCONT=1 IN LINE NO.2, THE LINES 2-5 CAN BE REPEATED AN ARBI-
TRARY NUMBER OF TIMES TO PLOT SYNTHETIC SEISMOGRAMS WITH DIFFE-
RENT AMPLITUDE SCALING, DIFFERENT FRAMES, ETC.
THE COMPUTATION TERMINATES WHEN MCONT=0 IN THE INPUT DATA LINE
NO.2.


                                                      
OUTPUT TO THE FILE LOU
----------------------

ALL THE ADDITIONAL INPUT DATA ARE STORED IN THE FILE LOU.
THE STORAGE OF OTHER DATA IS CONTROLLED BY THE PARAMETER IPRINT,
SEE INPUT DATA LINE NO.1.
FOR IPRINT=0: ALSO THE DATA LU4/1,LU4/2 AND LU4/3 ARE PRINTED.
FOR IPRINT=1: THE SAME AS FOR IPRINT=0. IN ADDITION, ONE LINE IS
PRINTED FOR EACH RECEIVER POSITION:
XX,SMAX,FACTOR,SFMAX
WHERE XX IS THE COORDINATE OF THE RECEIVER, SMAX 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.