]> git.uio.no Git - u/mrichter/AliRoot.git/blame - GEANT321/block/flkdt3.F
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Initialize fSepctrum = 0 in ctors
[u/mrichter/AliRoot.git] / GEANT321 / block / flkdt3.F
CommitLineData
fe4da5cc 1*
2* $Id$
3*
4* $Log$
5* Revision 1.1.1.1 1995/10/24 10:19:35 cernlib
6* Geant
7*
8*
9#include "geant321/pilot.h"
10*CMZ : 3.21/02 29/03/94 15.41.45 by S.Giani
11*-- Author :
12*$ CREATE BLKDT3.FOR
13*COPY BLKDT3
14*
15*=== blkdt3 ===========================================================*
16*== *
17 SUBROUTINE FLKDT3
18
19#include "geant321/dblprc.inc"
20#include "geant321/dimpar.inc"
21#include "geant321/iounit.inc"
22* * * reaction channel cross section data *
23#include "geant321/reac.inc"
24 DIMENSION
25 & UMOPI(92), UMOKC(68), UMOP(39), UMON(63), UMOK0(34),
26 & PLAPI(92), PLAKC(68), PLAP(39), PLAN(63), PLAK0(34),
27 & SPIKP1(315), SPIKPU(278), SPIKPV(372),
28 & SPIKPW(278), SPIKPX(372), SPIKP4(315),
29 & SPIKP5(187), SPIKP6(289),
30 & SKMPEL(102), SPIKP7(289), SKMNEL(68), SPIKP8(187),
31 & SPIKP9(143), SPIKP0(169), SPKPV(143),
32 & SAPPEL(105), SPIKPE(399), SAPNEL(84), SPIKPZ(273),
33 & SANPEL(84) , SPIKPF(273),
34 & SPKP15(187), SPKP16(272),
35 & NRKPI(164), NRKKC(132), NRKP(70), NRKN(116), NRKK0(54),
36 & NURELN(60)
37*
38 DIMENSION NRKLIN(532)
39 EQUIVALENCE (NRK(1,1), NRKLIN(1))
40 EQUIVALENCE ( UMO( 1), UMOPI(1)), ( UMO( 93), UMOKC(1))
41 EQUIVALENCE ( UMO(161), UMOP(1)), ( UMO(200), UMON(1))
42 EQUIVALENCE ( UMO(263), UMOK0(1))
43 EQUIVALENCE ( PLABF( 1), PLAPI(1)), ( PLABF( 93), PLAKC(1))
44 EQUIVALENCE ( PLABF(161), PLAP(1)), ( PLABF(200), PLAN(1))
45 EQUIVALENCE ( PLABF(263), PLAK0(1))
46 EQUIVALENCE ( WK( 1), SPIKP1(1)), ( WK( 316), SPIKPU(1))
47 EQUIVALENCE ( WK( 594), SPIKPV(1)), ( WK( 966), SPIKPW(1))
48 EQUIVALENCE ( WK(1244), SPIKPX(1)), ( WK(1616), SPIKP4(1))
49 EQUIVALENCE ( WK(1931), SPIKP5(1)), ( WK(2118), SPIKP6(1))
50 EQUIVALENCE ( WK(2407), SKMPEL(1)), ( WK(2509), SPIKP7(1))
51 EQUIVALENCE ( WK(2798), SKMNEL(1)), ( WK(2866), SPIKP8(1))
52 EQUIVALENCE ( WK(3053), SPIKP9(1)), ( WK(3196), SPIKP0(1))
53 EQUIVALENCE ( WK(3365), SPKPV(1)), ( WK(3508), SAPPEL(1))
54 EQUIVALENCE ( WK(3613), SPIKPE(1)), ( WK(4012), SAPNEL(1))
55 EQUIVALENCE ( WK(4096), SPIKPZ(1)), ( WK(4369), SANPEL(1))
56 EQUIVALENCE ( WK(4453), SPIKPF(1)), ( WK(4726), SPKP15(1))
57 EQUIVALENCE ( WK(4913), SPKP16(1))
58 EQUIVALENCE (NRK(1,1), NRKLIN(1))
59 EQUIVALENCE (NRKLIN( 1), NRKPI(1)), (NRKLIN( 165), NRKKC(1))
60 EQUIVALENCE (NRKLIN( 297), NRKP(1)), (NRKLIN( 367), NRKN(1))
61 EQUIVALENCE (NRKLIN( 483), NRKK0(1))
62 EQUIVALENCE (NURE(1,1), NURELN(1))
63*
64 DIMENSION UMOPI1(92), UMOKC1(68), UMOP1(39), UMON1(63),
65 & UMOK01(34), PLAPI1(92), PLAKC1(68), PLAP1(39), PLAN1(63),
66 & PLAK01(34), SIIN1(296),
67 & SPKP11(315), SPKPU1(278), CPKPV1(372),
68 & SPKPW1(278), SPKPX1(372), SPKP41(315),
69 & SPKP51(187), SPKP61(289),
70 & SKMPL1(102), SPKP71(289), SKMNL1(68), SPKP81(187),
71 & SPKP91(143), SPKP01(169), SPKPV1(143),
72 & SAPPL1(105), SPKPE1(399), SAPNL1(84), SPKPZ1(273),
73 & SANPL1(84) , SPKPF1(273),
74 & CPKP51(187), CPKP61(272),
75 & NRKPI1(164), NRKKC1(132), NRKP1(70), NRKN1(116),
76 & NRKK01(54), NURE1(60)
77*
78**** pi- p data *
79**** pi+ n data *
80 DATA PLAPI1/ 0.D0, .3D0, .5D0, .6D0, .7D0, .8D0, .9D0, .95D0,1.D0,
81 & 1.15D0, 1.3D0, 1.5D0, 1.6D0, 1.8D0, 2.D0, 2.3D0, 2.5D0, 2.8D0,
82 & 3.D0, 3.5D0, 4.D0, 0.D0, .285D0, .4D0, .45D0, .5D0, .6D0, .7D0,
83 & .75D0, .8D0, .85D0, .9D0, 1.D0, 1.15D0, 1.3D0, 1.5D0, 1.6D0,
84 & 1.8D0, 2.D0, 2.3D0, 2.5D0, 2.8D0, 3.D0, 3.5D0, 4.D0, 4.5D0, 0.D0,
85 & .285D0, .4D0, .45D0, .5D0, .6D0, .7D0, .75D0, .8D0, .85D0, .9D0,
86 & 1.D0, 1.15D0, 1.3D0, 1.5D0, 1.6D0, 1.8D0, 2.D0, 2.3D0, 2.5D0,
87 & 2.8D0, 3.D0, 3.5D0, 4.D0, 4.5D0, 0.D0, .3D0, .5D0, .6D0, .7D0,
88 & .8D0, .9D0, .95D0, 1.D0, 1.15D0, 1.3D0, 1.5D0, 1.6D0, 1.8D0,
89 & 2.D0, 2.3D0, 2.5D0, 2.8D0, 3.D0, 3.5D0, 4.D0 /
90 DATA PLAKC1/
91 & 0.D0, .58D0, .8D0, 1.01D0, 1.23D0, 1.45D0, 1.68D0, 1.94D0,
92 & 2.18D0, 2.42D0, 2.68D0, 2.96D0, 3.24D0,
93 & 3.51D0, 3.84D0, 4.16D0, 4.49D0,
94 & 0.D0, .58D0, .8D0, 1.01D0, 1.23D0, 1.45D0, 1.68D0, 1.94D0,
95 & 2.18D0, 2.42D0, 2.68D0, 2.96D0, 3.24D0,
96 & 3.51D0, 3.84D0, 4.16D0, 4.49D0,
97 & 0.D0, .58D0, .8D0, 1.01D0, 1.23D0, 1.45D0, 1.68D0, 1.94D0,
98 & 2.18D0, 2.42D0, 2.68D0, 2.96D0, 3.24D0,
99 & 3.51D0, 3.84D0, 4.16D0, 4.49D0,
100 & 0.D0, .58D0, .8D0, 1.01D0, 1.23D0, 1.45D0, 1.68D0, 1.94D0,
101 & 2.18D0, 2.42D0, 2.68D0, 2.96D0, 3.24D0,
102 & 3.51D0, 3.84D0, 4.16D0, 4.49D0/
103 DATA PLAK01/
104 & 0.D0, .58D0, .8D0, 1.01D0, 1.23D0, 1.45D0, 1.68D0, 1.94D0,
105 & 2.18D0, 2.42D0, 2.68D0, 2.96D0, 3.24D0,
106 & 3.51D0, 3.84D0, 4.16D0, 4.49D0,
107 & 0.D0, .58D0, .8D0, 1.01D0, 1.23D0, 1.45D0, 1.68D0, 1.94D0,
108 & 2.18D0, 2.42D0, 2.68D0, 2.96D0, 3.24D0,
109 & 3.51D0, 3.84D0, 4.16D0, 4.49D0/
110* pp pn np nn *
111 DATA PLAP1/
112 & 0.D0, 1.06D0, 1.34D0, 1.63D0, 1.92D0, 2.2D0, 2.5D0,2.8D0,3.1D0,
113 & 3.43D0, 3.75D0, 4.07D0, 4.43D0,
114 & 0.D0, 1.06D0, 1.34D0, 1.63D0, 1.92D0, 2.2D0, 2.5D0,2.8D0,3.1D0,
115 & 3.43D0, 3.75D0, 4.07D0, 4.43D0,
116 & 0.D0, 1.06D0, 1.34D0, 1.63D0, 1.92D0, 2.2D0, 2.5D0,2.8D0,3.1D0,
117 & 3.43D0, 3.75D0, 4.07D0, 4.43D0 /
118* app apn anp ann *
119 DATA PLAN1/
120 & 0.D0, 1.D-3, .1D0, .2D0, .3D0, .4D0, .5D0, .6D0,
121 & .74D0, 1.06D0, 1.34D0, 1.63D0, 1.92D0, 2.2D0, 2.5D0,2.8D0,3.1D0,
122 & 3.43D0, 3.75D0, 4.07D0, 4.43D0,
123 & 0.D0, 1.D-3, .1D0, .2D0, .3D0, .4D0, .5D0, .6D0,
124 & .74D0, 1.06D0, 1.34D0, 1.63D0, 1.92D0, 2.2D0, 2.5D0,2.8D0,3.1D0,
125 & 3.43D0, 3.75D0, 4.07D0, 4.43D0,
126 & 0.D0, 1.D-3, .1D0, .2D0, .3D0, .4D0, .5D0, .6D0,
127 & .74D0, 1.06D0, 1.34D0, 1.63D0, 1.92D0, 2.2D0, 2.5D0,2.8D0,3.1D0,
128 & 3.43D0, 3.75D0, 4.07D0, 4.43D0 /
129 DATA SIIN1/ 296*0.D0 /
130 DATA UMOPI1/ 1.08D0,1.233D0,1.302D0,1.369D0,1.496D0,
131 & 1.557D0,1.615D0,1.6435D0,
132 & 1.672D0,1.753D0,1.831D0,1.930D0,1.978D0,2.071D0,2.159D0,
133 & 2.286D0,2.366D0,2.482D0,2.56D0,
134 & 2.735D0,2.90D0,
135 & 1.08D0,1.222D0,1.302D0,1.3365D0,1.369D0,1.434D0,
136 & 1.496D0,1.527D0,1.557D0,
137 & 1.586D0,1.615D0,1.672D0,1.753D0,1.831D0,1.930D0,1.978D0,
138 & 2.071D0,2.159D0,2.286D0,2.366D0,
139 & 2.482D0,2.560D0,2.735D0,2.90D0,3.06D0,
140 & 1.08D0,1.222D0,1.302D0,1.3365D0,1.369D0,1.434D0,
141 & 1.496D0,1.527D0,1.557D0,
142 & 1.586D0,1.615D0,1.672D0,1.753D0,1.831D0,1.930D0,1.978D0,
143 & 2.071D0,2.159D0,2.286D0,2.366D0,
144 & 2.482D0,2.560D0,2.735D0,2.90D0,3.06D0,
145 & 1.08D0,1.233D0,1.302D0,1.369D0,1.496D0,
146 & 1.557D0,1.615D0,1.6435D0,
147 & 1.672D0,1.753D0,1.831D0,1.930D0,1.978D0,2.071D0,2.159D0,
148 & 2.286D0,2.366D0,2.482D0,2.56D0,
149 & 2.735D0, 2.90D0/
150 DATA UMOKC1/ 1.44D0,
151 & 1.598D0,1.7D0,1.8D0,1.9D0,2.0D0,2.1D0,2.2D0,2.3D0,2.4D0,2.5D0,
152 & 2.6D0,2.7D0,2.8D0,2.9D0,3.0D0,
153 & 3.1D0,1.44D0,
154 & 1.598D0,1.7D0,1.8D0,1.9D0,2.0D0,2.1D0,2.2D0,2.3D0,2.4D0,2.5D0,
155 & 2.6D0,2.7D0,2.8D0,2.9D0,3.0D0,
156 & 3.1D0,1.44D0,
157 & 1.598D0,1.7D0,1.8D0,1.9D0,2.0D0,2.1D0,2.2D0,2.3D0,2.4D0,2.5D0,
158 & 2.6D0,2.7D0,2.8D0,2.9D0,3.0D0,
159 & 3.1D0,1.44D0,
160 & 1.598D0,1.7D0,1.8D0,1.9D0,2.0D0,2.1D0,2.2D0,2.3D0,2.4D0,2.5D0,
161 & 2.6D0,2.7D0,2.8D0,2.9D0,3.0D0,
162 & 3.1D0/
163 DATA UMOK01/ 1.44D0,
164 & 1.598D0,1.7D0,1.8D0,1.9D0,2.0D0,2.1D0,2.2D0,2.3D0,2.4D0,2.5D0,
165 & 2.6D0,2.7D0,2.8D0,2.9D0,3.0D0,
166 & 3.1D0,1.44D0,
167 & 1.598D0,1.7D0,1.8D0,1.9D0,2.0D0,2.1D0,2.2D0,2.3D0,2.4D0,2.5D0,
168 & 2.6D0,2.7D0,2.8D0,2.9D0,3.0D0,
169 & 3.1D0/
170* pp pn np nn *
171 DATA UMOP1/
172 & 1.88D0,2.102D0,2.2D0,2.3D0,2.4D0,2.5D0,2.6D0,2.7D0,2.8D0,2.9D0,
173 & 3.D0,3.1D0,3.2D0,
174 & 1.88D0,2.102D0,2.2D0,2.3D0,2.4D0,2.5D0,2.6D0,2.7D0,2.8D0,2.9D0,
175 & 3.D0,3.1D0,3.2D0,
176 & 1.88D0,2.102D0,2.2D0,2.3D0,2.4D0,2.5D0,2.6D0,2.7D0,2.8D0,2.9D0,
177 & 3.D0,3.1D0,3.2D0/
178* app apn anp ann *
179 DATA UMON1 /
180 & 1.877D0,1.87701D0,1.879D0,1.887D0,1.9D0,1.917D0,1.938D0,1.962D0,
181 & 2.D0,2.102D0,2.2D0,2.3D0,2.4D0,2.5D0,2.6D0,2.7D0,2.8D0,2.9D0,
182 & 3.D0,3.1D0,3.2D0,
183 & 1.877D0,1.87701D0,1.879D0,1.887D0,1.9D0,1.917D0,1.938D0,1.962D0,
184 & 2.D0,2.102D0,2.2D0,2.3D0,2.4D0,2.5D0,2.6D0,2.7D0,2.8D0,2.9D0,
185 & 3.D0,3.1D0,3.2D0,
186 & 1.877D0,1.87701D0,1.879D0,1.887D0,1.9D0,1.917D0,1.938D0,1.962D0,
187 & 2.D0,2.102D0,2.2D0,2.3D0,2.4D0,2.5D0,2.6D0,2.7D0,2.8D0,2.9D0,
188 & 3.D0,3.1D0,3.2D0/
189**** reaction channel state particles *
190 DATA NRKPI1/ 13, 1, 15, 21, 81, 0, 13, 54, 23, 53, 13, 63, 13, 58,
191 & 23, 57, 13, 65, 1, 32, 53, 31, 54, 32, 53, 33, 53, 35, 63, 32,
192 & 13, 8, 23, 1, 17, 15, 21, 24, 22, 15, 82, 0, 61, 0, 13, 55, 23,
193 & 54, 14, 53, 13, 64, 23, 63, 13, 59, 23, 58, 14, 57, 13, 66, 23,
194 & 65, 1, 31, 8, 32, 1, 33, 1, 35, 54, 31, 55, 32, 54, 33, 53, 34,
195 & 54, 35, 14, 1, 23, 8, 17, 24, 20, 15, 22, 24, 83, 0, 62, 0, 14,
196 & 54, 23, 55, 13, 56, 14, 63, 23, 64, 14, 58, 23, 59, 13, 60, 14,
197 & 65, 23, 66, 8, 31, 1, 34, 8, 33, 8, 35, 55, 31, 54, 34, 55, 33,
198 & 56, 32, 55, 35, 14, 8, 24, 20, 84, 0, 14, 55, 23, 56, 14, 64, 14,
199 & 59, 23, 60, 14, 66, 8, 34, 56, 31, 55, 34, 56, 33, 56, 35, 64,34/
200 DATA NRKKC1/15, 1, 89, 0, 24, 53, 15, 54, 1, 36, 1, 40, 1, 44, 36,
201 & 63, 15, 63, 45, 53, 44, 54, 15, 8, 24, 1, 91, 0, 24, 54, 15, 55,
202 & 8, 36, 1, 37, 8, 40, 1, 41, 8, 44, 1, 45, 36, 64, 37, 63, 15, 64,
203 & 24, 63, 45, 54, 44, 55, 16, 1, 25, 8, 17, 23, 21, 14, 20,
204 & 13, 22, 23, 90, 0, 38, 1, 39, 8, 16, 54, 25, 55, 1, 42, 8, 43,
205 & 16, 63, 25, 64, 39, 64, 38, 63, 46, 54, 47, 55, 8, 47, 1, 46, 52,
206 & 0, 51, 0, 16, 8, 17, 14, 20, 23, 22, 14, 92, 0, 8, 38, 16, 55,
207 & 25, 56, 8, 42, 16, 64, 38, 64, 46, 55, 47, 56, 8, 46, 94, 0 /
208* *
209* k0 p k0 n ak0 p ak/ n *
210* *
211 DATA NRKK01/ 24, 8, 106, 0, 15, 56, 24, 55, 37, 8, 41, 8, 45, 8,
212 & 37, 64, 24, 64, 44, 56, 45, 55, 25, 1, 17, 13, 22, 13, 21, 23,
213 & 107, 0, 39, 1, 25, 54, 16, 53, 43, 1, 25, 63, 39, 63, 47, 54, 46,
214 & 53, 47, 1, 103, 0, 93, 0/
215* pp pn np nn *
216 DATA NRKP1/ 1, 1, 85, 0, 8, 53, 1, 54, 1, 63, 8, 57, 1, 58, 2*54,
217 & 53, 55, 63, 54, 64, 53, 1, 8, 86, 0, 8, 54, 1, 55, 8, 63, 1, 64,
218 & 8, 58, 1, 59, 64, 54, 63, 55, 54, 55, 53, 56, 77, 0, 2*8, 95, 0,
219 & 8, 55, 1, 56, 8, 64, 8, 59, 1, 60, 2*55, 54, 56, 64, 55, 63, 56 /
220* app apn anp ann *
221 DATA NRKN1/ 1, 2, 17, 18, 15, 16, 8, 9, 13, 14, 99, 0, 87, 0, 1,
222 & 68, 8, 69, 2, 54, 9, 55, 102, 0, 2, 63, 9, 64, 1, 75, 8, 76, 53,
223 & 67, 54, 68, 55, 69, 56, 70, 63, 68, 64, 69, 75, 54, 76, 55, 2, 8,
224 & 18, 20, 16, 24, 14, 23, 101, 0, 88, 0, 2, 55, 9, 56, 1, 67, 8,
225 & 68, 2, 64, 8, 75, 2, 59, 8, 72, 68, 55, 67, 54, 69, 56, 1, 9, 18,
226 & 21, 15, 25, 13, 23, 100, 0, 96, 0, 2, 53, 9, 54, 1, 69, 8, 70, 1,
227 & 76, 9, 63, 1, 73, 9, 58, 55, 70, 53, 68, 54, 69 /
228**** channel cross section *
229 DATA SPKP11/ 0.D0, 300.D0, 40.D0, 20.D0, 13.D0,8.5D0,8.D0, 9.5D0,
230 & 12.D0,14.D0,15.5D0,20.D0,17.D0,13.D0,10.D0,9.D0,8.5D0,8.D0,7.8D0,
231 & 7.3D0, 6.7D0, 9*0.D0,.23D0,.35D0,.7D0,.52D0,.4D0,.3D0,.2D0,.15D0,
232 & .13D0, .11D0, .09D0, .07D0, 0.D0, .033D0,.8D0,1.35D0,1.35D0,.5D0,
233 & 15*0.D0, 3*0.D0,.00D0,0.80D0,2.2D0,3.6D0,4.6D0,4.7D0,3.5D0,2.4D0,
234 &1.8D0,1.4D0,.75D0,.47D0,.25D0,.13D0,.08D0,6*0.D0,0.D0,1.2D0,3.3D0,
235 & 5.4D0,6.9D0,7.3D0,5.3D0,3.6D0,2.7D0,2.2D0,1.1D0,.73D0,.4D0,.22D0,
236 & .12D0,9*0.D0,.0D0,0.D0,2.0D0,4.4D0,6.8D0,9.9D0,7.9D0,6.0D0,3.8D0,
237 &2.5D0,2.D0,1.4D0,1.D0,.6D0,.35D0,10*0.D0,.25D0,.55D0,.75D0,1.25D0,
238 & 1.9D0,2.D0,1.8D0,1.5D0,1.25D0,1.D0,.8D0,6*0.D0,4*0.D0,.4D0,.85D0,
239 & 1.1D0, 1.85D0, 2.8D0, 3.D0,2.7D0,2.2D0,1.85D0,1.5D0,1.2D0,6*0.D0,
240 & 6*0.D0, .5D0, 1.2D0, 1.7D0, 3.4D0, 5.2D0, 6.4D0, 6.1D0, 5.6D0,
241 & 5.2D0, 6*0.D0, 2*0.D0, .0D0, 1.D0, 3.3D0, 5.2D0, 4.45D0, 3.6D0,
242 & 2.75D0, 1.9D0, 1.65D0, 1.3D0, .95D0, .6D0, .45D0, 6*0.D0, 3*0.D0,
243 & .0D0, .45D0, 1.4D0, 1.5D0, 1.1D0, .85D0, .5D0, .3D0, .2D0, .15D0,
244 & 8*0.D0, 5*0.D0, .0D0, .0D0, .6D0, .8D0, .95D0, .8D0, .7D0, .6D0,
245 & .5D0, .4D0, 6*0.D0, 5*0.D0, .0D0, .00D0, .85D0, 1.2D0, 1.4D0,
246 & 1.2D0, 1.05D0, .9D0, .7D0, .55D0, 6*0.D0, 5*0.D0, .0D0, .00D0,
247 & 1.D0, 1.5D0, 3.5D0, 4.15D0, 3.7D0, 2.7D0, 2.3D0, 1.75D0, 6*0.D0,
248 & 10*0.D0, .5D0, 2.0D0, 3.3D0, 5.4D0, 7.D0 /
249**** pi+ n data *
250 DATA SPKPU1/ 0.D0, 25.D0, 13.D0, 11.D0, 10.5D0, 14.D0, 20.D0,
251 & 20.D0, 16.D0, 14.D0, 19.D0, 28.D0, 17.5D0, 13.5D0, 12.D0, 10.5D0,
252 & 10.D0, 10.D0, 9.5D0, 9.D0, 8.D0, 7.5D0, 7.D0, 6.5D0, 6.D0, 0.D0,
253 & 48.D0, 19.D0, 15.D0, 11.5D0, 10.D0, 8.D0, 6.5D0, 5.5D0, 4.8D0,
254 & 4.2D0, 7.5D0, 3.4D0, 2.5D0, 2.5D0, 2.1D0, 1.4D0, 1.D0, .8D0,
255 & .6D0, .46D0, .3D0, .2D0, .15D0, .13D0, 11*0.D0, .95D0, .65D0,
256 & .48D0, .35D0, .2D0, .18D0, .17D0, .16D0, .15D0, .1D0, .09D0,
257 & .065D0, .05D0, .04D0, 12*0.D0, .2D0, .25D0, .25D0, .2D0, .1D0,
258 & .08D0, .06D0, .045D0, .03D0, .02D0, .01D0, .005D0, .003D0,
259 & 12*0.D0, .3D0, .24D0, .18D0, .15D0, .13D0, .12D0, .11D0, .1D0,
260 & .09D0, .08D0, .05D0, .04D0, .03D0, 0.D0, 0.16D0, .7D0, 1.3D0,
261 & 3.1D0, 4.5D0, 2.D0, 18*0.D0, 3*.0D0, 0.D0, 0.D0, 4.0D0, 11.D0,
262 & 11.4D0, 10.3D0, 7.5D0, 6.8D0, 4.75D0, 2.5D0, 1.5D0, .9D0, .55D0,
263 & .35D0, 13*0.D0, .1D0, .34D0, .5D0, .8D0, 1.1D0, 2.25D0, 3.3D0,
264 & 2.3D0, 1.6D0, .95D0, .45D0, .28D0, .15D0, 10*0.D0, 2*0.D0, .17D0,
265 & .64D0, 1.D0, 1.5D0, 2.1D0, 4.25D0, 6.2D0, 4.4D0, 3.D0, 1.8D0,
266 & .9D0, .53D0, .28D0, 10*0.D0, 2*0.D0, .25D0, .82D0,
267 & 1.3D0, 1.9D0, 2.8D0, 5.5D0 , 8.D0, 5.7D0, 3.9D0, 2.35D0, 1.15D0,
268 & .69D0, .37D0, 10*0.D0, 7*0.D0, .0D0, .34D0, 1.5D0, 3.47D0,
269 & 5.87D0, 6.23D0, 4.27D0, 2.6D0, 1.D0, .6D0, .3D0, .15D0, 6*0.D0/
270*
271 DATA CPKPV1/ 7*0.D0, .00D0, .16D0, .75D0, 1.73D0, 2.93D0, 3.12D0,
272 & 2.13D0, 1.3D0, .5D0, .3D0, .15D0, .08D0, 6*0.D0, 10*0.D0, .2D0,
273 & .6D0, .92D0, 2.4D0, 4.9D0, 6.25D0, 5.25D0, 3.5D0, 2.15D0, 1.4D0,
274 & 1.D0, .7D0, 13*0.D0, .13D0, .4D0, .62D0, 1.6D0, 3.27D0, 4.17D0,
275 & 3.5D0, 2.33D0, 1.43D0, .93D0, .66D0, .47D0, 13*0.D0, .07D0, .2D0,
276 & .31D0, .8D0, 1.63D0, 2.08D0, 1.75D0, 1.17D0, .72D0, .47D0, .34D0,
277 & .23D0, 17*0.D0, .33D0, 1.D0, 1.8D0, 2.67D0, 5.33D0, 6.D0, 5.53D0,
278 & 5.D0, 17*0.D0, .17D0, .5D0, .9D0, 1.83D0, 2.67D0, 3.0D0, 2.77D0,
279 & 2.5D0, 3*0.D0, 3*0.D0, 1.D0, 3.3D0, 2.8D0, 2.5D0, 2.3D0, 1.8D0,
280 & 1.5D0, 1.1D0, .8D0, .7D0, .55D0, .3D0, 10*0.D0, 9*0.D0, .1D0,
281 & .4D0, 1.D0, 1.4D0, 2.2D0, 2.5D0, 2.2D0, 1.65D0, 1.35D0, 1.1D0,
282 & .8D0, .6D0, .4D0, 12*0.D0, .15D0, .6D0, 1.5D0, 2.1D0, 3.3D0,
283 & 3.8D0, 3.3D0, 2.45D0, 2.05D0, 1.65D0, 1.2D0, .9D0, .6D0, 3*0.D0,
284 & 9*0.D0, .10D0, .2D0, .5D0, .7D0, 1.3D0, 1.55D0, 1.9D0, 1.8D0,
285 & 1.55D0, 1.35D0, 1.15D0, .95D0, .7D0, 13*0.D0, .2D0, .5D0, .7D0,
286 & 1.3D0, 1.55D0, 1.9D0, 1.8D0, 1.55D0, 1.35D0, 1.15D0, .95D0, .7D0,
287 & 17*0.D0, .2D0, .5D0, .85D0, 2.D0, 2.15D0, 2.05D0, 1.75D0, 1.D0,
288 & 17*0.D0, .13D0, .33D0, .57D0, 1.33D0, 1.43D0, 1.36D0, 1.17D0,
289 & .67D0, 17*0.D0, .07D0, .17D0, .28D0, .67D0, .72D0, .69D0, .58D0,
290 & .33D0,17*0.D0,.4D0, .7D0, 1.D0, 1.6D0, 1.8D0, 2.3D0,1.9D0,1.7D0 /
291**** pi- p data *
292 DATA SPKPW1/ 0.D0, 25.D0, 13.D0, 11.D0, 10.5D0, 14.D0, 2*20.D0,
293 & 16.D0, 14.D0, 19.D0, 28.D0, 17.5D0, 13.5D0, 12.D0, 10.5D0,
294 & 2*10.D0, 9.5D0, 9.D0, 8.D0, 7.5D0, 7.D0, 6.5D0, 6.D0, 0.D0,
295 & 48.D0, 19.D0, 15.D0, 11.5D0, 10.D0, 8.D0, 6.5D0, 5.5D0, 4.8D0,
296 & 4.2D0, 7.5D0, 3.4D0, 2*2.5D0, 2.1D0, 1.4D0, 1.D0, .8D0, .6D0,
297 & .46D0, .3D0, .2D0, .15D0, .13D0, 11*0.D0, .95D0, .65D0, .48D0,
298 & .35D0, .2D0, .18D0, .17D0, .16D0, .15D0, .1D0, .09D0, .065D0,
299 & .05D0, .04D0, 12*0.D0, .2D0, 2*.25D0, .2D0, .1D0, .08D0, .06D0,
300 & .045D0, .03D0, .02D0, .01D0, .005D0, .003D0, 12*0.D0, .3D0,
301 & .24D0, .18D0, .15D0, .13D0, .12D0, .11D0, .1D0, .09D0, .08D0,
302 & .05D0, .04D0, .03D0, 0.D0, 0.16D0, .7D0, 1.3D0, 3.1D0, 4.5D0,
303 & 2.D0, 23*0.D0, 4.0D0, 11.D0, 11.4D0, 10.3D0, 7.5D0, 6.8D0,
304 & 4.75D0, 2.5D0, 1.5D0, .9D0, .55D0, .35D0, 13*0.D0, .1D0, .34D0,
305 & .5D0, .8D0, 1.1D0, 2.25D0, 3.3D0, 2.3D0, 1.6D0, .95D0, .45D0,
306 & .28D0, .15D0, 12*0.D0, .17D0, .64D0, 1.D0, 1.5D0, 2.1D0, 4.25D0,
307 & 6.2D0, 4.4D0, 3.D0, 1.8D0, .9D0, .53D0, .28D0, 12*0.D0, .25D0,
308 & .82D0, 1.3D0, 1.9D0, 2.8D0, 5.5D0, 8.D0, 5.7D0, 3.9D0, 2.35D0,
309 & 1.15D0, .69D0, .37D0, 18*0.D0, .34D0, 1.5D0, 3.47D0, 5.87D0,
310 & 6.23D0, 4.27D0, 2.6D0, 1.D0, .6D0, .3D0, .15D0, 6*0.D0/
311*
312 DATA SPKPX1/ 8*0.D0, .16D0, .75D0, 1.73D0, 2.93D0, 3.12D0,
313 & 2.13D0, 1.3D0, .5D0, .3D0, .15D0, .08D0, 16*0.D0, .2D0, .6D0,
314 & .92D0, 2.4D0, 4.9D0, 6.25D0, 5.25D0, 3.5D0, 2.15D0, 1.4D0, 1.D0,
315 & .7D0, 13*0.D0, .13D0, .4D0, .62D0, 1.6D0, 3.27D0, 4.17D0, 3.5D0,
316 & 2.33D0, 1.43D0, .93D0, .66D0, .47D0, 13*0.D0, .07D0, .2D0, .31D0,
317 & .8D0, 1.63D0, 2.08D0, 1.75D0, 1.17D0, .72D0, .47D0, .34D0, .23D0,
318 & 17*0.D0, .33D0, 1.D0, 1.8D0, 2.67D0, 5.33D0, 6.D0, 5.53D0, 5.D0,
319 & 17*0.D0, .17D0, .5D0, .9D0, 1.83D0, 2.67D0, 3.0D0, 2.77D0, 2.5D0,
320 & 6*0.D0, 1.D0, 3.3D0, 2.8D0, 2.5D0, 2.3D0, 1.8D0, 1.5D0, 1.1D0,
321 & .8D0, .7D0, .55D0, .3D0, 19*0.D0, .1D0, .4D0, 1.D0, 1.4D0, 2.2D0,
322 & 2.5D0, 2.2D0, 1.65D0, 1.35D0, 1.1D0, .8D0, .6D0, .4D0, 12*0.D0,
323 & .15D0, .6D0, 1.5D0, 2.1D0, 3.3D0, 3.8D0, 3.3D0, 2.45D0, 2.05D0,
324 & 1.65D0, 1.2D0, .9D0, .6D0, 12*0.D0, .10D0, .2D0, .5D0, .7D0,
325 & 1.3D0, 1.55D0, 1.9D0, 1.8D0, 1.55D0, 1.35D0, 1.15D0, .95D0, .7D0,
326 & 13*0.D0, .2D0, .5D0, .7D0, 1.3D0, 1.55D0, 1.9D0, 1.8D0, 1.55D0,
327 & 1.35D0, 1.15D0, .95D0, .7D0, 17*0.D0, .2D0, .5D0, .85D0, 2.D0,
328 & 2.15D0, 2.05D0, 1.75D0, 1.D0, 17*0.D0, .13D0, .33D0, .57D0,
329 & 1.33D0, 1.43D0, 1.36D0, 1.17D0, .67D0, 17*0.D0, .07D0, .17D0,
330 & .28D0, .67D0, .72D0, .69D0, .58D0, .33D0, 17*0.D0, .4D0, .7D0,
331 & 1.D0, 1.6D0, 1.8D0, 2.3D0, 1.9D0, 1.7D0 /
332**** pi- n data *
333 DATA SPKP41 / 0.D0, 300.D0, 40.D0, 20.D0, 13.D0, 8.5D0, 8.D0,
334 & 9.5D0, 12.D0, 14.D0, 15.5D0, 20.D0, 17.D0, 13.D0, 10.D0, 9.D0,
335 & 8.5D0, 8.D0, 7.8D0, 7.3D0, 6.7D0, 9*0.D0, .23D0, .35D0, .7D0,
336 & .52D0, .4D0, .3D0, .2D0, .15D0, .13D0, .11D0, .09D0, .07D0, 0.D0,
337 & .033D0, .8D0, 2*1.35D0, .5D0, 19*0.D0, 0.8D0, 2.2D0, 3.6D0,
338 & 4.6D0, 4.7D0, 3.5D0, 2.4D0, 1.8D0, 1.4D0, .75D0, .47D0, .25D0,
339 & .13D0, .08D0, 7*0.D0, 1.2D0, 3.3D0, 5.4D0, 6.9D0, 7.3D0, 5.3D0,
340 & 3.6D0, 2.7D0, 2.2D0, 1.1D0, .73D0, .4D0, .22D0, .12D0, 11*0.D0,
341 & 2.0D0, 4.4D0, 6.8D0, 9.9D0, 7.9D0, 6.0D0, 3.8D0, 2.5D0, 2.D0,
342 & 1.4D0, 1.D0, .6D0, .35D0, 10*0.D0, .25D0, .55D0, .75D0, 1.25D0,
343 & 1.9D0, 2.D0, 1.8D0, 1.5D0, 1.25D0, 1.D0, .8D0, 10*0.D0, .4D0,
344 & .85D0, 1.1D0, 1.85D0, 2.8D0, 3.D0, 2.7D0, 2.2D0, 1.85D0, 1.5D0,
345 & 1.2D0, 12*0.D0, .5D0, 1.2D0, 1.7D0, 3.4D0, 5.2D0, 6.4D0, 6.1D0,
346 & 5.6D0, 5.2D0, 9*0.D0, 1.D0, 3.3D0, 5.2D0, 4.45D0, 3.6D0, 2.75D0,
347 & 1.9D0, 1.65D0, 1.3D0, .95D0, .6D0, .45D0, 10*0.D0, .45D0, 1.4D0,
348 & 1.5D0, 1.1D0, .85D0, .5D0, .3D0, .2D0, .15D0, 15*0.D0, .6D0,
349 & .8D0, .95D0, .8D0, .7D0, .6D0, .5D0, .4D0, 13*0.D0, .85D0, 1.2D0,
350 & 1.4D0, 1.2D0, 1.05D0, .9D0, .7D0, .55D0, 13*0.D0, 1.D0, 1.5D0,
351 & 3.5D0, 4.15D0, 3.7D0, 2.7D0, 2.3D0, 1.75D0, 16*0.D0, .5D0, 2.0D0,
352 & 3.3D0, 5.4D0, 7.D0 /
353**** k+ p data *
354 DATA SPKP51/ 0.D0, 20.D0, 14.D0, 12.D0, 11.5D0, 10.D0, 8.D0,
355 & 7.D0, 6.D0, 5.5D0, 5.3D0, 5.D0, 4.5D0, 4.4D0, 3.8D0, 3.D0, 2.8D0,
356 & 0.D0, .5D0, 1.15D0, 2.D0, 1.3D0, .8D0, .45D0, 13*0.D0, 0.9D0,
357 & 2.5D0, 3.D0, 2.5D0, 2.3D0, 2.D0, 1.7D0, 1.5D0, 1.2D0, .9D0, .6D0,
358 & .45D0, .21D0, .2D0, 3*0.D0, .9D0, 2.5D0, 3.D0, 2.5D0, 2.3D0,
359 & 2.D0, 1.7D0, 1.5D0, 1.2D0, .9D0, .6D0, .45D0, .21D0, .2D0,
360 & 4*0.D0, 1.D0, 2.1D0, 2.6D0, 2.3D0, 2.1D0, 1.8D0, 1.7D0, 1.4D0,
361 & 1.2D0, 1.05D0, .9D0, .66D0, .5D0, 7*0.D0, .3D0, 2*1.D0, .9D0,
362 & .7D0, .4D0, .3D0, .2D0, 11*0.D0, .1D0, 1.D0, 2.2D0, 3.5D0, 4.2D0,
363 & 4.55D0, 4.85D0, 4.9D0, 10*0.D0, .2D0, .7D0, 1.6D0, 2.5D0, 2.2D0,
364 & 1.71D0, 1.6D0, 6*0.D0, 1.4D0, 3.8D0, 5.D0, 4.7D0, 4.4D0, 4.D0,
365 & 3.5D0, 2.85D0, 2.35D0, 2.01D0, 1.8D0, 12*0.D0, .1D0, .8D0,2.05D0,
366 & 3.31D0, 3.5D0, 12*0.D0, .034D0, .2D0, .75D0, 1.04D0, 1.24D0 /
367**** k+ n data *
368 DATA SPKP61/ 0.D0, 6.D0, 11.D0, 13.D0, 6.D0, 5.D0, 3.D0, 2.2D0,
369 & 1.5D0, 1.2D0, 1.D0, .7D0, .6D0, .5D0, .45D0, .35D0, .3D0, 0.D0,
370 & 6.D0, 11.D0, 13.D0, 6.D0, 5.D0, 3.D0, 2.2D0, 1.5D0, 1.2D0, 1.D0,
371 & .7D0, .6D0, .5D0, .45D0, .35D0, .3D0, 0.D0, .5D0, 1.3D0, 2.8D0,
372 & 2.3D0, 1.6D0, .9D0, 13*0.D0, 0.9D0, 2.5D0, 3.D0, 2.5D0, 2.3D0,
373 & 2.D0, 1.7D0, 1.5D0,1.2D0,.9D0,.6D0,.45D0,.21D0,.2D0,3*0.D0,0.9D0,
374 & 2.5D0, 3.D0, 2.5D0, 2.3D0,2.D0,1.7D0,1.5D0,1.2D0,.9D0,.6D0,.45D0,
375 & .21D0, .2D0,4*0.D0,1.D0,2.1D0,2.6D0,2.3D0,2.D0,1.8D0,1.7D0,1.4D0,
376 & 1.2D0,1.15D0,.9D0,.66D0,.5D0,4*0.D0,1.D0,2.1D0,2.6D0,2.3D0,2.1D0,
377 & 1.8D0,1.7D0,1.4D0,1.2D0, 1.15D0, .9D0, .66D0, .5D0, 7*0.D0, .3D0,
378 & 2*1.D0, .9D0, .7D0, .4D0, .35D0, .2D0, 9*0.D0, .3D0, 2*1.D0,.9D0,
379 & .7D0, .4D0, .35D0, .2D0, 11*0.D0, .1D0, 1.D0, 2.4D0,3.5D0,4.25D0,
380 & 4.55D0, 4.85D0, 4.9D0, 9*0.D0, .1D0, 1.D0, 2.4D0, 3.5D0, 4.25D0,
381 & 4.55D0, 4.85D0, 4.9D0, 10*0.D0, .2D0, .7D0, 1.6D0, 2.5D0, 2.2D0,
382 & 1.71D0, 1.6D0, 10*0.D0, .2D0, .7D0, 1.6D0, 2.5D0, 2.2D0, 1.71D0,
383 & 1.6D0, 6*0.D0, 1.4D0, 3.8D0, 5.D0, 4.7D0,4.4D0,4.D0,3.5D0,2.85D0,
384 & 2.35D0, 2.01D0, 1.8D0, 6*0.D0, 1.4D0,3.8D0,5.D0,4.7D0,4.4D0,4.D0,
385 & 3.5D0,2.85D0,2.35D0,2.01D0,1.8D0,12*0.D0,.1D0,.8D0,2.05D0,3.31D0,
386 & 3.5D0, 12*0.D0, .034D0,.2D0,.75D0,1.04D0,1.24D0 /
387**** k- p data *
388 DATA SKMPL1/ 0.D0, 35.D0, 22.D0, 25.D0, 17.D0, 9.D0, 9.5D0, 8.D0,
389 & 7.D0, 6.5D0, 6.1D0, 5.D0, 4.8D0, 4.6D0, 4.45D0, 4.3D0, 4.2D0,
390 & 0.D0, 8.D0, 3.5D0, 8.D0, 3.D0, 1.9D0, 1.7D0, 1.D0, .9D0, .8D0,
391 & .75D0, .5D0, .42D0, .38D0, .34D0, .25D0, .2D0,
392 & 0.D0, 3.D0, 3.2D0, 3.5D0, 1.5D0, 1.4D0, 1.1D0, .6D0, .5D0,
393 & .35D0, .28D0, .25D0, .18D0, .12D0, .1D0, .08D0, .04D0,
394 & 0.D0, 8.5D0, 2.4D0, 1.7D0, 1.3D0, 1.3D0, 1.1D0, .5D0,
395 & .4D0, .4D0, .35D0, .3D0, .28D0, .2D0, .16D0, .13D0, .11D0,
396 & 0.D0, 7.D0, 4.8D0, 1.4D0, 1.9D0, .9D0, .4D0, .2D0, .13D0,
397 & .1D0, .08D0, .06D0, .04D0, .02D0, .015D0, .01D0, .01D0,
398 & 0.D0, 5.5D0, 1.D0, .8D0, .75D0, .32D0, .2D0, .1D0, .09D0,
399 & .08D0, .065D0, .05D0, .04D0, .022D0, .017D0, 2*.01D0/
400 DATA SPKP71 / 0.D0, .56D0, 1.46D0, 3.16D0, 2.01D0, 1.28D0, .74D0,
401 & 14*0.D0, 1.13D0, 2.61D0, 2.91D0, 2.58D0, 2.35D0, 2.02D0,
402 & 1.91D0, 1.57D0, 1.35D0, 1.29D0, 1.01D0, .74D0, .65D0, 4*0.D0,
403 & 1.13D0, 2.61D0, 2.91D0, 2.58D0, 2.35D0, 2.02D0, 1.91D0, 1.57D0,
404 & 1.35D0, 1.29D0, 1.01D0, .74D0, .65D0, 3*0.D0, 1.0D0, 3.03D0,
405 & 3.36D0, 2.8D0, 2.58D0, 2.24D0, 1.91D0, 1.68D0, 1.35D0, 1.01D0,
406 & .67D0, .5D0, .24D0, .23D0, 3*0.D0, 1.0D0, 3.03D0, 3.36D0, 2.8D0,
407 & 2.58D0, 2.24D0, 1.91D0, 1.68D0, 1.35D0, 1.01D0, .67D0, .5D0,
408 & .24D0, .23D0, 7*0.D0, .34D0, 1.12D0, 1.12D0, 1.01D0, .78D0,
409 & .45D0, .39D0, .22D0, .07D0, 0.D0, 7*0.D0, .34D0, 1.12D0, 1.12D0,
410 & 1.01D0, .78D0, .45D0, .39D0, .22D0, .07D0, 0.D0, 6*0.D0, 1.71D0,
411 & 4.26D0, 5.6D0, 5.57D0, 4.93D0, 4.48D0, 3.92D0, 3.19D0, 2.63D0,
412 & 2.25D0, 2.D0, 6*0.D0, 1.71D0, 4.26D0, 5.6D0, 5.57D0, 4.93D0,
413 & 4.48D0, 3.92D0, 3.19D0, 2.63D0, 2.25D0, 2.D0, 10*0.D0, .22D0,
414 & .8D0, .75D0, 1.D0, 1.3D0, 1.5D0, 1.3D0, 10*0.D0, .22D0, .8D0,
415 & .75D0, 1.D0, 1.3D0, 1.5D0, 1.3D0, 13*0.D0, .1D0, .3D0, .7D0,1.D0,
416 & 13*0.D0, .1D0, .3D0, .7D0, 1.D0, 9*0.D0, .11D0, 1.72D0, 2.69D0,
417 & 3.92D0, 4.76D0, 5.10D0, 5.44D0, 5.3D0, 9*0.D0, .11D0, 1.72D0,
418 & 2.69D0, 3.92D0, 4.76D0, 5.1D0, 5.44D0, 5.3D0, 5*0.D0,9.2D0,4.7D0,
419 & 1.9D0, 10*0.D0, 2.5D0, 15.D0, 21.5D0, 15.3D0, 3.D0,1.5D0,10*0.D0/
420***** k- n data *
421 DATA SKMNL1/0.D0, 4.D0, 9.5D0, 20.D0, 13.D0, 9.5D0, 6.D0, 4.4D0,
422 & 3.D0, 2.4D0, 2.D0, 1.4D0, 1.2D0, 1.D0, .9D0, .7D0, .6D0,
423 & 0.D0, 4.5D0, 6.D0, 5.D0, 2.5D0, 2.D0, 1.7D0, 2.1D0,
424 & 1.9D0, .9D0, .5D0, .3D0, .24D0, .2D0, .18D0, .1D0, .09D0,
425 & 0.D0, 1.8D0, 2.D0, 1.1D0, .9D0, .5D0, .5D0, .4D0, .4D0,
426 & .2D0, .1D0, .06D0, .05D0, .04D0, .03D0, .02D0, .02D0,
427 & 0.D0, 1.5D0, 2.D0, .9D0, 1.1D0, .4D0, .6D0, .7D0, .65D0,
428 & .3D0, .17D0, .1D0, .08D0, .07D0, .06D0, .04D0, .03D0/
429 DATA SPKP81/0.D0, .56D0, 1.29D0, 2.26D0, 1.01D0, .64D0, .37D0,
430 & 14*0.D0, 1.13D0, 2.61D0, 2.91D0, 2.58D0, 2.35D0, 2.02D0,
431 & 1.91D0, 1.57D0, 1.35D0, 1.29D0, 1.01D0, .74D0, .65D0,
432 & 3*0.D0, 1.D0, 3.03D0, 3.36D0, 2.8D0, 2.58D0, 2.24D0,
433 & 1.91D0, 1.68D0, 1.35D0, 1.01D0, .67D0, .5D0, .24D0, .23D0,
434 & 3*0.D0, 1.D0, 3.03D0, 3.36D0, 2.8D0, 2.58D0, 2.24D0,
435 & 1.91D0, 1.68D0, 1.35D0, 1.01D0, .67D0, .5D0, .24D0, .23D0,
436 & 7*0.D0, .34D0, 1.12D0, 1.12D0, 1.01D0, .78D0, .45D0,
437 & .39D0, .22D0, .07D0, 0.D0,
438 & 6*0.D0, 1.71D0, 4.26D0, 5.6D0, 5.57D0, 4.93D0,
439 & 4.48D0, 3.92D0, 3.19D0, 2.63D0, 2.25D0, 2.D0,
440 & 10*0.D0, .22D0, .8D0, .75D0, 1.D0, 1.3D0, 1.5D0, 1.3D0,
441 & 13*0.D0, .1D0, .3D0, .7D0, 1.D0,
442 & 13*0.D0, .1D0, .3D0, .7D0, 1.D0,
443 & 9*0.D0, .11D0, 1.72D0, 2.69D0, 3.92D0, 4.76D0,
444 & 5.10D0, 5.44D0, 5.3D0,
445 & 4*0.D0, 0.00D0, 9.2D0, 4.7D0, 1.9D0, 9*0.D0/
446***** p p data *
447 DATA SPKP91/ 0.D0, 24.D0, 25.D0, 27.D0, 23.D0, 21.D0, 20.D0,
448 & 19.D0, 17.D0, 15.5D0, 14.D0, 13.5D0, 13.D0,
449 & 0.D0, 3.6D0, 1.7D0, 10*0.D0,
450 & .0D0, 0.D0, 8.7D0, 17.7D0, 18.8D0, 15.9D0,
451 & 11.7D0, 8.D0, 6.D0, 5.3D0, 4.5D0, 3.9D0, 3.5D0,
452 & .0D0, .0D0, 2.8D0, 5.8D0, 6.2D0, 5.1D0, 3.8D0,
453 & 2.7D0, 2.1D0, 1.8D0, 1.5D0, 1.3D0, 1.1D0,
454 & 5*0.D0, 4.6D0, 10.2D0, 15.1D0,
455 & 16.9D0, 16.5D0, 11.D0, 5.5D0, 3.5D0,
456 & 10*0.D0, 4.3D0, 7.6D0, 9.D0,
457 & 10*0.D0, 1.7D0, 2.6D0, 3.D0,
458 & 6*0.D0, .3D0, .6D0, 1.D0, 1.6D0, 1.3D0, .8D0, .6D0,
459 & 6*0.D0, .7D0, 1.2D0, 1.8D0, 2.5D0, 1.8D0, 1.3D0,
460 & 1.2D0, 10*0.D0, .6D0, 1.4D0, 1.7D0,
461 & 10*0.D0, 1.9D0, 4.1D0, 5.2D0/
462***** p n data *
463 DATA SPKP01/ 0.D0, 24.D0, 25.D0, 27.D0, 23.D0, 21.D0, 20.D0,
464 & 19.D0, 17.D0, 15.5D0, 14.D0, 13.5D0, 13.D0,
465 & 0.D0, 1.8D0, .2D0, 12*0.D0,
466 & 3.2D0, 6.05D0, 9.9D0, 5.1D0,
467 & 3.8D0, 2.7D0, 1.9D0, 1.5D0, 1.4D0, 1.3D0, 1.1D0,
468 & 2*.0D0, 3.2D0, 6.05D0, 9.9D0, 5.1D0,
469 & 3.8D0, 2.7D0, 1.9D0, 1.5D0, 1.4D0, 1.3D0, 1.1D0,
470 & 5*0.D0, 4.6D0, 10.2D0, 15.1D0,
471 & 16.4D0, 15.2D0, 11.D0, 5.4D0, 3.5D0,
472 & 5*0.D0, 4.6D0, 10.2D0, 15.1D0,
473 & 16.4D0, 15.2D0, 11.D0, 5.4D0, 3.5D0,
474 & 10*0.D0, .7D0, 5.1D0, 8.D0,
475 & 10*0.D0, .7D0, 5.1D0, 8.D0,
476 & 10*.0D0, .3D0, 2.8D0, 4.7D0,
477 & 10*.0D0, .3D0, 2.8D0, 4.7D0,
478 & 7*0.D0, 1.2D0, 2.5D0, 3.5D0, 6.D0, 5.3D0, 2.9D0,
479 & 7*0.D0, 1.7D0, 3.6D0, 5.4D0, 9.D0, 7.6D0, 4.2D0,
480 & 5*0.D0, 7.7D0, 6.1D0, 2.9D0, 5*0.D0/
481* nn - data *
482* *
483 DATA SPKPV1/ 0.D0, 24.D0, 25.D0, 27.D0, 23.D0, 21.D0, 20.D0,
484 & 19.D0, 17.D0, 15.5D0, 14.D0, 13.5D0, 13.D0,
485 & 0.D0, 3.6D0, 1.7D0, 12*0.D0,
486 & 8.7D0, 17.7D0, 18.8D0, 15.9D0,
487 & 11.7D0, 8.D0, 6.D0, 5.3D0, 4.5D0, 3.9D0, 3.5D0,
488 & .0D0, .0D0, 2.8D0, 5.8D0, 6.2D0, 5.1D0, 3.8D0,
489 & 2.7D0, 2.1D0, 1.8D0, 1.5D0, 1.3D0, 1.1D0,
490 & 5*0.D0, 4.6D0, 10.2D0, 15.1D0, 16.9D0, 16.5D0,
491 & 11.D0, 5.5D0, 3.5D0,
492 & 10*0.D0, 4.3D0, 7.6D0, 9.D0,
493 & 10*0.D0, 1.7D0, 2.6D0, 3.D0,
494 & 6*0.D0, .3D0, .6D0, 1.D0, 1.6D0, 1.3D0, .8D0, .6D0,
495 & 6*0.D0, .7D0, 1.2D0, 1.8D0, 2.5D0, 1.8D0, 1.3D0,
496 & 1.2D0, 10*0.D0, .6D0, 1.4D0, 1.7D0,
497 & 10*0.D0, 1.9D0, 4.1D0, 5.2D0/
498**************** ap - p - data *
499 DATA SAPPL1/ 0.D0, 176.D0, 160.D0, 105.D0, 75.D0, 68.D0, 65.D0,
500 & 50.D0, 50.D0, 43.D0, 42.D0, 40.5D0, 35.D0, 30.D0, 28.D0,
501 & 25.D0, 22.D0, 21.D0, 20.D0, 18.D0, 17.D0, 11*0.D0,
502 & .05D0, .15D0, .18D0, .2D0, .2D0, .3D0, .4D0, .6D0, .7D0, .85D0,
503 & 0.D0, 1.D0, .9D0, .46D0, .3D0, .23D0, .18D0, .16D0, .14D0,
504 & .1D0, .08D0, .05D0, .02D0, .015D0, 4*.011D0, 3*.005D0,
505 & 0.D0, 55.D0, 50.D0, 25.D0, 15.D0, 15.D0, 14.D0, 12.D0,
506 & 10.D0, 7.D0, 6.D0, 4.D0, 3.3D0, 2.8D0, 2.4D0, 2.D0, 1.8D0,
507 & 1.55D0, 1.3D0, .95D0, .75D0,
508 & 0.D0, 3.3D0, 3.D0, 1.5D0, 1.D0, .7D0, .4D0, .35D0, .4D0,
509 & .25D0, .18D0, .08D0, .04D0, .03D0, .023D0, .016D0, .014D0,
510 & .01D0, .008D0, .006D0, .005D0/
511 DATA SPKPE1/0.D0, 215.D0, 193.D0, 170.D0, 148.D0, 113.D0, 97.D0,
512 & 84.D0, 78.D0, 68.D0, 64.D0, 61.D0, 46.D0, 36.D0, 31.3D0, 28.5D0,
513 & 25.7D0, 22.6D0, 21.4D0, 20.7D0, 19.9D0,
514 & 9*0.D0, 2.D0, 2.5D0, .2D0, 19*0.D0, .3D0, 1.4D0, 2.2D0, 1.2D0,
515 & 1.1D0, 1.D0, .8D0, .6D0, .5D0, .4D0, .3D0, 10*0.D0, .3D0, 1.4D0,
516 & 2.2D0, 1.2D0, 1.1D0, 1.D0, .8D0, .6D0, .5D0, .4D0, .3D0, 10*0.D0,
517 & .3D0, 1.4D0, 2.2D0, 1.2D0, 1.1D0, 1.D0, .8D0, .6D0, .5D0, .4D0,
518 & .3D0, 10*0.D0, .3D0, 1.4D0, 2.2D0, 1.2D0, 1.1D0, 1.D0, .8D0,
519 & .6D0, .5D0, .4D0, .3D0, 9*0.D0, .6D0, 2.5D0, 5.D0, 5.2D0, 5.1D0,
520 & 5.4D0, 5.8D0, 2.8D0, 2.1D0, 1.8D0, 1.6D0, 1.2D0, 13*0.D0, 1.3D0,
521 & 1.5D0, 2.D0, 2.5D0, 2.5D0, 2.3D0, 1.8D0, 1.4D0, 13*0.D0, 1.3D0,
522 & 1.5D0, 2.D0, 2.5D0, 2.5D0, 2.3D0, 1.8D0, 1.4D0, 13*0.D0, 1.3D0,
523 & 1.5D0, 2.D0, 2.5D0, 2.5D0, 2.3D0, 1.8D0, 1.4D0, 13*0.D0, 1.3D0,
524 & 1.5D0, 2.D0, 2.5D0, 2.5D0, 2.3D0, 1.8D0, 1.4D0, 14*0.D0, .2D0,
525 & .5D0, 1.1D0, 1.6D0, 1.4D0, 1.1D0, .9D0, 14*0.D0, .2D0, .5D0,
526 & 1.1D0, 1.6D0, 1.4D0, 1.1D0, .9D0, 14*0.D0, .2D0, .5D0, 1.1D0,
527 & 1.6D0, 1.4D0, 1.1D0, .9D0, 14*0.D0, .2D0, .5D0, 1.1D0, 1.6D0,
528 & 1.4D0, 1.1D0, .9D0, 17*0.D0, .3D0, 1.6D0, 2.6D0, 3.6D0, 17*0.D0,
529 & .3D0, 1.6D0, 2.6D0, 3.6D0, 17*0.D0, .3D0, 1.6D0, 2.6D0,
530 & 3.6D0, 17*0.D0, .3D0, 1.6D0, 2.6D0, 3.6D0 /
531**************** ap - n - data *
532 DATA SAPNL1/
533 & 0.D0, 176.D0, 160.D0, 105.D0, 75.D0, 68.D0, 65.D0,
534 & 50.D0, 50.D0, 43.D0, 42.D0, 40.5D0, 35.D0, 30.D0, 28.D0,
535 & 25.D0, 22.D0, 21.D0, 20.D0, 18.D0, 17.D0, 11*0.D0,
536 & .05D0, .15D0, .18D0, .2D0, .2D0, .3D0, .4D0, .6D0, .7D0,
537 & .85D0, 0.D0, 1.D0, .9D0, .46D0, .3D0, .23D0, .18D0, .16D0,
538 & .14D0, .1D0, .08D0, .05D0, .02D0, .015D0, 4*.011D0, 3*.005D0,
539 & 0.D0, 3.3D0, 3.D0, 1.5D0, 1.D0, .7D0, .4D0, .35D0, .4D0,
540 & .25D0, .18D0, .08D0, .04D0, .03D0, .023D0, .016D0, .014D0,
541 & .01D0, .008D0, .006D0, .005D0 /
542 DATA SPKPZ1/ 0.D0, 215.D0, 193.D0, 170.D0, 148.D0, 113.D0, 97.D0,
543 & 84.D0, 78.D0, 68.D0, 64.D0, 61.D0, 46.D0, 36.D0, 31.3D0, 28.5D0,
544 & 25.7D0, 22.6D0, 21.4D0, 20.7D0, 19.9D0, 9*0.D0, 2.4D0, .2D0,
545 & 20*0.D0, 1.8D0, 2.8D0, 3.6D0, 2.3D0, 1.8D0, 1.5D0, 1.3D0, 1.D0,
546 & .7D0, .5D0, .3D0, 10*0.D0, 1.8D0, 2.8D0, 3.6D0, 2.3D0, 1.8D0,
547 & 1.5D0, 1.3D0, 1.D0, .7D0, .5D0, .3D0, 10*0.D0, 1.8D0, 2.8D0,
548 & 3.6D0, 2.3D0, 1.8D0, 1.5D0, 1.3D0, 1.D0, .7D0, .5D0, .3D0,
549 & 10*0.D0, 1.8D0, 2.8D0, 3.6D0, 2.3D0, 1.8D0, 1.5D0, 1.3D0, 1.D0,
550 & .7D0, .5D0, .3D0, 13*0.D0, 5.2D0, 8.7D0, 11.4D0, 14.D0, 11.9D0,
551 & 7.6D0, 6.D0, 5.D0, 13*0.D0, 5.2D0, 8.7D0, 11.4D0, 14.D0, 11.9D0,
552 & 7.6D0, 6.D0, 5.D0, 18*0.D0, 1.D0, 4.9D0, 8.5D0, 18*0.D0, 1.D0,
553 & 4.9D0, 8.5D0, 15*0.D0, 1.9D0, 2.3D0, 4.D0, 6.5D0, 5.2D0, 3.4D0,
554 & 15*0.D0, 1.9D0, 2.3D0, 4.D0, 6.5D0, 5.2D0, 3.4D0, 15*0.D0, 1.9D0,
555 & 2.3D0, 4.D0, 6.5D0, 5.2D0, 3.4D0 /
556* *
557* *
558**************** an - p - data *
559* *
560 DATA SANPL1/
561 & 0.D0, 176.D0, 160.D0, 105.D0, 75.D0, 68.D0, 65.D0, 50.D0,
562 & 50.D0, 43.D0, 42.D0, 40.5D0, 35.D0, 30.D0, 28.D0,
563 & 25.D0, 22.D0, 21.D0, 20.D0, 18.D0, 17.D0, 11*0.D0, .05D0,
564 & .15D0, .18D0, .2D0, .2D0, .3D0, .4D0, .6D0, .7D0, .85D0,
565 & 0.D0, 1.D0, .9D0, .46D0, .3D0, .23D0, .18D0, .16D0, .14D0,
566 & .1D0, .08D0, .05D0, .02D0, .015D0, 4*.011D0, 3*.005D0,
567 & 0.D0, 3.3D0, 3.D0, 1.5D0, 1.D0, .7D0, .4D0, .35D0, .4D0, .25D0,
568 & .18D0, .08D0, .04D0, .03D0, .023D0, .016D0, .014D0,
569 & .01D0, .008D0, .006D0, .005D0 /
570 DATA SPKPF1/ 0.D0, 215.D0, 193.D0, 170.D0, 148.D0, 113.D0, 97.D0,
571 & 84.D0, 78.D0, 68.D0, 64.D0, 61.D0, 46.D0, 36.D0, 31.3D0, 28.5D0,
572 & 25.7D0, 22.6D0, 21.4D0, 20.7D0, 19.9D0, 9*0.D0, 2.4D0, .2D0,
573 & 20*0.D0, 1.8D0, 2.8D0, 3.6D0, 2.3D0, 1.8D0, 1.5D0, 1.3D0, 1.D0,
574 & .7D0, .5D0, .3D0, 10*0.D0, 1.8D0, 2.8D0, 3.6D0, 2.3D0, 1.8D0,
575 & 1.5D0, 1.3D0, 1.D0, .7D0, .5D0, .3D0, 10*0.D0, 1.8D0, 2.8D0,
576 & 3.6D0, 2.3D0, 1.8D0, 1.5D0, 1.3D0, 1.D0, .7D0, .5D0, .3D0,
577 & 10*0.D0, 1.8D0, 2.8D0, 3.6D0, 2.3D0, 1.8D0, 1.5D0, 1.3D0, 1.D0,
578 & .7D0, .5D0, .3D0, 13*0.D0, 5.2D0, 8.7D0, 11.4D0, 14.D0, 11.9D0,
579 & 7.6D0, 6.D0, 5.D0, 13*0.D0, 5.2D0, 8.7D0, 11.4D0, 14.D0, 11.9D0,
580 & 7.6D0, 6.D0, 5.D0, 18*0.D0, 1.D0, 4.9D0, 8.5D0, 18*0.D0, 1.D0,
581 & 4.9D0, 8.5D0, 15*0.D0, 1.9D0, 2.3D0, 4.D0, 6.5D0, 5.2D0, 3.4D0,
582 & 15*0.D0, 1.9D0, 2.3D0, 4.D0, 6.5D0, 5.2D0, 3.4D0, 15*0.D0, 1.9D0,
583 & 2.3D0, 4.D0, 6.5D0, 5.2D0, 3.4D0 /
584**** ko - n - data *
585 DATA CPKP51/0.D0, 20.D0, 14.D0, 12.D0, 11.5D0, 10.D0, 8.D0, 7.D0,
586 & 6.D0, 5.5D0, 5.3D0, 5.D0, 4.5D0, 4.4D0, 3.8D0, 3.D0, 2.8D0,
587 & 0.D0, .5D0, 1.15D0, 2.D0, 1.3D0, .8D0, .45D0, 10*0.D0,
588 & 3*0.D0, 0.9D0, 2.5D0, 3.D0, 2.5D0, 2.3D0, 2.D0, 1.7D0,
589 & 1.5D0, 1.2D0, .9D0, .6D0, .45D0, .21D0, .2D0,
590 & 3*0.D0, 0.9D0, 2.5D0, 3.D0, 2.5D0, 2.3D0, 2.D0, 1.7D0,
591 & 1.5D0, 1.2D0, .9D0, .6D0, .45D0, .21D0, .2D0,
592 & 4*0.D0, 1.D0, 2.1D0, 2.6D0, 2.3D0, 2.1D0, 1.8D0, 1.7D0,
593 & 1.4D0, 1.2D0, 1.05D0, .9D0, .66D0, .5D0,
594 & 7*0.D0, .3D0, 1.D0, 1.D0, .9D0, .7D0, .4D0, .30D0, .2D0,
595 & 11*0.D0, .1D0, 1.D0, 2.2D0, 3.5D0, 4.20D0, 4.55D0,
596 & 4.85D0, 4.9D0,
597 & 10*0.D0, .2D0, .7D0, 1.6D0, 2.5D0, 2.2D0, 1.71D0, 1.6D0,
598 & 6*0.D0, 1.4D0, 3.8D0, 5.D0, 4.7D0, 4.4D0, 4.D0, 3.5D0,
599 & 2.85D0, 2.35D0, 2.01D0, 1.8D0,
600 & 12*0.D0, .1D0, .8D0, 2.05D0, 3.31D0, 3.5D0,
601 & 12*0.D0, .034D0, .20D0, .75D0, 1.04D0, 1.24D0 /
602**** ako - p - data *
603 DATA CPKP61/ 0.D0,4.D0,9.5D0,20.D0,13.D0,9.5D0,6.D0,4.4D0,
604 & 3.D0,2.4D0,2.D0,1.4D0,1.2D0,1.D0,.9D0,.7D0,.6D0,0.D0,
605 & 4.5D0,6.D0,5.D0,2.5D0,2.D0,1.7D0,2.1D0,1.9D0,.9D0,.5D0,
606 & .3D0,.24D0,.2D0,.18D0,.1D0,.09D0,0.D0,1.8D0,2.D0,1.1D0,
607 & .9D0,.5D0,.5D0,.4D0,.4D0,.2D0,.1D0,.06D0,.05D0,.04D0,
608 & .03D0,.02D0,.02D0,0.D0,1.5D0,2.D0,.9D0,1.1D0,.4D0,.6D0,
609 & .7D0,.65D0,.3D0,.17D0,.1D0,.08D0,.07D0,.06D0,.04D0,
610 & .03D0,0.D0,.56D0,1.29D0,2.26D0,1.01D0,.64D0,.37D0,
611 & 14*0.D0,1.13D0,2.61D0,2.91D0,2.58D0,2.35D0,2.02D0,1.91D0,
612 & 1.57D0,1.35D0,1.29D0,1.01D0,.74D0,.65D0,3*0.D0,1.0D0,
613 & 3.03D0,3.36D0,2.8D0,2.58D0,2.24D0,1.91D0,1.68D0,1.35D0,
614 & 1.01D0,.67D0,.5D0,.24D0,.23D0,3*0.D0,1.0D0,3.03D0,3.36D0,
615 & 2.8D0,2.58D0,2.24D0,1.91D0,1.68D0,1.35D0,1.01D0,.67D0,
616 & .5D0,.24D0,.23D0,7*0.D0,.34D0,1.12D0,1.12D0,1.01D0,.78D0,
617 & .45D0,.39D0,.22D0,.07D0,7*0D0,1.71D0,4.26D0,5.6D0,5.57D0,4.93D0,
618 & 4.48D0,3.92D0,3.19D0,2.63D0,2.25D0,2.D0,10*0.D0,.22D0,.8D0,
619 & .75D0,1.D0,1.3D0,1.5D0,1.3D0,13*0.D0,.1D0,.3D0,.7D0,1.D0,
620 & 13*0.D0,.1D0,.3D0,.7D0,1.D0,9*0.D0,.11D0,1.72D0,2.69D0,3.92D0,
621 & 4.76D0,5.10D0,5.44D0,5.3D0,5*0.D0,9.2D0,4.7D0,1.9D0,9*0.D0,
622 & .0D0,2.5D0,15.D0,21.5D0,15.3D0,3.D0,1.5D0,10*0.D0 /
623 DATA NURE1/9, 12, 5*0, 10, 14, 3*0, 1, 3, 5, 7, 6*0, 2, 6, 16,
624 & 5*0, 10, 13, 5*0, 11, 12, 3*0, 2, 4, 6, 8, 6*0, 3, 15, 7, 5*0 /
625*
626 DO 10 J=1,92
627 UMOPI(J)=UMOPI1(J)
628 PLAPI(J)=PLAPI1(J)
629 10 CONTINUE
630 DO 20 J=1,68
631 UMOKC(J)=UMOKC1(J)
632 PLAKC(J)=PLAKC1(J)
633 20 CONTINUE
634 DO 30 J=1,39
635 UMOP(J)=UMOP1(J)
636 PLAP(J)=PLAP1(J)
637 30 CONTINUE
638 DO 40 J=1,63
639 UMON(J)=UMON1(J)
640 PLAN(J)=PLAN1(J)
641 40 CONTINUE
642 DO 50 J=1,34
643 UMOK0(J)=UMOK01(J)
644 PLAK0(J)=PLAK01(J)
645 50 CONTINUE
646 DO 60 J=1,296
647 SIIN(J)=SIIN1(J)
648 60 CONTINUE
649 DO 70 J=1,315
650 SPIKP1(J)=SPKP11(J)
651 SPIKP4(J)=SPKP41(J)
652 70 CONTINUE
653 DO 80 J=1,278
654 SPIKPU(J)=SPKPU1(J)
655 SPIKPW(J)=SPKPW1(J)
656 80 CONTINUE
657 DO 90 J=1,372
658 SPIKPV(J)=CPKPV1(J)
659 SPIKPX(J)=SPKPX1(J)
660 90 CONTINUE
661 DO 100 J=1,187
662 SPIKP5(J)=SPKP51(J)
663 SPIKP8(J)=SPKP81(J)
664 SPKP15(J)=CPKP51(J)
665 100 CONTINUE
666 DO 110 J=1,289
667 SPIKP6(J)=SPKP61(J)
668 110 CONTINUE
669 DO 120 J=1,102
670 SKMPEL(J)=SKMPL1(J)
671 120 CONTINUE
672 DO 130 J=1,289
673 SPIKP7(J)=SPKP71(J)
674 130 CONTINUE
675 DO 140 J=1,68
676 SKMNEL(J)=SKMNL1(J)
677 140 CONTINUE
678 DO 150 J=1,143
679 SPIKP9(J)=SPKP91(J)
680 SPKPV(J)=SPKPV1(J)
681 150 CONTINUE
682 DO 160 J=1,169
683 SPIKP0(J)=SPKP01(J)
684 160 CONTINUE
685 DO 170 J=1,105
686 SAPPEL(J)=SAPPL1(J)
687 170 CONTINUE
688 DO 180 J=1,399
689 SPIKPE(J)=SPKPE1(J)
690 180 CONTINUE
691 DO 190 J=1,84
692 SAPNEL(J)=SAPNL1(J)
693 SANPEL(J)=SANPL1(J)
694 190 CONTINUE
695 DO 200 J=1,273
696 SPIKPZ(J)=SPKPZ1(J)
697 SPIKPF(J)=SPKPF1(J)
698 200 CONTINUE
699 DO 210 J=1,272
700 SPKP16(J)=CPKP61(J)
701 210 CONTINUE
702 DO 220 J=1,164
703 NRKPI(J)=NRKPI1(J)
704 220 CONTINUE
705 DO 230 J=1,132
706 NRKKC(J)=NRKKC1(J)
707 230 CONTINUE
708 DO 240 J=1,70
709 NRKP(J)=NRKP1(J)
710 240 CONTINUE
711 DO 250 J=1,116
712 NRKN(J)=NRKN1(J)
713 250 CONTINUE
714 DO 270 J=1,54
715 NRKK0(J)=NRKK01(J)
716 270 CONTINUE
717 DO 280 J=1,60
718 NURELN(J)=NURE1(J)
719 280 CONTINUE
720
721
722*= end*block.blkdt3 *
723 END