C PROGRAM R A Y P L O T C *********************** C C PROGRAM RAYPLOT IS DESIGNED FOR PLOTTING OF RAY DIAGRAMS, C TRAVEL TIMES AND AMPLITUDES OF SEISMIC BODY WAVES FROM ANY FILE C WHICH HAS A FORM OF THE FILE LU GENERATED FOR THIS PURPOSE BY C PROGRAM SEIS88. C ALTERNATIVELY, PROGRAM RAYPLOT CAN BE USED FOR THE INSPECTION C OF THE CONTENT OF THE FILE LU, WHICH IS WRITTEN IN BINARY FORM. C THE FILE LU CAN BE GENERATED E.G. BY PROGRAM SEIS88. C C *************************************************************** C C A SHORT DESCRIPTION OF THE PROGRAM. C ********************************** C C RAY DIAGRAMS, TRAVEL TIMES AND AMPLITUDES ARE PLOTTED IN THREE C DIFFERENT FRAMES. FOR EACH ELEMENTARY WAVE, FOR WHICH COMPUTATIONS C WERE PERFORMED, IT IS POSSIBLE TO PLOT ALL THE PICTURES OR ONLY C SOME OF THEM. SEVERAL AMPLITUDE PICTURES MAY BE PLOTTED FOR ONE C ELEMENTARY WAVE, CORRESPONDING TO DIFFERENT COMPONENTS OR, IF C SLIGHT ABSORPTION IS CONSIDERED, TO DIFFERENT FREQUENCIES. THE C HORIZONTAL COORDINATE IN THE RAY DIAGRAM CORRESPONDS TO THE RANGE C IN THE MODEL, AND IT INCREASES FROM THE LEFT TO THE RIGHT. THE C VERTICAL COORDINATE IN THE RAY DIAGRAM PLOTS CORRESPONDS TO THE C DEPTH. IN TRAVEL TIME AND AMPLITUDE PLOTS, THE HORIZONTAL AXIS C CORRESPONDS TO THE RANGE WHEN THE RECEIVERS ARE SITUATED ALONG THE C EARTH'S SURFACE OR ALONG AN INTERFACE. FOR RECEIVERS SITUATED C ALONG A VERTICAL PROFILE, THE HORIZONTAL AXIS OF TRAVEL TIME AND C AMPLITUDE PLOTS CORRESPONDS TO THE DEPTH IN THE MODEL. THE VERTICAL C COORDINATE CORRESPONDS TO THE TRAVEL TIME OR REDUCED TRAVEL TIME C PLOTS AND DECADIC LOGARITHM OF AN AMPLITUDE FOR AMPLITUDE PLOTS. C ALL THE RAYS STORED IN LU OR ONLY SELECTED RAYS FROM IT CAN BE C PLOTTED. THE SELECTION CAN BE DONE IN TWO WAYS: C A) GLOBAL SELECTION, ACCORDING TO THE POSITION OF THE TERMI- C NATION POINT OF THE RAY. THE INDEX "IND" (SEE THE DESCRIPTION C OF PROGRAM SEIS88) IS USED FOR THIS SELECTION. IN THIS WAY, IT C IS POSSIBLE TO PLOT THE RAYS WITH THE TERMINATION POINTS ALONG C THE EARTH'S SURFACE, ALONG THE LEFT VERTICAL BORDER, RIGHT VERTI- C CAL BORDER, OR ALONG SOME INTERFACE. IT IS EVEN POSSIBLE TO SE- C LECT AND PLOT THE RAYS WHICH DO NOT TERMINATE PROPERLY IN COR- C RESPONDENCE WITH THE NUMERICAL CODE OF THE WAVE (IND=8, IND=9, C ETC.). GLOBAL SELECTION IS NOT APPLICABLE TO THE FILE LU GENE- C RATED BY PROGRAM SEIS88. THE FILE GENERATED BY PROGRAM SEIS88 C CONTAINS RAYS WHICH ARE ALREADY SELECTED. THEY HAVE TERMINATION C POINTS ALONG THE EARTH'S SURFACE (IND=3), AT A SPECIFIED SYS- C TEM OF RECEIVERS. C B) MANUAL SELECTION, USING SUCCESSIVE NUMBERS OF RAYS, SEE C INDEX NRAY IN INPUT DATA CARD NO.3. C THE TRAVEL TIMES OBTAINED EXPERIMENTALLY OR IN SOME OTHER C WAY CAN BE PLOTTED TOGETHER WITH COMPUTED TRAVEL TIMES (FOR C APPLICATIONS IN MATHEMATICAL MODELLING). THE SAME APPLIES C TO AMPLITUDES. AMPLITUDES OF RADIAL, TRANSVERSE AND VERTICAL C COMPONENTS OF DISPLACEMENT VECTOR CAN BE PLOTTED. C CALCOMP PLOTTING ROUTINES ARE REQUIRED. C C **************************************************************** C C INPUT DATA C ********** C C INPUT DATA CONSIST PARTIALLY OF THE DATA GENERATED IN THE PROGRAM C SEIS88 (OR ANY OTHER PROGRAM WHICH GENERATES A FILE OF THE SAME C STRUCTURE AS SEIS88), STORED IN THE BINARY FORM IN THE FILE LU, C AND PARTIALLY OF THE ADDITIONAL INPUT DATA CONTROLLING PLOTTING, C INTRODUCED DIRECTLY (CARD READER, TERMINAL). C C THE DATA STORED IN LU. C ********************* C FOR DETAILS, SEE THE DESCRIPTION OF THE CONTENTS OF THE FILE LU1 C IN PROGRAM SEIS88. C THE DATA ARE STORED IN THE FOLLOWING ORDER: C 1) ICONT,ITPR C 2) NINT, (NPNT(I),I=1,NINT) C 3) FOR I=1,2,...,NINT: C (A1(J,I),B1(J,I),C1(J,I),D1(J,I),X1(J,I),III(J,I),J=1,NC), C WHERE NC=NPNT(I)-1. C 4) X0,Z0,ROS,VPS,VSS C 5) N,IND C 6) (X(J),Y(J),J=1,N) C THE DATA 5,6 ARE REPEATED SUCCESSIVELY FOR ALL RAYS. AFTER C THE LAST RAY, N=0. C 7) NS C 8) (INDI(I),X(I),T(I),TAS(I),ANG(I),AX(I),AY(I),AZ(I),PHX(I),PHY(I), C PHZ(I),I=1,NS) C AFTER THIS, DATA FOR A NEW RAY DIAGRAM ARE STORED, STARTING C FROM 1). ICONT=0 IS AN INDICATION OF THE END OF THE FILE LU. C C THE ADDITIONAL INPUT DATA C ************************* C THE ADDITIONAL INPUT DATA CAN BE INTRODUCED BY A CARD READER C OR BY A TERMINAL. THE PLACES, WHERE THE DATA FROM LU ARE READ IN, C ARE DENOTED BY **LU/1, **LU/2, ETC. C C 1)ONE CARD C LU,ISHIFT,IPRINT FORMAT(16I5) C LU... THE NUMBER OF THE FILE IN WHICH THE RAY DIAGRAMS, C TRAVEL TIMES AND AMPLITUDES ARE STORED. C IF LU NOT SPECIFIED, LU=7. C ISHIFT... SHIFT BETWEEN INDIVIDUAL PICTURES, IN CM. C IF NOT SPECIFIED, ISHIFT=8. C IPRINT... CONTROLS THE PRINTING OF SOME DATA AND RESULTS C ON THE LINE PRINTER. C IPRINT=0: ONLY ADDITIONAL INPUT DATA ARE REPRO- C DUCED. C IPRINT=1: IN ADDITION TO IPRINT=0, ALSO THE TABLES C OF VALUES OF AMPLITUDES SHOWN IN THE PLOTS. C IPRINT=2: IN ADDITION TO IPRINT=0, ALSO THE TABLES C OF TRAVEL TIMES AND AMPLITUDES STORED IN THE FILE C LU. C IPRINT=3: IN ADDITION TO IPRINT=2, ALSO THE TABLES C OF RAYS. SEE MORE DETAILS ABOUT THE OUTPUT IN C THE END OF THE DESCRIPTION. C **LU/1 C C 2) ONE CARD C ITEXT FORMAT(20A4) C ARBITRARY ALPHANUMERIC TEXT WHICH IS SHOWN BELOW THE PLOTS. C C 3) ONE CARD, VARIOUS SWITCHES. C NTICX,NTICY,NTICT,NTICA,INDO,INDT,ICOMP,NRAY, C IBOUND,IRED,IRS,NDX,NDY FORMAT(16I5) C NTICX... NTICX.NE.0: THE NUMBER OF MARKED INTERVALS ON C THE HORIZONTAL RANGE AXIS BETWEEN TWO NEIGHBOUR- C HOOD TICS WITH CORRESPONDING COORDINATE VALUES. C NTICX=0: THE PLOTTING TERMINATES, INDEPENDENTLY C OF THE CONTENT OF LU. C NTICY... NTICY.NE.0: THE SAME AS NTICX, BUT FOR THE DEPTH C AXIS. C NTICY=0: NO RAY DIAGRAM WILL BE PLOTTED C NTICT... NTICT.NE.0: THE SAME AS NTICX, BUT FOR A TIME C AXIS C NTICT=0: NO PLOT OF TRAVEL TIMES. C NTICA... NTICA.NE.0: THE SAME AS NTICX, BUT FOR THE AMPLI- C TUDE AXIS. C NTICA=0: NO PLOT OF AMPLITUDES. C INDO... CONTROLS THE GLOBAL SELECTION OF RAYS IN THE RAY C DIAGRAM. IN THE RAY DIAGRAM, ONLY THE RAYS WITH C IND=INDO ARE PLOTTED. C IF INDO=0, ALL RAYS STORED AT LU ARE PLOTTED, IF C THEY ARE NOT EXCLUDED DURING THE MANUAL SELECTION C OF RAYS, SEE NRAY. C INDT... CONTROLS THE GLOBAL SELECTION OF TRAVEL TIMES C IN THE PLOT OF THE TRAVEL-TIMES. ONLY THE TRAVEL C TIMES FOR IND=INDT ARE PLOTTED. FOR INDT=1,2,3, C AMPLITUDES ARE ALSO PLOTTED. IF INDT NOT SPECI- C FIED, INDT=3. C NRAY... CONTROLS THE MANUAL SELECTION OF RAYS FOR PLOT- C TING. C FOR NRAY.NE.0, NRAY IS A NUMBER OF RAYS STORED C IN LU WHICH WILL NOT APPEAR IN THE RAY DIAGRAM. C THESE RAYS ARE SPECIFIED IN CARD NO.6. C NRAY=0: ALL RAYS ARE PLOTTED (WHEN, OF COURSE, C THEY ARE NOT EXCLUDED IN THE GLOBAL SELECTION). C IBOUND... IABS(IBOUND)... NUMBER OF POINTS USED IN THE PLOT C OF THE RAY DIAGRAMS TO APPROXIMATE INTERFACES. C IBOUND.LT.0... TO DISTINGUISH INTERFACES FROM C RAYS, EACH INTERFACE IS PLOTTED BY A HEAVY LINE C (3 SLIGHTLY SHIFTED SINGLE LINES). C IBOUND.GT.0... EACH INTERFACE IS PLOTTED BY A C SINGLE LINE. C IF NOT SPECIFIED, IBOUND=100 AUTOMATICALLY. C NOTE: THE FICTITIOUS PARTS OF INTERFACES AND THE C PARTS OF INTERFACES WHICH COINCIDE WITH THOSE C DIRECTLY ABOVE THEM ARE NOT PLOTTED. C IRED... IRED=0: NO REDUCTION OF TRAVEL TIMES C IRED=1: REDUCTION OF TRAVEL TIMES, THE REDUCTION C VELOCITY IS GIVEN IN CARD NO.7. C NOTE: FOR RECEIVERS SITUATED ALONG A VERTICAL PRO- C FILE (ITPR=1), AUTOMATICALLY IRED=0. C IRS... IRS+1: NUMBER OF THE POINT (THE SAME ON ALL RAYS), C FROM WHICH PLOT OF THE RAY SHOULD START. C FOR IRS=0, RAYS ARE PLOTTED RIGHT FROM THE SOURCE. C NDX,NDY... CONTROL THE PRECISION OF NUMBERS DESCRIBING THE CO- C ORDINATE AXES IN THE PLOTS. C ND.GT.0: THE NUMBER OF DIGITS TO THE RIGHT OF THE C DECIMAL POINT. C ND.EQ.0: ONLY INTEGER PORTIONS OF THE NUMBERS WITH C DECIMAL POINTS. C ND.LT.0: INTEGERS. C C 4) ONE CARD, DESCRIPTION OF THE RANGE AXIS C XMIN,XMAX,XLEN,DTICX,SC FORMAT(8F10.5) C XMIN,XMAX... THE MINIMUM AND MAXIMUM VALUES ON THE RANGE C AXIS (IN THE USERS LENGTH UNITS). C XLEN... LENGTH OF THE RANGE AXIS, IN CM. C DTICX... THE DISTANCE BETWEEN TWO NEIGHBOURING MARKS ON C THE RANGE 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 SC... CONTROLS THE SCALES OF TICS AND ALPHANUMERICAL C TEXTS. FOR SC=1.0, THE TICS ARE 0.15 CM LONG AND C COORDINATES AND TEXTS DESCRIBING THE PLOTS ARE C 0.4 AND 0.45 CM HIGH, RESPECTIVELY. FOR SC LARGER C (LESS) THAN 1.0, THE SCALES ARE PROPORTIONALLY C INCREASED (DECREASED). DEFAULT VALUE, SC=1.0. C C NOTE: THE ABOVE PARAMETERS SPECIFY THE HORIZONTAL AXIS IN C PLOTS OF RAY DIAGRAMS FOR ANY RECEIVER CONFIGURATION AND IN C PLOTS OF TRAVEL TIMES AND AMPLITUDES FOR RECEIVERS SITUATED C ALONG THE EARTH'S SURFACE OR ALONG AN INTERFACE. FOR RECEI- C VERS SITUATED ALONG A VERTICAL PROFILE, PARAMETERS FROM INPUT C DATA CARD NO.5 ARE AUTOMATICALLY USED FOR THE SPECIFICATION C OF HORIZONTAL AXES OF PLOTS OF TRAVEL TIMES AND AMPLITUDES. C C 5)ONE CARD, DESCRIPTION OF THE DEPTH AXIS C THIS CARD IS INCLUDED ONLY WHEN NTICY.NE.0. C YMIN,YMAX,YLEN,DTICY FORMAT(8F10.5) C THE SAME AS IN 4, BUT FOR THE DEPTH AXIS. C C **LU/2 C **LU/3 C **LU/4 C C 6) ONE CARD, MANUAL SELECTION OF RAYS C INCLUDED ONLY WHEN NRAY.NE.0. C (NR(I),I=1,NRAY) FORMAT(26I3) C NR(I)... SUCCESSIVE NUMBERS OF RAYS STORED IN LU WHICH C ARE NOT TO BE PLOTTED IN THE RAY DIAGRAM. C C **LU/5 C **LU/6 C **LU/7 C **LU/8 C C 7) ONE CARD, DESCRIPTION OF THE TRAVEL-TIME AXIS C INCLUDED ONLY WHEN NTICT.NE.0. C TMIN,TMAX,TLEN,DTICT,VRED FORMAT(8F10.5) C THE SAME MEANING AS IN 4, BUT FOR THE TRAVEL-TIME AXIS. C VRED IS THE REDUCTION VELOCITY. FOR IRED=1, ALL THE QUAN- C TITIES CORRESPOND TO THE REDUCED TRAVEL TIMES. C NOTE: IN CASE OF ITPR=1, I.E. FOR RECEIVERS SITUATED ALONG C A VERTICAL PROFILE), THE HORIZONTAL AXIS OF THE PLOT COR- C RESPONDS TO THE DEPTH, NOT TO THE RANGE. THE PARAMETERS FROM C INPUT DATA CARD NO.5 ARE AUTOMATICALLY USED FOR THE DESCRIP- C TION OF THE HORIZONTAL AXIS OF THE PLOT C C 8) ONE CARD C INCLUDED ONLY WHEN NTICY.NE.0. C NEXP FORMAT(26I3) C NUMBER OF POINTS OF THE ADDITIONAL TRAVEL-TIME CURVE WHICH C IS TO BE PLOTTED IN ONE FRAME TOGETHER WITH THE COMPUTED C TRAVEL-TIME CURVE. C NEXP=0: NO ADDITIONAL TRAVEL-TIME POINTS ARE PLOTTED C C 9) ONE CARD C INCLUDED ONLY WHEN NTICY.NE.0 AND NEXP.NE.0. C (X(I),T(I),I=1,NEXP) FORMAT(8F10.5) C RECEIVER COORDINATES AND CORRESPONDING TRAVEL TIMES IN AN C ADDITIONAL TRAVEL-TIME CURVE. THE TIMES SHOULD BE IN A NON- C RE-DUCED FORM. IN THE PLOT, THEY WILL BE REDUCED IN THE C SAME WAY AS THE COMPUTED TIMES. C C 10) ONE CARD, DESCRIPTION OF THE AMPLITUDE AXIS C INCLUDED ONLY WHEN NTICA.NE.0. C AMIN,AMAX,ALEN,DTICA,FREQ,KABS,ICOMP,MSOUR FORMAT(5F10.5,2I5) C AMIN,AMAX,ALEN,DTICA... THE SAME MEANING AS IN 4, BUT FOR C THE AMPLITUDE AXIS (DECADIC LOGARITHMES OF AMPLI- C TUDES). C IF(ALEN.LT..00001), NO AMPLITUDE-DISTANCE CURVE C WILL BE PLOTTED. C FREQ... FREQUENCY FOR WHICH AMPLITUDES ARE EVALUATED C (IN CASE OF ABSORPTION, KABS.NE.0). C KABS... KABS=0... NO ABSORPTION IS CONSIDERED. C KABS=1... NON-CAUSAL OR CAUSAL ABSORPTION C (FUTTERMAN'S MODEL) IS CONSIDERED. C ICOMP... ICOMP=0: AMPLITUDES OF VERTICAL COMPONENT. C ICOMP=1: AMPLITUDES OF RADIAL COMPONENT. C ICOMP=2: AMPLITUDES OF TRANSVERSE COMPONENT. C MSOUR... CONTROLS THE CHOICE OF THE TYPE OF THE POINT C SOURCE. C MSOUR=0: UNIT ISOTROPIC RADIATION PATTERN. C MSOUR=1: SINGLE FORCE. C MSOUR=2: DOUBLE COUPLE. C MSOUR=3: EXPLOSIVE (IMPLOSIVE) SOURCE. C C NOTE1: ABSORPTION EFFECTS CAN ONLY BE CONSIDERED IF GLOBAL C ABSORPTION FACTOR 'TAS', SEE DATA STORED IN LU/8, IS NON- C ZERO. OTHERWISE KABS HAS ONLY A FORMAL MEANING. C C NOTE2: ABSORPTION AFFECTS ONLY AMPLITUDES STORED IN FILE LU, C NOT THE AMPLITUDES FORMING THE ADDITONAL AMPLITUDE CURVES C READ IN IN THE INPUT DATA CARD NO.13. C C NOTE3: THE NOTE FROM THE DESCRIPTION OF THE INPUT DATA CARD C NO.7 HOLDS ALSO FOR AMPLITUDES. C C 11) ONE CARD, SPECIFICATION OF PARAMETERS OF THE SOURCE C INCLUDED ONLY WHEN MSOUR.NE.0. C IPAR(1),...,IPAR(4),PAR(1),...PAR(6) FORMAT(4I5,6F10.5) C FOR THE MEANING OF THE ABOVE PARAMETERS SEE DESCRIPTION OF C THE INPUT DATA CARD NO.7 IN PROGRAM SYNTPL IN THIS PACKAGE. C C 12) ONE CARD C INCLUDED ONLY WHEN NTICA.NE.0. C NEXP FORMAT(26I3) C NUMBER OF POINTS IN THE ADDITIONAL AMPLITUDE CURVE WHICH C IS TO BE PLOTTED IN THE SAME FRAME AS THE COMPUTED AMPLI- C TUDE CURVE. FOR NEXP=0, NO ADDITIONAL AMPLITUDE CURVE IS C PLOTTED. C C 12) ONE CARD C INCLUDED ONLY WHEN NTICA.NE.0 AND NEXP.NE.0. C (X(I),Y(I),I=1,NEXP) FORMAT(8F10.5) C RECEIVER COORDINATES AND CORRESPONDING AMPLITUDES OF THE C ADDITIONAL AMPLITUDE CURVE. REAL AMPLITUDES, NOT LOGA- C RITHMES OF AMPLITUDES, SHOULD BE SPECIFIED. C C C **************************************************************** C C TERMINATION OF COMPUTATIONS C *************************** C C THE INPUT DATA CARDS NO.10-13 CAN BE REPEATED FOR THE GIVEN C ELEMENTARY WAVE ANY NUMBER TIMES (TO PERFORM PLOTTING OF AMPLITUDE- C DISTANCE CURVES OF THE ELEMENTARY WAVE E.G. FOR DIFFERENT FREQUEN- C CIES, IF ABSORPTION IS CONSIDERED, OR FOR DIFFERENT COMPONENTS, OR C FOR DIFFERENT SPECIFICATIONS OF THE SOURCE). THE LAST CARD SHOULD C BE A BLANK CARD (TO GUARANTEE ALEN.LT..00001). AFTER THE BLANK CARD, C A NEW SYSTEM OF CARDS 2-13 FOR A NEW ELEMENTARY WAVE STORED IN THE C FILE LU MAY FOLLOW. C C C AS SOON AS ICONT=0 IS FOUND IN LU, THE PLOTTING FROM THE FILE LU C TERMINATES. THE PLOTTING FROM THE FILE LU MAY BE ALSO TERMINATED C IF WE PUT NTICX=0 IN THE CARD NO.3, EVEN IF LU STILL CONTAINS C SOME OTHER DATA. IN BOTH CASES, THE FILE IS CLOSED AND THE COM- C PUTATIONS TERMINATE. C C ****************************************************************** C C THE OUTPUT ON A LINE PRINTER C **************************** C C THE OUTPUT IS CONTROLLED BY THE PARAMETER IPRINT, SEE INPUT DATA C CARD NO 1. C FOR IPRINT=0, ALL THE ADDITIONAL DATA SHOWN ABOVE AND THE DATA C LU/1, LU/4, LU/5 AND LU/7 ARE PRINTED. C FOR IPRINT=1, THE DATA FOR IPRINT=0 AND ALSO THE RECEIVER COORDI- C NATES AND VALUES OF CORRESPONDING AMPLITUDES WHICH APPEAR IN THE C PLOTS ARE PRINTED. C FOR IPRINT=2, THE DATA FOR IPRINT=0 AND ALSO THE DATA LU/8 ARE C PRINTED. C FOR IPRINT=3, THE DATA FOR IPRINT=2 AND ALSO ALL THE DATA FROM C LU ARE PRINTED, EVEN THE COORDINATES OF ALL RAYS. C C **************************************************************** C