C C INCLUDE 'rp2d.inc'
C ------------------------------------------------------------------
INTEGER NEWMAX
PARAMETER (NEWMAX=100)
INTEGER ITER,NEW,INDOLD,IND(NEWMAX),IREOLD,IRE(NEWMAX)
INTEGER IRAOLD,IRA(NEWMAX)
REAL AOLD,A(NEWMAX),DAOLD,DA(NEWMAX)
REAL XOLD,X(NEWMAX),XAOLD,XA(NEWMAX)
COMMON/RP2DC/ITER,NEW,INDOLD,IND,IREOLD,IRE,IRAOLD,IRA,
* AOLD,A,DAOLD,DA,XOLD,X,XAOLD,XA
SAVE /RP2DC/
C ------------------------------------------------------------------
C
C ITER... Number of iterations when searching for the last
C boundary-value ray.
C NEW... Number of computed rays in the current cut (i.e.
C corresponding to the last value of B) to the right (i.e.
C for greater values of A) from the interval under study.
C The interval under study is A = AOLD to A(NEW). NEW=1
C when beginning to study a new basic interval. NEW=NEW+1
C if the interval under study is divided by a new ray.
C AOLD=A(NEW) and then NEW=NEW-1 if we change to the next
C interval.
C INDOLD..History of the left-end ray of the interval under study.
C IND... Ray histories of the rays to the right from the interval
C under study.
C IREOLD,IRE... Indices of the receivers for two-point rays,
C otherwise zeros.
C IRAOLD,IRA... Indices of the rays.
C AOLD... Value of a corresponding to the left-end ray of the
C interval under study.
C A... Values of a corresponding to the computed rays situated to
C the right of the interval under study.
C DAOLD...Element of a corresponding to the ray A=AOLD.
C DA... DA(I) is the element of a corresponding to the ray A=A(I).
C XOLD... Value of X-function corresponding to the ray A=AOLD.
C Defined only for successful rays, i.e. ray crossing the
C surface ISRFR at least IABS(IPOINT) times.
C X... X(I) is the value of X-function corresponding to the ray
C A=A(I). Defined only for successful rays.
C XAOLD...Derivative of X-function with respect to a corresponding
C to the ray A=AOLD. Defined only for successful rays.
C XA... XA(I) is the derivative of X-function with respect to A
C corresponding to the ray A=A(I). Defined only for
C successful rays.
C
C Common block /RP2DC/ is included in external procedures RP2D,
C RP2DM, and RP2DG of 'rp2d.for'.
C
C Date: 2001, October 29
C Coded by Ludek Klimes
C
C=======================================================================
C