*****************************************
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) XPRF,YPRF,0.,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.