This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / GEANT321 / gphys / gshlin.F
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
13 C.
14 C.    ******************************************************************
15 C.    *                                                                *
16 C     *       Initialize atomic shells' potentials                     *
17 C.    *                                                                *
18 C.    * Potentials according to Rev.Mod.Phys. vol. 39 nr. 1 125 (1967) *
19 C.    * Xe upadetd according to Grishin et al. NIM A309, 476 (1991)    *
20 C.    *                                                                *
21 C.    *    ==>Called by : GPHINI                                       *
22 C.    *       Author    J. Chwastowski                                 *
23 C.    *                                                                *
24 C.    ******************************************************************
25 C.
26 C  Shells
27 c  See alse Handbook of Phys. & Chemistry
28 c
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