***************************************** Description of the program A N R A Y P L ***************************************** Program ANRAYPL is designed for plotting of ray diagrams (projections of rays into vertical and horizontal planes), travel times and amplitudes of seismic body waves from any file which has formatting of data like the file LU1 generated for this purpose by program ANRAY. Short description of the program ANRAYPL ---------------------------------------- Ray diagrams, travel times and amplitudes are plotted in separate frames. For each elementary wave, for which computations were performed, it is possible to plot all frames or only some of them. Two ray diagram frames may be plotted for each elementary wave, containing projections of rays into horizontal and vertical planes. Three amplitude frames may be plotted for each elementary wave, corresponding to three different components of the polarization vector: vertical, radial and transverse. The horizontal axis in the ray diagram frames coresponds to the range. Vertical axis corresponds either to the depth (in case of projections of rays into a vertical plane) or to the distance transverse to range (in case of projections of rays into a horizontal plane). The z-axis is positive downwards. The horizontal axis in the frames of travel times and ray amplitudes corresponds to the range (in case of the profile with receivers situated along surface or along an interface) or to the depth (in case of vertical profile with receivers). The z-axis, corresponding to travel times or amplitudes, is positive upwards. Note: Contours of the interfaces in ray diagrams are not plotted. The travel times and amplitudes of vertical, radial or transverse components obtained experimentally or in some other way can be plotted together with computed travel times. The following CALCOMP plotting routines are required: PLOTS, PLOT, NUMBER and SYMBOL. Description of input and output data ------------------------------------ Input data consist partially of the data obtained in the program ANRAY, which are stored in the file LU1, and partially of the additional data specified by the user in the file LIN. Output data describing the computations are stored in the file LOU. Specification of the files LIN, LOU and LU1 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 LU1 is read from the file LIN, see below. The data stored at LU1 ---------------------- The data are stored in the file LU1 in the formatted form. For details see the description of the content of the file LU1 in program ANRAY. 1) ICONT,NDST,ILOC FORMAT(26I3) 2) VP,VS,RO FORMAT(8F10.5) 3) NPN,NPN,NPN FORMAT(26I3) 4) APN,APN,APN,APN,APN FORMAT(5E15.5) 5) XSOUR,YSOUR,ZSOUR,PROF FORMAT(8F10.5) 6) (DST(I),I=1,NDST) FORMAT(8F10.5) 7) N,II FORMAT(26I3) 8) (AY(2,I),AY(3,I),AY(4,I),I=1,N) FORMAT(3E15.5) 9) INDEX FORMAT(26I3) 10)(INDI(I),XCOOR(I),ZCOOR(I),TIME(I),(AMPX(J,I),AMPY(J,I), AMPZ(J,I),J=1,2),AZM(I),DCL(I),I=1,INDEX) FORMAT(I5,3E15.5/5x,6e15.5/5x,8e15.5) The sequence of data 1-10 repeats for each considered elementary wave. Indication of the end of information in LU1 is ICONT=0 in line 1. The additional input data in the file LIN ----------------------------------------- These data are specified by the user. They control plotting of ray, travel time and amplitude diagrams. The places, where the data from LU1 are read in are denoted by **LU1/1, **LU1/2, etc. 1) ISHIFTX,ISHIFTY,IPRINT,LU1 FORMAT(16I5) ISHIFTX,ISHIFTY... shifts between individual pictures in x and y directions, in cm. If not specified, the default value is ISHIFT=8. IPRINT... controls the output of some data and results into the file output file LOU. IPRINT=0: Only additional input data are reproduced. IPRINT=1: In addition to the previous data, the tables of travel times and amplitudes are given. IPRINT=2: In addition, the tables of rays are printed. LU1... the number of the file in which ray diagrams, travel times and amplitudes are stored. It has no meaning if routine SERV is used. **LU1/1 2) TEXT FORMAT(20A4) Arbitrary alphanumerical text which is shown below the plots. 3) NTICX,NTICY,NTICH,NTICT,NTICA,IPROF,NRAY,IBOUND,IRED,IRS, NDX,NDZ,NDH FORMAT(16I5) NTICX... NTICX.NE.0: The number of marked intervals on the x axis between two neighbouring tics with corresponding coordinate values. NTICX=0: Indication of termination of the plotting, independently of the content of LU1. NTICY... NTICY.NE.0: The same as NTICX, but for the depth axis in the plot of the diagram of ray projections into the vertical plane containing the profile. NTICY=0: No ray diagram is to be plotted. NTICH... NTICH.NE.0: The same as NTICX, but for the transverse axis in the plot of the diagram of ray projections into the horizontal plane. NTICH=0: No ray diagram is to be plotted. NTICT... NTICT.NE.0: The same as NTICX, but for a time axis. NTICT=0: No travel times are to be plotted. NTICA... NTICA.NE.0: The same as NTICX, but for the amplitude axis. NTICA=0: No amplitudes are to be plotted. IPROF... dummy parameter in this version. NRAY... controls the manual selection of rays for plotting. NRAY.NE.0: NRAY is a number of rays stored in LU1 which are not to appear in the ray diagram. These rays are specified in card No.6. NRAY=0: Rays to all receivers along the considered profile are plotted. IBOUND... dummy parameter in this version. IRED... IRED=0: No reduction of travel times. IRED=1: Reduction of travel times; the reduction velocity is given in card No.7. IRS... IRS=0: Rays are plotted right from the source. IRS.NE.0: IRS+1 is the number of the point on any ray from which plot of the ray is to start. NDX,NDY,NDH... control the precision of numbers describing the coordinate axes in the plots. NDX corresponds to range axis. The same value NDY is used for ray, travel time and amplitude axes. NDH corresponds to the transverse axis in the plots of horizontal projections of rays. 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. 4) XMIN,XMAX,XLEN,DTICX,SC,AROT FORMAT(8F10.5) XMIN,XMAX... the minimum and maximum values along the range axis (in the user length units, e.g. km). Remember, the epicenter is automatically considered to be situated at x=0. XLEN... length of the x axis, in cm. DTICX... the distance between two neighbouring tics on the x axis which are denoted by the corresponding coordinate values (in the user length units). DTICX.GT.0.0: Tic marks starting from XMIN and appearing at the subsequent points XMIN+DTICX, XMIN+2.0*DTICX,... DTICX.LT.0.0: Tic marks start and continue to be plotted from the first integer multiple of ABS(DTICX) greater than XMIN. SC... scaling factor, controls the scales of tics and alphanumeric texts. For SC=1.0, the tics are 0.15 cm long and coordinates and text describing the plot are 0.4 and 0.45 cm high, respectively. AROT... an angle by which the horizontal coordinate axes of the recording coordinate system are to be rotated. For AROT=0, the radial axis is parallel to the horizontal projection of the line connecting source with receivers. Transverse axis is perpendicular to it. For AROT.NE.0, the axes are rotated anticlockwise in the horizontal plane by the angle AROT. 5) YMIN,YMAX,YLEN,DTICY,HMIN,HMAX,HLEN,DTICH FORMAT(8F10.5) The same as in card No.4, but for the depth axis in the frame for the diagram of ray projections into the vertical plane (quantities with Y) and for the diagram of ray projections into the horizontal plane (quantities with H). **LU1/2 **LU1/3 **LU1/4 **LU1/5 6) This card is included only when NRAY.NE.0. (NR(I),I=1,NRAY) FORMAT(26I3) NR(I)... successive numbers of rays, which are not to be plotted in the ray diagrams. Rays are numbered by the indices of receivers to which rays arrive, see the field DST in the data stored in LU1, line 6. **LU1/6 **LU1/7 **LU1/8 **LU1/9 7) This card is included only when NTICT.NE.0. TMIN,TMAX,TLEN,DTICT,VRED FORMAT(8F10.5) The same meaning as in the case of the x axis, but for travel time axis. VRED is the reduction velocity. For IRED=1, all the above input quantities correspond to the reduced travel times. The default value of VRED is VRED=6.0. 8) This card is included only when NTICT.NE.0. NEXP FORMAT(26I3) NEXP... NEXP.NE.0: Number of points of the additional travel-time curve which is to be plotted in one frame together with the computed travel-time curve. NEXP=0: No additional travel-time points are to be plotted. 9) This card is included only when NTICT.NE.0 and NEXP.NE.0. (X(I),T(I),I=1,NEXP) FORMAT(8F10.5) X(I),T(I)... x coordinates and corresponding travel times in an additional travel-time curve. The times should be in a non-reduced form. In the plot, they will be reduced in the same way as the computed travel times. 10) This card is included only when NTICA.NE.0. AMIN,AMAX,ALEN,DTICA,ICOMP,MSOUR FORMAT(4F10.5,4I5) AMIN,AMAX,ALEN,DTICA... the same meaning as in the case of the x axis, but for the amplitude axis (decadic loga- rithms of amplitudes). ALEN=0.: Indication of the end of plotting amplitude-distance curves for the considered elementary wave. ICOMP... ICOMP=0: Amplitudes of the vertical component. ICOMP=1: Amplitudes of the radial component. ICOMP=2: Amplitudes of the transverse component. MSOUR... MSOUR=0: Unit isotropic radiation pattern. MSOUR=1: Single force. MSOUR=2: Double couple. MSOUR=3: Explosive (implosive) source. Note: MSOUR.NE.0 can be used only for the source situated in an isotropic layer. 11) Specification of parameters of the radiation pattern. This card is included only when MSOUR.NE.0. IPAR(1),...,IPAR(4),PAR(1),...,PAR(6) FORMAT(4I5,6F10.5) MSOUR=1 - single force: IPAR(1)-IPAR(4)... free parameters. PAR(1)... force azimuth - the angle, in radians, made by the force vector and the positive x axis. PAR(2)... magnitude of the force. PAR(3)... force declination - the angle, in radians, made by the force vector and the horizontal plane. PAR(4)-PAR(6)... free parameters. MSOUR=2 - double couple: IPAR(1)-IPAR(4)... free parameters. PAR(1)... dip angle, in radians. PAR(2)... source moment. PAR(3)... strike angle, in radians. PAR(4)... rake angle, inradians. PAR(5)-PAR(6)... free parameters. Note: The formulae for the double couple follow the definition in Aki and Richards (1980). MSOUR=3 - explosive (implosive) source: IPAR(1)=1 (-1)... for explosive (implosive) source. IPAR(2)-IPAR(4)... free parameters. PAR(1)... magnitude of the source. PAR(2)-PAR(6)... free parameters. 12) This card is included only when NTICA.NE.0. NEXP FORMAT(26I3) NEXP... NEXP.NE.0: Number of points in the additional amplitude-distance curve which is to be plotted in the same frame as the computed amplitude-distance curve. NEXP=0: No additional amplitude-distance curve is to be plotted. 13) This card is included only when NTICA.NE.0 and NEXP.NE.0. (X(I),Y(I),I=1,NEXP) FORMAT(8F10.5) X(I),Y(I)... x coordinates and coresponding amplitudes of the additional amplitude-distance curve. Real amplitudes, not logarithms of amplitudes should be specified. Termination of computations: ---------------------------- Input data cards No.10,11,12 and 13 can be repeated for the given elementary wave several times (to perform plotting of amplitude- distance curves of the elementary wave corresponding to different components of the polarization vector, or for different types and orientations of point sources). The last card should be a blank card (to guarantee ALEN=0. After the blank card, a new system of cards 2-13 for a new elementary wave stored in the file LU1 may follow. As soon as ICONT=0 is found in LU1, plotting from the file LU1 terminates. The plotting from LU1 may be also terminated if we put NTICX=0 in the card No.3, even if LU1 still contains some other data. In both cases, the file is closed and the computations terminate.