--- DarthVader/CalorimeterLevel2/src/elio.for 2006/12/15 14:41:54 1.3 +++ DarthVader/CalorimeterLevel2/src/elio.for 2009/08/04 14:01:19 1.7 @@ -31,11 +31,11 @@ C REAL LUNGH REAL RPIANO(2) - REAL BAR(2,NPLA) - REAL NQQ(2*NPLA) + REAL BAR(2,NPLAV) + REAL NQQ(2*NPLAV) REAL RIL(NCHA/2) REAL TG(2) - INTEGER IBAR(2,NPLA) + INTEGER IBAR(2,NPLAV) INTEGER NTOT(2) COMMON / ANGOLO / BAR, IBAR SAVE / ANGOLO / @@ -68,9 +68,17 @@ C DO J = 1,NPLA IF (MOD(J,2).EQ.0) THEN - SHIFT = -0.5 + IF (REVERSE.EQ.0) THEN + SHIFT = -0.5 + ELSE + SHIFT = +0.5 + ENDIF ELSE - SHIFT = 0.5 + IF (REVERSE.EQ.0) THEN + SHIFT = +0.5 + ELSE + SHIFT = -0.5 + ENDIF ENDIF C DO K = 1,NGIR @@ -153,7 +161,8 @@ QQQQ = QQQ / FLOAT(NPIA) ! Truncated mean . ENDIF C - 400 CONTINUE +c 400 CONTINUE + CONTINUE RETURN END