C     INCLUDE 'rparc.inc'
C     ------------------------------------------------------------------
      REAL G1,G2,X1,X2,X1G1,X2G1,X1G2,X2G2,GG11,GG12,GG22
     *             G1,G2,X1,X2,X1G1,X2G1,X1G2,X2G2,GG11,GG12,GG22
      SAVE /RPARC/
C     ..................................................................
      PARAMETER (MRPARH=10000)
      SAVE /RPARH/
C     ------------------------------------------------------------------
C     LURPAR..Logical unit number of the external input device
C             containing the input data.
C     IRAY0...Difference IRAY0=IRAY-JRAY between index of a ray within
C             the elementary wave, and the index within one shooting
C             domain.
C     JRAY... Reference index of the ray within one shooting domain
C             (not within one elementary wave).
C     JTYPE...Reserved for future use.
C     G1,G2...Normalized shooting parameters.  Here, normalized means
C             that the shooting domain is (0,1)*(0,1) in the normalized
C             shooting parameters.
C             For one-parametric shooting:
C               G1,G2: multiplied by ANUM and BNUM, respectively.
C               G1: the first, shooting parameter.
C               G2: the second, fixed parameter.
C     X1,X2...Values of the X-functions at the point of intersection
C             with the reference surface (e.g., two selected coordinates
C             of the point of intersection).
C     X1G1,X2G1,X1G2,X2G2... Derivatives of the X-functions with respect
C             to the normalized shooting parameters.
C     GG11,GG12,GG22... Metric tensor in ray-parameter domain measuring
C             the distance between rays.
C     JPOINT..The number of crosssections of the computed ray with the
C             surface ISRFR.  Because, for the structural interfaces,
C             the positive or negative side of the surface is specified,
C             either zero or two crosssections correspond to the
C             reflection from the surface ISRFR.
C     KMAH... KMAH index (number of caustic points along the ray) at the
C             reference surface.  Zero if the ray does not reach the
C             reference surface.
C     INDMIN..Absolute value of the ray-history index minus the
C             location of the corresponding ray history in the array
C             IRPARH of the common block /RPARH/.
C     NRPARH..Index in the array IRPARH of the last element of the
C             history of the computed ray, see below.
C     IRPARH..Array containing different histories of rays.  Under the
C             ray history we understand here the block-surface code (see
C             the file 'code.for') describing the behaviour of
C             individual rays, complemented by the indication IEND
C             specifying the reason of the termination of the
C             computation of the ray (see the subroutine RAY2 of the
C             file 'ray.for').
C             IRPARH(0)... Number of different ray histories encountered
C               yet.
C             IRPARH(I),I=1,...,IRPARH(0)... Indices in the array IRPARH
C               of the last elements of individual ray histories.
C             IRPARH(I),I=IRPARH(0)+1,...,IRPARH(IRPARH(0))...
C               Individual encountered ray histories.
C             IRPARH(I),I=IRPARH(IRPARH(0)),...,NRPARH... History of the
C               computed ray.  After finishing computation of the ray,
C               it is compared with the recorded histories.  If it
C               appears to be a new history, it is recorded, otherwise
C               it is forgotten.
C             If the array IRPARH is too small to keep all different ray
C             histories corresponding to the computed elementary wave,
C             the oldest histories are forgotten.
C     Common block /RPARC/ is included in external procedures RPAR1,
C     RPAR2, RPAR32, and RPAR4 of 'rpar.for', in 'scro.for'.
C     Common block /RPARH/ is included in external procedures RPAR2,
C     RPAR31, RPAR32, IHIST, and SRPARH of 'rpar.for'.
C Date: 2003, April 3
C Coded by Ludek Klimes