C
C INCLUDE 'initc.inc' C ------------------------------------------------------------------ INTEGER MYI PARAMETER (MYI=29) INTEGER MSRFCA PARAMETER (MSRFCA=128) INTEGER ISB1I,ICB1I REAL YLI(6),YI(MYI),FSRFCA(MSRFCA) COMMON/INITC/ISB1I,ICB1I,YLI,YI,FSRFCA SAVE /INITC/ C ------------------------------------------------------------------ C ISB1I,ICB1I... Indices of a simple and a complex blocks in which C the initial point of the ray is situated, see C C.R.T.6.1. C YLI... Array containing the values of the quantities YL(1)-YL(6), C see C C.R.T.5.5.4, C describing the local properties of the model C at the initial point of the ray, see C C.R.T.6.1. C They must not be changed outside the subroutine INIT2. C C Description of YL. C C YI... Array containing the following quantities describing the C properties of the rays and of the travel-time field, see C C.R.T.6.1. C C Following quantities, describing the properties defined for C a single ray, are defined in subroutine INIT2 of file 'init.for'. C INIT2 C YI(1)...Initial travel time. C YI(2)...Initial imaginary part of the complex travel time. C YI(3)-YI(5)... Coordinates of the initial point of the ray. C YI(6)-YI(8)... Covariant components of the initial slowness C vector. C YI(9)-YI(11)... Covariant components of the first basis vector of C the ray-centred coordinate system at the initial point of C the ray (perpendicular to the slowness vector C YI(6)-YI(8)). C YI(12),YI(16) QR11,QR12 C YI(13),YI(17) QR21,QR22 C YI(14),YI(18) PR11,PR12 C YI(15),YI(19)... PR21,PR22 C Elements of the ray geometrical spreading matrix QR, and C of the matrix PR (see C.R.T., C eq.(5.13)) C at the initial point of the ray. C YI(20),YI(21)... Take-off parameters of the ray. C C Following quantities YI(22)-YI(29), describing the properties C of the discrete system of computed rays in the vicinity of the C computed ray, are defined in subroutine RPAR4 of file 'rpar.for'. C RPAR4 C YI(22)..Area of the element of the ray-parameter surface, C corresponding to the ray, see C.R.T., C eq.(6.1). C YI(23),YI(24),YI(25)... Components 11, 12, 22 of the symmetric C matrix inverse to the specific moment of the element of C the ray-parameter surface corresponding to the ray, see C eq.(6.2). C C Additional quantities related to the shooting algorithm: C YI(26),YI(27)... Normalized take-off parameters of the ray, both C taking the values between 0 and 1. C YI(28),YI(29)... For a successful ray, values of the C X1 and X2 C functions parametrizing the reference surface. C Otherwise zeros. C C The index of the last allocated numeric unit of array FSRFCA is C named MSRFCA. Dimension MSRFCA may be adjusted if necessary. C C Common block /INITC/ is included in external procedures INIT1 and C INIT2 of 'init.for', in OUTP of 'raycb.for', in 'rpar.for', C in 'writ.for', in 'scro.for', and may be included in any other C subroutine. C C Date: 2004, December 7 C Coded by Ludek Klimes C C======================================================================= C