******************************************** 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 LU3, 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. DESCRIPTION OF INPUT AND OUTPUT DATA ------------------------------------ INPUT DATA CONSIST PARTIALLY OF THE DATA GENERATED BY THE PROGRAM SYNTAN AND STORED IN A FORMATTED FORM IN THE FILE LU3, 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. SPECIFICATION OF THE FILES LIN, LOU AND LU3 CAN BE MADE THROUGH THE ROUTINE SERV, WHICH IS PART OF THIS PROGRAM PACKAGE. IF THE ROUTINE SERV IS NOT USED (PUT 'C' IN FRONT OF 'CALL SERV' IN THE BEGINNING OF THE MAIN PROGRAM), THE FILES LIN AND LOU ARE AUTOMATICALLY SPECIFIED LIN=5 LOU=6, THE FILE WITH PLOT HAS NUMBER 7, AND THE NUMBER LU3 IS READ FROM THE FILE LIN, SEE BELOW. THE DATA STORED IN LU3 ---------------------- IN THE FILE LU3, THE COMPUTED RAY SYNTHETIC SEISMOGRAMS AT INDI- VIDUAL 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 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 LU3 ARE READ IN AFTER ALL THE ADDITIONAL INPUT DATA ARE READ IN. 1) ONE CARD LU3,ISHIFTX,ISHIFTY,IPRINT FORMAT(16I5) LU3... THE NUMBER OF THE FILE IN WHICH THE SYNTHETIC SEISMOGRAMS ARE STORED. ISHIFTX,ISHIFTY... THE X- AND Y-SHIFTS WITH RESPECT TO THE ORIGIN OF COORDINATES, IN CM. 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 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 CARD 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 CARD, 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 CARD, 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 CARD. 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 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, RES- PECTIVELY. TERMINATION OF COMPUTATIONS --------------------------- IF MCONT=1 IN CARD NO.2, THE CARDS 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 CARD 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 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.