********************************************
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
X-AXIS CORRESPONDS TO THE X-AXIS IN THE MODEL FOR THE RECEIVERS
SITUATED ALONG THE EARTH'S SURFACE OR ALONG AN INTERFACE, AND TO
THE Z-AXIS FOR THE RECEIVERS SITUATED ALONG A VERTICAL PROFILE.
THE Y-AXIS OF THE FRAME CORRESPONDS TO THE TRAVEL TIME OR THE
REDUCED TRAVEL TIME. IT IS POSSIBLE TO USE VARIOUS FORMS OF AMP-
LITUDE SCALING ALONG THE PROFILE OF RECEIVERS, I.E. ALONG THE X-
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.
   CALCOMP PLOTTING ROUTINES PLOTS, PLOT, NUMBER AND SYMBOL
ARE REQUIRED.

***
POSITIVE POLARITY OF IMPULSES IS TO THE LEFT OF VERTICAL TIME AXIS.

THE COMPONENTS OF THE DISPLACEMENT VECTOR INTO THE X AND Y
COORDINATES OF THE MODEL COORDINATE SYSTEM ARE FORMALLY CALLED
"RADIAL" AND "TRANSVERSE" COMPONENTS, RESPECTIVELY. THEY ARE
TRULY RADIAL AND TRANSVERSE IF THE PROFILE OF RECEIVERS LIES
ALONG THE X-AXIS. BY USE OF THE ANGLE AROT, SEE PROGRAMS ANRAYPL
AND FRESAN, HORIZONTAL COMPONENTS CAN BE TRANSFORMED INTO
ROTATED COORDINATE SYSTEM.
***

                                                     
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.out'
EXAMPLE OF THE MAIN INPUT DATA:
    'seispl.sch' /

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 CONTROLLING PLOTTING,
SPECIFIED BY THE USER AND 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.

                                                      
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(17A4)
   ARBITRARY ALPHANUMERIC TEXT DESCRIBING THE COMPUTATIONS. THIS
   TEXT WILL APPEAR UNDER THE PLOTS. IT WAS READ IN IN THE PRO-
   GRAM SEIS88.
 2) MDIST,MRED,MCOMP,ITPR,VRED,RSTEP,XSOUR,DT  FORMAT(4I5,4F10.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: RADIAL COMPONENT.
                MCOMP=2: TRANSVERSE COMPONENT.
      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...  X-COORDINATE OF THE SOURCE, (XSOUR=0.).
      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.9,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(20I4)
 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 SHOULD BE INTRODUCED DIRECTLY (CARDS, TERMI-
 NAL). THEY CONTROL THE PLOTTING OF SYNTHETIC SEISMOGRAMS. THE
 DATA FROM LU4 ARE READ IN AFTER ALL THE ADDITIONAL INPUT DATA
 ARE READ IN.

1) ONE LINE
   IPRINT,XSHIFT,YSHIFTY                     FORMAT(I5,2F10.5)
        IPRINT...CONTROLS STORAGE OF SOME DATA IN THE FILE LOU.
                 SEE DETAILS IN THE SECTION ON OUTPUT TABLES.
        XSHIFT,YSHIFT... THE X- AND Y-SHIFTS WITH RESPECT TO
                 THE ORIGIN OF COORDINATES, IN CM.

2) ONE LINE, VARIOUS SWITCHES
  MCONT,MEPIC,NTICX,NTICY,NDX,NDY          FORMAT(16I5)
         MCONT... MCONT=0: THE COMPUTATION TERMINATES.
                  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
                  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.
         NTICY... THE SAME AS NTICY, BUT FOR THE TRAVEL-TIME AXIS.
         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.
                  ND.LT.0: INTEGERS.

3) ONE LINE, SELECTION OF RECEIVER POSITIONS.
  INCLUDED ONLY WHEN MEPIC.NE.0.
     NEPIC,(IEP(I),I=1,NEPIC)              FORMAT(16I5)
         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    FORMAT(8F10.5)
        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,S               FORMAT(8F10.5)
    WHEN EPICS NOT SPECIFIED, EPICS=10.
    WHEN B1 NOT SPECIFIED, B1=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 THE SCALES OF TICS AND ALPHANUMERICAL TEXTS.
       FOR SC=1.0, 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


 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.