C PROGRAM S E I S P L O T C ************************* C C PROGRAM SEISPLOT IS DESIGNED FOR PLOTTING OF SYNTHETIC SEISMO- C GRAMS STORED IN THE FILE LU3, GENERATED IN THE PROGRAM SYNTPL. C IT CAN BE ALSO USED TO PLOT OTHER SYNTHETIC SEISMOGRAMS WHICH C ARE STORED IN THE SAME FORM AS THOSE GENERATED IN PROGRAM SYNTPL C (SEE A SHORT INFORMATION IN THE FOLLOWING). C C ***************************************************************** C C A SHORT DESCRIPTION OF THE PROGRAM C ********************************** C C SYNTHETIC SEISMOGRAMS ARE PLOTTED IN A FRAME, IN WHICH THE C X-AXIS CORRESPONDS TO THE X-AXIS IN THE MODEL FOR THE RECEIVERS C SITUATED ALONG THE EARTH'S SURFACE OR ALONG AN INTERFACE, AND TO C THE Z-AXIS FOR THE RECEIVERS SITUATED ALONG A VERTICAL PROFILE. C THE Y-AXIS OF THE FRAME CORRESPONDS TO THE TRAVEL TIME OR THE C REDUCED TRAVEL TIME. IT IS POSSIBLE TO USE VARIOUS FORMS OF AMP- C LITUDE SCALING ALONG THE PROFILE OF RECEIVERS, I.E. ALONG THE X- C AXIS OF THE FRAME). THE SCALING OF SYNTHETIC SEISMOGRAMS ALONG C TIME AXIS IS NOT USED, BUT MAY BE SIMPLY INTRODUCED. A SELECTION C OF RECEIVERS AT WHICH THE SYNTHETIC SEISMOGRAMS SHOULD BE COMPUTED C CAN BE PERFORMED. C CALCOMP PLOTTING ROUTINES ARE REQUIRED. C C ****************************************************************** C C INPUT DATA C ********** C C INPUT DATA CONSIST PARTIALLY OF THE DATA GENERATED BY A PROGRAM C SYNTPL AND STORED IN A FORMATTED FORM IN THE FILE LU3, AND PAR- C TIALLY OF THE ADDITIONAL INPUT DATA CONTROLLING PLOTTING, INTRO- C DUCED DIRECTLY (CARD READER, TERMINAL). C C DATA STORED IN LU3 C ****************** C IN THE FILE LU3, THE COMPUTED RAY SYNTHETIC SEISMOGRAMS AT INDI- C VIDUAL RECEIVER POSITIONS ARE STORED, TOGETHER WITH SOME OTHER C RELEVANT IMPORTANT INFORMATION. C THE DATA IN LU3 ARE STORED IN THE FOLLOWING SUCCESSION: C 1) MTEXT FORMAT(17A4) C ARBITRARY ALPHANUMERIC TEXT DESCRIBING THE COMPUTATIONS. THIS C TEXT WILL APPEAR UNDER THE PLOTS. IT WAS READ IN IN THE PRO- C GRAM SEIS88. C 2) MDIST,MRED,MCOMP,ITPR,VRED,RSTEP,XSOUR,DT FORMAT(4I5,4F10.5) C MDIST... THE NUMBER OF RECEIVER POSITIONS. C MRED... MRED=0:NON REDUCED TRAVEL TIMES, C MRED=1: REDUCED TRAVEL TIMES. C MCOMP... MCOMP=0: VERTICAL COMPONENT. C MCOMP=1: RADIAL COMPONENT. C MCOMP=2: TRANSVERSE COMPONENT. C ITPR... ITPR=0: RECEIVERS ALONG THE EARTH'S SURFACE. C ITPR=1: RECEIVERS ALONG A VERTICAL PROFILE. C ITPR.GT.100: RECEIVERS ALONG AN INTERFACE. C VRED... REDUCTION VELOCITY C RSTEP... AVERAGE DIFFERENCE BETWEEN X-COORDINATES OF C NEIGHBOURING RECEIVER POSITIONS. C XSOUR... X-COORDINATE OF THE SOURCE. C DT... TIME STEP IN SYNTHETIC SEISMOGRAM. C 3) XMX,SMAXIM FORMAT(22X,F10.5,9X,E15.9) C XMX... COORDINATE OF OF THE RECEIVER AT WHICH THE C MAXIMUM AMPLITUDE 'SMAXIM' WAS RECORDED. C SMAXIM... THE MAXIMUM AMPLITUDE OF ALL TRACES. C THE FOLLOWING CARDS 4,5 ARE SUCCESSIVELY REPEATED FOR ALL C RECEIVER POSITIONS: C 4) XX,SMAX,TMIN,NPTS FORMAT(F10.5,E15.9,F10.5,I5) C XX... COORDINATE OF THE RECEIVER (X-COORDINATE WHEN C RECEIVERS ARE SITUATED ALONG THE EARTH'S SURFACE C OR ALONG AN INTERFACE AND Z-COORDINATE WHEN RE- C CEIVERS ARE SITUATED ALONG A VERTICAL PROFILE. C SMAX... MAXIMUM AMPLITUDE IN THE SYNTHETIC SEISMOGRAM C AT XX. C TMIN... THE TIME CORRESPONDING TO THE FIRST NON-ZERO C POINT IN THE SYNTHETIC SEISMOGRAM. C NPTS... NUMBER OF POINTS IN SYNTHETIC SEISMOGRAM. C 5) SYNTHETIC SEISMOGRAM FORMAT(20I4) C THE MAXIMUM LENGTH OF ONE SYNTHETIC SEISMOGRAM IS 3001 POINTS. C ONLY THE PARTS OF THE SEISMOGRAMS SITUATED BETWEEN THE FIRST C AND THE LAST NON-ZERO ARRIVAL ARE STORED. EACH SEISMOGRAM IS C STORED IN A NORMALIZED FORM, SEE THE FORMAT, WITH MAXIMUM C AMPLITUDE EQUAL 999. C C ADDITIONAL INPUT DATA C ********************* C C THE FOLLOWING DATA SHOULD BE INTRODUCED DIRECTLY (CARDS, TERMI- C NAL). THEY CONTROL THE PLOTTING OF SYNTHETIC SEISMOGRAMS. THE C DATA FROM LU3 ARE READ IN AFTER ALL THE ADDITIONAL INPUT DATA C ARE READ IN. C C 1) ONE CARD C LU3,ISHIFT,IPRINT FORMAT(16I5) C LU3... THE NUMBER OF THE FILE IN WHICH THE SYNTHETIC C SEISMOGRAMS ARE STORED. C ISHIFT...THE SHIFT BETWEEN INDIVIDUAL PICTURES, IN CM. C IPRINT...CONTROLS THE PRINTOUT OF SOME DATA ON THE LINE C PRINTER. SEE DETAILS IN THE SECTION ON OUTPUT C TABLES. C C 2) ONE CARD, VARIOUS SWITCHES C MCONT,MEPIC,NTICX,NTICY,NDX,NDY FORMAT(16I5) C MCONT... MCONT=0: THE COMPUTATION TERMINATES. C MCONT=1: A NEW PLOT OF SYNTHETIC SEISMOGRAMS C WILL BE MADE FROM THE SAME FILE AS THE PRECE- C DING PLOT. C MCONT=-1: A NEW PLOT OF SYNTHETIC SEISMOGRAMS C WILL BE MADE FROM A NEW FILE. GO TO THE INPUT C DATA CARD NO.1. C MEPIC... CONTROLS THE SELECTION OF RECEIVER POSITIONS C MEPIC=0: NO SELECTION C MEPIC=1: SELECTION APPLIED, SEE CARD NO.3 C NTICX... THE NUMBER OF MARKED INTERVALS ALONG THE X-AXIS C BETWEEN TWO ADJOINING TICS DENOTED BY COORDI- C NATE VALUES. C NTICY... THE SAME AS NTICY, BUT FOR THE TRAVEL-TIME AXIS. C NDX,NDY... CONTROL THE PRECISION OF NUMBERS DESCRIBING C THE COORDINATE AXES IN THE PLOTS. C ND.GT.0: THE NUMBER OF DIGITS TO THE RIGHT OF C THE DECIMAL POINT. C ND=0: ONLY INTEGER PORTIONS OF THE NUMBERS WITH C DECIMAL POINTS. C ND.LT.0: INTEGERS. C C 3) ONE CARD, SELECTION OF RECEIVER POSITIONS. C INCLUDED ONLY WHEN MEPIC.NE.0. C NEPIC,(IEP(I),I=1,NEPIC) FORMAT(16I5) C NEPIC... NUMBER OF SELECTED RECEIVER POSITIONS AT WHICH C SYNTHETIC SEISMOGRAMS ARE TO BE CONSTRUCTED. C IEP(1),IEP(2),...,IEP(NEPIC)... SEQUENTIAL NUMBERS OF C SELECTED RECEIVER POSITIONS (FROM THE LEFT TO C THE RIGHT). C NOTE: SIMILAR SELECTION CAN BE DONE ALREADY IN THE PROGRAM C SYNTPL. THE SELECTION HERE IS MADE AMONG THOSE RECEIVERS WHICH C REMAINED AFTER THE SELECTION IN SYNTPL. C C 4) ONE CARD, DESCRIPTION OF AXES IN THE PLOT C XMIN,XMAX,XLEN,DTICX,YMIN,YMAX,YLEN,DTICY... FORMAT(8F10.5) C XMIN,XMAX...THE MINIMUM AND MAXIMUM VALUES AT THE X-AXIS C (IN THE USERS LENGTH UNITS) C XLEN... LENGTH OF THE X-AXIS (IN CM) C DTICX... THE DISTANCE BETWEEN TWO NEIGHBOURING MARKS ON C THE X-AXIS WHICH ARE DENOTED BY CORRESPONDING C COORDINATE VALUES (IN THE USERS LENGTH UNITS). C DTICX.GT.0.: TIC MARKS STARTING FROM XMIN AND C APPEARING AT THE SUBSEQUENT POINTS XMIN+DTICX, C XMIN+2.*DTICX,... C DTICX.LT.0.: TIC MARKS START AND CONTINUED TO C BE PLOTTED FROM THE FIRST INTEGER MULTIPLE OF C ABS(DTICX) GREATER THAN XMIN. C YMIN,YMAX,YLEN,DTICY... THE SAME FOR THE TIME AXIS. FOR C REDUCED TIME AXIS (MRED=1), THE QUANTITIES C CORRESPOND TO THE REDUCED TRAVEL TIMES C C C 5) ONE CARD. PARAMETERS OF THE AMPLITUDE SCALING ALONG THE X-AXIS C AMP,B1,EPICS,EPS,SC FORMAT(8F10.5) C WHEN EPICS NOT SPECIFIED, EPICS=10. C WHEN B1 NOT SPECIFIED, B1=1. C THE PARAMETERS AMP,B1,EPICS,EPS CONTROL THE SCALING OF SEIS- C MOGRAMS. THE SCALING FACTOR OF SYNTHETIC SEISMOGRAMS MAY GENE- C RALLY DEPEND ON THE RECEIVER POSITION. DENOTE THE COORDINATE C OF THE RECEIVER BY XX, AND THE AVERAGE DISTANCE BETWEEN NEIGH- C BOURING RECEIVERS BY DDX (I.E., THE AVERAGE DISTANCE BETWEEN C INDIVIDUAL TRACES, IN PLOT UNITS). THEN THE FOLLOWING FIVE C SYSTEMS OF SCALING CAN BE USED: C C 1) FACTOR=B1*DDX/SMAXI, C WHEN AMP=0. (THE WHOLE CARD 5 MAY BE BLANK IN THIS CASE). C NORMALIZATION OF MAXIMUM AMPLITUDES IN ALL TRACES: C FOR B1=1, THE MAXIMUM AMPLITUDES IN ALL TRACES EQUAL TO C THE AVERAGE DISTANCE BETWEEN INDIVIDUAL TRACES. C 2) FACTOR=B1*DDX/SMAXIM, C WHEN AMP=-1 AND EPS=0. C THE WHOLE CARD 5 WITH THE EXCEPTION OF AMP MAY BE BLANK IN C THIS CASE. C REAL AMPLITUDES: FOR B1=1, THE MAXIMUM AMPLITUDE IN THE C WHOLE PLOT EQUALS TO THE DISTANCE BETWEEN INDIVIDUAL TRACES C 3) FACTOR=B1, C WHEN AMP=1 AND EPS=0. C REAL AMPLITUDES: MANUAL DETERMINATION OF THE SCALING FACTOR C BY SPECIFYING B1. C 4) FACTOR=B1*DDX*((ABS(XX-XSOUR)/EPICS)**EPS)/SMAXIM C WHEN AMP=-1. AND EPS.GT.0. C POWER SCALING OF AMPLITUDES: WHEN EPICS CORRESPONDS TO THE C EPICENTRAL DISTANCE OF THE TRACE WITH THE MAXIMUM AMPLITUDE C AND B1=1, THEN THE MAXIMUM AMPLITUDE OF THE TRACE AT EPICS C EQUALS TO THE AVERAGE DISTANCE BETWEEN INDIVIDUAL TRACES. C 5) FACTOR=B1*(ABS(XX-XSOUR)/EPICS)**EPS C WHEN AMP=1. AND EPS.GT.0. C POWER SCALING OF AMPLITUDES: FULLY MANUAL, ALL THE PARA- C METERS MUST BE SPECIFIED. C NOTE THAT THE SCALING SYSTEMS 4 AND 5 GIVE FOR XX=XSOUR TRACES C WITH ZERO AMPLITUDES. C SC... CONTROLS THE SCALES OF TICS AND ALPHANUMERICAL TEXTS. C FOR SC=1.0, THE TICS ARE 0.15 CM LONG AND COORDINATES AND C TEXT DESCRIBING THE PLOTS ARE 0.4 AND 0.45 CM HIGH, RES- C PECTIVELY. C C C **************************************************************** C C TERMINATION OF COMPUTATIONS C *************************** C C IF MCONT=1 IN CARD NO.2, THE CARDS 2-5 CAN BE REPEATED AN ARBI- C TRARY NUMBER OF TIMES TO PLOT SYNTHETIC SEISMOGRAMS WITH DIFFE- C RENT AMPLITUDE SCALING, DIFFERENT FRAMES, ETC. C IF MCONT=-1 IN CARD NO.2, THE FILE LU3 IS CLOSED, AND A NUMBER C OF A NEW FILE IS READ IN FROM INPUT DATA CARD NO.1. THE FILE IS C OPENED AND THE COMPUTATIONS PROCEED AS DESCRIBED ABOVE. NOTE C THAT THIS POSSIBILITY CAN BE USED ONLY IN OPERATING SYSTEMS C ALLOWING SPECIFICATION OF THE NUMBER OF A FILE OUTSIDE THE SOUR- C CE PROGRAM (E.G., IBM OPERATING SYSTEM). C THE COMPUTATION TERMINATES WHEN MCONT=0 IN THE INPUT DATA CARD C NO.2. C C ****************************************************************** C C OUTPUT ON A LINE PRINTER C ************************ C C ALL THE ADDITIONAL INPUT DATA ARE REPRODUCED ON THE LINE PRINTER. C THE PRINTOUT OF OTHER DATA IS CONTROLLED BY THE PARAMETER IPRINT, C SEE INPUT DATA CARD NO.1. C FOR IPRINT=0: ALSO THE DATA LU3/1,LU3/2 AND LU3/3 ARE PRINTED. C FOR IPRINT=1: THE SAME AS FOR IPRINT=0. IN ADDITION, ONE LINE IS C PRINTED FOR EACH RECEIVER POSITION: C XX,SMAX,FACTOR,SFMAX C WHERE XX IS THE COORDINATE OF THE RECEIVER, SMAX IS THE MAXIMUM C AMPLITUDE IN THE CORRESPONDING SEISMOGRAM, FACTOR IS THE AMPLITUDE C SCALING FACTOR CORRESPONDING TO THIS RECEIVER, AND SFMAX IS THE C MAXIMUM AMPLITUDE IN THE SYNTHETIC SEISMOGRAM AFTER SCALING. C C ***************************************************************** C C