Bugfix in AliPoints2Memory
[u/mrichter/AliRoot.git] / GEANT321 / gphys / gshlin.F
CommitLineData
fe4da5cc 1*
2* $Id$
3*
4* $Log$
5* Revision 1.1.1.1 1995/10/24 10:21:33 cernlib
6* Geant
7*
8*
9#include "geant321/pilot.h"
10*CMZ : 3.21/02 29/03/94 15.41.23 by S.Giani
11*-- Author :
12 SUBROUTINE GSHLIN
13C.
14C. ******************************************************************
15C. * *
16C * Initialize atomic shells' potentials *
17C. * *
18C. * Potentials according to Rev.Mod.Phys. vol. 39 nr. 1 125 (1967) *
19C. * Xe upadetd according to Grishin et al. NIM A309, 476 (1991) *
20C. * *
21C. * ==>Called by : GPHINI *
22C. * Author J. Chwastowski *
23C. * *
24C. ******************************************************************
25C.
26C Shells
27c See alse Handbook of Phys. & Chemistry
28c
29 IMPLICIT NONE
30#include "geant321/gcpmxz.inc"
31#include "geant321/gcshpt.inc"
32 INTEGER NS,N1,NP,I
33 REAL TMP
34 DIMENSION NS(MAXSHL),N1(MAXELZ),NP(MAXELZ)
35 DIMENSION TMP(LENGTH)
36 DATA NS / 1,3*5,3*17,2*21,35,2*33,
37 + 2*41,2*58,3*49,2*73,3*82 /
38 DATA N1( 1),NP( 1) / 1, 1 /
39 DATA TMP( 1) / 13.59811 /
40 DATA N1( 2),NP( 2) / 2, 1 /
41 DATA TMP( 2) / 24.58678 /
42 DATA N1( 3),NP( 3) / 3, 1 /
43 DATA TMP( 3) / 54.75000 /
44 DATA N1( 4),NP( 4) / 4, 1 /
45 DATA TMP( 4) / 111.00000 /
46 DATA N1( 5),NP( 5) / 5, 4 /
47 DATA (TMP(I),I= 5, 8) /
48 + 188.00, 12.58, 4.70, 4.70 /
49 DATA N1( 6),NP( 6) / 9, 4 /
50 DATA (TMP(I),I= 9, 12) /
51 + 283.80, 17.56, 6.40, 6.40 /
52 DATA N1( 7),NP( 7) / 13, 4 /
53 DATA (TMP(I),I= 13, 16) /
54 + 401.60, 23.10, 9.20, 9.20 /
55 DATA N1( 8),NP( 8) / 17, 4 /
56 DATA (TMP(I),I= 17, 20) /
57 + 532.00, 23.70, 7.10, 7.10 /
58 DATA N1( 9),NP( 9) / 21, 4 /
59 DATA (TMP(I),I= 21, 24) /
60 + 685.40, 31.00, 8.60, 8.60 /
61 DATA N1( 10),NP( 10) / 25, 4 /
62 DATA (TMP(I),I= 25, 28) /
63 + 866.90, 45.00, 18.30, 18.30 /
64 DATA N1( 11),NP( 11) / 29, 4 /
65 DATA (TMP(I),I= 29, 32) /
66 + 1072.10, 63.30, 31.10, 31.10 /
67 DATA N1( 12),NP( 12) / 33, 4 /
68 DATA (TMP(I),I= 33, 36) /
69 + 1305.00, 89.40, 51.40, 51.40 /
70 DATA N1( 13),NP( 13) / 37, 4 /
71 DATA (TMP(I),I= 37, 40) /
72 + 1559.60, 117.70, 73.10, 73.10 /
73 DATA N1( 14),NP( 14) / 41, 4 /
74 DATA (TMP(I),I= 41, 44) /
75 + 1838.90, 148.70, 99.20, 99.20 /
76 DATA N1( 15),NP( 15) / 45, 4 /
77 DATA (TMP(I),I= 45, 48) /
78 + 2145.50, 189.30, 132.20, 132.20 /
79 DATA N1( 16),NP( 16) / 49, 4 /
80 DATA (TMP(I),I= 49, 52) /
81 + 2472.00, 229.20, 164.80, 164.80 /
82 DATA N1( 17),NP( 17) / 53, 7 /
83 DATA (TMP(I),I= 53, 59) /
84 + 2822.40, 270.20, 201.60, 200.00
85 + , 17.50, 6.80, 6.80 /
86 DATA N1( 18),NP( 18) / 60, 7 /
87 DATA (TMP(I),I= 60, 66) /
88 + 3202.90, 320.00, 247.30, 245.20
89 + , 25.30, 12.40, 12.40 /
90 DATA N1( 19),NP( 19) / 67, 7 /
91 DATA (TMP(I),I= 67, 73) /
92 + 3607.40, 377.10, 296.30, 293.60
93 + , 33.90, 17.80, 17.80 /
94 DATA N1( 20),NP( 20) / 74, 7 /
95 DATA (TMP(I),I= 74, 80) /
96 + 4038.10, 437.80, 350.00, 346.40
97 + , 43.70, 25.40, 25.40 /
98 DATA N1( 21),NP( 21) / 81, 9 /
99 DATA (TMP(I),I= 81, 89) /
100 + 4492.80, 500.40, 406.70, 402.20
101 + , 53.80, 32.30, 32.30, 6.60, 6.60 /
102 DATA N1( 22),NP( 22) / 90, 9 /
103 DATA (TMP(I),I= 90, 98) /
104 + 4966.40, 563.70, 461.50, 455.50
105 + , 60.30, 34.60, 34.60, 3.70, 3.70 /
106 DATA N1( 23),NP( 23) / 99, 9 /
107 DATA (TMP(I),I= 99, 107) /
108 + 5465.10, 628.20, 520.50, 512.90
109 + , 66.50, 37.80, 37.80, 2.20, 2.20 /
110 DATA N1( 24),NP( 24) / 108, 9 /
111 DATA (TMP(I),I= 108, 116) /
112 + 5989.20, 694.60, 583.70, 574.50
113 + , 74.10, 42.50, 42.50, 2.30, 2.30 /
114 DATA N1( 25),NP( 25) / 117, 9 /
115 DATA (TMP(I),I= 117, 125) /
116 + 6539.00, 769.00, 651.40, 640.30
117 + , 83.90, 48.60, 48.60, 3.30, 3.30 /
118 DATA N1( 26),NP( 26) / 126, 9 /
119 DATA (TMP(I),I= 126, 134) /
120 + 7112.00, 846.10, 721.10, 708.10
121 + , 92.90, 54.00, 54.00, 3.60, 3.60 /
122 DATA N1( 27),NP( 27) / 135, 9 /
123 DATA (TMP(I),I= 135, 143) /
124 + 7708.90, 925.60, 793.60, 778.60
125 + , 100.70, 59.50, 59.50, 2.90, 2.90 /
126 DATA N1( 28),NP( 28) / 144, 9 /
127 DATA (TMP(I),I= 144, 152) /
128 + 8332.80, 1008.10, 871.90, 854.70
129 + , 111.80, 68.10, 68.10, 3.60, 3.60 /
130 DATA N1( 29),NP( 29) / 153, 9 /
131 DATA (TMP(I),I= 153, 161) /
132 + 8978.90, 1096.60, 951.00, 931.10
133 + , 119.80, 73.60, 73.60, 1.60, 1.60 /
134 DATA N1( 30),NP( 30) / 162, 9 /
135 DATA (TMP(I),I= 162, 170) /
136 + 9658.60, 1193.60, 1042.80, 1019.70
137 + , 135.90, 86.60, 86.60, 8.10, 8.10 /
138 DATA N1( 31),NP( 31) / 171, 9 /
139 DATA (TMP(I),I= 171, 179) /
140 + 10367.10, 1297.70, 1142.30, 1115.40
141 + , 158.10, 106.80, 102.90, 17.40, 17.40 /
142 DATA N1( 32),NP( 32) / 180, 9 /
143 DATA (TMP(I),I= 180, 188) /
144 + 11103.10, 1414.30, 1247.80, 1216.70
145 + , 180.00, 127.90, 120.80, 28.70, 28.70 /
146 DATA N1( 33),NP( 33) / 189, 12 /
147 DATA (TMP(I),I= 189, 200) /
148 + 11866.70, 1526.50, 1358.60, 1323.10
149 + , 203.50, 146.40, 140.50, 41.20
150 + , 41.20, 0.00, 2.50, 2.50 /
151 DATA N1( 34),NP( 34) / 201, 12 /
152 DATA (TMP(I),I= 201, 212) /
153 + 12657.80, 1653.90, 1476.20, 1435.80
154 + , 231.50, 168.20, 161.90, 56.70
155 + , 56.70, 0.00, 5.60, 5.60 /
156 DATA N1( 35),NP( 35) / 213, 12 /
157 DATA (TMP(I),I= 213, 224) /
158 + 13473.70, 1782.00, 1596.00, 1549.90
159 + , 256.50, 189.30, 181.50, 70.10
160 + , 69.00, 27.30, 5.20, 4.60 /
161 DATA N1( 36),NP( 36) / 225, 12 /
162 DATA (TMP(I),I= 225, 236) /
163 + 14325.60, 1921.00, 1727.20, 1674.90
164 + , 292.80, 222.70, 213.80, 95.00
165 + , 93.80, 24.00, 10.60, 10.60 /
166 DATA N1( 37),NP( 37) / 237, 12 /
167 DATA (TMP(I),I= 237, 248) /
168 + 15199.70, 2065.10, 1863.90, 1804.40
169 + , 322.10, 247.40, 238.50, 111.80
170 + , 110.30, 29.30, 14.80, 14.00 /
171 DATA N1( 38),NP( 38) / 249, 12 /
172 DATA (TMP(I),I= 249, 260) /
173 + 16104.60, 2216.30, 2006.80, 1939.60
174 + , 357.50, 279.80, 269.10, 135.00
175 + , 133.10, 37.70, 19.90, 19.90 /
176 DATA N1( 39),NP( 39) / 261, 12 /
177 DATA (TMP(I),I= 261, 272) /
178 + 17038.40, 2372.50, 2155.50, 2080.00
179 + , 393.60, 312.40, 300.30, 159.60
180 + , 157.40, 45.40, 25.60, 25.60 /
181 DATA N1( 40),NP( 40) / 273, 12 /
182 DATA (TMP(I),I= 273, 284) /
183 + 17997.60, 2531.60, 2306.70, 2222.30
184 + , 430.30, 344.20, 350.50, 182.40
185 + , 180.00, 51.30, 28.70, 28.70 /
186 DATA N1( 41),NP( 41) / 285, 14 /
187 DATA (TMP(I),I= 285, 298) /
188 + 18985.60, 2697.70, 2464.70, 2370.50
189 + , 468.40, 378.40, 383.00, 207.40
190 + , 204.60, 58.10, 33.90, 33.90
191 + , 3.20, 3.20 /
192 DATA N1( 42),NP( 42) / 299, 14 /
193 DATA (TMP(I),I= 299, 312) /
194 + 19999.50, 2865.50, 2625.10, 2520.20
195 + , 504.60, 409.70, 392.30, 230.30
196 + , 227.00, 61.80, 34.80, 34.80
197 + , 1.80, 1.80 /
198 DATA N1( 43),NP( 43) / 313, 12 /
199 DATA (TMP(I),I= 313, 324) /
200 + 21044.00, 3042.50, 2793.20, 2676.90
201 + , 0.00, 444.90, 425.00, 256.40
202 + , 252.90, 0.00, 38.90, 38.90 /
203 DATA N1( 44),NP( 44) / 325, 14 /
204 DATA (TMP(I),I= 325, 338) /
205 + 22117.20, 3224.00, 2966.90, 2837.90
206 + , 585.00, 482.80, 460.60, 283.60
207 + , 279.40, 74.90, 43.10, 43.10
208 + , 2.00, 2.00 /
209 DATA N1( 45),NP( 45) / 339, 14 /
210 DATA (TMP(I),I= 339, 352) /
211 + 23219.90, 3411.90, 3146.10, 3003.80
212 + , 627.10, 521.00, 496.20, 311.70
213 + , 307.00, 81.00, 47.90, 47.90
214 + , 2.50, 2.50 /
215 DATA N1( 46),NP( 46) / 353, 14 /
216 DATA (TMP(I),I= 353, 366) /
217 + 24350.30, 3604.30, 3330.30, 3173.30
218 + , 669.90, 559.10, 531.50, 340.00
219 + , 334.70, 86.40, 51.10, 51.10
220 + , 1.50, 1.50 /
221 DATA N1( 47),NP( 47) / 367, 14 /
222 DATA (TMP(I),I= 367, 380) /
223 + 25514.00, 3805.80, 3523.70, 3351.10
224 + , 717.50, 602.40, 571.40, 372.80
225 + , 366.70, 95.20, 62.60, 55.90
226 + , 3.30, 3.30 /
227 DATA N1( 48),NP( 48) / 381, 14 /
228 DATA (TMP(I),I= 381, 394) /
229 + 26711.20, 4018.00, 3727.00, 3537.50
230 + , 770.20, 650.70, 616.50, 410.50
231 + , 403.70, 107.60, 66.90, 66.90
232 + , 9.30, 9.30 /
233 DATA N1( 49),NP( 49) / 395, 19 /
234 DATA (TMP(I),I= 395, 413) /
235 + 27939.90, 4237.50, 3938.00, 3730.10
236 + , 825.60, 702.20, 664.30, 450.80
237 + , 443.10, 121.90, 77.40, 77.40
238 + , 16.20, 16.20, 0.00, 0.00
239 + , 0.10, 0.80, 0.80 /
240 DATA N1( 50),NP( 50) / 414, 19 /
241 DATA (TMP(I),I= 414, 432) /
242 + 29200.10, 4464.70, 4156.10, 3928.80
243 + , 883.80, 756.40, 714.40, 493.30
244 + , 484.80, 136.50, 88.60, 88.60
245 + , 23.90, 23.90, 0.00, 0.00
246 + , 0.90, 1.10, 1.10 /
247 DATA N1( 51),NP( 51) / 433, 19 /
248 DATA (TMP(I),I= 433, 451) /
249 + 30491.20, 4698.30, 4380.40, 4132.20
250 + , 943.70, 811.90, 765.60, 536.90
251 + , 527.50, 152.00, 98.40, 98.40
252 + , 31.40, 31.40, 0.00, 0.00
253 + , 6.70, 2.10, 2.10 /
254 DATA N1( 52),NP( 52) / 452, 19 /
255 DATA (TMP(I),I= 452, 470) /
256 + 31813.80, 4939.20, 4612.00, 4341.40
257 + , 1006.00, 869.70, 818.70, 582.50
258 + , 572.10, 168.30, 110.20, 110.20
259 + , 39.80, 39.80, 0.00, 0.00
260 + , 11.60, 2.30, 2.30 /
261 DATA N1( 53),NP( 53) / 471, 19 /
262 DATA (TMP(I),I= 471, 489) /
263 + 33169.40, 5188.10, 4852.10, 4557.10
264 + , 1072.10, 930.50, 874.60, 631.30
265 + , 619.40, 186.40, 122.70, 122.70
266 + , 49.60, 49.60, 0.00, 0.00
267 + , 13.60, 3.30, 3.30 /
268 DATA N1( 54),NP( 54) / 490, 19 /
269 DATA (TMP(I),I= 490, 508) /
270 + 34561.40, 5452.80, 5103.70, 4782.20
271 + , 1145.70, 999.00, 937.00, 689.00
272 + , 672.30, 213.20, 157.00, 146.70
273 + , 69.50, 67.50, 0.00, 0.00
274 + , 23.39, 13.43, 12.13 /
275 DATA N1( 55),NP( 55) / 509, 19 /
276 DATA (TMP(I),I= 509, 527) /
277 + 35984.60, 5714.30, 5359.40, 5011.90
278 + , 1217.10, 1065.00, 997.60, 739.50
279 + , 725.50, 230.80, 172.30, 161.60
280 + , 78.80, 76.50, 0.00, 0.00
281 + , 22.70, 13.10, 11.40 /
282 DATA N1( 56),NP( 56) / 528, 19 /
283 DATA (TMP(I),I= 528, 546) /
284 + 37440.60, 5988.80, 5623.60, 5247.00
285 + , 1292.80, 1136.70, 1062.20, 796.10
286 + , 780.70, 253.00, 191.80, 179.70
287 + , 92.50, 89.90, 0.00, 0.00
288 + , 39.10, 16.60, 14.60 /
289 DATA N1( 57),NP( 57) / 547, 19 /
290 DATA (TMP(I),I= 547, 565) /
291 + 38924.60, 6266.30, 5890.60, 5482.70
292 + , 1361.30, 1204.40, 1123.40, 848.50
293 + , 831.70, 270.40, 205.80, 191.40
294 + , 98.90, 98.90, 0.00, 0.00
295 + , 32.30, 14.40, 14.40 /
296 DATA N1( 58),NP( 58) / 566, 19 /
297 DATA (TMP(I),I= 566, 584) /
298 + 40443.00, 6548.80, 6164.20, 5723.40
299 + , 1434.60, 1272.80, 1185.40, 901.30
300 + , 883.30, 289.60, 223.30, 207.20
301 + , 110.00, 110.00, 0.10, 0.10
302 + , 37.80, 19.80, 19.80 /
303 DATA N1( 59),NP( 59) / 585, 19 /
304 DATA (TMP(I),I= 585, 603) /
305 + 41990.60, 6834.80, 6440.40, 5964.30
306 + , 1511.00, 1337.40, 1242.20, 951.10
307 + , 931.00, 304.50, 236.30, 217.60
308 + , 113.20, 113.20, 2.00, 2.00
309 + , 37.40, 22.30, 22.30 /
310 DATA N1( 60),NP( 60) / 604, 19 /
311 DATA (TMP(I),I= 604, 622) /
312 + 43568.90, 7126.00, 6721.50, 6207.90
313 + , 1575.30, 1402.80, 1297.40, 999.90
314 + , 977.70, 315.20, 243.30, 224.60
315 + , 117.50, 117.50, 1.50, 1.50
316 + , 37.50, 21.20, 21.80 /
317 DATA N1( 61),NP( 61) / 623, 14 /
318 DATA (TMP(I),I= 623, 636) /
319 + 45184.00, 7427.90, 7012.80, 6459.30
320 + , 1650.00, 1471.40, 1356.90, 1051.50
321 + , 1026.90, 0.00, 242.00, 242.00
322 + , 120.40, 120.40 /
323 DATA N1( 62),NP( 62) / 637, 19 /
324 DATA (TMP(I),I= 637, 655) /
325 + 46834.20, 7736.80, 7311.80, 6716.20
326 + , 1722.80, 1540.70, 1419.80, 1106.00
327 + , 1080.20, 345.70, 265.60, 247.40
328 + , 129.00, 129.00, 5.50, 5.50
329 + , 37.40, 21.30, 21.30 /
330 DATA N1( 63),NP( 63) / 656, 19 /
331 DATA (TMP(I),I= 656, 674) /
332 + 48519.00, 8052.00, 7617.10, 6976.90
333 + , 1800.00, 1613.90, 1480.60, 1160.60
334 + , 1130.90, 360.20, 283.90, 256.60
335 + , 133.20, 133.20, 0.00, 0.00
336 + , 31.80, 22.00, 22.00 /
337 DATA N1( 64),NP( 64) / 675, 19 /
338 DATA (TMP(I),I= 675, 693) /
339 + 50239.10, 8375.60, 7930.30, 7242.80
340 + , 1880.80, 1688.30, 1544.00, 1217.20
341 + , 1185.20, 375.80, 288.50, 270.90
342 + , 140.50, 140.50, 0.10, 0.10
343 + , 36.10, 20.30, 20.30 /
344 DATA N1( 65),NP( 65) / 694, 19 /
345 DATA (TMP(I),I= 694, 712) /
346 + 51995.70, 8708.00, 8251.60, 7514.00
347 + , 1967.50, 1767.70, 1611.30, 1275.00
348 + , 1241.20, 397.90, 310.20, 285.00
349 + , 147.00, 147.00, 2.60, 2.60
350 + , 39.00, 25.40, 25.40 /
351 DATA N1( 66),NP( 66) / 713, 19 /
352 DATA (TMP(I),I= 713, 731) /
353 + 53788.50, 9045.80, 8580.60, 7790.10
354 + , 2046.80, 1841.80, 1675.60, 1332.50
355 + , 1294.90, 416.30, 331.80, 292.90
356 + , 154.20, 154.20, 4.20, 4.20
357 + , 62.90, 26.30, 26.30 /
358 DATA N1( 67),NP( 67) / 732, 19 /
359 DATA (TMP(I),I= 732, 750) /
360 + 55617.70, 9394.20, 8917.80, 8071.10
361 + , 2128.30, 1922.80, 1741.20, 1391.50
362 + , 1351.40, 435.70, 343.50, 306.60
363 + , 161.00, 161.00, 3.70, 3.70
364 + , 51.20, 20.30, 20.30 /
365 DATA N1( 68),NP( 68) / 751, 19 /
366 DATA (TMP(I),I= 751, 769) /
367 + 57485.50, 9751.30, 9264.30, 8357.90
368 + , 2206.50, 2005.80, 1811.80, 1453.30
369 + , 1409.30, 449.10, 366.20, 320.00
370 + , 176.70, 167.60, 4.30, 4.30
371 + , 59.80, 29.40, 29.40 /
372 DATA N1( 69),NP( 69) / 770, 19 /
373 DATA (TMP(I),I= 770, 788) /
374 + 59389.60, 10115.70, 9616.90, 8648.00
375 + , 2306.80, 2089.80, 1884.50, 1514.60
376 + , 1467.70, 471.70, 385.90, 336.60
377 + , 179.60, 179.60, 5.30, 5.30
378 + , 53.20, 32.30, 32.30 /
379 DATA N1( 70),NP( 70) / 789, 19 /
380 DATA (TMP(I),I= 789, 807) /
381 + 61332.30, 10486.40, 9978.20, 8943.60
382 + , 2398.10, 2173.00, 1949.80, 1576.30
383 + , 1527.80, 487.20, 396.70, 343.50
384 + , 198.10, 184.90, 6.30, 6.30
385 + , 54.10, 23.40, 23.40 /
386 DATA N1( 71),NP( 71) / 808, 19 /
387 DATA (TMP(I),I= 808, 826) /
388 + 63313.80, 10870.40, 10348.60, 9244.10
389 + , 2491.20, 2263.50, 2023.60, 1639.40
390 + , 1588.50, 506.20, 410.10, 359.30
391 + , 204.80, 195.00, 6.90, 6.90
392 + , 56.80, 28.00, 28.00 /
393 DATA N1( 72),NP( 72) / 827, 19 /
394 DATA (TMP(I),I= 827, 845) /
395 + 65350.80, 11270.70, 10739.40, 9560.70
396 + , 2600.90, 2365.40, 2107.60, 1716.40
397 + , 1661.70, 538.10, 437.00, 380.40
398 + , 223.80, 213.70, 17.10, 17.10
399 + , 64.90, 38.10, 30.60 /
400 DATA N1( 73),NP( 73) / 846, 21 /
401 DATA (TMP(I),I= 846, 866) /
402 + 67416.40, 11681.50, 11136.10, 9881.10
403 + , 2708.00, 2468.70, 2194.00, 1793.20
404 + , 1735.10, 565.50, 464.80, 404.50
405 + , 241.30, 229.30, 25.00, 25.00
406 + , 71.10, 44.90, 36.40, 5.70, 5.70 /
407 DATA N1( 74),NP( 74) / 867, 21 /
408 DATA (TMP(I),I= 867, 887) /
409 + 69525.00, 12099.80, 11544.00, 10206.80
410 + , 2819.60, 2574.90, 2281.00, 1871.60
411 + , 1809.20, 595.00, 491.60, 425.30
412 + , 258.80, 245.40, 36.50, 33.60
413 + , 77.10, 46.80, 35.60, 6.10, 6.10 /
414 DATA N1( 75),NP( 75) / 888, 21 /
415 DATA (TMP(I),I= 888, 908) /
416 + 71676.40, 12526.70, 11958.70, 10535.30
417 + , 2931.70, 2681.60, 2367.30, 1948.90
418 + , 1882.90, 625.00, 517.90, 444.40
419 + , 273.70, 260.20, 40.60, 40.60
420 + , 82.80, 45.60, 34.60, 3.50, 3.50 /
421 DATA N1( 76),NP( 76) / 909, 19 /
422 DATA (TMP(I),I= 909, 927) /
423 + 73870.80, 12968.00, 12385.00, 10870.90
424 + , 3048.50, 2792.20, 2457.20, 2030.80
425 + , 1960.10, 654.30, 546.50, 468.20
426 + , 289.40, 272.80, 46.30, 43.40
427 + , 83.70, 58.00, 45.40 /
428 DATA N1( 77),NP( 77) / 928, 21 /
429 DATA (TMP(I),I= 928, 948) /
430 + 76111.00, 13418.50, 12824.10, 11215.20
431 + , 3173.70, 2908.70, 2550.70, 2116.10
432 + , 2040.40, 690.10, 577.10, 494.30
433 + , 311.40, 294.90, 63.40, 60.50
434 + , 95.20, 63.00, 50.50, 1.80, 1.80 /
435 DATA N1( 78),NP( 78) / 949, 21 /
436 DATA (TMP(I),I= 949, 969) /
437 + 78394.80, 13879.90, 13272.60, 11563.70
438 + , 3296.00, 3026.50, 2645.40, 2201.90
439 + , 2121.60, 722.00, 609.20, 519.00
440 + , 330.80, 313.30, 74.30, 71.10
441 + , 101.70, 65.30, 51.70, 2.20, 2.20 /
442 DATA N1( 79),NP( 79) / 970, 20 /
443 DATA (TMP(I),I= 970, 989) /
444 + 80724.90, 14352.80, 13733.60, 11918.70
445 + , 3424.90, 3147.80, 2743.00, 2291.10
446 + , 2205.70, 758.80, 643.70, 545.40
447 + , 352.00, 333.90, 86.40, 82.80
448 + , 107.80, 71.70, 53.70, 2.50 /
449 DATA N1( 80),NP( 80) / 990, 21 /
450 DATA (TMP(I),I= 990,1010) /
451 + 83102.30, 14839.30, 14208.70, 12283.90
452 + , 3561.60, 3278.50, 2847.10, 2384.90
453 + , 2294.90, 800.30, 676.90, 571.00
454 + , 378.30, 359.80, 102.20, 98.50
455 + , 120.30, 80.50, 57.60, 6.40, 6.40 /
456 DATA N1( 81),NP( 81) / 1011, 21 /
457 DATA (TMP(I),I=1011,1031) /
458 + 85530.40, 15346.70, 14697.90, 12657.50
459 + , 3704.10, 3425.70, 2956.60, 2485.10
460 + , 2389.30, 845.50, 721.30, 609.00
461 + , 406.60, 386.20, 122.80, 118.50
462 + , 136.30, 99.60, 75.40, 15.30, 13.10 /
463 DATA N1( 82),NP( 82) / 1032, 24 /
464 DATA (TMP(I),I=1032,1055) /
465 + 88004.50, 15860.80, 15200.00, 13035.20
466 + , 3850.70, 3554.20, 3066.40, 2585.60
467 + , 2484.00, 893.60, 763.90, 644.50
468 + , 435.20, 412.90, 142.90, 138.10
469 + , 147.30, 104.80, 86.00, 21.80
470 + , 19.20, 3.10, 0.70, 0.70 /
471 DATA N1( 83),NP( 83) / 1056, 24 /
472 DATA (TMP(I),I=1056,1079) /
473 + 90525.90, 16387.50, 15711.10, 13418.60
474 + , 3999.10, 3696.30, 3176.90, 2687.60
475 + , 2579.60, 938.20, 805.30, 678.90
476 + , 463.60, 440.00, 161.90, 157.40
477 + , 159.30, 116.80, 92.80, 26.50
478 + , 24.40, 0.00, 2.70, 2.70 /
479 DATA N1( 84),NP( 84) / 1080, 21 /
480 DATA (TMP(I),I=1080,1100) /
481 + 93105.00, 16939.30, 16244.30, 13813.80
482 + , 4149.40, 3854.10, 3301.90, 2798.00
483 + , 2683.00, 995.30, 851.00, 705.00
484 + , 500.20, 473.40, 0.00, 0.00
485 + , 0.00, 0.00, 0.00, 31.40, 31.40 /
486 DATA N1( 85),NP( 85) / 1101, 14 /
487 DATA (TMP(I),I=1101,1114) /
488 + 95729.90, 17493.00, 16784.70, 14213.50
489 + , 4317.00, 4008.00, 3426.00, 2908.70
490 + , 2786.70, 1042.00, 886.00, 740.00
491 + , 533.20, 506.00 /
492 DATA N1( 86),NP( 86) / 1115, 14 /
493 DATA (TMP(I),I=1115,1128) /
494 + 98404.00, 18049.00, 17337.10, 14619.40
495 + , 4482.00, 4159.00, 3538.00, 3021.50
496 + , 2892.40, 1097.00, 929.00, 768.00
497 + , 566.60, 541.00 /
498 DATA N1( 87),NP( 87) / 1129, 14 /
499 DATA (TMP(I),I=1129,1142) /
500 + 101137.00, 18639.00, 17906.50, 15031.20
501 + , 4652.00, 4327.00, 3663.00, 3136.20
502 + , 2999.90, 1153.00, 980.00, 810.00
503 + , 603.30, 577.00 /
504 DATA N1( 88),NP( 88) / 1143, 24 /
505 DATA (TMP(I),I=1143,1166) /
506 + 103921.90, 19236.70, 18484.30, 15444.40
507 + , 4822.00, 4489.50, 3791.80, 3248.40
508 + , 3104.90, 1208.40, 1057.60, 879.10
509 + , 635.90, 602.70, 298.90, 298.90
510 + , 254.40, 200.40, 152.80, 67.20
511 + , 67.20, 43.50, 18.80, 18.80 /
512 DATA N1( 89),NP( 89) / 1167, 14 /
513 DATA (TMP(I),I=1167,1180) /
514 + 106755.30, 19840.00, 19083.20, 15871.00
515 + , 5002.00, 4656.00, 3909.00, 3370.20
516 + , 3219.00, 1269.00, 1080.00, 890.00
517 + , 674.90, 644.00 /
518 DATA N1( 90),NP( 90) / 1181, 24 /
519 DATA (TMP(I),I=1181,1204) /
520 + 109650.90, 20472.10, 19693.20, 16300.30
521 + , 5182.30, 4830.40, 4046.10, 3490.80
522 + , 3332.00, 1329.50, 1168.20, 967.30
523 + , 714.10, 676.40, 344.40, 335.20
524 + , 290.20, 229.40, 181.80, 94.30
525 + , 87.90, 59.50, 49.00, 43.00 /
526 DATA N1( 91),NP( 91) / 1205, 21 /
527 DATA (TMP(I),I=1205,1225) /
528 + 112601.40, 21104.60, 20313.70, 16733.10
529 + , 5366.90, 5000.90, 4173.80, 3611.20
530 + , 3441.80, 1387.10, 1224.30, 1006.70
531 + , 743.40, 708.20, 371.20, 359.50
532 + , 309.60, 222.90, 222.90, 94.10, 94.10 /
533 DATA N1( 92),NP( 92) / 1226, 24 /
534 DATA (TMP(I),I=1226,1249) /
535 + 115606.10, 21757.40, 20947.60, 17166.30
536 + , 5548.00, 5182.20, 4303.40, 3727.60
537 + , 3551.70, 1440.80, 1272.60, 1044.90
538 + , 780.40, 737.70, 391.30, 380.90
539 + , 323.70, 259.30, 195.10, 105.00
540 + , 96.30, 70.70, 42.30, 32.30 /
541 DATA N1( 93),NP( 93) / 1250, 21 /
542 DATA (TMP(I),I=1250,1270) /
543 + 118678.00, 22426.80, 21600.50, 17610.00
544 + , 5723.20, 5366.20, 4434.70, 3850.30
545 + , 3665.80, 1500.70, 1327.70, 1086.80
546 + , 815.90, 770.30, 415.00, 404.40
547 + , 0.00, 283.40, 206.10, 109.30, 101.30 /
548 DATA N1( 94),NP( 94) / 1271, 21 /
549 DATA (TMP(I),I=1271,1291) /
550 + 121818.00, 23097.20, 22266.20, 18056.80
551 + , 5932.90, 5541.20, 4556.60, 3972.60
552 + , 3778.10, 1558.60, 1372.10, 1114.80
553 + , 848.90, 801.40, 445.80, 432.40
554 + , 351.90, 274.10, 206.50, 116.00, 105.40 /
555 DATA N1( 95),NP( 95) / 1292, 21 /
556 DATA (TMP(I),I=1292,1312) /
557 + 125027.00, 23772.90, 22944.00, 18504.10
558 + , 6120.50, 5710.20, 4667.00, 4092.10
559 + , 3886.90, 1617.10, 1411.80, 1135.70
560 + , 878.70, 827.60, 0.00, 0.00
561 + , 0.00, 0.00, 0.00, 115.80, 103.30 /
562 DATA N1( 96),NP( 96) / 1313, 21 /
563 DATA (TMP(I),I=1313,1333) /
564 + 128220.00, 24460.00, 23779.00, 18930.00
565 + , 6288.00, 5895.00, 4797.00, 4227.00
566 + , 3971.00, 1643.00, 1440.00, 1154.00
567 + , 0.00, 853.00, 0.00, 0.00
568 + , 385.00, 0.00, 0.00, 0.00, 100.00 /
569 DATA N1( 97),NP( 97) / 1334, 21 /
570 DATA (TMP(I),I=1334,1354) /
571 + 131590.00, 25275.00, 24385.00, 19452.00
572 + , 6556.00, 6147.00, 4977.00, 4366.00
573 + , 4132.00, 1755.00, 1554.00, 1235.00
574 + , 0.00, 877.00, 0.00, 0.00
575 + , 398.00, 0.00, 0.00, 0.00, 100.00 /
576 DATA N1( 98),NP( 98) / 1355, 21 /
577 DATA (TMP(I),I=1355,1375) /
578 + 135960.00, 26110.00, 25250.00, 19930.00
579 + , 6754.00, 6359.00, 5109.00, 4497.00
580 + , 4253.00, 1799.00, 1616.00, 1279.00
581 + , 0.00, 902.00, 0.00, 0.00
582 + , 419.00, 0.00, 0.00, 0.00, 100.00 /
583 DATA N1( 99),NP( 99) / 1376, 17 /
584 DATA (TMP(I),I=1376,1392) /
585 + 139490.00, 26900.00, 26020.00, 20410.00
586 + , 6977.00, 6574.00, 5252.00, 4630.00
587 + , 4374.00, 1868.00, 1680.00, 1321.00
588 + , 0.00, 927.00, 0.00, 0.00, 435.00 /
589 DATA N1(100),NP(100) / 1393, 17 /
590 DATA (TMP(I),I=1393,1409) /
591 + 143090.00, 27700.00, 26810.00, 20900.00
592 + , 7205.00, 6793.00, 5397.00, 4766.00
593 + , 4498.00, 1937.00, 1747.00, 1366.00
594 + , 0.00, 952.00, 0.00, 0.00, 454.00 /
595 DO 10 I = 1,LENGTH
596 ESHELL(I) = TMP(I)
597 10 CONTINUE
598 DO 20 I = 1,MAXSHL
599 NSHLST(I) = NS(I)
600 20 CONTINUE
601 DO 30 I = 1,MAXELZ
602 N1ST(I) = N1(I)
603 NSHLLS(I) = NP(I)
604 30 CONTINUE
605 END