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 TG(2) |
REAL TG(2) |
38 |
INTEGER IBAR(2,NPLA) |
INTEGER IBAR(2,NPLAV) |
39 |
INTEGER NTOT(2) |
INTEGER NTOT(2) |
40 |
COMMON / ANGOLO / BAR, IBAR |
COMMON / ANGOLO / BAR, IBAR |
41 |
SAVE / ANGOLO / |
SAVE / ANGOLO / |
68 |
C |
C |
69 |
DO J = 1,NPLA |
DO J = 1,NPLA |
70 |
IF (MOD(J,2).EQ.0) THEN |
IF (MOD(J,2).EQ.0) THEN |
71 |
SHIFT = -0.5 |
IF (REVERSE.EQ.0) THEN |
72 |
|
SHIFT = -0.5 |
73 |
|
ELSE |
74 |
|
SHIFT = +0.5 |
75 |
|
ENDIF |
76 |
ELSE |
ELSE |
77 |
SHIFT = 0.5 |
IF (REVERSE.EQ.0) THEN |
78 |
|
SHIFT = +0.5 |
79 |
|
ELSE |
80 |
|
SHIFT = -0.5 |
81 |
|
ENDIF |
82 |
ENDIF |
ENDIF |
83 |
C |
C |
84 |
DO K = 1,NGIR |
DO K = 1,NGIR |
99 |
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 |
100 |
IF (NGIR.EQ.0) GOTO 5 |
IF (NGIR.EQ.0) GOTO 5 |
101 |
C IF (NGIR.NE.1) GOTO 5 |
C IF (NGIR.NE.1) GOTO 5 |
102 |
IPP = 0. |
IPP = 0 |
103 |
DO L = 1,NGIR |
DO L = 1,NGIR |
104 |
CALL LASTRISCIA(CLUS(M,J,L),IN) |
CALL LASTRISCIA(CLUS(M,J,L),IN) |
105 |
RIL(L) = FLOAT(ABS(NN - IN)) |
RIL(L) = FLOAT(ABS(NN - IN)) |
161 |
QQQQ = QQQ / FLOAT(NPIA) ! Truncated mean . |
QQQQ = QQQ / FLOAT(NPIA) ! Truncated mean . |
162 |
ENDIF |
ENDIF |
163 |
C |
C |
164 |
400 CONTINUE |
c 400 CONTINUE |
165 |
|
CONTINUE |
166 |
|
|
167 |
RETURN |
RETURN |
168 |
END |
END |