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

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

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

revision 1.1 by mocchiut, Fri May 19 13:15:50 2006 UTC revision 1.3 by mocchiut, Thu Apr 10 11:41:48 2008 UTC
# Line 9  C Line 9  C
9        INCLUDE 'INTEST.TXT'        INCLUDE 'INTEST.TXT'
10  C  C
11        INTEGER A,B, INPIA        INTEGER A,B, INPIA
12        INTEGER IBAR(2,NPLA)        INTEGER IBAR(2,NPLAV)
13        REAL BAR(2,NPLA)        REAL BAR(2,NPLAV)
14        REAL RQT1, RQT2        REAL RQT1, RQT2
15        COMMON / ANGOLO / BAR, IBAR        COMMON / ANGOLO / BAR, IBAR
16        SAVE / ANGOLO /        SAVE / ANGOLO /
# Line 29  C     Line 29  C    
29  C      C    
30  C     END EM  C     END EM
31  C      C    
32           IF (A.LE.2) A = 3           IF (A.LE.2) A = 2
33           IF (A.GE.(NCHA-1)) A = NCHA - 2           IF (A.GE.NCHA) A = NCHA - 1
34           DO J = A-1,A+1           DO J = A-1,A+1
35              IF (DEXY(1,I,J).GE.EMIN) RQT1 = RQT1 + DEXY(1,I,J)              IF (DEXY(1,I,J).GE.EMIN) RQT1 = RQT1 + DEXY(1,I,J)
36           ENDDO           ENDDO
37  C  C
38           DO J=1,A-2           IF ((A-2).GE.1) THEN
39              PXY = DEXY(1,I,J)              DO J=1,A-2
40              IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY                 PXY = DEXY(1,I,J)
41           ENDDO                 IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY
42  C              ENDDO
43           DO J=A+2,NCHA           ENDIF
44              PXY = DEXY(1,I,J)  C
45              IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY           IF ((A+2).LE.NCHA) THEN
46           ENDDO              DO J=A+2,NCHA
47                   PXY = DEXY(1,I,J)
48                   IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY
49                ENDDO
50             ENDIF
51  C  C
52   600     CONTINUE   600     CONTINUE
53  C  C
# Line 55  C Line 59  C
59  C  C
60  C     END EM  C     END EM
61  C  C
62           IF (B.LE.2) B = 3           IF (B.LE.2) B = 2
63           IF (B.GE.(NCHA-1)) B = NCHA - 2           IF (B.GE.NCHA) B = NCHA - 1
64           DO J = B-1,B+1           DO J = B-1,B+1
65              IF (DEXY(2,I,J).GE.EMIN) RQT1 = RQT1 + DEXY(2,I,J)              IF (DEXY(2,I,J).GE.EMIN) RQT1 = RQT1 + DEXY(2,I,J)
66           ENDDO           ENDDO
67  C      C    
68           DO J=1,B-2           IF ((B-2).GE.1) THEN
69              PXY = DEXY(2,I,J)              DO J=1,B-2
70              IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY                 PXY = DEXY(2,I,J)
71           ENDDO                 IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY
72  C                  ENDDO
73           DO J=B+2,NCHA           ENDIF
74              PXY = DEXY(2,I,J)  C    
75              IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY           IF ((B+2).LE.NCHA) THEN
76           ENDDO              DO J=B+2,NCHA
77                   PXY = DEXY(2,I,J)
78                   IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY
79                ENDDO
80             ENDIF
81  C  C
82   700     CONTINUE   700     CONTINUE
83        ENDDO        ENDDO
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.23