/[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.2 by mocchiut, Fri Jul 20 08:24:54 2007 UTC revision 1.3 by mocchiut, Thu Apr 10 11:41:48 2008 UTC
# 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                ENDDO
73             ENDIF
74  C      C    
75           DO J=B+2,NCHA           IF ((B+2).LE.NCHA) THEN
76              PXY = DEXY(2,I,J)              DO J=B+2,NCHA
77              IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY                 PXY = DEXY(2,I,J)
78           ENDDO                 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.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.23