GEANT321/gdraw/gdtr6.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:20:29  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.27  by  S.Giani
  11 *-- Author :
  12       SUBROUTINE GDTR6(JMAX)
  13 C.
  14 C.    ******************************************************************
  15 C.    *                                                                *
  16 C.    *       Build all optimized tree                                 *
  17 C.    *                                                                *
  18 C.    *       JMAX = total number of nodes (output)                    *
  19 C.    *                                                                *
  20 C.    *    ==>Called by : GDTR0                                        *
  21 C.    *       Author : P.Zanarini   *********                          *
  22 C.    *                                                                *
  23 C.    ******************************************************************
  24 C.
  25 #include "geant321/gcdraw.inc"
  26 #include "geant321/gcbank.inc"
  27 C.
  28 C.    ------------------------------------------------------------------
  29 C.
  30       DO 10 J=1,NUMNOD
  31          IQ(JSCA+J)=0
  32    10 CONTINUE
  33       JNEW=0
  34       DO 40 J=1,NUMNOD
  35          IF (IQ(JDUP+J).EQ.0) THEN
  36             JX=IQ(JXON+J)
  37             IF (JX.NE.0) THEN
  38                IQ(JDUP+JX)=0
  39                JB=IQ(JBRO+JX)
  40    20          IF (JB.EQ.0) GO TO 30
  41                IQ(JDUP+JB)=0
  42                JB=IQ(JBRO+JB)
  43                GO TO 20
  44    30          CONTINUE
  45             ENDIF
  46          ELSE
  47             JNEW=JNEW+1
  48             IQ(JSCA+J)=JNEW
  49          ENDIF
  50    40 CONTINUE
  51       JMAX=JNEW
  52       IQ(JNAM1+1)=IQ(JNAM+1)
  53       IQ(JMOT1+1)=0
  54       IQ(JBRO1+1)=0
  55       DO 50 J=1,NUMNOD
  56          IF (IQ(JSCA+J).NE.0) CALL GDTR5(J)
  57    50 CONTINUE
  58       END