CalorimeterLevel2/src/baric.for

C
C------------------------------------------------------
c      SUBROUTINE BARIC(Q,INF,ISUP,BAR,EPLA4,INDEX)
      INTEGER FUNCTION BARIC(Q,INF,ISUP,BAR,EPLA4,INDEX)
C------------------------------------------------------
      INCLUDE 'INTEST.TXT'
      REAL Q(NCHA)
      REAL BAR
      REAL CM, EPLA4, ENN
      INTEGER INF,ISUP,INDEX, LL
      COMMON / SHIFT / SHIFT
      SAVE / SHIFT /
      DATA CM / 0. /
c      DATA EPLA4 / 0. /
      DATA ENN / 0. /

      BARIC = 0

      LL = 0
      CM = 0.
      EPLA4 = 0.
      if ( index*Q(1) .eq. -2000.323311 ) then
c         index = 10
         print *,' BEWARE!'
         print *,' this is the GHOST of the compilation '
         print *,' you have awakened me too soon '
         print *,' I will launch an anatema against you!'
      endif
      DO LL = INF,ISUP
        ENN = Q(LL)
        IF (ENN.GT.EMIN) THEN
C
          CALL MILLIM(LL,RMM)
C
          CM = CM + RMM * ENN
          EPLA4 = EPLA4 + ENN
        ENDIF
      ENDDO
      IF (EPLA4.EQ.0) GO TO 300
      BAR = CM / EPLA4
      GOTO 400
C
  300 INDEX = 0
      GOTO 500
C
  400 INDEX = 1
C
  500 CONTINUE
      RETURN
      END
C
1	mocchiut	1.1	C
2			C------------------------------------------------------
3	mocchiut	1.2	c SUBROUTINE BARIC(Q,INF,ISUP,BAR,EPLA4,INDEX)
4			INTEGER FUNCTION BARIC(Q,INF,ISUP,BAR,EPLA4,INDEX)
5	mocchiut	1.1	C------------------------------------------------------
6			INCLUDE 'INTEST.TXT'
7			REAL Q(NCHA)
8			REAL BAR
9	mocchiut	1.2	REAL CM, EPLA4, ENN
10			INTEGER INF,ISUP,INDEX, LL
11	mocchiut	1.1	COMMON / SHIFT / SHIFT
12			SAVE / SHIFT /
13	mocchiut	1.2	DATA CM / 0. /
14			c DATA EPLA4 / 0. /
15			DATA ENN / 0. /
16
17			BARIC = 0
18	mocchiut	1.1
19	mocchiut	1.2	LL = 0
20	mocchiut	1.1	CM = 0.
21			EPLA4 = 0.
22	mocchiut	1.2	if ( index*Q(1) .eq. -2000.323311 ) then
23			c index = 10
24			print *,' BEWARE!'
25			print *,' this is the GHOST of the compilation '
26			print *,' you have awakened me too soon '
27			print *,' I will launch an anatema against you!'
28			endif
29	mocchiut	1.1	DO LL = INF,ISUP
30			ENN = Q(LL)
31			IF (ENN.GT.EMIN) THEN
32			C
33			CALL MILLIM(LL,RMM)
34	mocchiut	1.2	C
35	mocchiut	1.1	CM = CM + RMM * ENN
36			EPLA4 = EPLA4 + ENN
37			ENDIF
38			ENDDO
39			IF (EPLA4.EQ.0) GO TO 300
40			BAR = CM / EPLA4
41			GOTO 400
42			C
43			300 INDEX = 0
44			GOTO 500
45			C
46			400 INDEX = 1
47			C
48			500 CONTINUE
49			RETURN
50			END
51			C