Bugfix in AliPoints2Memory
[u/mrichter/AliRoot.git] / GEANT321 / ggeom / ggnlev.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1995/10/24 10:20:49  cernlib
6 * Geant
7 *
8 *
9 #include "geant321/pilot.h"
10 *CMZ :  3.21/02 29/03/94  15.41.28  by  S.Giani
11 *-- Author :
12       SUBROUTINE GGNLEV
13 C.
14 C.    ******************************************************************
15 C.    *                                                                *
16 C.    *    SUBR. GGNLEV                                                *
17 C.    *                                                                *
18 C.    *   Scans the JVOLUM structure to find the maximum Number of     *
19 C.    *   LEVels and maximum number of volumes in each level.top       *
20 C.    *   Stores corresponding information in bank next to JVOLUM      *
21 C.    *   mother bank.                                                 *
22 C.    *                                                                *
23 C.    *   Called by : GGCLOS                                           *
24 C.    *   Author    : S.Banerjee                                       *
25 C.    *                                                                *
26 C.    ******************************************************************
27 C.
28 #include "geant321/gcbank.inc"
29       PARAMETER (NLVMAX=15)
30       INTEGER ILINK(NLVMAX), LJVOM(NLVMAX), NLINK(NLVMAX)
31       INTEGER NDVLP(NLVMAX), MAXVO(NLVMAX)
32 C.
33 C.    ------------------------------------------------------------------
34 *
35       CALL VZERO (MAXVO, NLVMAX)
36       NLEV = 1
37       NLVT = NLEV
38       LJVOM(NLEV) = LQ(JVOLUM-1)
39       ILINK(NLEV) = 1
40       MAXVO(NLEV) = 1
41 *
42 * *** Loop over volumes going down the tree
43 *
44    10 IF (NLEV.GT.NLVT) NLVT = NLEV
45 *
46 *  ** Check if current volume is a source of local development
47 *
48       IF (LQ(LJVOM(NLEV)).GT.0) THEN
49          NDVLP(NLEV) = NLEV
50       ELSE IF (NLEV.GT.1) THEN
51          NDVLP(NLEV) = NDVLP(NLEV-1)
52       ELSE
53          NDVLP(NLEV) = NLVMAX
54       ENDIF
55       NIN = Q(LJVOM(NLEV)+3)
56       IF (NIN.GT.0) THEN
57 *
58 *   *    Current volume has contents defined by position
59 *
60          NLINK(NLEV) = NIN
61          IF (NIN.GT.MAXVO(NLEV+1)) MAXVO(NLEV+1) = NIN
62          IF (NLINK(NLEV).EQ.0) GO TO 30
63          IN   = NLINK(NLEV)
64          NLINK(NLEV) = IN - 1
65          JIN  = LQ(LJVOM(NLEV)-IN)
66          IVOT = Q(JIN+2)
67          NLEV = NLEV + 1
68          LJVOM(NLEV) = LQ(JVOLUM-IVOT)
69          ILINK(NLEV) = IN
70          GO TO 10
71 *
72       ELSE IF (NIN.LT.0) THEN
73 *
74 *   *    Current volume is divided
75 *
76          JDIV = LQ(LJVOM(NLEV)-1)
77          IVOT = Q(JDIV+2)
78          IF (NLEV.GT.NDVLP(NLEV)) THEN
79             JPAR = LQ(LJVOM(NDVLP(NLEV)))
80             DO 15 ILEV = NDVLP(NLEV), NLEV-1
81                IF (IQ(JPAR+1).EQ.0) THEN
82                   JPAR = LQ(JPAR-ILINK(ILEV+1))
83                   IF (JPAR.EQ.0) GO TO 20
84                ELSE
85                   JPAR = LQ(JPAR-1)
86                ENDIF
87                IF (ILEV.EQ.NLEV-1) THEN
88                   NDIV = IQ(JPAR+1)
89                   GO TO 25
90                ENDIF
91    15       CONTINUE
92          ENDIF
93    20    NDIV = Q(JDIV+3)
94          IF (NDIV.LT.0) NDIV = -NDIV
95    25    NLINK(NLEV) = 0
96          IF (NDIV.GT.MAXVO(NLEV+1)) MAXVO(NLEV+1) = NDIV
97          NLEV = NLEV + 1
98          LJVOM(NLEV) = LQ(JVOLUM-IVOT)
99          ILINK(NLEV) = 1
100          GO TO 10
101       ENDIF
102 *
103 *  ** No more content in this volume
104 *   * Go up the tree
105 *
106    30 NLEV = NLEV - 1
107       IF (NLEV.GT.0) THEN
108          IF (NLINK(NLEV).EQ.0) GO TO 30
109          IN   = NLINK(NLEV)
110          NLINK(NLEV) = IN - 1
111          JIN  = LQ(LJVOM(NLEV)-IN)
112          IVOT = Q(JIN+2)
113          NLEV = NLEV + 1
114          LJVOM(NLEV) = LQ(JVOLUM-IVOT)
115          ILINK(NLEV) = IN
116          GO TO 10
117       ENDIF
118 *
119 * *** Scan is complete; now store the information in next-of-same-type
120 *  ** link of JVOLUM
121 *
122       JVOLX = LQ(JVOLUM)
123       IF (JVOLX.GT.0) CALL MZDROP (IXCONS, JVOLX, ' ')
124       CALL MZBOOK (IXCONS, JVOLX, JVOLUM, 0, 'NLEV', 0, 0, NLVT, 2, 0)
125       CALL UCOPY (MAXVO, IQ(JVOLX+1), NLVT)
126 #if defined(CERNLIB_DEBUGG)
127       IF (IDEBUG.NE.0) THEN
128          WRITE (CHMAIL, 2001) NLVT
129          CALL GMAIL (1, 0)
130          IF (NLVT.GT.0) THEN
131             WRITE (CHMAIL, 2002) (MAXVO(I), I = 1, NLVT)
132             CALL GMAIL (0, 0)
133          ENDIF
134       ENDIF
135  2001 FORMAT (' GGNLEV : Maximum number of levels = ',I4,
136      +        ' -- maximum number of volumes in each level')
137  2002 FORMAT (5X,20I6)
138 #endif
139 *                                                             END GGNLEV
140       END