31 |
C |
C |
32 |
REAL LUNGH |
REAL LUNGH |
33 |
REAL RPIANO(2) |
REAL RPIANO(2) |
34 |
REAL BAR(2,NPLA) |
REAL BAR(2,NPLAV) |
35 |
REAL NQQ(2*NPLA) |
REAL NQQ(2*NPLAV) |
36 |
REAL RIL(NCHA/2) |
REAL RIL(NCHA/2) |
37 |
REAL REL(NCHA/2) |
REAL REL(NCHA/2) |
38 |
REAL TG(2) |
REAL TG(2) |
39 |
INTEGER IBAR(2,NPLA) |
INTEGER IBAR(2,NPLAV) |
40 |
INTEGER NTOT(2),LH |
INTEGER NTOT(2),LH |
41 |
COMMON / ANGOLO / BAR, IBAR |
COMMON / ANGOLO / BAR, IBAR |
42 |
SAVE / ANGOLO / |
SAVE / ANGOLO / |
68 |
RPIANO(M) = 0. |
RPIANO(M) = 0. |
69 |
C |
C |
70 |
DO J = 1,NPLA |
DO J = 1,NPLA |
71 |
|
C |
72 |
IF (MOD(J,2).EQ.0) THEN |
IF (MOD(J,2).EQ.0) THEN |
73 |
SHIFT = -0.5 |
IF (REVERSE.EQ.0) THEN |
74 |
|
SHIFT = -0.5 |
75 |
|
ELSE |
76 |
|
SHIFT = +0.5 |
77 |
|
ENDIF |
78 |
ELSE |
ELSE |
79 |
SHIFT = 0.5 |
IF (REVERSE.EQ.0) THEN |
80 |
|
SHIFT = +0.5 |
81 |
|
ELSE |
82 |
|
SHIFT = -0.5 |
83 |
|
ENDIF |
84 |
ENDIF |
ENDIF |
85 |
C |
C |
86 |
DO K = 1,NGIR |
DO K = 1,NGIR |
101 |
NGIR = NCL(M,J) ! number of clusters in layer J of the M View |
NGIR = NCL(M,J) ! number of clusters in layer J of the M View |
102 |
IF (NGIR.EQ.0) GOTO 5 |
IF (NGIR.EQ.0) GOTO 5 |
103 |
C IF (NGIR.NE.1) GOTO 5 |
C IF (NGIR.NE.1) GOTO 5 |
104 |
IPP = 0. |
IPP = 0 |
105 |
DO L = 1,NGIR |
DO L = 1,NGIR |
106 |
CALL LASTRISCIA(CLUS(M,J,L),IN) |
CALL LASTRISCIA(CLUS(M,J,L),IN) |
107 |
RIL(L) = FLOAT(ABS(NN - IN)) |
RIL(L) = FLOAT(ABS(NN - IN)) |
166 |
QQQQ = QQQ / FLOAT(NPIA) ! Truncated mean . |
QQQQ = QQQ / FLOAT(NPIA) ! Truncated mean . |
167 |
ENDIF |
ENDIF |
168 |
C |
C |
169 |
400 CONTINUE |
c 400 CONTINUE |
170 |
|
CONTINUE |
171 |
|
|
172 |
RETURN |
RETURN |
173 |
END |
END |