Parent Directory
|
Revision Log
|
Patch
| revision 1.1 by mocchiut, Tue Jan 23 11:52:27 2007 UTC | revision 1.5 by mocchiut, Tue Aug 4 14:01:22 2009 UTC | |
|---|---|---|
| # | Line 31 C--------------------------------------- | Line 31 C--------------------------------------- |
| 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 / |
| # | Line 68 C | Line 68 C |
| 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 |
| # | Line 157 C | Line 166 C |
| 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 |
|
||||||||
| ViewVC Help | |
| Powered by ViewVC 1.1.23 |