/[PAMELA software]/DarthVader/CalorimeterLevel2/src/cluster.for
ViewVC logotype

Diff of /DarthVader/CalorimeterLevel2/src/cluster.for

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1.1.1 by mocchiut, Fri May 19 13:15:50 2006 UTC revision 1.5 by mocchiut, Thu Jan 16 15:29:12 2014 UTC
# Line 30  C--------------------------------------- Line 30  C---------------------------------------
30        SAVE / SHIFT /        SAVE / SHIFT /
31  C  C
32        I = 0        I = 0
33        DO J = 1,NPLA        DO J = 1,NPLAV
34          DO N = 1,2          DO N = 1,2
35            DO M = 1,NCHA            DO M = 1,NCHA
36              CLUS(N,J,M) = 0.              CLUS(N,J,M) = 0.
# Line 45  C Line 45  C
45              QQQ(M) = 0.              QQQ(M) = 0.
46              IF (DEXY(N,J,M).GT.EMIN) THEN              IF (DEXY(N,J,M).GT.EMIN) THEN
47                QQQ(M) = DEXY(N,J,M)                QQQ(M) = DEXY(N,J,M)
48    c              print *,' N ',N,' J ',J,' M ',M,' QQQ ',QQQ(M)
49              ENDIF              ENDIF
50            ENDDO            ENDDO
51  C  C
52            IGIRO = 0.            IGIRO = 0
53  C  C
54     50     CONTINUE     50     CONTINUE
55  C  C
# Line 74  C Line 75  C
75            IA = MSTR + 1            IA = MSTR + 1
76  C  C
77            NN = 0            NN = 0
78            IF (M.EQ.2) NN = 1  c          IF (N.EQ.2) NN = 1
79            IF (MOD(I,2).EQ.NN) THEN            IF (MOD(J,2).EQ.NN) THEN
80               SHIFT = +0.5               IF (REVERSE.EQ.0) THEN
81                    SHIFT = -0.5
82                 ELSE
83                    SHIFT = +0.5
84                 ENDIF
85            ELSE            ELSE
86               SHIFT = -0.5               IF (REVERSE.EQ.0) THEN
87                    SHIFT = +0.5
88                 ELSE
89                    SHIFT = -0.5
90                 ENDIF          
91            ENDIF            ENDIF
92            CALL BARIC(QQQ,IA,IB,BAR,ENE,IERR1)            CALL BARIC(QQQ,IA,IB,BAR,ENE,IERR1)
93            IGIRO = IGIRO + 1            IGIRO = IGIRO + 1
94            IGI = IGIRO + NCHA / 2            IGI = IGIRO + NCHA / 2
95    c          print *,'IGIRO ',IGIRO,' IGI ',IGI,' NCHA ',NCHA
96    c          PRINT *,' ENE ',ENE,' BAR ',BAR
97            CLUS(N,J,IGIRO) = BAR            CLUS(N,J,IGIRO) = BAR
98            CLUS(N,J,IGI) = ENE            CLUS(N,J,IGI) = ENE
99  C  C
Legend:
Removed from v.1.1.1.1  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.23