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 |
IF (REVERSE.EQ.0) THEN |
IF (REVERSE.EQ.0) THEN |
|
SHIFT = +0.5 |
|
|
ELSE |
|
72 |
SHIFT = -0.5 |
SHIFT = -0.5 |
73 |
|
ELSE |
74 |
|
SHIFT = +0.5 |
75 |
ENDIF |
ENDIF |
76 |
ELSE |
ELSE |
77 |
IF (REVERSE.EQ.0) THEN |
IF (REVERSE.EQ.0) THEN |
|
SHIFT = -0.5 |
|
|
ELSE |
|
78 |
SHIFT = +0.5 |
SHIFT = +0.5 |
79 |
|
ELSE |
80 |
|
SHIFT = -0.5 |
81 |
ENDIF |
ENDIF |
82 |
ENDIF |
ENDIF |
83 |
C |
C |
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 |