/[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.5 by mocchiut, Tue Aug 4 14:01:20 2009 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    c         print *,' I ',I,' A ',RQT1
53   600     CONTINUE   600     CONTINUE
54  C  C
55           B = IBAR(2,I)           B = IBAR(2,I)
# Line 55  C Line 60  C
60  C  C
61  C     END EM  C     END EM
62  C  C
63           IF (B.LE.2) B = 3           IF (B.LE.2) B = 2
64           IF (B.GE.(NCHA-1)) B = NCHA - 2           IF (B.GE.NCHA) B = NCHA - 1
65           DO J = B-1,B+1           DO J = B-1,B+1
66              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)
67    c            print *,' ==> J ',j,' dexy ',DEXY(2,I,J)
68           ENDDO           ENDDO
69  C      C    
70           DO J=1,B-2           IF ((B-2).GE.1) THEN
71              PXY = DEXY(2,I,J)              DO J=1,B-2
72              IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY                 PXY = DEXY(2,I,J)
73           ENDDO                 IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY
74  C                  ENDDO
75           DO J=B+2,NCHA           ENDIF
76              PXY = DEXY(2,I,J)  C    
77              IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY           IF ((B+2).LE.NCHA) THEN
78           ENDDO              DO J=B+2,NCHA
79                   PXY = DEXY(2,I,J)
80                   IF (PXY.GE.EMIN) RQT2 = RQT2 + PXY
81                ENDDO
82             ENDIF
83  C  C
84    c         print *,' I ',I,' B ',RQT1
85   700     CONTINUE   700     CONTINUE
86        ENDDO        ENDDO
87  C      C    
88  C  C
89   800  RETURN  c 800  RETURN
90          RETURN
91        END        END
92    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.23