73 |
CALL VZERO(YY,NPLAV) |
CALL VZERO(YY,NPLAV) |
74 |
VARFIT(1) = 0. |
VARFIT(1) = 0. |
75 |
VARFIT(2) = 0. |
VARFIT(2) = 0. |
76 |
NPFIT(1) = 0. |
NPFIT(1) = 0 |
77 |
NPFIT(2) = 0. |
NPFIT(2) = 0 |
78 |
VARCFIT(1) = 0. |
VARCFIT(1) = 0. |
79 |
VARCFIT(2) = 0. |
VARCFIT(2) = 0. |
80 |
NPCFIT(1) = 0. |
NPCFIT(1) = 0 |
81 |
NPCFIT(2) = 0. |
NPCFIT(2) = 0 |
82 |
ELEN = 0. |
ELEN = 0. |
83 |
SELEN = 0. |
SELEN = 0. |
84 |
NSTRIP = 0. |
NSTRIP = 0. |
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 |