PROGRAM   S E I S P L O T
*************************

  PROGRAM SEISPLOT IS DESIGNED FOR PLOTTING OF SYNTHETIC SEISMOGRAMS
STORED IN THE FILE LU3, GENERATED IN THE PROGRAM SYNTPL. IT CAN BE
ALSO USED TO PLOT OTHER SYNTHETIC SEISMOGRAMS WHICH ARE STORED IN
THE SAME FORM AS THOSE GENERATED IN PROGRAM SYNTPL (SEE A SHORT
INFORMATION IN THE FOLLOWING).

*****************************************************************

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
AMPLITUDE 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 ARE REQUIRED.
  THE PROGRAM GENERATES A POSTSCRIPT FILE WITH THE DESIRED PLOT.

******************************************************************

INPUT DATA
**********

INPUT DATA CONSIST PARTIALLY OF THE DATA GENERATED BY A PROGRAM
SYNTPL AND STORED IN A FORMATTED FORM IN THE FILE LU3, AND
PARTIALLY OF THE ADDITIONAL INPUT DATA CONTROLLING PLOTTING.
SOME OF THE ADDITIONAL INPUT DATA ARE READ IN FROM STANDARD INPUT
BY LIST-DIRECTED INPUT (FREE FORMAT). SOME OF THE DATA MAY BE ALSO
SPECIFIED IN THE FORM OF THE SEP PARAMETERS,
SEE ALSO THE SEP PARAMETERS CONTROLLING THE FORM OF THE OUTPUT
POSTSCRIPT FILE.


DATA STORED IN LU3
******************

  IN THE FILE LU3, THE COMPUTED RAY SYNTHETIC SEISMOGRAMS AT
INDIVIDUAL RECEIVER POSITIONS ARE STORED, TOGETHER WITH SOME OTHER
RELEVANT IMPORTANT INFORMATION.
  THE DATA IN LU3 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 PROGRAM
   SEIS.
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.
      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
                RECEIVERS 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.

ADDITIONAL INPUT DATA
*********************

  THE FOLLOWING DATA SHOULD BE INTRODUCED DIRECTLY (CARDS, TERMINAL).
THEY CONTROL THE PLOTTING OF SYNTHETIC SEISMOGRAMS. THE DATA FROM LU3
ARE READ IN AFTER ALL THE ADDITIONAL INPUT DATA ARE READ IN.

1) ONE CARD
   ISHIFT,IPRINT
      ISHIFT...THE SHIFT BETWEEN INDIVIDUAL PICTURES, IN CM.
      IPRINT...CONTROLS THE PRINTOUT OF SOME DATA ON THE LINE
               PRINTER. SEE DETAILS IN THE SECTION ON OUTPUT
               TABLES.

2) ONE CARD, VARIOUS SWITCHES
   MCONT,MEPIC,NTICX,NTICY,NDX,NDY
      MCONT..MCONT=0: THE COMPUTATION TERMINATES.
             MCONT=1: A NEW PLOT OF SYNTHETIC SEISMOGRAMS WILL BE
                MADE FROM THE SAME FILE AS THE PRECEDING PLOT.
             MCONT=-1: A NEW PLOT OF SYNTHETIC SEISMOGRAMS WILL BE
                MADE FROM A NEW FILE. GO TO THE INPUT DATA CARD NO.1.
      MEPIC...  CONTROLS THE SELECTION OF RECEIVER POSITIONS
             MEPIC=0: NO SELECTION
             MEPIC=1: SELECTION APPLIED, SEE CARD NO.3
      NTICX...  THE NUMBER OF MARKED INTERVALS ALONG THE X-AXIS
                BETWEEN TWO ADJOINING TICS DENOTED BY COORDINATE
                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 CARD, 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 CARD, 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 CARD. PARAMETERS OF THE AMPLITUDE SCALING ALONG THE X-AXIS
   AMP,B1,EPICS,EPS,SC
   WHEN EPICS NOT SPECIFIED, EPICS=10.
   WHEN B1 NOT SPECIFIED, B1=1.
   THE PARAMETERS AMP,B1,EPICS,EPS CONTROL THE SCALING OF SEISMOGRAMS.
   THE SCALING FACTOR OF SYNTHETIC SEISMOGRAMS MAY GENERALLY DEPEND
   ON THE RECEIVER POSITION. DENOTE THE COORDINATE OF THE RECEIVER
   BY XX, AND THE AVERAGE DISTANCE BETWEEN NEIGHBOURING 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 CARD 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 CARD 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,
             RESPECTIVELY.

****************************************************************

TERMINATION OF COMPUTATIONS
***************************

  IF MCONT=1 IN CARD NO.2, THE CARDS 2-5 CAN BE REPEATED AN ARBITRARY
NUMBER OF TIMES TO PLOT SYNTHETIC SEISMOGRAMS WITH DIFFERENT AMPLITUDE
SCALING, DIFFERENT FRAMES, ETC.
  IF MCONT=-1 IN CARD NO.2, THE FILE LU3 IS CLOSED, AND A NUMBER
OF A NEW FILE IS READ IN FROM INPUT DATA CARD NO.1. THE FILE IS
OPENED AND THE COMPUTATIONS PROCEED AS DESCRIBED ABOVE. NOTE
THAT THIS POSSIBILITY CAN BE USED ONLY IN OPERATING SYSTEMS
ALLOWING SPECIFICATION OF THE NUMBER OF A FILE OUTSIDE THE SOURCE
PROGRAM (E.G., IBM OPERATING SYSTEM).
  THE COMPUTATION TERMINATES WHEN MCONT=0 IN THE INPUT DATA CARD NO.2.

******************************************************************

OUTPUT ON A LINE PRINTER
************************

  ALL THE ADDITIONAL INPUT DATA ARE REPRODUCED ON THE LINE PRINTER.
THE PRINTOUT OF OTHER DATA IS CONTROLLED BY THE PARAMETER IPRINT,
SEE INPUT DATA CARD NO.1.
  FOR IPRINT=0: ALSO THE DATA LU3/1,LU3/2 AND LU3/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.

*****************************************************************