| 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 |