[u/mrichter/AliRoot.git] / GEANT321 / erdecks / evolio.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1996/03/06 15:37:35  mclareni
* Add geane321 source directories
*
*
#include "geant321/pilot.h"
*CMZ :  3.21/02 29/03/94  15.41.49  by  S.Giani
*-- Author :
      SUBROUTINE EVOLIO (NVLIN,NAMIN,NUMIN,NVLOUT,NAMOUT,NUMOUT)
*
C.    ******************************************************************
C.    *                                                                *
C.    *    Called with INWVOL=2 this routine save the volume's tree    *
C.    *                                                                *
C.    *    Called with INWVOL=1 this routine returns                   *
C.    *             - the list of volumes (and number) left            *
C.    *             - the list of volumes (and number) entered         *
C.    *                                                                *
C.    *    ==>Called by : ERTRGO                                       *
C.    *       Author    M.Maire  *********                             *
C.    *                                                                *
C.    ******************************************************************
C.
#include "geant321/gcvolu.inc"
#include "geant321/gctrak.inc"
*
      PARAMETER (NVMX=15)
      COMMON/OLDTRE/NLVOLD,NAMOLD(NVMX),NUMOLD(NVMX),IOFLAG(NVMX),KONLY
*
      DIMENSION NAMIN(NVMX),NUMIN(NVMX),NAMOUT(NVMX),NUMOUT(NVMX)
*
*
      NVLOUT = 0
      NVLIN  = 0
      IF(INWVOL.EQ.0) RETURN
      IF(NLEVEL.EQ.0) RETURN
*
* *** Save the tree when leaving a volume
*
      IF(INWVOL.EQ.2) THEN
         NLVOLD = NLEVEL
         KONLY  = 1
         DO 1 ILV =1,NLEVEL
           NAMOLD(ILV) = NAMES (ILV)
           NUMOLD(ILV) = NUMBER(ILV)
           IOFLAG(ILV) = 2
           IF(GONLY(ILV).NE.1.) KONLY = 0
    1    CONTINUE
         RETURN
      ENDIF
*
* *** Entering in a new volume : compare the trees
*
      IF(INWVOL.NE.1) RETURN
*
* *** Simplest case : old tree is 'ONLY'
*
      IF(KONLY.EQ.1) THEN
        NLMIN = MIN(NLEVEL,NLVOLD)
        NL    = 0
   10   NL    = NL + 1
        IF(NL.GT.NLMIN) GO TO 12
        IF((NAMES (NL).EQ.NAMOLD(NL)).AND.
     *     (NUMBER(NL).EQ.NUMOLD(NL))) GO TO 10
*
   12   CONTINUE
        IF(NL.LE.NLVOLD) THEN
          DO 14 IL = NL,NLVOLD
            NVLOUT = NVLOUT + 1
            NAMOUT(NVLOUT) = NAMOLD(IL)
            NUMOUT(NVLOUT) = NUMOLD(IL)
   14     CONTINUE
        ENDIF
*
        IF(NL.LE.NLEVEL) THEN
          DO 16 IL = NL,NLEVEL
            NVLIN  = NVLIN  + 1
            NAMIN (NVLIN ) = NAMES (IL)
            NUMIN (NVLIN ) = NUMBER(IL)
   16     CONTINUE
        ENDIF
*
        RETURN
      ENDIF
*
* *** Special case for 'NOT ONLY' volumes
*
      IF(KONLY.EQ.0) THEN
         DO 24 IL1 = 1,NLEVEL
           DO 22 IL2 = 1,NLVOLD
              IF((NAMES (IL1).EQ.NAMOLD(IL2)).AND.
     *           (NUMBER(IL1).EQ.NUMOLD(IL2))) THEN
                   IOFLAG(IL2) = 0
                   GO TO 24
              ENDIF
   22      CONTINUE
           NVLIN  = NVLIN  + 1
           NAMIN (NVLIN ) = NAMES (IL1)
           NUMIN (NVLIN ) = NUMBER(IL1)
   24    CONTINUE
*
         DO 26 IL2 = 1,NLVOLD
           IF(IOFLAG(IL2).EQ.0) GO TO 26
           NVLOUT = NVLOUT + 1
           NAMOUT(NVLOUT) = NAMOLD(IL2)
           NUMOUT(NVLOUT) = NUMOLD(IL2)
   26    CONTINUE
*
      ENDIF
*
      END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.1 1996/03/06 15:37:35 mclareni
	6	* Add geane321 source directories
	7	*
	8	*
	9	#include "geant321/pilot.h"
	10	*CMZ : 3.21/02 29/03/94 15.41.49 by S.Giani
	11	*-- Author :
	12	SUBROUTINE EVOLIO (NVLIN,NAMIN,NUMIN,NVLOUT,NAMOUT,NUMOUT)
	13	*
	14	C. ******************************************************************
	15	C. * *
	16	C. * Called with INWVOL=2 this routine save the volume's tree *
	17	C. * *
	18	C. * Called with INWVOL=1 this routine returns *
	19	C. * - the list of volumes (and number) left *
	20	C. * - the list of volumes (and number) entered *
	21	C. * *
	22	C. * ==>Called by : ERTRGO *
	23	C. * Author M.Maire ********* *
	24	C. * *
	25	C. ******************************************************************
	26	C.
	27	#include "geant321/gcvolu.inc"
	28	#include "geant321/gctrak.inc"
	29	*
	30	PARAMETER (NVMX=15)
	31	COMMON/OLDTRE/NLVOLD,NAMOLD(NVMX),NUMOLD(NVMX),IOFLAG(NVMX),KONLY
	32	*
	33	DIMENSION NAMIN(NVMX),NUMIN(NVMX),NAMOUT(NVMX),NUMOUT(NVMX)
	34	*
	35	*
	36	NVLOUT = 0
	37	NVLIN = 0
	38	IF(INWVOL.EQ.0) RETURN
	39	IF(NLEVEL.EQ.0) RETURN
	40	*
	41	* *** Save the tree when leaving a volume
	42	*
	43	IF(INWVOL.EQ.2) THEN
	44	NLVOLD = NLEVEL
	45	KONLY = 1
	46	DO 1 ILV =1,NLEVEL
	47	NAMOLD(ILV) = NAMES (ILV)
	48	NUMOLD(ILV) = NUMBER(ILV)
	49	IOFLAG(ILV) = 2
	50	IF(GONLY(ILV).NE.1.) KONLY = 0
	51	1 CONTINUE
	52	RETURN
	53	ENDIF
	54	*
	55	* *** Entering in a new volume : compare the trees
	56	*
	57	IF(INWVOL.NE.1) RETURN
	58	*
	59	* *** Simplest case : old tree is 'ONLY'
	60	*
	61	IF(KONLY.EQ.1) THEN
	62	NLMIN = MIN(NLEVEL,NLVOLD)
	63	NL = 0
	64	10 NL = NL + 1
65	IF(NL.GT.NLMIN) GO TO 12
66	IF((NAMES (NL).EQ.NAMOLD(NL)).AND.
67	* (NUMBER(NL).EQ.NUMOLD(NL))) GO TO 10
68	*
69	12 CONTINUE
70	IF(NL.LE.NLVOLD) THEN
71	DO 14 IL = NL,NLVOLD
72	NVLOUT = NVLOUT + 1
73	NAMOUT(NVLOUT) = NAMOLD(IL)
74	NUMOUT(NVLOUT) = NUMOLD(IL)
75	14 CONTINUE
76	ENDIF
77	*
78	IF(NL.LE.NLEVEL) THEN
79	DO 16 IL = NL,NLEVEL
80	NVLIN = NVLIN + 1
81	NAMIN (NVLIN ) = NAMES (IL)
82	NUMIN (NVLIN ) = NUMBER(IL)
83	16 CONTINUE
84	ENDIF
85	*
86	RETURN
87	ENDIF
88	*
89	* *** Special case for 'NOT ONLY' volumes
90	*
91	IF(KONLY.EQ.0) THEN
92	DO 24 IL1 = 1,NLEVEL
93	DO 22 IL2 = 1,NLVOLD
94	IF((NAMES (IL1).EQ.NAMOLD(IL2)).AND.
95	* (NUMBER(IL1).EQ.NUMOLD(IL2))) THEN
96	IOFLAG(IL2) = 0
97	GO TO 24
98	ENDIF
99	22 CONTINUE
100	NVLIN = NVLIN + 1
101	NAMIN (NVLIN ) = NAMES (IL1)
102	NUMIN (NVLIN ) = NUMBER(IL1)
103	24 CONTINUE
104	*
105	DO 26 IL2 = 1,NLVOLD
106	IF(IOFLAG(IL2).EQ.0) GO TO 26
107	NVLOUT = NVLOUT + 1
108	NAMOUT(NVLOUT) = NAMOLD(IL2)
109	NUMOUT(NVLOUT) = NUMOLD(IL2)
110	26 CONTINUE
111	*
112	ENDIF
113	*
114	END