C GUIDE TO THE PROGRAM PACKAGE SEIS88 C *********************************** C C C C 1. DESCRIPTION OF THE PACKAGE C ***************************** C C THE PROGRAM PACKAGE IS DESIGNED FOR THE NUMERICAL MODELLING OF C SEISMIC WAVE FIELDS IN 2-D LATERALLY VARYING LAYERED STRUCTURES C BY THE RAY METHOD. C THE PACKAGE CONSISTS OF SIX PROGRAMS: THE BASIC PROGRAM, ALSO C CALLED SEIS88 AND FOUR AUXILIARY PROGRAMS, SMOOTH, RAYPLOT, SYNTPL C SEISPLOT AND POLARPLOT. C C THE BASIC PART OF ALL COMPUTATIONS IS PERFORMED IN THE PROGRAM C SEIS88. IT CAN BE SUBDIVIDED INTO THREE PARTS: SEIS, MODEL AND C RAY. C IN THE PART 'SEIS', THE DATA CONCERNING THE COMPUTED WAVE FIELD C ARE READ IN. TWO-POINT RAY TRACING AND GENERATION OF ELEMENTARY C WAVES ARE CONTROLLED FROM THIS PART OF THE PROGRAM. C IN THE PART 'MODEL', THE DATA CONCERNING THE MODEL OF MEDIUM C ARE READ IN AND CORRESPONDING APPROXIMATION OF INTERFACES AND VE- C LOCITY DISTRIBUTION IS PERFORMED. THIS PART OF THE PROGRAM SEIS88 C IS INTERCHANGEABLE. IN THE PRESENT VERSION, THE FOLLOWING TWO C 'MODEL' PARTS OF SEIS88 ARE AVAILABLE (THEY DIFFER BY VELOCITY C APPROXIMATION WITHIN INDIVIDUAL LAYERS): C I) BICUBIC SPLINE INTERPOLATION C II) LINEAR INTERPOLATION BETWEEN ISOVELOCITY INTERFACES C IN ALL SPECIFIED CASES, INTERFACES ARE APPROXIMATED BY INTER- C POLATING CUBIC SPLINES. C IN THE PART 'RAY', STANDARD RAY TRACING, DYNAMIC RAY TRACING C AND EVALUATION OF AMPLITUDES ARE PERFORMED. THIS PART AND ANY OF C 'MODEL' PARTS CAN BE ALSO USED IN OTHER PROGRAMS BASED ON RAY C TRACING. C C THE PROGRAM SEIS88 MAY USE AS INPUT DATA THE DATA GENERATED IN C THE PROGRAM SMOOTH AND STORED IN THE FORMATTED FORM IN THE FILE LU, C (SEE BELOW). THE INPUT DATA CAN BE ALSO READ IN DIRECTLY FROM TER- C MINAL OR CARDS. THE PROGRAM SEIS88 GENERATES TWO FILES: FILE LU1 C (UNFORMATTED) WITH DATA FOR PLOTTING RAY DIAGRAMS, TRAVEL TIMES C AND RAY AMPLITUDES AND FILE LU2 (FORMATTED OR UNFORMATTED) WITH C DATA FOR COMPUTATION OF RAY SYNTHETIC SEISMOGRAMS. C PROGRAM SEIS88 CONSISTS OF THE FOLLOWING PARTS: C SEIS.FOR, RAY.FOR, MODIS.FOR (OR MODSPL.FOR) AND SERV.FOR. C THESE PARTS MUST BE LINKED TOGETHER BEFORE RUNNING THE C PROGRAM. NONE OF THESE PARTS CAN WORK INDEPENDENTLY. C C THE PROGRAM SMOOTH MAY BE USED FOR THE PREPARATION OF A C SMOOTH VELOCITY INPUT DATA FOR THE PROGRAM SEIS88 IF BICUBIC C SPLINE INTERPOLATION IS USED IN SEIS88 TO APPROXIMATE THE C VELOCITY DISTRIBUTION, SEE I) ABOVE. THE PROGRAM PRODUCES C ONE FILE, FILE LU WITH THE VELOCITY INPUT DATA IN THE FORMAT- C TED FORM AS REQUIRED BY SEIS88 (SEE INPUT DATA CARD SUB 4, C BELOW). C PROGRAM SMOOTH CONSISTS OF THE FOLLOWING PARTS: C SMOOTH.FOR, SERV.FOR. C C THE PROGRAM RAYPLOT PERFORMS PLOTTING OF RESULTS STORED C IN THE FILE LU1, GENERATED BY SEIS88. IT PRODUCES PLOTS C OF RAY DIAGRAMS, TIME-DISTANCE AND AMPLITUDE-DISTANCE CURVES. C CALCOMP PLOTTING ROUTINES ARE REQUIRED. C PROGRAM RAYPLOT CONSISTS OF THE FOLLOWING PARTS: C RAYPL.FOR, SOURCE.FOR, BORDER.FOR, SERV.FOR. C C THE PROGRAM SYNTPL COMPUTES RAY SYNTHETIC SEISMOGRAMS FROM C THE DATA STORED IN FILE LU2, GENERATED IN SEIS88. THE PRO- C GRAM GENERATES TWO FILES: FILE LU3, WHICH CONTAINS SAMPLED C SYNTHETIC SEISMOGRAMS (IN A FORMATTED AND NORMALIZED FORM), C WHICH MIGHT BE FURTHER PROCESSED AND PLOTTED IN PROGRAM SEIS- C PLOT; UNFORMATTED FILE LU4, WHICH CAN BE USED FOR PLOTTING C OF SEISMOGRAMS USING STANDARD SEISMIC ROUTINES SUCH AS ROUTINE C TRPLOT IN UNIX OPERATING SYSTEM IN SEP STANFORD. C PROGRAM SYNTPL CONSISTS OF THE FOLLOWING PARTS: C SYNTPL.FOR, SOURCE.FOR, SERV.FOR. C C THE PROGRAM SEISPLOT CAN BE USED TO PLOT THE SYNTHETIC C SEISMOGRAMS STORED IN FILE LU3, GENERATED BY THE PROGRAM C SYNTPL, IN A WAY COMMON IN SEISMOLOGY (EACH SEISMOGRAM IS C PLOTTED AS A SINGLE CURVE). CALCOMP PLOTTING ROUTINES ARE C REQUIRED. C PROGRAM SEISPLOT CONSISTS OF THE FOLLOWING PARTS: C SEISPL.FOR, BORDER.FOR, SERV.FOR. C C THE PROGRAM POLARPLOT CAN BE USED TO PLOT THE SYNTHETIC C PARTICLE MOTION DIAGRAMS FROM DATA STORED IN FILES LU3, GENE- C RATED BY THE PROGRAM SYNTPL. CALCOMP PLOTTING ROUTINES ARE C REQUIRED. C PROGRAM POLARPLOT CONSISTS OF THE FOLLOWING PARTS: C POLAR.FOR, BORDER.FOR, SERV.FOR. C C NOTE: ALL THE PROGRAMS ARE WRITTEN IN FORTRAN IV. IN THE BEGIN- C NING OF EACH PROGRAM ROUTINE SERV IS CALLED WHICH ENABLES INTER- C ACTIVE OPENING OF REQUIRED INPUT AND OUTPUT FILE AND ALSO CON- C TROLS INITIALIZATION OF PLOTTING. C C C 2. CONTENT OF THE DISC C ********************** C C THIS DISC CONTAINS 4 KINDS OF FILES: SOURCE PROGRAMS, TEXT C FILES DESCRIBING THE SOURCE PROGRAMS, INPUT DATA FOR TEST C EXAMPLE AND OUTPUT FILES FOR THE TEST EXAMPLES (UNFORMATTED C FILES ARE NOT INCLUDED). THE FILES ARE ORDERED IN THE FOLLOWING C WAY: C C I) FILES WITH SOURCE PROGRAMS (OR PARTS OF SOURCE PROGRAMS): C SEIS.FOR, MODIS.FOR, MODSPL.FOR, RAY.FOR, SMOOTH.FOR, C RAYPL.FOR, SYNTPL.FOR, SEISPL.FOR, POLAR.FOR, BORDER.FOR, C SOURCE.FOR SERV.FOR. C C II) TEXT FILES: C SEIS.TXT, SMOOTH.TXT, RAYPL.TXT, SYNTPL.TXT, SEISPL.TXT, C POLAR.TXT. C C III) INPUT FILES FOR TEST EXAMPLE: C SEIS.UNC, RAYPL.UNC, SYNT.UNC, SEISPL.UNC, POLAR.UNC, C C IV) OUTPUT FILES FOR TEST EXAMPLES: C THE FILES ARE ORDERED AS FOLLOWS: C LU1.UNC, LU2.UNC, LU3.UNC. C C C GOOD LUCK! C C ************************************************************** C