| 202 |
C IF (M.EQ.2) NN = 1 |
C IF (M.EQ.2) NN = 1 |
| 203 |
IF (MOD(I,2).EQ.NN) THEN |
IF (MOD(I,2).EQ.NN) THEN |
| 204 |
IF (REVERSE.EQ.0) THEN |
IF (REVERSE.EQ.0) THEN |
|
SHIFT = +0.5 |
|
|
ELSE |
|
| 205 |
SHIFT = -0.5 |
SHIFT = -0.5 |
| 206 |
|
ELSE |
| 207 |
|
SHIFT = +0.5 |
| 208 |
ENDIF |
ENDIF |
| 209 |
ELSE |
ELSE |
| 210 |
IF (REVERSE.EQ.0) THEN |
IF (REVERSE.EQ.0) THEN |
|
SHIFT = -0.5 |
|
|
ELSE |
|
| 211 |
SHIFT = +0.5 |
SHIFT = +0.5 |
| 212 |
|
ELSE |
| 213 |
|
SHIFT = -0.5 |
| 214 |
ENDIF |
ENDIF |
| 215 |
ENDIF |
ENDIF |
| 216 |
C |
C |
| 247 |
C CALCULATE ELEN AND SELEN |
C CALCULATE ELEN AND SELEN |
| 248 |
C |
C |
| 249 |
c ELEN = 40.82 * QTOT * 1.06 * 1E-4 |
c ELEN = 40.82 * QTOT * 1.06 * 1E-4 |
| 250 |
ELEN = 260. * QTOT |
C ELEN = 260. * QTOT |
| 251 |
|
ELEN = QTOT / 260. |
| 252 |
IF (ELEN.GT.0.) THEN |
IF (ELEN.GT.0.) THEN |
| 253 |
SELEN = ELEN * (0.01183 + 0.121 / SQRT(ELEN)) |
SELEN = ELEN * (0.01183 + 0.121 / SQRT(ELEN)) |
| 254 |
ELSE |
ELSE |
| 257 |
C |
C |
| 258 |
9696 CONTINUE |
9696 CONTINUE |
| 259 |
C |
C |
| 260 |
45 continue |
c 45 continue |
| 261 |
|
|
| 262 |
50 continue |
c 50 continue |
| 263 |
|
|
| 264 |
return |
return |
| 265 |
END |
END |