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 |