COSMOS v7.655  COSMOSv7655
(AirShowerMC)
jamdat.f
Go to the documentation of this file.
1 !***********************************************************************
2 !***********************************************************************
3 ! *
4 ! PART 1: *
5 ! *
6 ! List of subprograms in rough order of relevance with main purpose *
7 ! (S = subroutine, F = function, B = block data, E = entry) *
8 ! *
9 ! b jamdata to give default parameters for jam *
10 ! b jydata to contain all default values for pythia *
11 ! *
12 !***********************************************************************
13 !***********************************************************************
14 
15  block data jamdata
16 
17 !...Purpose: to give default values to switches and parameters
18  implicit double precision(a-h, o-z)
19 #include "jam2.inc"
20 
21 !...Set defaults in the mstc array
22 !... 1 2 3 4 5 6 7 8 9 10
23  data (mstc(i),i=1,100)/
24  $ 19780503, 1, 1, 2, 1, 0, 0, 0, 0, 1,
25  1 1, 30, 1, 30, 0, 0, 0, 0, 27, 0,
26  2 0, 9*0,
27  3 0, 0, 0, 0, 0, 1, 2, 3, 4, 8,
28  4 1, 1, 1, 2, 0, 1, 4*0,
29  5 3, 5, 1, 1, 0, 1, 4*0,
30  6 3, 22, 3, 3, 1, 0, 2, 2, 2, 2,
31  7 1, 3, 2, 1, 1, 0, 0, 0, 0, 0,
32  8 1, 3, 2, 0, 2, 0, 0, 0, 0, 0,
33  9 10*0/
34 
35 !... 1 2 3 4 5 6 7 8 9 10
36  data (mstc(i),i=101,200)/
37  $ 0,1,8*0,
38  1 10*0,
39  2 10*0,
40  3 10*0,
41  4 10*0,
42  5 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
43  6 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
44  7 20*0,
45  9 1, 090, 2005,08,14,5*0/
46 
47 
48 !...Set defaults in the parc array
49 !.. 1 2 3 4 5 6 7 8 9 10
50  data (parc(i),i=1,100)/
51  $ -1.d0,100.0d0, 0.0d0, 0.0d0, -1.d0, 2.d0, 1.0d0, 1.0d0,
52  & 0.0d0, 0.0d0,
53  1 0.54d0,1.124d0, 0.8d0, 0.8d0, 2.0d0, 0.0d0, 0.0d0, 0.0d0,
54  & 0.0d0, 0.0d0,
55  2 0.168d0,-16.0d0,0.001439767d0,.9396d0,.9383d0,.135d0,.1396d0,
56  & .93895d0,.1373d0,0.d0,
57  3 5.0d0, 55.0d0, 55.0d0,200.0d0,200.0d0,150.0d0,350.0d0, 0.05d0,
58  & 0.0d0, 0.0d0,
59  4 0.001d0, 0.5d0, 0.4d0, 2.5d0, 0.7d0, 2.d0, 0.36d0, 0.3d0,
60  & 0.2d0, 0.d0,
61  5 2.0d0, 3.5d0, 2.0d0, 1.0d0, 1.0d0, 4.0d0, 0.1d0, 0.d0,
62  & 0.d0, 0.d0,
63  6 4.6d0, 2.8d0, 1.8d0, 1.09d0, 0.7d0, 0.36d0, 0.1d0, 1.4d0,
64  & 0.05d0, 1.0d0,
65  7 10.d0, 0.71d0, 0.1d0, 0.1d0, 1.0d0, 0.d0, 0.0d0, 0.0d0,
66  & 0.d0, 0.0d0,
67  8 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0,
68  & 0.d0, 0.d0,
69  9 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0,
70  & 0.d0, 0.d0/
71 
72 
73  data (parc(i),i=101,200)/
74  $ 10*0.d0,
75  1 0.d0, 9*0.d0,
76  2 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0,
77  & 0.d0, 0.d0,
78  3 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0,
79  & 0.d0, 0.d0,
80  4 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0,
81  & 0.d0, 0.d0,
82  5 2.1d0, 0.3d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0,
83  & 0.d0, 0.d0,
84  6 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0,
85  & 0.d0, 0.d0,
86  7 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0,
87  & 0.d0, 0.d0,
88  8 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0,
89  & 0.d0, 0.d0,
90  9 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0, 0.d0,
91  & 0.d0, 0.d0/
92 
93 !...Elementary particles.
94  data id_quark/1/,id_lept/2/,id_exc/3/,id_boson/4/
95  $ ,id_diq/5/
96  $ ,id_tec/6/
97  $ ,id_susy/7/
98  $ ,id_special/6/
99 
100 !...Light Mesons.
101  data id_pi/101/,id_light1/100/
102  $ id_light0/120/
103  $ ,id_str/130/
104  $ ,id_charm/140/
105  $ ,id_bott/150/
106  $ ,id_cc/160/
107  $ ,id_bb/170/
108  $ ,id_mdiff/199/
109 
110 !...Baryons
111  data id_nucl/11/,id_nucls/12/,id_delt/13/,id_delts/14/
112  $ ,id_lamb/21/,id_lambs/22/
113  $ ,id_sigm/31/,id_sigms/32/
114  $ ,id_xi/41/,id_xis/42/
115  $ ,id_omega/51/
116  $ ,id_charmb/61/,id_bottb/72/
117  $ ,id_bdiff/99/
118 
119 
120  data fname/ '0','JAMRUN%DAT','JAMINFO%DAT','JAMMULTI%DAT'
121  $ ,'JAMDISP%DAT', ' ', ' ', ' '/
122 
123  data check/10*' '/
124 
125  end
126 
127 !*********************************************************************
128 !*********************************************************************
129 !...PYDATA
130 !...Default values for switches and parameters,
131 !...and particle, decay and process data.
132 
133  block data pjdata
134 
135 !...Double precision and integer declarations.
136  implicit double precision(a-h, o-z)
137 !...Commonblocks.
138 #include "jam2.inc"
139  common/jydatr/mrpy(6),rrpy(100)
140  common/pjsubs/msel,mselpd,msub(500),kfin(2,-40:40),ckin(200)
141  common/pjpars/mstp(200),parp(200),msti(200),pari(200)
142  common/pjint1/mint(400),vint(400)
143  common/pjint2/iset(500),kfpr(500,2),coef(500,20),icol(40,4,2)
144  common/pjint3/xsfx(2,-40:40),isig(1000,3),sigh(1000)
145  common/pjint4/mwid(500),wids(500,5)
146  common/pjint5/ngenpd,ngen(0:500,3),xsec(0:500,3)
147  common/pjint6/proc(0:500)
148  character proc*28
149  common/pjint7/sigt(0:6,0:6,0:5)
150  common/pjmssm/imss(0:99),rmss(0:99)
151  common/pjssmt/zmix(4,4),umix(2,2),vmix(2,2),smz(4),smw(2),
152  &sfmix(16,4)
153 ! common/pybins/ihist(4),indx(1000),bin(20000)
154  save /jydatr/,/pjsubs/,
155  &/pjpars/,/pjint1/,/pjint2/,/pjint3/,/pjint4/,/pjint5/,
156  &/pjint6/,/pjint7/,/pjmssm/,/pjssmt/
157 
158 !...PYDAT1, containing status codes and most parameters.
159  data mstu/
160  & 0, 0, 0, 1000,10000, 500, 4000, 0, 0, 2,
161  1 6, 1, 1, 0, 1, 1, 0, 0, 0, 0,
162  2 2, 10, 0, 0, 1, 10, 0, 0, 0, 0,
163  3 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
164  4 2, 2, 1, 4, 2, 1, 1, 0, 0, 0,
165  5 25, 24, 0, 1, 0, 0, 0, 0, 0, 0,
166  6 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
167  7 30*0,
168  1 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
169  2 1, 5, 3, 5, 0, 0, 0, 0, 0, 0,
170  & 80*0/
171  data paru/
172  & 3.141592653589793d0, 6.283185307179586d0,
173  & 0.197327d0, 5.06773d0, 0.389380d0, 2.56819d0, 4*0d0,
174  1 0.001d0, 0.09d0, 0.01d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0,
175  2 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0,
176  3 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0,
177  4 2.0d0, 1.0d0, 0.25d0, 2.5d0, 0.05d0,
178  4 0d0, 0d0, 0.0001d0, 0d0, 0d0,
179  5 2.5d0,1.5d0,7.0d0,1.0d0,0.5d0,2.0d0,3.2d0, 0d0, 0d0, 0d0,
180  6 40*0d0,
181  & 0.00729735d0, 0.232d0, 0.007764d0, 1.0d0, 1.16639d-5,
182  & 0d0, 0d0, 0d0, 0d0, 0d0,
183  1 0.20d0, 0.25d0, 1.0d0, 4.0d0, 10d0, 0d0, 0d0, 0d0, 0d0, 0d0,
184  2 -0.693d0, -1.0d0, 0.387d0, 1.0d0, -0.08d0,
185  2 -1.0d0, 1.0d0, 1.0d0, 1.0d0, 0d0,
186  3 1.0d0,-1.0d0, 1.0d0,-1.0d0, 1.0d0, 0d0, 0d0, 0d0, 0d0, 0d0,
187  4 5.0d0, 1.0d0, 1.0d0, 0d0, 1.0d0, 1.0d0, 0d0, 0d0, 0d0, 0d0,
188  5 1.0d0, 0d0, 0d0, 0d0, 1000d0, 1.0d0, 1.0d0, 1.0d0, 1.0d0,0d0,
189  6 1.0d0, 1.0d0, 1.0d0, 1.0d0, 1.0d0, 0d0, 0d0, 0d0, 0d0, 0d0,
190  7 1.0d0, 1.0d0, 1.0d0, 1.0d0, 1.0d0, 1.0d0, 1.0d0, 0d0,0d0,0d0,
191  8 1.0d0, 1.0d0, 1.0d0, 0.0d0, 0.0d0, 1.0d0, 1.0d0, 0d0,0d0,0d0,
192  9 0d0, 0d0, 0d0, 0d0, 1.0d0, 0d0, 0d0, 0d0, 0d0, 0d0/
193  data mstj/
194  & 1, 3, 0, 0, 0, 0, 0, 0, 0, 2,
195  1 4, 2, 0, 1, 0, 0, 0, 0, 0, 0,
196  2 2, 1, 1, 2, 1, 2, 2, 0, 0, 0,
197  3 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
198  4 2, 2, 4, 2, 5, 3, 3, 0, 0, 3,
199  5 0, 3, 0, 0, 0, 0, 0, 0, 0, 0,
200  6 40*0,
201  & 5, 2, 7, 5, 1, 1, 0, 2, 0, 2,
202  1 0, 0, 0, 0, 1, 1, 0, 0, 0, 0,
203  2 80*0/
204  data parj/
205  & 0.10d0, 0.30d0, 0.40d0, 0.05d0, 0.50d0,
206  & 0.50d0, 0.50d0, 0.6d0, 1.2d0, 0.6d0,
207  1 0.50d0,0.60d0,0.75d0,0.5d0,0.2d0, 0d0, 0d0, 1.0d0, 1.0d0, 0d0,
208  2 0.36d0, 1.0d0,0.01d0, 2.0d0,1.0d0,0.4d0,0.5d0, 0d0, 0d0, 0d0,
209  3 0.10d0, 1.0d0, 0.8d0, 1.5d0,0d0,2.0d0,0.2d0,2.5d0,0.6d0,0d0,
210  4 0.3d0, 0.58d0, 0.5d0, 0.9d0,0.5d0,1.0d0,1.0d0,1.0d0,0d0,0d0,
211  5 0.77d0, 0.77d0, 0.77d0, -0.05d0, -0.005d0,
212  5 -0.00001d0, -0.00001d0, -0.00001d0, 1.0d0, 0d0,
213  6 4.5d0, 0.7d0, 0d0,0.003d0, 0.5d0, 0.5d0, 0d0, 0d0, 0d0, 0d0,
214  7 10d0, 1000d0, 100d0, 1000d0, 0d0, 0.7d0,10d0, 0d0, 0d0, 0d0,
215  8 0.29d0, 1.0d0, 1.0d0, 0d0, 10d0, 10d0, 0d0, 0d0, 0d0, 0d0,
216  9 0.02d0, 1.0d0, 0.2d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0,
217  & 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0,
218  1 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0,
219  2 1.0d0, 0.25d0,91.187d0,2.489d0, 0.01d0,
220  2 2.0d0, 1.0d0, 0.25d0,0.002d0, 0d0,
221  3 0d0, 0d0, 0d0, 0d0, 0.01d0, 0.99d0, 0d0, 0d0, 0.2d0, 0d0,
222  4 60*0d0/
223 
224 
225 !...PYDATR, with initial values for the random number generator.
226  data mrpy/19780503,0,0,97,33,0/
227 
228 !...Default values for allowed processes and kinematics constraints.
229  data msel/1/
230  data msub/500*0/
231  data ((kfin(i,j),j=-40,40),i=1,2)/16*0,4*1,4*0,6*1,5*0,5*1,0,
232  &5*1,5*0,6*1,4*0,4*1,16*0,16*0,4*1,4*0,6*1,5*0,5*1,0,5*1,5*0,
233  &6*1,4*0,4*1,16*0/
234  data ckin/
235  & 2.0d0, -1.0d0, 0.0d0, -1.0d0, 1.0d0,
236  & 1.0d0, -10d0, 10d0, -10d0, 10d0,
237  1 -10d0, 10d0, -10d0, 10d0, -10d0,
238  1 10d0, -1.0d0, 1.0d0, -1.0d0, 1.0d0,
239  2 0.0d0, 1.0d0, 0.0d0, 1.0d0, -1.0d0,
240  2 1.0d0, -1.0d0, 1.0d0, 0d0, 0d0,
241  3 2.0d0, -1.0d0, 0d0, 0d0, 0.0d0,
242  3 -1.0d0, 0.0d0, -1.0d0, 4.0d0, -1.0d0,
243  4 12.0d0, -1.0d0, 12.0d0, -1.0d0, 12.0d0,
244  4 -1.0d0, 12.0d0, -1.0d0, 0d0, 0d0,
245  5 0.0d0, -1.0d0, 0.0d0, -1.0d0, 0.0d0,
246  5 -1.0d0, 0d0, 0d0, 0d0, 0d0,
247  6 140*0d0/
248 
249 !...Default values for main switches and parameters. Reset information.
250  data (mstp(i),i=1,100)/
251  & 3, 1, 2, 0, 0, 0, 0, 0, 0, 0,
252  1 1, 0, 1, 0, 5, 0, 0, 0, 0, 0,
253  2 1, 0, 1, 0, 0, 0, 0, 0, 0, 1,
254  3 1, 2, 0, 1, 0, 2, 1, 5, 2, 0,
255  4 1, 1, 3, 7, 3, 1, 1, 0, 1, 0,
256  5 4, 1, 3, 1, 5, 1, 1, 6, 1, 7,
257  6 1, 3, 2, 2, 1, 1, 2, 0, 0, 0,
258  7 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
259  8 1, 1, 100, 0, 0, 0, 0, 0, 0, 0,
260  9 1, 4, 1, 2, 0, 0, 0, 0, 0, 0/
261  data (mstp(i),i=101,200)/
262  & 3, 1, 0, 0, 0, 0, 0, 0, 0, 0,
263  1 1, 1, 1, 0, 0, 0, 0, 0, 0, 0,
264  2 0, 1, 2, 1, 1, 50, 0, 0, 10, 0,
265  3 0, 4, 0, 1, 0, 0, 0, 0, 0, 0,
266  4 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
267  5 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
268  6 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
269  7 0, 2, 0, 0, 0, 0, 0, 0, 0, 0,
270  8 6, 115, 1998, 01, 27, 0, 0, 0, 0, 0,
271  9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
272  data (parp(i),i=1,100)/
273  & 0.25d0, 10d0, 8*0d0,
274  1 0d0, 0d0, 1.0d0, 0.01d0, 0.6d0, 1.0d0, 1.0d0, 3*0d0,
275  2 10*0d0,
276  3 1.5d0,2.0d0,0.075d0,1.0d0,0.2d0,0d0,2.0d0,0.70d0,0.006d0,0d0,
277  4 0.02d0,2.0d0,0.10d0,1000d0,2054d0, 123d0, 246d0, 50d0, 2*0d0,
278  5 1.0d0, 9*0d0,
279  6 0.25d0, 1.0d0,0.25d0, 1.0d0, 2.0d0,1d-3, 4.0d0,1d-3,2*0d0,
280  7 4.0d0, 0.25d0, 8*0d0,
281  8 1.40d0,1.55d0,0.5d0, 0.2d0,0.33d0,0.66d0, 0.7d0, 0.5d0,2*0d0,
282  9 0.44d0,0.20d0,2.0d0,1.0d0,0d0,3.0d0,1.0d0,0.75d0,0.44d0,2.0d0/
283  data (parp(i),i=101,200)/
284  & 0.5d0, 0.28d0, 1.0d0, 0.8d0, 6*0d0,
285  1 2.0d0, 3*0d0, 1.5d0, 0.5d0, 0.6d0, 2.5d0, 2.0d0, 1.0d0,
286  2 1.0d0, 0.4d0, 8*0d0,
287  3 0.01d0, 9*0d0,
288  4 0.33333d0, 82d0, 1d0, 4d0, 200d0, 5*0d0,
289  5 0d0, 0d0, 0d0, 0d0, 6*0d0,
290  6 2.20d0, 23.6d0, 18.4d0, 11.5d0, 6*0d0,
291  7 0d0, 0d0, 0d0, 1.0d0, 6*0d0,
292  8 20*0d0/
293  data msti/200*0/
294  data pari/200*0d0/
295  data mint/400*0/
296  data vint/400*0d0/
297 
298 !...Constants for the generation of the various processes.
299  data (iset(i),i=1,100)/
300  & 1, 1, 1, -1, 3, -1, -1, 3, -2, 2,
301  1 2, 2, 2, 2, 2, 2, -1, 2, 2, 2,
302  2 -1, 2, 2, 2, 2, 2, -1, 2, 2, 2,
303  3 2, -1, 2, 2, 2, 2, -1, -1, -1, -1,
304  4 -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
305  5 -1, -1, 2, 2, -1, -1, -1, 2, -1, -1,
306  6 -1, -1, -1, -1, -1, -1, -1, 2, 2, 2,
307  7 4, 4, 4, -1, -1, 4, 4, -1, -1, 2,
308  8 2, 2, 2, 2, 2, 2, 2, 2, 2, -2,
309  9 0, 0, 0, 0, 0, 9, -2, -2, -2, -2/
310  data (iset(i),i=101,200)/
311  & -1, 1, 1, -2, -2, 2, 2, 2, -2, 2,
312  1 2, 2, 2, 2, 2, -1, -1, -1, -2, -2,
313  2 5, 5, 5, 5, -2, -2, -2, -2, -2, -2,
314  3 -1, -2, -2, -2, -2, -2, -2, -2, -2, -2,
315  4 1, 1, 1, 1, 1, -2, 1, 1, 1, -2,
316  5 1, 1, 1, -2, -2, 1, 1, 1, -2, -2,
317  6 2, 2, 2, 2, 2, 2, 2, 2, -2, -2,
318  7 2, 2, 5, 5, -2, 2, 2, 5, 5, -2,
319  8 5, 5, -2, -2, -2, 5, 5, -2, -2, -2,
320  9 1, 1, 1, 2, -2, -2, -2, -2, -2, -2/
321  data (iset(i),i=201,300)/
322  & 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
323  1 2, 2, 2, 2, -2, 2, 2, 2, 2, 2,
324  2 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
325  3 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
326  4 2, 2, 2, 2, -1, 2, 2, 2, 2, 2,
327  5 2, 2, 2, 2, -1, 2, -1, 2, 2, -2,
328  6 2, 2, 2, 2, 2, -1, -1, -1, -1, -1,
329  7 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
330  8 -2, -2, -2, -2, -2, -2, -2, -2, -2, -2,
331  9 -2, -2, -2, -2, -2, -2, -2, -2, -2, -2/
332  data (iset(i),i=301,500)/200*-2/
333  data ((kfpr(i,j),j=1,2),i=1,50)/
334  & 23, 0, 24, 0, 25, 0, 24, 0, 25, 0,
335  & 24, 0, 23, 0, 25, 0, 0, 0, 0, 0,
336  1 0, 0, 0, 0, 21, 21, 21, 22, 21, 23,
337  1 21, 24, 21, 25, 22, 22, 22, 23, 22, 24,
338  2 22, 25, 23, 23, 23, 24, 23, 25, 24, 24,
339  2 24, 25, 25, 25, 0, 21, 0, 22, 0, 23,
340  3 0, 24, 0, 25, 0, 21, 0, 22, 0, 23,
341  3 0, 24, 0, 25, 0, 21, 0, 22, 0, 23,
342  4 0, 24, 0, 25, 0, 21, 0, 22, 0, 23,
343  4 0, 24, 0, 25, 0, 21, 0, 22, 0, 23/
344  data ((kfpr(i,j),j=1,2),i=51,100)/
345  5 0, 24, 0, 25, 0, 0, 0, 0, 0, 0,
346  5 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
347  6 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
348  6 0, 0, 0, 0, 21, 21, 24, 24, 23, 24,
349  7 23, 23, 24, 24, 23, 24, 23, 25, 22, 22,
350  7 23, 23, 24, 24, 24, 25, 25, 25, 0, 211,
351  8 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
352  8 443, 21,10441, 21,20443, 21, 445, 21, 0, 0,
353  9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
354  9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
355  data ((kfpr(i,j),j=1,2),i=101,150)/
356  & 23, 0, 25, 0, 25, 0, 0, 0, 0, 0,
357  & 443, 22, 443, 21, 443, 22, 0, 0, 22, 25,
358  1 21, 25, 0, 25, 21, 25, 22, 22, 21, 22,
359  1 22, 23, 23, 23, 24, 24, 0, 0, 0, 0,
360  2 25, 6, 25, 6, 25, 0, 25, 0, 0, 0,
361  2 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
362  3 23, 5, 0, 0, 0, 0, 0, 0, 0, 0,
363  3 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
364  4 32, 0, 34, 0, 37, 0, 40, 0, 39, 0,
365  4 0, 0, 4000001, 0, 4000002, 0, 38, 0, 0, 0/
366  data ((kfpr(i,j),j=1,2),i=151,200)/
367  5 35, 0, 35, 0, 35, 0, 0, 0, 0, 0,
368  5 36, 0, 36, 0, 36, 0, 0, 0, 0, 0,
369  6 6, 37, 39, 0, 39, 39, 39, 39, 11, 0,
370  6 11, 0, 0, 4000001, 0, 4000002, 0, 0, 0, 0,
371  7 23, 35, 24, 35, 35, 0, 35, 0, 0, 0,
372  7 23, 36, 24, 36, 36, 0, 36, 0, 0, 0,
373  8 35, 6, 35, 6, 0, 0, 0, 0, 0, 0,
374  8 36, 6, 36, 6, 0, 0, 0, 0, 0, 0,
375  9 54, 0, 55, 0, 56, 0, 11, 0, 0, 0,
376  9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
377  data ((kfpr(i,j),j=1,2),i=201,240)/
378  & 1000011, 1000011, 2000011, 2000011, 1000011,
379  & 2000011, 1000013, 1000013, 2000013, 2000013,
380  & 1000013, 2000013, 1000015, 1000015, 2000015,
381  & 2000015, 1000015, 2000015, 1000011, 1000012,
382  1 1000015, 1000016, 2000015, 1000016, 1000012,
383  1 1000012, 1000016, 1000016, 0, 0,
384  1 1000022, 1000022, 1000023, 1000023, 1000025,
385  1 1000025, 1000035, 1000035, 1000022, 1000023,
386  2 1000022, 1000025, 1000022, 1000035, 1000023,
387  2 1000025, 1000023, 1000035, 1000025, 1000035,
388  2 1000024, 1000024, 1000037, 1000037, 1000024,
389  2 1000037, 1000022, 1000024, 1000023, 1000024,
390  3 1000025, 1000024, 1000035, 1000024, 1000022,
391  3 1000037, 1000023, 1000037, 1000025, 1000037,
392  3 1000035, 1000037, 1000021, 1000022, 1000021,
393  3 1000023, 1000021, 1000025, 1000021, 1000035/
394  data ((kfpr(i,j),j=1,2),i=241,280)/
395  4 1000021, 1000024, 1000021, 1000037, 1000021,
396  4 1000021, 1000021, 1000021, 0, 0,
397  4 1000002, 1000022, 2000002, 1000022, 1000002,
398  4 1000023, 2000002, 1000023, 1000002, 1000025,
399  5 2000002, 1000025, 1000002, 1000035, 2000002,
400  5 1000035, 1000001, 1000024, 2000005, 1000024,
401  5 1000001, 1000037, 2000005, 1000037, 1000002,
402  5 1000021, 2000002, 1000021, 0, 0,
403  6 1000006, 1000006, 2000006, 2000006, 1000006,
404  6 2000006, 1000006, 1000006, 2000006, 2000006,
405  6 0, 0, 0, 0, 0,
406  6 0, 0, 0, 0, 0,
407  7 1000002, 1000002, 2000002, 2000002, 1000002,
408  7 2000002, 1000002, 1000002, 2000002, 2000002,
409  7 1000002, 2000002, 1000002, 1000002, 2000002,
410  7 2000002, 1000002, 1000002, 2000002, 2000002/
411  data ((kfpr(i,j),j=1,2),i=281,500)/440*0/
412  data coef/10000*0d0/
413  data (((icol(i,j,k),k=1,2),j=1,4),i=1,40)/
414  &4,0,3,0,2,0,1,0,3,0,4,0,1,0,2,0,2,0,0,1,4,0,0,3,3,0,0,4,1,0,0,2,
415  &3,0,0,4,1,4,3,2,4,0,0,3,4,2,1,3,2,0,4,1,4,0,2,3,4,0,3,4,2,0,1,2,
416  &3,2,1,0,1,4,3,0,4,3,3,0,2,1,1,0,3,2,1,4,1,0,0,2,2,4,3,1,2,0,0,1,
417  &3,2,1,4,1,4,3,2,4,2,1,3,4,2,1,3,3,4,4,3,1,2,2,1,2,0,3,1,2,0,0,0,
418  &4,2,1,0,0,0,1,0,3,0,0,3,1,2,0,0,4,0,0,4,0,0,1,2,2,0,0,1,4,4,3,3,
419  &2,2,1,1,4,4,3,3,3,3,4,4,1,1,2,2,3,2,1,3,1,2,0,0,4,2,1,4,0,0,1,2,
420  &4,0,0,0,4,0,1,3,0,0,3,0,2,4,3,0,3,4,0,0,1,0,0,1,0,0,3,4,2,0,0,2,
421  &3,0,0,0,1,0,0,0,0,0,3,0,2,0,0,0,2,0,3,1,2,0,0,0,3,2,1,0,1,0,0,0,
422  &4,4,3,3,2,2,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
423  &0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0/
424 
425 
426 !...Character constants: name of processes.
427  data proc(0)/ 'All included subprocesses '/
428  data (proc(i),i=1,20)/
429  &'f + fbar -> gamma*/Z0 ', 'f + fbar'' -> W+/- ',
430  &'f + fbar -> h0 ', 'gamma + W+/- -> W+/- ',
431  &'Z0 + Z0 -> h0 ', 'Z0 + W+/- -> W+/- ',
432  &' ', 'W+ + W- -> h0 ',
433  &' ', 'f + f'' -> f + f'' (QFD) ',
434  1'f + f'' -> f + f'' (QCD) ','f + fbar -> f'' + fbar'' ',
435  1'f + fbar -> g + g ', 'f + fbar -> g + gamma ',
436  1'f + fbar -> g + Z0 ', 'f + fbar'' -> g + W+/- ',
437  1'f + fbar -> g + h0 ', 'f + fbar -> gamma + gamma ',
438  1'f + fbar -> gamma + Z0 ', 'f + fbar'' -> gamma + W+/- '/
439  data (proc(i),i=21,40)/
440  2'f + fbar -> gamma + h0 ', 'f + fbar -> Z0 + Z0 ',
441  2'f + fbar'' -> Z0 + W+/- ', 'f + fbar -> Z0 + h0 ',
442  2'f + fbar -> W+ + W- ', 'f + fbar'' -> W+/- + h0 ',
443  2'f + fbar -> h0 + h0 ', 'f + g -> f + g ',
444  2'f + g -> f + gamma ', 'f + g -> f + Z0 ',
445  3'f + g -> f'' + W+/- ', 'f + g -> f + h0 ',
446  3'f + gamma -> f + g ', 'f + gamma -> f + gamma ',
447  3'f + gamma -> f + Z0 ', 'f + gamma -> f'' + W+/- ',
448  3'f + gamma -> f + h0 ', 'f + Z0 -> f + g ',
449  3'f + Z0 -> f + gamma ', 'f + Z0 -> f + Z0 '/
450  data (proc(i),i=41,60)/
451  4'f + Z0 -> f'' + W+/- ', 'f + Z0 -> f + h0 ',
452  4'f + W+/- -> f'' + g ', 'f + W+/- -> f'' + gamma ',
453  4'f + W+/- -> f'' + Z0 ', 'f + W+/- -> f'' + W+/- ',
454  4'f + W+/- -> f'' + h0 ', 'f + h0 -> f + g ',
455  4'f + h0 -> f + gamma ', 'f + h0 -> f + Z0 ',
456  5'f + h0 -> f'' + W+/- ', 'f + h0 -> f + h0 ',
457  5'g + g -> f + fbar ', 'g + gamma -> f + fbar ',
458  5'g + Z0 -> f + fbar ', 'g + W+/- -> f + fbar'' ',
459  5'g + h0 -> f + fbar ', 'gamma + gamma -> f + fbar ',
460  5'gamma + Z0 -> f + fbar ', 'gamma + W+/- -> f + fbar'' '/
461  data (proc(i),i=61,80)/
462  6'gamma + h0 -> f + fbar ', 'Z0 + Z0 -> f + fbar ',
463  6'Z0 + W+/- -> f + fbar'' ', 'Z0 + h0 -> f + fbar ',
464  6'W+ + W- -> f + fbar ', 'W+/- + h0 -> f + fbar'' ',
465  6'h0 + h0 -> f + fbar ', 'g + g -> g + g ',
466  6'gamma + gamma -> W+ + W- ', 'gamma + W+/- -> Z0 + W+/- ',
467  7'Z0 + Z0 -> Z0 + Z0 ', 'Z0 + Z0 -> W+ + W- ',
468  7'Z0 + W+/- -> Z0 + W+/- ', 'Z0 + Z0 -> Z0 + h0 ',
469  7'W+ + W- -> gamma + gamma ', 'W+ + W- -> Z0 + Z0 ',
470  7'W+/- + W+/- -> W+/- + W+/- ', 'W+/- + h0 -> W+/- + h0 ',
471  7'h0 + h0 -> h0 + h0 ', 'q + gamma -> q'' + pi+/- '/
472  data (proc(i),i=81,100)/
473  8'q + qbar -> Q + Qbar, mass ', 'g + g -> Q + Qbar, massive ',
474  8'f + q -> f'' + Q, massive ', 'g + gamma -> Q + Qbar, mass ',
475  8'gamma + gamma -> F + Fbar, m', 'g + g -> J/Psi + g ',
476  8'g + g -> chi_0c + g ', 'g + g -> chi_1c + g ',
477  8'g + g -> chi_2c + g ', ' ',
478  9'Elastic scattering ', 'Single diffractive (XB) ',
479  9'Single diffractive (AX) ', 'Double diffractive ',
480  9'Low-pT scattering ', 'Semihard QCD 2 -> 2 ',
481  9' ', ' ',
482  9' ', ' '/
483  data (proc(i),i=101,120)/
484  &'g + g -> gamma*/Z0 ', 'g + g -> h0 ',
485  &'gamma + gamma -> h0 ', ' ',
486  &' ', 'g + g -> J/Psi + gamma ',
487  &'gamma + g -> J/Psi + g ', 'gamma+gamma -> J/Psi + gamma',
488  &' ', 'f + fbar -> gamma + h0 ',
489  1'f + fbar -> g + h0 ', 'q + g -> q + h0 ',
490  1'g + g -> g + h0 ', 'g + g -> gamma + gamma ',
491  1'g + g -> g + gamma ', 'g + g -> gamma + Z0 ',
492  1'g + g -> Z0 + Z0 ', 'g + g -> W+ + W- ',
493  1' ', ' '/
494  data (proc(i),i=121,140)/
495  2'g + g -> Q + Qbar + h0 ', 'q + qbar -> Q + Qbar + h0 ',
496  2'f + f'' -> f + f'' + h0 ',
497  2'f + f'' -> f" + f"'' + h0 ',
498  2' ', ' ',
499  2' ', ' ',
500  2' ', ' ',
501  3'g + g -> Z0 + q + qbar ', ' ',
502  3' ', ' ',
503  3' ', ' ',
504  3' ', ' ',
505  3' ', ' '/
506  data (proc(i),i=141,160)/
507  4'f + fbar -> gamma*/Z0/Z''0 ', 'f + fbar'' -> W''+/- ',
508  4'f + fbar'' -> H+/- ', 'f + fbar'' -> R ',
509  4'q + l -> LQ ', ' ',
510  4'd + g -> d* ', 'u + g -> u* ',
511  4'g + g -> eta_techni ', ' ',
512  5'f + fbar -> H0 ', 'g + g -> H0 ',
513  5'gamma + gamma -> H0 ', ' ',
514  5' ', 'f + fbar -> A0 ',
515  5'g + g -> A0 ', 'gamma + gamma -> A0 ',
516  5' ', ' '/
517  data (proc(i),i=161,180)/
518  6'f + g -> f'' + H+/- ', 'q + g -> LQ + lbar ',
519  6'g + g -> LQ + LQbar ', 'q + qbar -> LQ + LQbar ',
520  6'f + fbar -> f'' + fbar'' (g/Z)',
521  6'f +fbar'' -> f" + fbar"'' (W) ',
522  6'q + q'' -> q" + d* ', 'q + q'' -> q" + u* ',
523  6' ', ' ',
524  7'f + fbar -> Z0 + H0 ', 'f + fbar'' -> W+/- + H0 ',
525  7'f + f'' -> f + f'' + H0 ',
526  7'f + f'' -> f" + f"'' + H0 ',
527  7' ', 'f + fbar -> Z0 + A0 ',
528  7'f + fbar'' -> W+/- + A0 ',
529  7'f + f'' -> f + f'' + A0 ',
530  7'f + f'' -> f" + f"'' + A0 ',
531  7' '/
532  data (proc(i),i=181,200)/
533  8'g + g -> Q + Qbar + H0 ', 'q + qbar -> Q + Qbar + H0 ',
534  8' ', ' ',
535  8' ', 'g + g -> Q + Qbar + A0 ',
536  8'q + qbar -> Q + Qbar + A0 ', ' ',
537  8' ', ' ',
538  9'f + fbar -> rho_tech0 ', 'f + f'' -> rho_tech+/- ',
539  9'f + fbar -> omega_tech0 ', 'f+fbar -> f''+fbar'' (technic)',
540  9' ', ' ',
541  9' ', ' ',
542  9' ', ' '/
543  data (proc(i),i=201,220)/
544  &'f + fbar -> ~e_L + ~e_Lbar ', 'f + fbar -> ~e_R + ~e_Rbar ',
545  &'f + fbar -> ~e_R + ~e_Lbar ', 'f + fbar -> ~mu_L + ~mu_Lbar',
546  &'f + fbar -> ~mu_R + ~mu_Rbar', 'f + fbar -> ~mu_L + ~mu_Rbar',
547  &'f+fbar -> ~tau_1 + ~tau_1bar', 'f+fbar -> ~tau_2 + ~tau_2bar',
548  &'f+fbar -> ~tau_1 + ~tau_2bar', 'q + qbar'' -> ~l_L + ~nulbar ',
549  1'q+qbar''-> ~tau_1 + ~nutaubar', 'q+qbar''-> ~tau_2 + ~nutaubar',
550  1'f + fbar -> ~nul + ~nulbar ', 'f+fbar -> ~nutau + ~nutaubar',
551  1' ', 'f + fbar -> ~chi1 + ~chi1 ',
552  1'f + fbar -> ~chi2 + ~chi2 ', 'f + fbar -> ~chi3 + ~chi3 ',
553  1'f + fbar -> ~chi4 + ~chi4 ', 'f + fbar -> ~chi1 + ~chi2 '/
554  data (proc(i),i=221,240)/
555  2'f + fbar -> ~chi1 + ~chi3 ', 'f + fbar -> ~chi1 + ~chi4 ',
556  2'f + fbar -> ~chi2 + ~chi3 ', 'f + fbar -> ~chi2 + ~chi4 ',
557  2'f + fbar -> ~chi3 + ~chi4 ', 'f+fbar -> ~chi+-1 + ~chi-+1 ',
558  2'f+fbar -> ~chi+-2 + ~chi-+2 ', 'f+fbar -> ~chi+-1 + ~chi-+2 ',
559  2'q + qbar'' -> ~chi1 + ~chi+-1', 'q + qbar'' -> ~chi2 + ~chi+-1',
560  3'q + qbar'' -> ~chi3 + ~chi+-1', 'q + qbar'' -> ~chi4 + ~chi+-1',
561  3'q + qbar'' -> ~chi1 + ~chi+-2', 'q + qbar'' -> ~chi2 + ~chi+-2',
562  3'q + qbar'' -> ~chi3 + ~chi+-2', 'q + qbar'' -> ~chi4 + ~chi+-2',
563  3'q + qbar -> ~chi1 + ~g ', 'q + qbar -> ~chi2 + ~g ',
564  3'q + qbar -> ~chi3 + ~g ', 'q + qbar -> ~chi4 + ~g '/
565  data (proc(i),i=241,260)/
566  4'q + qbar'' -> ~chi+-1 + ~g ', 'q + qbar'' -> ~chi+-2 + ~g ',
567  4'q + qbar -> ~g + ~g ', 'g + g -> ~g + ~g ',
568  4' ', 'qj + g -> ~qj_L + ~chi1 ',
569  4'qj + g -> ~qj_R + ~chi1 ', 'qj + g -> ~qj_L + ~chi2 ',
570  4'qj + g -> ~qj_R + ~chi2 ', 'qj + g -> ~qj_L + ~chi3 ',
571  5'qj + g -> ~qj_R + ~chi3 ', 'qj + g -> ~qj_L + ~chi4 ',
572  5'qj + g -> ~qj_R + ~chi4 ', 'qj + g -> ~qk_L + ~chi+-1 ',
573  5'qj + g -> ~qk_R + ~chi+-1 ', 'qj + g -> ~qk_L + ~chi+-2 ',
574  5'qj + g -> ~qk_R + ~chi+-2 ', 'qj + g -> ~qj_L + ~g ',
575  5'qj + g -> ~qj_R + ~g ', ' '/
576  data (proc(i),i=261,280)/
577  6'f + fbar -> ~t_1 + ~t_1bar ', 'f + fbar -> ~t_2 + ~t_2bar ',
578  6'f + fbar -> ~t_1 + ~t_2bar ', 'g + g -> ~t_1 + ~t_1bar ',
579  6'g + g -> ~t_2 + ~t_2bar ', ' ',
580  6' ', ' ',
581  6' ', ' ',
582  7'qi + qj -> ~qi_L + ~qj_L ', 'qi + qj -> ~qi_R + ~qj_R ',
583  7'qi + qj -> ~qi_L + ~qj_R ', 'qi+qjbar -> ~qi_L + ~qj_Lbar',
584  7'qi+qjbar -> ~qi_R + ~qj_Rbar', 'qi+qjbar -> ~qi_L + ~qj_Rbar',
585  7'f + fbar -> ~qi_L + ~qi_Lbar', 'f + fbar -> ~qi_R + ~qi_Rbar',
586  7'g + g -> ~qi_L + ~qi_Lbar ', 'g + g -> ~qi_R + ~qi_Rbar '/
587  data (proc(i),i=281,500)/220*' '/
588 
589 !...Cross sections and slope offsets.
590  data sigt/294*0d0/
591 
592 !...Supersymmetry switches and parameters.
593  data imss/0,
594  & 0, 0, 0, 1, 0, 0, 0, 1, 0, 0,
595  1 89*0/
596  data rmss/0d0,
597  & 80d0,160d0,500d0,800d0,2d0,250d0,200d0,800d0,700d0,800d0,
598  1 700d0,500d0,250d0,200d0,800d0,400d0,0d0,0.1d0,850d0,0.041d0,
599  2 1d0,800d0,1d4,1d4,1d4,0d0,0d0,24d17,2*0d0,
600  3 69*0d0/
601 
602 !...Data for histogramming routines.
603 ! data ihist/1000,20000,55,1/
604 ! data indx/1000*0/
605 
606  data parf/
607  & 0.5d0,0.25d0, 0.5d0,0.25d0, 1d0, 0.5d0, 0d0, 0d0, 0d0, 0d0,
608  1 0.5d0, 0d0, 0.5d0, 0d0, 1d0, 1d0, 0d0, 0d0, 0d0, 0d0,
609  2 0.5d0, 0d0, 0.5d0, 0d0, 1d0, 1d0, 0d0, 0d0, 0d0, 0d0,
610  3 0.5d0, 0d0, 0.5d0, 0d0, 1d0, 1d0, 0d0, 0d0, 0d0, 0d0,
611  4 0.5d0, 0d0, 0.5d0, 0d0, 1d0, 1d0, 0d0, 0d0, 0d0, 0d0,
612  5 0.5d0, 0d0, 0.5d0, 0d0, 1d0, 1d0, 0d0, 0d0, 0d0, 0d0,
613  6 0.75d0, 0.5d0, 0d0,0.1667d0,0.0833d0,0.1667d0,0d0,0d0,0d0, 0d0,
614  7 0d0, 0d0, 1d0,0.3333d0,0.6667d0,0.3333d0,0d0,0d0,0d0, 0d0,
615  8 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0,
616  9 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0,
617  & 0.325d0,0.325d0,0.5d0,1.6d0, 5.0d0, 0d0, 0d0, 0d0, 0d0, 0d0,
618  1 0d0,0.11d0,0.16d0,0.048d0,0.50d0,0.45d0,0.55d0,0.60d0,0d0,0d0,
619  2 0.2d0, 0.1d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0,
620  3 60*0d0,
621  4 0.2d0, 0.5d0, 8*0d0,
622  5 1800*0d0/
623  data ((vckm(i,j),j=1,4),i=1,4)/
624  & 0.95113d0, 0.04884d0, 0.00003d0, 0.00000d0,
625  & 0.04884d0, 0.94940d0, 0.00176d0, 0.00000d0,
626  & 0.00003d0, 0.00176d0, 0.99821d0, 0.00000d0,
627  & 0.00000d0, 0.00000d0, 0.00000d0, 1.00000d0/
628 
629 ! include 'jydat.f'
630 
631 !...PYDAT2, with particle data and flavour treatment parameters.
632  data (kchg(i,1),i= 1, 500)/-1,2,-1,2,-1,2,-1,2,2*0,-3,0,-3,0,
633  &-3,0,-3,6*0,3,9*0,3,2*0,3,0,-1,12*0,3,2*0,3,28*0,2,-1,15*0,3,4*0,
634  &3,-2,2*1,4,2*-2,2*1,-2,2*1,3*4,2*1,2*-2,2*1,2*-2,2*1,-2,0,3,0,3,
635  &0,3,0,3,0,3,0,3,0,3,0,3,0,3,0,3,22*0,3,0,3,0,3,0,3,0,3,0,3,0,3,0,
636  &3,0,3,0,3,0,4*3,3*0,3*3,4*0,2*3,2*0,2*3,2*0,3,0,3,5*0,3*3,3*0,
637  &3*3,5*0,2*3,2*0,2*3,3*0,3,0,3,3*0,3,0,3,0,3,0,3,0,3,0,3,0,3,0,3,
638  &0,3,0,3,0,3,-3,0,3,6,-3,0,3,6,-3,0,3,6,-3,0,3,6,-3,0,3,6,-3,0,3,
639  &6,-3,0,3,6,-3,0,3,6,-3,0,3,6,-3,0,3,6,13*0,-3,0,3,-3,0,3,-3,0,3,
640  &-3,0,3,-3,0,3,-3,0,3,-3,0,3,-3,0,3,-3,0,3,-3,0,-3,0,-3,0,-3,0,-3,
641  &0,-3,0,-3,2*0,3,0,2*3,2*6,3,2*0,2*3,2*0,2*3,2*6,2*3,6,2*-3,0,-3,
642  &3*0,2*3,0,3,2*-3,2*0,2*-3,3*0,2*3,2*0,2*3,2*-3,2*0,2*-3,2*0,-3,
643  &-1,2,-1,2,-1,2,-3,0,-3,0,-3,4*0,3,2*0,3,0,-1,2,-1,2,-1,2,-3,0,-3,
644  &0,-3,0,-1,2,-3,42*0/
645  data (kchg(i,2),i= 1, 500)/8*1,12*0,2,16*0,2,1,67*0,25*-1,
646  &292*0,6*1,6*0,2,7*0,6*1,6*0,2*1,43*0/
647  data (kchg(i,3),i= 1, 500)/8*1,2*0,8*1,5*0,1,9*0,1,2*0,1,0,2*1,
648  &11*0,1,2*0,1,26*0,1,0,2*1,15*0,1,3*0,27*1,0,1,0,1,0,1,0,1,0,1,0,
649  &1,0,1,0,1,0,1,0,1,21*0,31*1,4*0,6*1,2*0,3*1,2*0,12*1,3*0,8*1,2*0,
650  &4*1,2*0,184*1,3*0,1,2*0,1,0,16*1,41*0/
651  data (kchg(i,4),i= 1, 281)/1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,
652  &16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,
653  &37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,
654  &58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,
655  &79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,
656  &100,210,110,330,220,2110,2210,1103,2101,2103,2203,3101,3103,3201,
657  &3203,3303,4101,4103,4201,4203,4403,4301,4303,5101,5103,5201,5203,
658  &5301,5303,5401,5403,5503,111,211,113,213,10111,10211,20113,20213,
659  &115,215,10113,10213,20111,20211,10115,10215,30113,30213,40113,
660  &40213,10220,221,223,225,331,333,335,10221,10223,10331,10333,
661  &20221,20223,20333,30221,30223,30333,50223,60223,130,310,311,321,
662  &313,323,10313,10323,315,325,10311,10321,20313,20323,30313,30323,
663  &40313,40323,317,327,10315,10325,20315,20325,411,413,415,421,423,
664  &425,431,433,435,440,441,443,445,10411,10413,10421,10423,10431,
665  &10433,10441,10443,20413,20423,20433,20443,100443,511,513,515,521,
666  &523,525,531,533,535,541,543,545,551,553,555,10511,10513,10521,
667  &10523,10531,10533,10541,10543,10551,10553,20513,20523,20533,
668  &20543,20553,100553,2112,2212,12112,12212,1214,2124,22112,22212,
669  &32112,32212,2116,2216,12116,12216,21214,22124,42112,42212,31214,
670  &32124,11218,12128,1114,2114,2214,2224,31114,32114,32214,32224/
671  data (kchg(i,4),i= 282, 500)/1112,1212,2122,2222,11114,12114,
672  &12214,12224,11112,11212,12122,12222,1116,1216,2126,2226,21112,
673  &21212,22122,22222,21114,22114,22214,22224,11116,11216,12126,
674  &12226,1118,2118,2218,2228,3122,13122,3124,23122,33122,13124,
675  &43122,53122,3126,13126,23124,3128,23126,3112,3212,3222,3114,3214,
676  &3224,13112,13212,13222,13114,13214,13224,23112,23212,23222,3116,
677  &3216,3226,13116,13216,13226,23114,23214,23224,3118,3218,3228,
678  &3312,3322,3314,3324,13312,13322,13314,13324,23312,23322,3316,
679  &3326,3334,4112,4114,4122,4132,4212,4214,4222,4224,4232,4312,4314,
680  &4322,4324,4332,4334,4412,4414,4422,4424,4432,4434,4444,5112,5114,
681  &5122,5132,5142,5212,5214,5222,5224,5232,5242,5312,5314,5322,5324,
682  &5332,5334,5342,5412,5414,5422,5424,5432,5434,5442,5444,5512,5514,
683  &5522,5524,5532,5534,5542,5544,5554,1000001,1000002,1000003,
684  &1000004,1000005,1000006,1000011,1000012,1000013,1000014,1000015,
685  &1000016,1000021,1000022,1000023,1000024,1000025,1000035,1000037,
686  &1000039,2000001,2000002,2000003,2000004,2000005,2000006,2000011,
687  &2000012,2000013,2000014,2000015,2000016,4000001,4000002,4000011,
688  &4000012,41*0/
689  data (kchg(i,5),i= 1, 500)/8*1,2*0,8*2,2*0,5*4,2*0,2*4,2*0,6*4,
690  &3*6,10*0,6*6,44*0,4*199,2*99,25*5,2*101,18*100,19*120,24*130,
691  &9*140,4*160,6*140,2*160,3*140,2*160,12*150,3*170,8*150,2*170,
692  &4*150,2*170,2*11,20*12,4*13,36*14,21,12*22,3*31,24*32,2*41,10*42,
693  &51,22*61,35*72,32*7,4*3,41*0/
694  data (kchg(i,6),i= 1, 500)/8*1,96*0,2*3,25*2,120*0,172*3,6*1,
695  &14*0,6*1,6*0,2*1,43*0/
696  data (kchg(i,7),i= 1, 500)/2*0,-1,107*0,4*-1,-2,5*0,2*-1,4*0,
697  &2*-1,42*0,24*1,6*0,3*1,8*0,2*1,4*0,1,8*0,3*-1,10*0,2*-1,6*0,-1,
698  &65*0,40*-1,12*-2,-3,3*0,-1,4*0,5*-1,2*-2,4*0,2*-1,4*0,-1,5*0,-1,
699  &0,4*-1,2*-2,-1,4*0,2*-1,6*0,2*-1,5*0,-1,19*0,-1,54*0/
700  data (pmas(i,1),i= 1, 236)/0.0099d0,0.0056d0,0.199d0,1.35d0,
701  &5d0,175d0,2*400d0,2*0d0,0.00051d0,0d0,0.10566d0,0d0,1.777d0,0d0,
702  &400d0,5*0d0,91.187d0,80.33d0,80d0,6*0d0,500d0,900d0,500d0,
703  &3*300d0,350d0,200d0,5000d0,10*0d0,3*100d0,3*200d0,26*0d0,1d0,2d0,
704  &5d0,21*0d0,0.77133d0,0.57933d0,2*0.77133d0,0.80473d0,0.92953d0,
705  &0.80473d0,0.92953d0,1.09361d0,1.96908d0,2.00808d0,1.96908d0,
706  &2.00808d0,3.27531d0,2.15432d0,2.17967d0,5.38897d0,5.40145d0,
707  &5.38897d0,5.40145d0,5.56725d0,5.57536d0,6.67143d0,6.67397d0,
708  &10.07354d0,0.13498d0,0.13957d0,0.7685d0,0.7669d0,2*0.9835d0,
709  &2*1.23d0,2*1.318d0,2*1.231d0,2*1.3d0,2*1.67d0,2*1.465d0,2*1.7d0,
710  &0.55d0,0.54745d0,0.78194d0,1.275d0,0.95777d0,1.0194d0,1.525d0,
711  &0.98d0,1.17d0,2*1.4d0,1.295d0,1.282d0,1.427d0,1.3d0,1.419d0,
712  &1.68d0,1.512d0,1.662d0,3*0.49767d0,0.4936d0,0.8961d0,0.8916d0,
713  &2*1.29d0,1.432d0,1.425d0,2*1.429d0,2*1.402d0,2*1.4124d0,
714  &2*1.714d0,2*1.77d0,2*1.773d0,2*1.816d0,1.8693d0,2.01d0,2.46d0,
715  &1.8645d0,2.0067d0,2.46d0,1.9685d0,2.1124d0,2.5735d0,0d0,2.9798d0,
716  &3.09688d0,3.5562d0,2.272d0,2.424d0,2.272d0,2.424d0,2.5d0,2.536d0,
717  &3.4151d0,3.46d0,2*2.372d0,2.56d0,3.5106d0,3.686d0,5.2792d0,
718  &5.3248d0,5.83d0,5.2789d0,5.3248d0,5.83d0,5.3693d0,5.4163d0,
719  &6.07d0,6.594d0,6.602d0,7.35d0,9.4d0,9.4603d0,9.9132d0,5.68d0/
720  data (pmas(i,1),i= 237, 500)/5.73d0,5.68d0,5.73d0,5.92d0,5.97d0,
721  &7.25d0,7.3d0,9.8598d0,9.875d0,2*5.78d0,6.02d0,7.3d0,9.8919d0,
722  &10.0233d0,0.93957d0,0.93827d0,2*1.44d0,2*1.52d0,2*1.535d0,
723  &2*1.65d0,2*1.675d0,2*1.68d0,2*1.7d0,2*1.71d0,2*1.72d0,2*1.99d0,
724  &1.234d0,1.233d0,1.232d0,1.231d0,4*1.6d0,4*1.62d0,4*1.7d0,4*1.9d0,
725  &4*1.905d0,4*1.91d0,4*1.92d0,4*1.93d0,4*1.95d0,1.11568d0,1.405d0,
726  &1.52d0,1.6d0,1.67d0,1.69d0,1.8d0,1.81d0,1.82d0,1.83d0,1.89d0,
727  &2.1d0,2.11d0,1.19744d0,1.19255d0,1.18937d0,1.3872d0,1.3837d0,
728  &1.3828d0,3*1.66d0,3*1.67d0,3*1.75d0,3*1.775d0,3*1.915d0,3*1.94d0,
729  &3*2.03d0,1.3213d0,1.3149d0,1.535d0,1.5318d0,2*1.69d0,2*1.82d0,
730  &2*1.95d0,2*2.03d0,1.67245d0,2.4521d0,2.5d0,2.2849d0,2.4703d0,
731  &2.4535d0,2.5d0,2.4529d0,2.5d0,2.4656d0,2.55d0,2.63d0,2.55d0,
732  &2.63d0,2.704d0,2.8d0,3.59798d0,3.65648d0,3.59798d0,3.65648d0,
733  &3.78663d0,3.82466d0,4.91594d0,5.8d0,5.81d0,5.641d0,5.84d0,
734  &7.00575d0,5.8d0,5.81d0,5.8d0,5.81d0,5.84d0,7.00575d0,5.96d0,
735  &5.97d0,5.96d0,5.97d0,6.12d0,6.13d0,7.19099d0,7.03724d0,7.0485d0,
736  &7.03724d0,7.0485d0,7.21101d0,7.219d0,8.30945d0,8.31325d0,
737  &10.42272d0,10.44144d0,10.42272d0,10.44144d0,10.60209d0,
738  &10.61426d0,11.70767d0,11.71147d0,15.11061d0,32*500d0,4*400d0,
739  &41*0d0/
740  data (pmas(i,2),i= 1, 500)/5*0d0,1.4d0,16*0d0,2.47833d0,
741  &2.069d0,0.00295d0,6*0d0,14.67788d0,0d0,16.79392d0,8.45231d0,
742  &4.93534d0,5.80468d0,19.1898d0,0.39162d0,417.35284d0,93*0d0,
743  &0.151d0,0.149d0,2*0.06d0,2*0.4d0,2*0.107d0,2*0.142d0,2*0.4d0,
744  &2*0.24d0,2*0.31d0,2*0.235d0,0.5d0,0d0,0.00843d0,0.185d0,0.0002d0,
745  &0.00443d0,0.076d0,0.04d0,0.36d0,0.25d0,0.08d0,0.053d0,0.025d0,
746  &0.053d0,0.2d0,0.174d0,0.15d0,0.035d0,0.28d0,4*0d0,0.0505d0,
747  &0.0498d0,2*0.09d0,0.109d0,0.098d0,2*0.287d0,2*0.174d0,2*0.227d0,
748  &2*0.323d0,2*0.164d0,2*0.186d0,2*0.276d0,2*0d0,0.023d0,2*0d0,
749  &0.023d0,2*0d0,0.015d0,0d0,0.0013d0,0d0,0.002d0,0.05d0,0.02d0,
750  &0.05d0,0.02d0,0.05d0,0d0,0.014d0,0.01d0,3*0.05d0,0.0009d0,3*0d0,
751  &0.02d0,2*0d0,0.02d0,2*0d0,0.02d0,2*0d0,0.02d0,3*0d0,8*0.05d0,0d0,
752  &0.01d0,4*0.05d0,4*0d0,2*0.35d0,2*0.12d0,6*0.15d0,2*0.13d0,
753  &4*0.1d0,2*0.15d0,2*0.35d0,4*0.12d0,4*0.35d0,4*0.15d0,4*0.3d0,
754  &4*0.2d0,4*0.35d0,4*0.25d0,4*0.2d0,4*0.35d0,4*0.3d0,0d0,0.05d0,
755  &0.0156d0,0.15d0,0.035d0,0.065d0,0.3d0,0.15d0,0.08d0,0.095d0,
756  &0.1d0,2*0.2d0,3*0d0,0.0394d0,0.036d0,0.0358d0,3*0.1d0,3*0.06d0,
757  &3*0.09d0,6*0.12d0,3*0.22d0,3*0.18d0,2*0d0,0.0099d0,0.0091d0,
758  &2*0.05d0,2*0.024d0,2*0.06d0,2*0.02d0,58*0d0,19*1d0,0d0,7*1d0,0d0,
759  &1d0,0d0,1d0,0d0,2.60511d0,2.60839d0,0.42904d0,0.41921d0,41*0d0/
760  data (pmas(i,3),i= 1, 500)/5*0d0,14d0,16*0d0,24.78326d0,
761  &20.69d0,0.02954d0,6*0d0,146.77876d0,0d0,167.93924d0,84.52308d0,
762  &49.35344d0,58.04675d0,191.89803d0,3.91624d0,4173.52832d0,93*0d0,
763  &2*0.4d0,2*0.05d0,2*0.3d0,4*0.25d0,2*0.2d0,2*0.25d0,2*0.2d0,
764  &2*0.3d0,0.24d0,0d0,0.1d0,0.17d0,0.002d0,0.015d0,3*0.2d0,0.35d0,
765  &0.001d0,2*0.05d0,0.02d0,0.3d0,3*0.2d0,0.28d0,4*0d0,2*0.2d0,
766  &0.005d0,0.01d0,2*0.12d0,2*0.4d0,2*0.3d0,2*0.227d0,2*0.323d0,
767  &6*0.2d0,2*0d0,0.12d0,2*0d0,0.12d0,2*0d0,0.05d0,0d0,0.005d0,0d0,
768  &0.01d0,0.1d0,0.08d0,0.1d0,0.08d0,0.1d0,0d0,0.05d0,0.02d0,2*0.1d0,
769  &0.03d0,0.001d0,3*0d0,0.05d0,2*0d0,0.05d0,2*0d0,0.05d0,2*0d0,
770  &0.05d0,3*0d0,6*0.1d0,0.05d0,0.1d0,0d0,0.02d0,4*0.1d0,4*0d0,
771  &2*0.35d0,2*0.12d0,6*0.15d0,2*0.13d0,4*0.1d0,2*0.15d0,2*0.35d0,
772  &4*0.14d0,4*0.35d0,4*0.15d0,4*0.3d0,4*0.2d0,4*0.35d0,4*0.25d0,
773  &4*0.2d0,4*0.35d0,4*0.3d0,0d0,0.05d0,0.062d0,0.15d0,0.035d0,
774  &0.065d0,0.3d0,0.15d0,0.08d0,0.095d0,0.1d0,2*0.2d0,3*0d0,0.04d0,
775  &2*0.035d0,3*0.236d0,3*0.06d0,3*0.09d0,6*0.12d0,3*0.22d0,3*0.18d0,
776  &2*0d0,2*0.05d0,2*0.143d0,2*0.077d0,2*0.103d0,2*0.02d0,58*0d0,
777  &19*10d0,0.00001d0,7*10d0,0.00001d0,10d0,0.00001d0,10d0,0.00001d0,
778  &26.05109d0,26.08388d0,4.29043d0,4.19206d0,41*0d0/
779  data (pmas(i,4),i= 1, 500)/12*0d0,658654d0,0d0,0.0872d0,68*0d0,
780  &0.1d0,0.387d0,46*0d0,0.00003d0,7804.5d0,37*0d0,15500d0,26.762d0,
781  &0d0,3709d0,20*0d0,0.317d0,2*0d0,0.1244d0,2*0d0,0.14d0,19*0d0,
782  &0.468d0,2*0d0,0.462d0,2*0d0,0.483d0,2*0d0,0.15d0,83*0d0,78.88d0,
783  &12*0d0,44.34d0,0d0,23.96d0,24*0d0,49.1d0,87.1d0,10*0d0,24.6d0,
784  &2*0d0,0.0618d0,0.029d0,4*0d0,0.106d0,4*0d0,0.019d0,0d0,7*0.1d0,
785  &2*0d0,0.342d0,2*0.387d0,4*0d0,2*0.387d0,4*0d0,0.387d0,0d0,
786  &18*0.387d0,36*0d0,41*0d0/
787 
788 !...PYDAT3, with particle decay parameters and data.
789  data (mdcy(i,1),i= 1, 500)/5*0,3*1,6*0,1,0,1,5*0,3*1,6*0,1,0,
790  &7*1,10*0,2*1,0,3*1,26*0,3*1,48*0,19*1,0,17*1,4*0,20*1,0,2*1,0,
791  &2*1,0,2*1,2*0,15*1,0,2*1,0,2*1,0,24*1,2*0,60*1,0,12*1,3*0,24*1,
792  &2*0,10*1,2*0,1,3*0,1,0,1,2*0,1,0,1,0,1,0,1,0,1,0,2*1,0,1,4*0,1,0,
793  &1,3*0,1,0,1,0,1,2*0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,21*1,0,7*1,0,1,0,
794  &1,0,4*1,41*0/
795  data (mdcy(i,2),i= 1, 366)/1,9,17,25,33,41,54,64,2*0,74,78,80,
796  &85,87,141,143,148,2*0,151,160,172,188,208,6*0,287,0,309,332,414,
797  &494,521,524,525,10*0,534,539,0,544,564,588,26*0,606,607,611,46*0,
798  &620,622,624,629,631,633,635,637,639,648,655,656,657,660,663,670,
799  &677,683,689,690,691,694,702,708,717,723,734,740,745,748,755,759,
800  &763,777,781,786,789,791,793,796,805,807,809,815,818,821,826,831,
801  &841,852,854,856,858,860,866,872,878,884,893,902,909,916,923,930,
802  &993,996,1004,1065,1067,1075,1108,1110,0,1114,1115,1118,1120,1122,
803  &1124,1126,1128,1130,1132,1134,1135,1137,1139,1141,1143,1154,1190,
804  &1191,1199,1235,1236,1244,1283,1284,1288,1319,1320,1324,1325,1334,
805  &1336,1338,1340,1342,1344,1346,1348,1350,1352,1354,1355,1357,1359,
806  &1361,1363,1365,2*0,1379,1387,1395,1404,1413,1422,1431,1441,1451,
807  &1461,1471,1479,1487,1496,1505,1517,1529,1542,1555,1569,1583,1584,
808  &1587,1590,1591,1596,1605,1614,1619,1624,1633,1642,1647,1652,1661,
809  &1670,1675,1680,1689,1698,1703,1708,1717,1726,1731,1736,1745,1754,
810  &1759,1764,1773,1782,1787,1789,1793,1797,1799,1805,1816,1827,1833,
811  &1837,1840,1848,1856,1865,1874,1884,1894,1905,1914,1925,1939,1950,
812  &1952,1953,1956,1959,1962,1965,1972,1981,1988,1995,2004,2011,2019,
813  &2028,2036,2042,2050,2056,2066,2080,2090,2100,2113,2123,2133,2147,
814  &2157,2162,2165,2167,2169,2176,2183,2190,2197,2204,2211,2214,2217/
815  data (mdcy(i,2),i= 367, 500)/2221,2222,2223,2299,2303,2304,2305,
816  &2306,2307,2311,2312,2314,2315,2317,2321,2322,2326,2330,2334,2338,
817  &2342,2346,2350,2351,2352,2369,2378,2387,2388,2389,2390,2391,2400,
818  &2409,2410,2411,2412,2413,2422,2423,2432,2441,2450,2459,2468,2477,
819  &2486,2495,2504,2513,2522,2531,2540,2549,2558,2567,2576,2585,2597,
820  &2609,2621,2633,2645,2658,2669,2680,2691,2702,2713,2724,2785,2790,
821  &2892,2948,3066,3200,0,3271,3287,3303,3319,3335,3351,3367,0,3382,
822  &0,3397,0,3412,3416,3420,3423,41*0/
823  data (mdcy(i,3),i= 1, 500)/5*8,13,2*10,2*0,4,2,5,2,54,2,5,3,
824  &2*0,9,12,16,20,79,6*0,22,0,23,82,80,27,3,1,9,10*0,2*5,0,20,24,18,
825  &26*0,1,4,9,46*0,2*2,5,5*2,9,7,2*1,2*3,2*7,2*6,2*1,3,8,6,9,6,11,6,
826  &5,3,7,2*4,14,4,5,3,2*2,3,9,2*2,6,2*3,2*5,10,11,4*2,4*6,2*9,4*7,
827  &63,3,8,61,2,8,33,2,4,0,1,3,8*2,1,4*2,11,36,1,8,36,1,8,39,1,4,31,
828  &1,4,1,9,10*2,1,5*2,14,2*0,2*8,4*9,4*10,2*8,2*9,2*12,2*13,2*14,1,
829  &2*3,1,5,2*9,2*5,2*9,2*5,2*9,2*5,2*9,2*5,2*9,2*5,2*9,2*5,2*9,5,2,
830  &2*4,2,6,2*11,6,4,3,2*8,2*9,2*10,11,9,11,14,11,2,1,4*3,7,9,2*7,9,
831  &7,8,9,8,6,8,6,10,14,2*10,13,2*10,14,10,5,3,2*2,6*7,2*3,4,2*1,76,
832  &4,4*1,4,1,2,1,2,4,1,7*4,2*1,17,2*9,4*1,2*9,4*1,9,1,18*9,5*12,13,
833  &6*11,61,5,102,56,118,134,71,0,6*16,15,0,15,0,15,0,2*4,3,2,41*0/
834  data (mdme(i,1),i= 1,4000)/6*1,-1,7*1,-1,7*1,-1,7*1,-1,7*1,-1,
835  &7*1,-1,1,-1,12*1,2*-1,8*1,2*-1,73*1,-1,2*1,-1,6*1,2*-1,7*1,2*-1,
836  &3*1,-1,6*1,2*-1,6*1,2*-1,3*1,-1,3*1,-1,3*1,5*-1,3*1,-1,85*1,2*-1,
837  &6*1,8*-1,3*1,-1,3*1,-1,3*1,5*-1,3*1,4*-1,197*1,2*-1,2*1,-1,20*1,
838  &2*-1,6*1,2*-1,7*1,-1,3*1,-1,3*1,5*-1,3*1,-1,1,-1,6*1,2*-1,6*1,
839  &2*-1,2819*1,576*0/
840  data (mdme(i,2),i= 1,4000)/43*102,4*0,102,0,4*53,3*102,4*0,102,
841  &2*0,3*102,4*0,102,2*0,6*102,42,6*102,2*42,2*0,8*41,2*0,36*41,
842  &8*102,0,102,0,102,2*0,21*102,8*32,8*0,16*32,21*0,62*53,8*32,14*0,
843  &16*32,27*0,62*53,18*0,62*53,9*0,18*53,3*32,0,6*32,3*0,2*32,3*0,
844  &2*32,7*0,8*32,12*0,16*32,6*0,8*32,8*0,12,2*42,2*11,9*42,0,2,2*0,
845  &3,4*0,3,68*0,2,3*0,1,0,3,18*0,2*3,3*0,1,69*0,4*42,11*0,2*42,2*3,
846  &0,2*3,110*0,22*42,41*0,2*3,9*0,16*42,45*0,3,10*0,10*42,20*0,2*13,
847  &6*0,12,2*0,12,0,12,13*0,2*12,7*0,12,2*0,12,8*0,14*42,16*0,48,
848  &3*13,2*42,9*0,14*42,16*0,48,3*13,2*42,9*0,14*42,19*0,48,3*13,
849  &2*42,6*0,2*11,28*42,5*0,32,3*0,4*32,2*4,0,32,17*0,2*32,9*0,32,
850  &3*0,4*32,2*4,849*0,14*42,52*0,10*13,2*42,2*11,4*0,2*42,2*11,6*0,
851  &2*42,2*11,0,2*42,2*11,2*42,2*11,2*42,2*11,2*42,2*11,2*42,2*11,
852  &2*42,2*11,2*42,2*11,2*0,3*42,8*0,48,3*13,20*42,4*0,18*42,4*0,
853  &9*42,0,162*42,828*53,588*0/
854  data (mdme(i,3),i= 1,4000)/623*0,1,4*0,1,9*0,2*1,2,0,3*2,1,2,
855  &2*1,2,0,2,1,2,2*0,2*1,0,2*1,0,2,2*1,2,2*1,2*2,2*1,2,2*1,2,1,5*0,
856  &1,20*0,1,3*0,2*2,2*0,4*2,7*0,2*1,8*0,1,6*2,23*0,3*1,49*0,2*1,0,
857  &2*1,11*0,2*2,2*1,2*0,3*1,3*2,2*1,2*0,3*1,2,1,10*0,2*1,4*0,2*1,
858  &2*0,2*1,4*0,2*1,4*0,4*2,3*3,2*1,6*2,3*3,2*0,2*1,0,2*1,2*0,2*1,0,
859  &2*1,2*0,2*1,0,2*1,2*0,2*1,0,2*1,63*0,2*1,0,2*2,2*1,65*0,1,0,2*2,
860  &2*1,38*0,1,2*2,2*1,5*0,1,71*0,2*2,2*1,41*0,2*2,2*1,44*0,2*2,2*1,
861  &32*0,2*2,2*1,55*0,16*1,3*2,3*0,6*2,3*0,3*2,3*0,3*2,6*0,3*2,6*0,
862  &3*2,7*0,3*2,4*0,20*2,3*3,5*1,3*3,5*1,6*2,2*0,7*2,2*0,2,50*1,28*3,
863  &3*1,0,2*1,0,29*1,2,0,2*2,0,2*2,2*0,3*2,2*0,2*2,2*0,3*2,2*0,2,0,
864  &2*2,0,2,3*0,3*2,5*0,4*2,5*0,3*2,3*0,2,2*0,3*2,4*0,5*2,4*0,5*2,
865  &2*0,3*2,3,3*2,3*3,5*2,4*3,5*2,3*3,3*2,3,56*1,12*2,34*3,7*0,8*2,
866  &8*1,9*0,6*2,11*0,2*2,10*1,6*3,3*1,2*3,9*2,5*1,3*3,3*1,12*4,2*2,
867  &5*3,3*1,3*3,6*0,32*1,4*2,0,7*2,2*0,6*2,0,2*2,25*0,20*2,6*3,1,
868  &11*3,2*1,10*3,1,3*3,5*2,0,2,3*0,6*2,2*0,2,4*0,5*2,0,2,3*0,6*3,2,
869  &12*3,2,10*3,2,3*3,8*0,18*1,14*2,14*0,6*2,3*1,2*0,1,81*0,1,0,1,
870  &5*0,1,2*0,1,35*0,1,36*0,1,0,1,241*0,1,3*0,1,7*0,1,3*0,2*1,6*0,2,
871  &1,47*0,1,10*0,2*1,206*0,2*1,10*0,2*1,29*0,2*1,44*0,8*1,16*0,8*1,
872  &40*0,2*1,60*0,8*1,16*0,8*1,36*0,3*1,5*0,1,17*0,2,1,2*0,4*2,4*1,
873  &4*0,4*2,4*1,69*0,1,3*0,1,11*0,1,3*0,2*1,10*0,2,33*0,1,600*0/
874  data (brat(i) ,i= 1, 418)/43*0d0,0.00003d0,0.00177d0,0.9982d0,
875  &33*0d0,1d0,6*0d0,0.1783d0,0.1735d0,0.1131d0,0.2494d0,0.003d0,
876  &0.09d0,0.0027d0,0.01d0,0.0014d0,0.0012d0,2*0.00025d0,0.0071d0,
877  &0.012d0,0.0004d0,0.00075d0,0.00006d0,2*0.00078d0,0.0034d0,0.08d0,
878  &0.011d0,0.0191d0,0.00006d0,0.005d0,0.0133d0,0.0067d0,0.0005d0,
879  &0.0035d0,0.0006d0,0.0015d0,0.00021d0,0.0002d0,0.00075d0,0.0001d0,
880  &0.0002d0,0.0011d0,3*0.0002d0,0.00022d0,0.0004d0,0.0001d0,
881  &2*0.00205d0,2*0.00069d0,0.00025d0,0.00051d0,0.00025d0,35*0d0,
882  &0.15403d0,0.11945d0,0.15402d0,0.11931d0,0.15215d0,3*0d0,
883  &0.03357d0,0.0668d0,0.03357d0,0.0668d0,0.0335d0,0.0668d0,2*0d0,
884  &0.32139d0,0.0165d0,2*0d0,0.0165d0,0.32067d0,2*0d0,0.00001d0,
885  &0.00059d0,6*0d0,2*0.10814d0,0.10806d0,3*0d0,0.00031d0,0.04438d0,
886  &0.88031d0,4*0d0,0.0002d0,0.05531d0,0d0,0.01838d0,0.00071d0,0d0,
887  &0.00009d0,0.00032d0,62*0d0,0.14449d0,0.11223d0,0.14449d0,
888  &0.11223d0,0.14443d0,0.05782d0,2*0d0,0.03172d0,0.06305d0,
889  &0.03172d0,0.06305d0,0.03172d0,0.06305d0,8*0d0,0.24928d0,0.0128d0,
890  &0.00001d0,0d0,0.0128d0,0.24882d0,0.00039d0,0d0,0.00001d0,
891  &0.00046d0,0.22153d0,5*0d0,2*0.08464d0,0.08463d0,7*0d0,0.00005d0,
892  &0.00097d0,5*0d0,0.00007d0,0d0,0.00049d0,0.00001d0,0.00006d0,
893  &0.30591d0,0.68863d0,0d0,0.0038d0,66*0d0,0.00008d0,0.00167d0/
894  data (brat(i) ,i= 419, 733)/5*0d0,0.00013d0,0d0,0.00294d0,
895  &0.00001d0,3*0d0,0.99517d0,63*0d0,0.00002d0,0.07231d0,2*0d0,
896  &0.00001d0,0.00269d0,0d0,0.92497d0,18*0d0,0.0024d0,0.99483d0,
897  &0.00278d0,1d0,3*0.21511d0,0.21478d0,2*0d0,2*0.06995d0,2*0d0,1d0,
898  &3*0d0,0.95d0,0.05d0,3*0d0,4*0.25d0,16*0d0,4*0.25d0,20*0d0,1d0,
899  &17*0d0,1d0,2*0.08d0,0.76d0,0.08d0,2*0.105d0,0.04d0,0.5d0,0.08d0,
900  &0.14d0,0.01d0,0.015d0,0.005d0,0.988d0,0.012d0,0.999877d0,
901  &0.000123d0,0.998739d0,0.00079d0,0.00038d0,0.000046d0,0.000045d0,
902  &0.99955d0,0.00045d0,0.78d0,0.22d0,0.78d0,0.22d0,4*0.5d0,
903  &2*0.34725d0,0.144d0,0.104d0,0.0245d0,2*0.01225d0,0.0028d0,
904  &0.0057d0,2*0.34725d0,0.144d0,0.104d0,0.049d0,0.0028d0,0.0057d0,
905  &2*1d0,2*0.25d0,0.5d0,2*0.25d0,0.5d0,0.562d0,2*0.155d0,0.08d0,
906  &2*0.021d0,0.006d0,0.562d0,2*0.155d0,0.08d0,2*0.021d0,0.006d0,
907  &0.5d0,0.14333d0,0.28667d0,0.04d0,0.02d0,0.01d0,0.5d0,0.14333d0,
908  &0.28667d0,0.04d0,0.02d0,0.01d0,2*1d0,3*0.33333d0,0.3923d0,
909  &0.321d0,0.2317d0,0.0478d0,0.0049d0,0.0013d0,0.0003d0,0.0007d0,
910  &0.89d0,0.08693d0,0.0221d0,0.00083d0,2*0.00007d0,0.564d0,0.282d0,
911  &0.072d0,0.028d0,0.023d0,2*0.0115d0,0.005d0,0.003d0,0.437d0,
912  &0.208d0,0.302d0,0.0302d0,0.0212d0,0.0016d0,0.48947d0,0.34d0,
913  &3*0.043d0,0.027d0,0.0126d0,0.0013d0,0.0003d0,0.00025d0,0.00008d0/
914  data (brat(i) ,i= 734, 905)/0.444d0,2*0.222d0,0.104d0,2*0.004d0,
915  &0.52d0,0.26d0,0.11d0,2*0.055d0,0.333d0,0.334d0,0.333d0,0.36d0,
916  &0.18d0,0.03d0,2*0.015d0,2*0.2d0,4*0.25d0,3*0.16666d0,0.5d0,
917  &3*0.146d0,3*0.05d0,0.15d0,2*0.05d0,4*0.024d0,0.066d0,4*0.25d0,
918  &0.624d0,0.312d0,0.032d0,2*0.016d0,3*0.33333d0,4*0.5d0,
919  &3*0.33333d0,0.2112d0,0.1256d0,2*0.1939d0,2*0.1359d0,0.002d0,
920  &0.001d0,0.0006d0,0.6861d0,0.3139d0,2*0.5d0,0.6352d0,0.2116d0,
921  &0.0559d0,0.0173d0,0.0482d0,0.0318d0,0.665d0,0.333d0,0.002d0,
922  &0.666d0,0.333d0,0.001d0,0.28d0,0.14d0,0.313d0,0.157d0,0.11d0,
923  &0.28d0,0.14d0,0.313d0,0.157d0,0.11d0,0.333d0,0.166d0,0.168d0,
924  &0.084d0,0.087d0,0.043d0,0.059d0,2*0.029d0,0.002d0,0.332d0,
925  &0.166d0,0.168d0,0.084d0,0.086d0,0.043d0,0.059d0,2*0.029d0,
926  &2*0.002d0,0.667d0,0.333d0,0.667d0,0.333d0,0.667d0,0.333d0,
927  &0.667d0,0.333d0,0.576d0,0.288d0,0.044d0,0.022d0,0.04667d0,
928  &0.02333d0,0.576d0,0.288d0,0.044d0,0.022d0,0.04667d0,0.02333d0,
929  &0.258d0,0.129d0,0.20934d0,0.10466d0,0.1993d0,0.0997d0,0.258d0,
930  &0.129d0,0.20934d0,0.10466d0,0.1993d0,0.0997d0,0.3d0,0.15d0,
931  &0.182d0,0.091d0,0.12867d0,0.06433d0,0.08d0,0.00267d0,0.00133d0,
932  &0.3d0,0.15d0,0.182d0,0.091d0,0.12867d0,0.06433d0,0.08d0,
933  &0.00267d0,0.00133d0,0.53334d0,0.26666d0,0.03334d0,0.01666d0/
934  data (brat(i) ,i= 906,1078)/3*0.05d0,0.53334d0,0.26666d0,
935  &0.03334d0,0.01666d0,3*0.05d0,0.53334d0,0.26666d0,0.03334d0,
936  &0.01666d0,3*0.05d0,0.53334d0,0.26666d0,0.03334d0,0.01666d0,
937  &3*0.05d0,0.07d0,0.065d0,2*0.005d0,2*0.011d0,5*0.001d0,0.07d0,
938  &0.065d0,2*0.005d0,2*0.011d0,5*0.001d0,0.026d0,0.019d0,0.066d0,
939  &0.041d0,0.045d0,0.076d0,0.0073d0,2*0.0047d0,0.026d0,0.001d0,
940  &0.0006d0,0.0066d0,0.005d0,2*0.003d0,2*0.0006d0,2*0.001d0,0.006d0,
941  &0.005d0,0.012d0,0.0057d0,0.067d0,0.008d0,0.0022d0,0.027d0,
942  &0.004d0,0.019d0,0.012d0,0.002d0,0.009d0,0.0218d0,0.001d0,0.022d0,
943  &0.087d0,0.001d0,0.0019d0,0.0015d0,0.0028d0,0.683d0,0.306d0,
944  &0.011d0,0.3d0,0.15d0,0.16d0,0.08d0,0.13d0,0.06d0,0.08d0,0.04d0,
945  &0.034d0,0.027d0,2*0.002d0,2*0.004d0,2*0.002d0,0.034d0,0.027d0,
946  &2*0.002d0,2*0.004d0,2*0.002d0,0.0365d0,0.045d0,0.073d0,0.062d0,
947  &3*0.021d0,0.0061d0,0.015d0,0.025d0,0.0088d0,0.074d0,0.0109d0,
948  &0.0041d0,0.002d0,0.0035d0,0.0011d0,0.001d0,0.0027d0,2*0.0016d0,
949  &0.0018d0,0.011d0,0.0063d0,0.0052d0,0.018d0,0.016d0,0.0034d0,
950  &0.0036d0,0.0009d0,0.0006d0,0.015d0,0.0923d0,0.018d0,0.022d0,
951  &0.0077d0,0.009d0,0.0075d0,0.024d0,0.0085d0,0.067d0,0.0511d0,
952  &0.017d0,0.0004d0,0.0028d0,0.619d0,0.381d0,0.3d0,0.15d0,0.16d0,
953  &0.08d0,0.13d0,0.06d0,0.08d0,0.04d0,0.01d0,2*0.02d0,0.03d0/
954  data (brat(i) ,i=1079,1249)/2*0.005d0,2*0.02d0,0.03d0,2*0.005d0,
955  &0.015d0,0.037d0,0.028d0,0.079d0,0.095d0,0.052d0,0.0078d0,
956  &4*0.001d0,0.028d0,0.033d0,0.026d0,0.05d0,0.01d0,4*0.005d0,0.25d0,
957  &0.0952d0,0.94d0,0.06d0,2*0.4d0,2*0.1d0,1d0,0.0602d0,0.0601d0,
958  &0.8797d0,0.135d0,0.865d0,0.667d0,0.333d0,0.667d0,0.333d0,0.667d0,
959  &0.333d0,0.667d0,0.333d0,4*0.5d0,0.007d0,0.993d0,1d0,0.667d0,
960  &0.333d0,0.667d0,0.333d0,2*0.5d0,0.273d0,0.727d0,2*0.0083d0,
961  &0.1866d0,0.324d0,0.184d0,0.027d0,0.001d0,0.093d0,0.087d0,0.078d0,
962  &0.0028d0,0.02d0,0.055d0,2*0.005d0,0.008d0,0.012d0,0.02d0,0.055d0,
963  &2*0.005d0,0.008d0,0.012d0,0.01d0,0.03d0,0.0035d0,0.011d0,
964  &0.0055d0,0.0042d0,0.009d0,0.018d0,0.015d0,0.0185d0,0.0135d0,
965  &0.025d0,0.0004d0,0.0007d0,0.0008d0,0.0014d0,0.0019d0,0.0025d0,
966  &0.4291d0,0.08d0,0.07d0,0.02d0,0.015d0,0.005d0,1d0,0.3d0,0.15d0,
967  &0.16d0,0.08d0,0.13d0,0.06d0,0.08d0,0.04d0,0.02d0,0.055d0,
968  &2*0.005d0,0.008d0,0.012d0,0.02d0,0.055d0,2*0.005d0,0.008d0,
969  &0.012d0,0.01d0,0.03d0,0.0035d0,0.011d0,0.0055d0,0.0042d0,0.009d0,
970  &0.018d0,0.015d0,0.0185d0,0.0135d0,0.025d0,0.0004d0,0.0007d0,
971  &0.0008d0,0.0014d0,0.0019d0,0.0025d0,0.4291d0,0.08d0,0.07d0,
972  &0.02d0,0.015d0,0.005d0,1d0,0.3d0,0.15d0,0.16d0,0.08d0,0.13d0,
973  &0.06d0,0.08d0,0.04d0,0.02d0,0.055d0,2*0.005d0,0.008d0,0.012d0/
974  data (brat(i) ,i=1250,1416)/0.02d0,0.055d0,2*0.005d0,0.008d0,
975  &0.012d0,0.01d0,0.03d0,0.0035d0,0.011d0,0.0055d0,0.0042d0,0.009d0,
976  &0.018d0,0.015d0,0.0185d0,0.0135d0,0.025d0,2*0.0002d0,0.0007d0,
977  &2*0.0004d0,0.0014d0,0.001d0,0.0009d0,0.0025d0,0.4291d0,0.08d0,
978  &0.07d0,0.02d0,0.015d0,0.005d0,1d0,2*0.3d0,2*0.2d0,0.047d0,
979  &0.122d0,0.006d0,0.012d0,0.035d0,0.012d0,0.035d0,0.003d0,0.007d0,
980  &0.15d0,0.037d0,0.008d0,0.002d0,0.05d0,0.015d0,0.003d0,0.001d0,
981  &0.014d0,0.042d0,0.014d0,0.042d0,0.24d0,0.065d0,0.012d0,0.003d0,
982  &0.001d0,0.002d0,0.001d0,0.002d0,0.014d0,0.003d0,1d0,2*0.3d0,
983  &2*0.2d0,1d0,0.0252d0,0.0248d0,0.0267d0,0.015d0,0.045d0,0.015d0,
984  &0.045d0,0.7743d0,0.029d0,0.22d0,0.78d0,0.667d0,0.333d0,0.667d0,
985  &0.333d0,0.667d0,0.333d0,0.667d0,0.333d0,8*0.5d0,0.02d0,0.98d0,
986  &1d0,0.667d0,0.333d0,0.667d0,0.333d0,4*0.5d0,0.35d0,0.65d0,
987  &3*0.014d0,0.008d0,0.024d0,0.008d0,0.024d0,0.425d0,0.02d0,0.185d0,
988  &0.088d0,0.043d0,0.067d0,0.066d0,0.05d0,0.21667d0,0.43333d0,
989  &0.125d0,0.08333d0,0.04167d0,0.01667d0,0.03333d0,0.05d0,0.21667d0,
990  &0.43333d0,0.125d0,0.08333d0,0.04167d0,0.01667d0,0.03333d0,
991  &0.049d0,0.18333d0,0.36667d0,0.1125d0,0.075d0,0.0375d0,0.05833d0,
992  &0.11667d0,0.001d0,0.049d0,0.18333d0,0.36667d0,0.1125d0,0.075d0,
993  &0.0375d0,0.05833d0,0.11667d0,0.001d0,0.05d0,0.15d0,0.3d0,0.025d0/
994  data (brat(i) ,i=1417,1558)/0.01667d0,0.00833d0,0.01667d0,
995  &0.03333d0,0.4d0,0.05d0,0.15d0,0.3d0,0.025d0,0.01667d0,0.00833d0,
996  &0.01667d0,0.03333d0,0.4d0,0.05d0,0.23333d0,0.46667d0,0.025d0,
997  &0.01667d0,0.00833d0,0.04333d0,0.08667d0,0.01d0,0.06d0,0.05d0,
998  &0.23333d0,0.46667d0,0.025d0,0.01667d0,0.00833d0,0.04333d0,
999  &0.08667d0,0.01d0,0.06d0,0.009d0,0.11333d0,0.22667d0,0.275d0,
1000  &0.18333d0,0.09167d0,0.03d0,0.06d0,0.01d0,0.001d0,0.009d0,
1001  &0.11333d0,0.22667d0,0.275d0,0.18333d0,0.09167d0,0.03d0,0.06d0,
1002  &0.01d0,0.001d0,0.15d0,0.21667d0,0.43333d0,0.05d0,0.03333d0,
1003  &0.01667d0,0.03333d0,0.06667d0,0.15d0,0.21667d0,0.43333d0,0.05d0,
1004  &0.03333d0,0.01667d0,0.03333d0,0.06667d0,0.399d0,0.03333d0,
1005  &0.06667d0,0.1995d0,0.133d0,0.0665d0,0.03333d0,0.06667d0,0.002d0,
1006  &0.399d0,0.03333d0,0.06667d0,0.1995d0,0.133d0,0.0665d0,0.03333d0,
1007  &0.06667d0,0.002d0,0.015d0,0.05d0,0.1d0,0.075d0,0.05d0,0.025d0,
1008  &0.04167d0,0.08333d0,0.3d0,0.2d0,0.02d0,0.04d0,0.015d0,0.05d0,
1009  &0.1d0,0.075d0,0.05d0,0.025d0,0.04167d0,0.08333d0,0.3d0,0.2d0,
1010  &0.02d0,0.04d0,0.1d0,0.05d0,0.1d0,0.05d0,0.03333d0,0.01667d0,
1011  &0.06667d0,0.13333d0,0.04d0,0.05d0,0.02d0,0.04d0,0.3d0,0.1d0,
1012  &0.05d0,0.1d0,0.05d0,0.03333d0,0.01667d0,0.06667d0,0.13333d0,
1013  &0.04d0,0.05d0,0.02d0,0.04d0,0.3d0,0.03d0,0.05d0,0.1d0,0.05d0/
1014  data (brat(i) ,i=1559,1702)/0.03333d0,0.01667d0,0.06667d0,
1015  &0.13333d0,0.04d0,0.05d0,0.03333d0,0.06667d0,0.3d0,2*0.03d0,
1016  &0.05d0,0.1d0,0.05d0,0.03333d0,0.01667d0,0.06667d0,0.13333d0,
1017  &0.04d0,0.05d0,0.03333d0,0.06667d0,0.3d0,0.03d0,1d0,0.331d0,
1018  &0.663d0,0.006d0,0.663d0,0.331d0,0.006d0,1d0,0.15d0,0.175d0,
1019  &0.33d0,0.22d0,0.125d0,0.1d0,0.05d0,0.11667d0,0.05833d0,0.03667d0,
1020  &0.29333d0,0.22d0,0.08333d0,0.04167d0,0.1d0,0.05d0,0.11667d0,
1021  &0.05833d0,0.03667d0,0.29333d0,0.22d0,0.08333d0,0.04167d0,0.15d0,
1022  &0.175d0,0.33d0,0.22d0,0.125d0,0.05d0,0.25d0,0.3d0,2*0.2d0,
1023  &0.03333d0,0.01667d0,0.16667d0,0.08333d0,0.03333d0,0.26667d0,
1024  &0.2d0,0.13333d0,0.06667d0,0.03333d0,0.01667d0,0.16667d0,
1025  &0.08333d0,0.03333d0,0.26667d0,0.2d0,0.13333d0,0.06667d0,0.05d0,
1026  &0.25d0,0.3d0,3*0.2d0,0.24d0,0.16d0,0.3984d0,0.0016d0,0.13333d0,
1027  &0.06667d0,0.02667d0,0.21333d0,0.16d0,0.2656d0,0.1328d0,0.00107d0,
1028  &0.00053d0,0.13333d0,0.06667d0,0.02667d0,0.21333d0,0.16d0,
1029  &0.2656d0,0.1328d0,0.00107d0,0.00053d0,0.2d0,0.24d0,0.16d0,
1030  &0.3984d0,0.0016d0,0.4d0,0.1d0,0.03d0,0.02d0,0.45d0,0.26667d0,
1031  &0.13333d0,0.06667d0,0.03333d0,0.00333d0,0.02667d0,0.02d0,0.3d0,
1032  &0.15d0,0.26667d0,0.13333d0,0.06667d0,0.03333d0,0.00333d0,
1033  &0.02667d0,0.02d0,0.3d0,0.15d0,0.4d0,0.1d0,0.03d0,0.02d0,0.45d0/
1034  data (brat(i) ,i=1703,1866)/0.1d0,0.08999d0,0.05999d0,0.75d0,
1035  &0.00002d0,0.06667d0,0.03333d0,0.01d0,0.07999d0,0.05999d0,0.5d0,
1036  &0.25d0,2*0.00001d0,0.06667d0,0.03333d0,0.01d0,0.07999d0,
1037  &0.05999d0,0.5d0,0.25d0,2*0.00001d0,0.1d0,0.08999d0,0.05999d0,
1038  &0.75d0,0.00002d0,0.6d0,0.225d0,0.015d0,0.01d0,0.15d0,0.4d0,0.2d0,
1039  &0.15d0,0.075d0,0.00167d0,0.01333d0,0.01d0,0.1d0,0.05d0,0.4d0,
1040  &0.2d0,0.15d0,0.075d0,0.00167d0,0.01333d0,0.01d0,0.1d0,0.05d0,
1041  &0.6d0,0.225d0,0.015d0,0.01d0,0.15d0,0.37d0,0.2d0,0.24d0,0.16d0,
1042  &0.03d0,0.24667d0,0.12333d0,0.13333d0,0.06667d0,0.02667d0,
1043  &0.21333d0,0.16d0,0.02d0,0.01d0,0.24667d0,0.12333d0,0.13333d0,
1044  &0.06667d0,0.02667d0,0.21333d0,0.16d0,0.02d0,0.01d0,0.37d0,0.2d0,
1045  &0.24d0,0.16d0,0.03d0,0.15d0,0.85d0,0.1d0,0.05d0,0.56667d0,
1046  &0.28333d0,0.1d0,0.05d0,0.56667d0,0.28333d0,0.15d0,0.85d0,0.19d0,
1047  &0.4d0,0.18d0,0.12d0,0.1d0,0.01d0,0.12667d0,0.06333d0,0.26667d0,
1048  &0.13333d0,0.02d0,0.16d0,0.12d0,0.06667d0,0.03333d0,0.00667d0,
1049  &0.00333d0,0.12667d0,0.06333d0,0.26667d0,0.13333d0,0.02d0,0.16d0,
1050  &0.12d0,0.06667d0,0.03333d0,0.00667d0,0.00333d0,0.19d0,0.4d0,
1051  &0.18d0,0.12d0,0.1d0,0.01d0,0.639d0,0.358d0,0.002d0,0.001d0,
1052  &3*0.33333d0,2*0.23d0,3*0.14d0,3*0.04d0,2*0.1d0,3*0.2d0,
1053  &3*0.06667d0,2*0.1d0,3*0.11667d0,0.35d0,3*0.03333d0,2*0.125d0/
1054  data (brat(i) ,i=1867,2097)/3*0.1d0,0.25d0,3*0.06667d0,
1055  &2*0.175d0,3*0.01667d0,3*0.1d0,2*0.15d0,2*0.125d0,3*0.08667d0,
1056  &3*0.00333d0,2*0.24d0,2*0.3d0,3*0.03333d0,0.06d0,3*0.04d0,
1057  &2*0.06d0,2*0.0275d0,3*0.16667d0,0.02d0,3*0.14167d0,2*0.15d0,
1058  &3*0.01667d0,3*0.05d0,0.35d0,2*0.075d0,2*0.125d0,3*0.01667d0,
1059  &0.03d0,2*0.015d0,3*0.09333d0,0.051d0,2*0.1545d0,2*0.025d0,
1060  &3*0.03333d0,3*0.04167d0,0.125d0,2*0.3d0,0.999d0,0.001d0,1d0,
1061  &0.516d0,0.483d0,0.001d0,0.88d0,2*0.06d0,0.88d0,2*0.06d0,0.88d0,
1062  &2*0.06d0,0.2d0,0.35d0,2*0.175d0,0.0375d0,0.0125d0,0.05d0,2*0.1d0,
1063  &0.35d0,2*0.175d0,4*0.025d0,0.2d0,0.35d0,2*0.175d0,0.0375d0,
1064  &0.0125d0,0.05d0,0.1d0,0.15d0,2*0.225d0,3*0.1d0,2*0.05d0,0.15d0,
1065  &2*0.225d0,2*0.05d0,3*0.1d0,0.15d0,2*0.225d0,3*0.1d0,0.2d0,0.05d0,
1066  &2*0.04d0,0.52d0,0.1d0,0.0375d0,0.0125d0,2*0.1d0,0.05d0,2*0.04d0,
1067  &0.52d0,0.1d0,2*0.025d0,0.2d0,0.05d0,2*0.04d0,0.52d0,0.1d0,
1068  &0.0375d0,0.0125d0,0.4d0,0.2d0,2*0.02d0,0.13d0,0.23d0,3*0.2d0,
1069  &2*0.02d0,2*0.065d0,0.23d0,0.4d0,0.2d0,2*0.02d0,0.13d0,0.23d0,
1070  &0.05d0,0.2d0,2*0.015d0,0.11d0,0.09d0,0.15d0,0.1125d0,0.0375d0,
1071  &0.22d0,2*0.025d0,0.2d0,2*0.015d0,0.11d0,2*0.045d0,4*0.075d0,
1072  &2*0.11d0,0.05d0,0.2d0,2*0.015d0,0.11d0,0.09d0,0.15d0,0.1125d0,
1073  &0.0375d0,0.22d0,2*0.05d0,2*0.01d0,0.1d0,0.18d0,0.2d0,0.15d0/
1074  data (brat(i) ,i=2098,2295)/0.05d0,0.2d0,2*0.025d0,0.05d0,
1075  &2*0.01d0,0.1d0,2*0.09d0,0.2d0,4*0.1d0,2*0.05d0,2*0.01d0,0.1d0,
1076  &0.18d0,0.2d0,0.15d0,0.05d0,3*0.2d0,2*0.03d0,0.06d0,0.07d0,0.13d0,
1077  &0.15d0,0.05d0,0.08d0,2*0.1d0,0.2d0,4*0.03d0,2*0.035d0,0.13d0,
1078  &2*0.1d0,2*0.04d0,2*0.2d0,2*0.03d0,0.06d0,0.07d0,0.13d0,0.15d0,
1079  &0.05d0,0.08d0,0.9988d0,0.0001d0,0.0006d0,0.0004d0,0.0001d0,
1080  &0.9954d0,0.0011d0,0.0035d0,0.667d0,2*0.333d0,0.667d0,0.03333d0,
1081  &0.06667d0,2*0.2d0,0.4d0,0.03333d0,0.06667d0,0.03333d0,0.06667d0,
1082  &2*0.2d0,0.4d0,0.03333d0,0.06667d0,0.03333d0,0.06667d0,0.3d0,
1083  &0.1d0,0.2d0,0.1d0,0.2d0,0.03333d0,0.06667d0,0.3d0,0.1d0,0.2d0,
1084  &0.1d0,0.2d0,0.1d0,0.2d0,0.3d0,0.1d0,0.2d0,0.03333d0,0.06667d0,
1085  &0.1d0,0.2d0,0.3d0,0.1d0,0.2d0,0.03333d0,0.06667d0,0.2d0,
1086  &0.26667d0,0.53333d0,0.2d0,0.26667d0,0.53333d0,0.676d0,0.234d0,
1087  &0.085d0,0.005d0,2*1d0,0.018d0,2*0.005d0,0.003d0,0.002d0,
1088  &2*0.006d0,0.018d0,2*0.005d0,0.003d0,0.002d0,2*0.006d0,0.0066d0,
1089  &0.025d0,0.016d0,0.0088d0,2*0.005d0,0.0058d0,0.005d0,0.0055d0,
1090  &4*0.004d0,2*0.002d0,2*0.004d0,0.003d0,0.002d0,2*0.003d0,
1091  &3*0.002d0,2*0.001d0,0.002d0,2*0.001d0,2*0.002d0,0.0013d0,
1092  &0.0018d0,5*0.001d0,4*0.003d0,2*0.005d0,2*0.002d0,2*0.001d0,
1093  &2*0.002d0,2*0.001d0,0.2432d0,0.057d0,2*0.035d0,0.15d0,2*0.075d0/
1094  data (brat(i) ,i=2296,2496)/0.03d0,2*0.015d0,2*0.08d0,0.76d0,
1095  &0.08d0,4*1d0,2*0.08d0,0.76d0,0.08d0,1d0,2*0.5d0,1d0,2*0.5d0,
1096  &2*0.08d0,0.76d0,0.08d0,1d0,2*0.08d0,0.76d0,3*0.08d0,0.76d0,
1097  &3*0.08d0,0.76d0,3*0.08d0,0.76d0,3*0.08d0,0.76d0,3*0.08d0,0.76d0,
1098  &3*0.08d0,0.76d0,0.08d0,2*1d0,2*0.105d0,0.04d0,0.0077d0,0.02d0,
1099  &0.0235d0,0.0285d0,0.0435d0,0.0011d0,0.0022d0,0.0044d0,0.4291d0,
1100  &0.08d0,0.07d0,0.02d0,0.015d0,0.005d0,2*0.105d0,0.04d0,0.5d0,
1101  &0.08d0,0.14d0,0.01d0,0.015d0,0.005d0,2*0.105d0,0.04d0,0.5d0,
1102  &0.08d0,0.14d0,0.01d0,0.015d0,0.005d0,4*1d0,2*0.105d0,0.04d0,
1103  &0.5d0,0.08d0,0.14d0,0.01d0,0.015d0,0.005d0,2*0.105d0,0.04d0,
1104  &0.5d0,0.08d0,0.14d0,0.01d0,0.015d0,0.005d0,4*1d0,2*0.105d0,
1105  &0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,0.015d0,0.005d0,1d0,2*0.105d0,
1106  &0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,0.015d0,0.005d0,2*0.105d0,
1107  &0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,0.015d0,0.005d0,2*0.105d0,
1108  &0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,0.015d0,0.005d0,2*0.105d0,
1109  &0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,0.015d0,0.005d0,2*0.105d0,
1110  &0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,0.015d0,0.005d0,2*0.105d0,
1111  &0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,0.015d0,0.005d0,2*0.105d0,
1112  &0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,0.015d0,0.005d0,2*0.105d0,
1113  &0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,0.015d0,0.005d0,2*0.105d0/
1114  data (brat(i) ,i=2497,4000)/0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,
1115  &0.015d0,0.005d0,2*0.105d0,0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,
1116  &0.015d0,0.005d0,2*0.105d0,0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,
1117  &0.015d0,0.005d0,2*0.105d0,0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,
1118  &0.015d0,0.005d0,2*0.105d0,0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,
1119  &0.015d0,0.005d0,2*0.105d0,0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,
1120  &0.015d0,0.005d0,2*0.105d0,0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,
1121  &0.015d0,0.005d0,2*0.105d0,0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,
1122  &0.015d0,0.005d0,2*0.105d0,0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,
1123  &0.015d0,0.005d0,2*0.105d0,0.04d0,0.5d0,0.08d0,0.14d0,0.01d0,
1124  &0.015d0,0.005d0,827*0d0,0.8516d0,0.00539d0,0.04483d0,0.09819d0,
1125  &0.85053d0,0.02152d0,0.02989d0,0.09806d0,0.29439d0,0.10943d0,
1126  &0.59618d0,0.38983d0,0.61017d0,576*0d0/
1127  data (kfdp(i,1),i= 1, 375)/21,22,23,4*-24,25,21,22,23,4*24,25,
1128  &21,22,23,4*-24,25,21,22,23,4*24,25,21,22,23,4*-24,25,21,22,23,
1129  &4*24,25,37,1000022,1000023,1000025,1000035,21,22,23,4*-24,25,
1130  &2*-37,21,22,23,4*24,25,2*37,22,23,-24,25,23,24,-12,22,23,-24,25,
1131  &23,24,-12,-14,48*16,22,23,-24,25,23,24,22,23,-24,25,-37,23,24,37,
1132  &1,2,3,4,5,6,7,8,21,1,2,3,4,5,6,7,8,11,13,15,17,1,2,3,4,5,6,7,8,
1133  &11,12,13,14,15,16,17,18,4*-1,4*-3,4*-5,4*-7,-11,-13,-15,-17,1,2,
1134  &3,4,5,6,7,8,11,13,15,17,21,2*22,23,24,1000022,2*1000023,
1135  &3*1000025,4*1000035,2*1000024,2*1000037,1000001,2000001,1000001,
1136  &-1000001,1000002,2000002,1000002,-1000002,1000003,2000003,
1137  &1000003,-1000003,1000004,2000004,1000004,-1000004,1000005,
1138  &2000005,1000005,-1000005,1000006,2000006,1000006,-1000006,
1139  &1000011,2000011,1000011,-1000011,1000012,2000012,1000012,
1140  &-1000012,1000013,2000013,1000013,-1000013,1000014,2000014,
1141  &1000014,-1000014,1000015,2000015,1000015,-1000015,1000016,
1142  &2000016,1000016,-1000016,1,2,3,4,5,6,7,8,11,12,13,14,15,16,17,18,
1143  &24,37,2*23,25,35,4*-1,4*-3,4*-5,4*-7,-11,-13,-15,-17,3*24,1,2,3,
1144  &4,5,6,7,8,11,13,15,17,21,2*22,23,24,23,25,36,1000022,2*1000023,
1145  &3*1000025,4*1000035,2*1000024,2*1000037,1000001,2000001,1000001,
1146  &-1000001,1000002,2000002,1000002,-1000002,1000003,2000003/
1147  data (kfdp(i,1),i= 376, 606)/1000003,-1000003,1000004,2000004,
1148  &1000004,-1000004,1000005,2000005,1000005,-1000005,1000006,
1149  &2000006,1000006,-1000006,1000011,2000011,1000011,-1000011,
1150  &1000012,2000012,1000012,-1000012,1000013,2000013,1000013,
1151  &-1000013,1000014,2000014,1000014,-1000014,1000015,2000015,
1152  &1000015,-1000015,1000016,2000016,1000016,-1000016,1,2,3,4,5,6,7,
1153  &8,11,13,15,17,21,2*22,23,24,23,1000022,2*1000023,3*1000025,
1154  &4*1000035,2*1000024,2*1000037,1000001,2000001,1000001,-1000001,
1155  &1000002,2000002,1000002,-1000002,1000003,2000003,1000003,
1156  &-1000003,1000004,2000004,1000004,-1000004,1000005,2000005,
1157  &1000005,-1000005,1000006,2000006,1000006,-1000006,1000011,
1158  &2000011,1000011,-1000011,1000012,2000012,1000012,-1000012,
1159  &1000013,2000013,1000013,-1000013,1000014,2000014,1000014,
1160  &-1000014,1000015,2000015,1000015,-1000015,1000016,2000016,
1161  &1000016,-1000016,-1,-3,-5,-7,-11,-13,-15,-17,24,2*1000022,
1162  &2*1000023,2*1000025,2*1000035,1000006,2000006,1000006,2000006,
1163  &-1000001,-1000003,-1000011,-1000013,-1000015,-2000015,5,6,21,2,1,
1164  &2,3,4,5,6,11,13,15,4,5,11,13,15,2*4,-11,-13,-15,2*24,2*52,1,2,3,
1165  &4,5,6,7,8,11,12,13,14,15,16,17,18,2*24,2*52,4*-1,4*-3,4*-5,4*-7,
1166  &-11,-13,-15,-17,22,23,1,2,3,4,5,6,7,8,11,12,13,14,15,16,17,18,82/
1167  data (kfdp(i,1),i= 607, 956)/-11,-13,2*2,-12,-14,-16,2*-2,2*-4,
1168  &-2,-4,2*22,-13,-11,211,111,221,13,11,2*211,221,321,221,321,213,
1169  &-213,113,2*213,-213,221,223,321,130,310,111,331,213,113,221,223,
1170  &321,211,331,2*223,213,-213,10220,213,113,10220,225,213,-213,
1171  &30221,321,-321,221,225,213,113,30221,321,-311,221,3*211,221,223,
1172  &333,3*211,221,223,333,113,213,211,111,-211,22,111,211,2*22,211,
1173  &22,111,211,22,211,221,111,11,211,111,2*211,321,130,310,221,111,
1174  &211,111,3*22,111,321,130,-213,113,213,211,22,111,11,13,211,321,
1175  &130,310,221,211,111,211,111,321,130,310,213,113,-213,211,111,321,
1176  &130,310,2*211,313,-313,323,-323,10111,10211,-10211,221,10211,
1177  &10111,-10211,2*221,213,2*113,-213,2*321,2*311,113,313,-313,323,
1178  &-323,211,111,-321,130,310,213,-213,113,321,-321,323,-323,213,
1179  &-213,113,111,211,-12,12,-14,14,211,111,22,211,111,130,310,-13,
1180  &3*211,12,14,321,3*311,3*321,311,323,313,2*311,321,313,323,2*321,
1181  &311,323,313,323,313,321,4*311,321,313,323,313,323,311,5*321,
1182  &2*311,321,323,2*313,2*323,313,321,311,321,311,313,323,311,321,
1183  &311,2*321,311,321,311,323,313,311,321,311,321,313,323,321,311,
1184  &323,313,321,2*311,325,315,311,321,313,323,311,2*321,315,2*325,
1185  &315,323,313,3*311,315,325,313,323,3*321,325,315,323,313,3*311,
1186  &315,325,313,323,3*321,11*-11,11*-13,-311,-313,-311,-313,-20313/
1187  data (kfdp(i,1),i= 957,1368)/2*-311,-313,-311,-313,2*111,2*221,
1188  &2*331,2*113,2*223,2*333,-311,-313,2*-321,211,-311,-321,333,-311,
1189  &-313,-321,211,2*-321,2*-311,-321,211,113,421,2*411,421,411,423,
1190  &413,423,413,421,411,8*-11,8*-13,-321,-323,-321,-323,-311,2*-313,
1191  &-311,-313,2*-311,-321,-10323,-321,-323,-321,-311,2*-313,211,111,
1192  &333,3*-321,-311,-313,-321,-313,310,333,211,2*-321,-311,-313,-311,
1193  &211,-321,3*-311,211,113,321,2*421,411,421,413,423,413,423,411,
1194  &421,-15,5*-11,5*-13,221,331,333,221,331,333,10221,211,213,211,
1195  &213,321,323,321,323,2212,221,331,333,221,2*2,2*431,421,411,423,
1196  &413,82,11,13,82,443,82,421,411,423,413,411,421,413,423,411,421,
1197  &423,413,443,2*82,423,2*413,2*423,413,443,82,11,13,82,4*443,10441,
1198  &20443,445,441,6*12,6*14,2*16,3*-411,3*-413,2*-411,2*-413,2*441,
1199  &2*443,2*20443,2*2,2*4,2,4,511,521,511,523,513,523,513,521,511,
1200  &6*12,6*14,2*16,3*-421,3*-423,2*-421,2*-423,2*441,2*443,2*20443,
1201  &2*2,2*4,2,4,521,511,521,513,523,513,523,511,521,6*12,6*14,2*16,
1202  &3*-431,3*-433,2*-431,2*-433,3*441,3*443,3*20443,2*2,2*4,2,4,531,
1203  &521,511,523,513,16,2*4,2*12,2*14,2*16,4*2,4*4,2*-11,2*-13,2*-1,
1204  &2*-3,2*-11,2*-13,2*-1,541,511,521,513,523,21,11,13,15,1,2,3,4,21,
1205  &22,553,21,521,511,523,513,511,521,513,523,521,511,523,513,511,
1206  &521,513,523,553,2*21,523,2*513,2*523,2*513,523,553,21,11,13,15,1/
1207  data (kfdp(i,1),i=1369,1633)/2,3,4,21,22,2*553,10551,20553,555,
1208  &2*2112,2212,1114,2114,2214,2112,3*2212,2112,2224,2214,2114,2212,
1209  &3*2112,2212,1114,2114,2214,2112,2212,2112,2*2212,2112,2224,2214,
1210  &2114,2212,2112,2212,2*2112,2212,1114,2114,2214,2112,2212,2112,
1211  &2*2212,2112,2224,2214,2114,2212,2112,2212,2*2112,2212,1114,2114,
1212  &2214,2112,2212,2112,3122,2*2212,2112,2224,2214,2114,2212,2112,
1213  &2212,3122,2*2112,2212,1114,2114,2214,2112,2212,2112,3122,2*2212,
1214  &2112,2224,2214,2114,2212,2112,2212,3122,2*2112,2212,1114,2114,
1215  &2214,2112,3*2212,2112,2224,2214,2114,2212,3*2112,2212,1114,2114,
1216  &2214,2112,2212,3122,2*2212,2112,2224,2214,2114,2212,2112,3122,
1217  &2*2112,2212,1114,2114,2214,2112,2212,2112,3122,3212,3112,2*2212,
1218  &2112,2224,2214,2114,2212,2112,2212,3122,3212,3222,2*2112,2212,
1219  &1114,2114,2214,2112,2212,2112,3122,3212,3112,2112,2*2212,2112,
1220  &2224,2214,2114,2212,2112,2212,3122,3212,3222,2212,2*2112,2212,
1221  &1114,2114,2214,2112,2212,2112,3122,3212,3112,2*2112,2*2212,2112,
1222  &2224,2214,2114,2212,2112,2212,3122,3212,3222,2*2212,2112,2212,
1223  &2*2112,2212,2112,2*2212,12112,2112,1114,2114,2112,12112,12212,
1224  &2112,2212,2114,2214,1114,2112,2212,12212,12112,2212,2112,2214,
1225  &2114,2224,2212,2112,12212,2212,2224,2214,2212,12112,2112,1114,
1226  &2114,2112,12112,12212,2112,2212,2114,2214,1114,2112,2212,12212/
1227  data (kfdp(i,1),i=1634,1879)/12112,2212,2112,2214,2114,2224,2212,
1228  &2112,12212,2212,2224,2214,2212,2112,1114,2114,2112,3112,2112,
1229  &2212,2114,2214,1114,2112,2212,3212,3112,2212,2112,2214,2114,2224,
1230  &2212,2112,3212,3222,2212,2224,2214,2212,3222,12112,2112,1114,
1231  &2114,2112,12112,12212,2112,2212,2114,2214,1114,2112,2212,12212,
1232  &12112,2212,2112,2214,2114,2224,2212,2112,12212,2212,2224,2214,
1233  &2212,2112,1114,2114,2112,3112,2112,2212,2114,2214,1114,2112,2212,
1234  &3212,3112,2212,2112,2214,2114,2224,2212,2112,3212,3222,2212,2224,
1235  &2214,2212,3222,12112,2112,1114,2114,2112,12112,12212,2112,2212,
1236  &2114,2214,1114,2112,2212,12212,12112,2212,2112,2214,2114,2224,
1237  &2212,2112,12212,2212,2224,2214,2212,12112,2112,1114,2114,3112,
1238  &12112,12212,2112,2212,2114,2214,1114,3212,3112,12212,12112,2212,
1239  &2112,2214,2114,2224,3212,3222,12212,2212,2224,2214,3222,3*2112,
1240  &2212,2112,2*2212,2112,2212,2112,2*2212,12116,2112,1114,2114,2112,
1241  &3112,12116,12216,2112,2212,2114,2214,1114,2112,2212,3212,3112,
1242  &12216,12116,2212,2112,2214,2114,2224,2212,2112,3212,3222,12216,
1243  &2212,2224,2214,2212,3222,2212,2*2112,-12,3212,3222,3112,2112,
1244  &2212,3212,3222,3112,3214,3224,3114,2112,2212,3212,3222,3112,3214,
1245  &3224,3114,2112,2212,3212,3222,3112,3122,3214,3224,3114,2112,2212,
1246  &3212,3222,3112,3122,3214,3224,3114,2112,2212,3212,3222,3112,3214/
1247  data (kfdp(i,1),i=1880,2129)/3224,3114,2112,2212,2112,2212,3212,
1248  &3222,3112,3214,3224,3114,2112,2212,2112,2212,3212,3222,3112,3122,
1249  &3214,3224,3114,2112,2212,2112,2212,3212,3222,3112,3122,3214,3224,
1250  &3114,2112,2212,3212,3222,3112,3214,3224,3114,3122,2112,2212,2112,
1251  &2212,3212,3222,3112,3122,3322,3312,3214,3224,3114,3122,2112,2212,
1252  &2112,2212,3212,3222,3112,3214,3224,3114,3122,2112,2212,2112,-12,
1253  &3122,2212,2112,2212,3122,3212,3112,3122,3222,3112,3122,3222,3212,
1254  &2112,3122,3212,3112,1114,2114,3314,2112,2212,3122,3222,3112,2114,
1255  &2214,3324,3314,2212,3122,3222,3212,2224,2214,3324,2112,3122,3212,
1256  &3112,3214,13122,3124,2112,2212,3122,3222,3112,3224,3114,13122,
1257  &3124,2212,3122,3222,3212,3224,13122,3124,2112,3122,3212,2*3112,
1258  &3124,1114,2114,2112,2212,3122,3222,3112,3212,3124,2114,2214,2212,
1259  &3122,3222,3212,3222,3124,2224,2214,2112,3122,3212,3112,3214,3124,
1260  &2112,2212,3122,3222,3112,3224,3114,3124,2212,3122,3222,3212,3224,
1261  &3124,2112,3122,3212,2*3112,3312,3214,1114,2114,3314,2112,2212,
1262  &3122,3222,3112,3212,3322,3312,3224,3114,2114,2214,3324,3314,2212,
1263  &3122,3222,3212,3222,3322,3224,2224,2214,3324,2112,3122,3212,
1264  &2*3112,3214,13122,1114,2114,3314,2112,2212,3122,3222,3112,3212,
1265  &3224,3114,13122,2114,2214,3324,3314,2212,3122,3222,3212,3222,
1266  &3224,13122,2224,2214,3324,2112,3122,3212,3112,3312,3214,3124/
1267  data (kfdp(i,1),i=2130,2504)/1114,2114,3314,2112,2212,3122,3222,
1268  &3112,3322,3312,3224,3114,3124,2114,2214,3324,3314,2212,3122,3222,
1269  &3212,3322,3224,3124,2224,2214,3324,3122,3112,-12,-14,-12,2*3122,
1270  &3212,3322,3312,3322,2*3312,3322,3122,3212,3112,3314,3324,3322,
1271  &3312,3122,3212,3222,3324,3314,3312,3322,3122,3212,3112,3314,3324,
1272  &3322,3312,3122,3212,3222,3324,3314,3312,3322,3122,3212,3112,3314,
1273  &3324,3322,3312,3122,3212,3222,3324,3314,3122,3212,3112,3122,3212,
1274  &3222,3122,3322,3312,-12,2*4122,7*-11,7*-13,2*2224,2*2212,2*2214,
1275  &2*3122,2*3212,2*3214,5*3222,4*3224,2*3322,3324,2*2224,7*2212,
1276  &5*2214,2*2112,2*2114,2*3122,2*3212,2*3214,2*3222,2*3224,4*2,3,
1277  &2*2,1,2*2,-11,-13,2*2,4*4122,-11,-13,2*2,3*4132,3*4232,-11,-13,
1278  &2*2,4332,-11,-13,2*2,-11,-13,2*2,-11,-13,2*2,-11,-13,2*2,-11,-13,
1279  &2*2,-11,-13,2*2,-11,-13,2*2,2*5122,-12,-14,-16,5*4122,441,443,
1280  &20443,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,
1281  &2*-2,2*-4,-2,-4,4*5122,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,
1282  &2*-2,2*-4,-2,-4,2*5132,2*5232,-12,-14,-16,2*-2,2*-4,-2,-4,5332,
1283  &-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,
1284  &-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,
1285  &2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,
1286  &-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12/
1287  data (kfdp(i,1),i=2505,2704)/-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,
1288  &2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,
1289  &2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,
1290  &-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,
1291  &-14,-16,2*-2,2*-4,-2,-4,1000039,-1000024,-1000037,1000022,
1292  &1000023,1000025,1000035,1000002,2000002,1000002,2000002,1000021,
1293  &1000039,1000024,1000037,1000022,1000023,1000025,1000035,1000001,
1294  &2000001,1000001,2000001,1000021,1000039,-1000024,-1000037,
1295  &1000022,1000023,1000025,1000035,1000004,2000004,1000004,2000004,
1296  &1000021,1000039,1000024,1000037,1000022,1000023,1000025,1000035,
1297  &1000003,2000003,1000003,2000003,1000021,1000039,-1000024,
1298  &-1000037,1000022,1000023,1000025,1000035,1000006,2000006,1000006,
1299  &2000006,1000021,1000039,1000024,1000037,1000022,1000023,1000025,
1300  &1000035,1000005,2000005,1000005,2000005,1000021,1000022,1000039,
1301  &-1000024,-1000037,1000022,1000023,1000025,1000035,1000012,
1302  &2000012,1000012,2000012,1000039,1000024,1000037,1000022,1000023,
1303  &1000025,1000035,1000011,2000011,1000011,2000011,1000039,-1000024,
1304  &-1000037,1000022,1000023,1000025,1000035,1000014,2000014,1000014,
1305  &2000014,1000039,1000024,1000037,1000022,1000023,1000025,1000035,
1306  &1000013,2000013,1000013,2000013,1000039,-1000024,-1000037/
1307  data (kfdp(i,1),i=2705,2886)/1000022,1000023,1000025,1000035,
1308  &1000016,2000016,1000016,2000016,1000039,1000024,1000037,1000022,
1309  &1000023,1000025,1000035,1000015,2000015,1000015,2000015,1000039,
1310  &1000001,-1000001,2000001,-2000001,1000002,-1000002,2000002,
1311  &-2000002,1000003,-1000003,2000003,-2000003,1000004,-1000004,
1312  &2000004,-2000004,1000005,-1000005,2000005,-2000005,1000006,
1313  &-1000006,2000006,-2000006,6*1000022,6*1000023,6*1000025,
1314  &6*1000035,1000024,-1000024,1000024,-1000024,1000024,-1000024,
1315  &1000037,-1000037,1000037,-1000037,1000037,-1000037,10*1000039,
1316  &16*1000022,1000024,-1000024,1000024,-1000024,1000024,-1000024,
1317  &1000024,-1000024,1000024,-1000024,1000024,-1000024,1000037,
1318  &-1000037,1000037,-1000037,1000037,-1000037,1000037,-1000037,
1319  &1000037,-1000037,1000037,-1000037,1000024,-1000024,1000037,
1320  &-1000037,1000001,-1000001,2000001,-2000001,1000002,-1000002,
1321  &2000002,-2000002,1000003,-1000003,2000003,-2000003,1000004,
1322  &-1000004,2000004,-2000004,1000005,-1000005,2000005,-2000005,
1323  &1000006,-1000006,2000006,-2000006,1000011,-1000011,2000011,
1324  &-2000011,1000012,-1000012,2000012,-2000012,1000013,-1000013,
1325  &2000013,-2000013,1000014,-1000014,2000014,-2000014,1000015,
1326  &-1000015,2000015,-2000015,1000016,-1000016,2000016,-2000016/
1327  data (kfdp(i,1),i=2887,3133)/5*1000021,2*1000039,6*1000022,
1328  &6*1000023,6*1000025,6*1000035,1000022,1000023,1000025,1000035,
1329  &1000002,2000002,-1000001,-2000001,1000004,2000004,-1000003,
1330  &-2000003,1000006,2000006,-1000005,-2000005,1000012,2000012,
1331  &-1000011,-2000011,1000014,2000014,-1000013,-2000013,1000016,
1332  &2000016,-1000015,-2000015,2*1000021,5*1000039,16*1000022,
1333  &16*1000023,1000024,-1000024,1000024,-1000024,1000024,-1000024,
1334  &1000024,-1000024,1000024,-1000024,1000024,-1000024,1000037,
1335  &-1000037,1000037,-1000037,1000037,-1000037,1000037,-1000037,
1336  &1000037,-1000037,1000037,-1000037,1000024,-1000024,1000037,
1337  &-1000037,1000001,-1000001,2000001,-2000001,1000002,-1000002,
1338  &2000002,-2000002,1000003,-1000003,2000003,-2000003,1000004,
1339  &-1000004,2000004,-2000004,1000005,-1000005,2000005,-2000005,
1340  &1000006,-1000006,2000006,-2000006,1000011,-1000011,2000011,
1341  &-2000011,1000012,-1000012,2000012,-2000012,1000013,-1000013,
1342  &2000013,-2000013,1000014,-1000014,2000014,-2000014,1000015,
1343  &-1000015,2000015,-2000015,1000016,-1000016,2000016,-2000016,
1344  &5*1000021,5*1000039,16*1000022,16*1000023,16*1000025,1000024,
1345  &-1000024,1000024,-1000024,1000024,-1000024,1000024,-1000024,
1346  &1000024,-1000024,1000024,-1000024,1000037,-1000037,1000037/
1347  data (kfdp(i,1),i=3134,3319)/-1000037,1000037,-1000037,1000037,
1348  &-1000037,1000037,-1000037,1000037,-1000037,1000024,-1000024,
1349  &1000037,-1000037,1000001,-1000001,2000001,-2000001,1000002,
1350  &-1000002,2000002,-2000002,1000003,-1000003,2000003,-2000003,
1351  &1000004,-1000004,2000004,-2000004,1000005,-1000005,2000005,
1352  &-2000005,1000006,-1000006,2000006,-2000006,1000011,-1000011,
1353  &2000011,-2000011,1000012,-1000012,2000012,-2000012,1000013,
1354  &-1000013,2000013,-2000013,1000014,-1000014,2000014,-2000014,
1355  &1000015,-1000015,2000015,-2000015,1000016,-1000016,2000016,
1356  &-2000016,5*1000021,2*1000039,15*1000024,6*1000022,6*1000023,
1357  &6*1000025,6*1000035,1000022,1000023,1000025,1000035,1000002,
1358  &2000002,-1000001,-2000001,1000004,2000004,-1000003,-2000003,
1359  &1000006,2000006,-1000005,-2000005,1000012,2000012,-1000011,
1360  &-2000011,1000014,2000014,-1000013,-2000013,1000016,2000016,
1361  &-1000015,-2000015,2*1000021,1000039,-1000024,-1000037,1000022,
1362  &1000023,1000025,1000035,4*1000001,1000002,2000002,1000002,
1363  &2000002,1000021,1000039,1000024,1000037,1000022,1000023,1000025,
1364  &1000035,4*1000002,1000001,2000001,1000001,2000001,1000021,
1365  &1000039,-1000024,-1000037,1000022,1000023,1000025,1000035,
1366  &4*1000003,1000004,2000004,1000004,2000004,1000021,1000039/
1367  data (kfdp(i,1),i=3320,4000)/1000024,1000037,1000022,1000023,
1368  &1000025,1000035,4*1000004,1000003,2000003,1000003,2000003,
1369  &1000021,1000039,-1000024,-1000037,1000022,1000023,1000025,
1370  &1000035,4*1000005,1000006,2000006,1000006,2000006,1000021,
1371  &1000039,1000024,1000037,1000022,1000023,1000025,1000035,
1372  &4*1000006,1000005,2000005,1000005,2000005,1000021,1000039,
1373  &-1000024,-1000037,1000022,1000023,1000025,1000035,4*1000011,
1374  &1000012,2000012,1000012,2000012,1000039,-1000024,-1000037,
1375  &1000022,1000023,1000025,1000035,4*1000013,1000014,2000014,
1376  &1000014,2000014,1000039,-1000024,-1000037,1000022,1000023,
1377  &1000025,1000035,4*1000015,1000016,2000016,1000016,2000016,21,22,
1378  &23,-24,21,22,23,24,22,23,-24,23,24,576*0/
1379  data (kfdp(i,2),i= 1, 337)/3*1,2,4,6,8,1,3*2,1,3,5,7,2,3*3,2,4,
1380  &6,8,3,3*4,1,3,5,7,4,3*5,2,4,6,8,5,3*6,1,3,5,7,6,5,4*1000006,3*7,
1381  &2,4,6,8,7,4,6,3*8,1,3,5,7,8,5,7,2*11,12,11,12,2*11,2*13,14,13,14,
1382  &13,11,13,-211,-213,-211,-213,-211,-213,-211,-213,2*-211,-321,
1383  &-323,-321,2*-323,3*-321,4*-211,-213,-211,-213,-211,-213,-211,
1384  &-213,-211,-213,3*-211,-213,4*-211,-323,-321,2*-211,2*-321,3*-211,
1385  &2*15,16,15,16,15,2*17,18,17,2*18,2*17,-1,-2,-3,-4,-5,-6,-7,-8,21,
1386  &-1,-2,-3,-4,-5,-6,-7,-8,-11,-13,-15,-17,-1,-2,-3,-4,-5,-6,-7,-8,
1387  &-11,-12,-13,-14,-15,-16,-17,-18,2,4,6,8,2,4,6,8,2,4,6,8,2,4,6,8,
1388  &12,14,16,18,-1,-2,-3,-4,-5,-6,-7,-8,-11,-13,-15,-17,21,22,2*23,
1389  &-24,2*1000022,1000023,1000022,1000023,1000025,1000022,1000023,
1390  &1000025,1000035,-1000024,-1000037,-1000024,-1000037,-1000001,
1391  &2*-2000001,2000001,-1000002,2*-2000002,2000002,-1000003,
1392  &2*-2000003,2000003,-1000004,2*-2000004,2000004,-1000005,
1393  &2*-2000005,2000005,-1000006,2*-2000006,2000006,-1000011,
1394  &2*-2000011,2000011,-1000012,2*-2000012,2000012,-1000013,
1395  &2*-2000013,2000013,-1000014,2*-2000014,2000014,-1000015,
1396  &2*-2000015,2000015,-1000016,2*-2000016,2000016,-1,-2,-3,-4,-5,-6,
1397  &-7,-8,-11,-12,-13,-14,-15,-16,-17,-18,-24,-37,22,25,2*36,2,4,6,8,
1398  &2,4,6,8,2,4,6,8,2,4,6,8,12,14,16,18,23,22,25,-1,-2,-3,-4,-5,-6/
1399  data (kfdp(i,2),i= 338, 524)/-7,-8,-11,-13,-15,-17,21,22,2*23,
1400  &-24,2*25,36,2*1000022,1000023,1000022,1000023,1000025,1000022,
1401  &1000023,1000025,1000035,-1000024,-1000037,-1000024,-1000037,
1402  &-1000001,2*-2000001,2000001,-1000002,2*-2000002,2000002,-1000003,
1403  &2*-2000003,2000003,-1000004,2*-2000004,2000004,-1000005,
1404  &2*-2000005,2000005,-1000006,2*-2000006,2000006,-1000011,
1405  &2*-2000011,2000011,-1000012,2*-2000012,2000012,-1000013,
1406  &2*-2000013,2000013,-1000014,2*-2000014,2000014,-1000015,
1407  &2*-2000015,2000015,-1000016,2*-2000016,2000016,-1,-2,-3,-4,-5,-6,
1408  &-7,-8,-11,-13,-15,-17,21,22,2*23,-24,25,2*1000022,1000023,
1409  &1000022,1000023,1000025,1000022,1000023,1000025,1000035,-1000024,
1410  &-1000037,-1000024,-1000037,-1000001,2*-2000001,2000001,-1000002,
1411  &2*-2000002,2000002,-1000003,2*-2000003,2000003,-1000004,
1412  &2*-2000004,2000004,-1000005,2*-2000005,2000005,-1000006,
1413  &2*-2000006,2000006,-1000011,2*-2000011,2000011,-1000012,
1414  &2*-2000012,2000012,-1000013,2*-2000013,2000013,-1000014,
1415  &2*-2000014,2000014,-1000015,2*-2000015,2000015,-1000016,
1416  &2*-2000016,2000016,2,4,6,8,12,14,16,18,25,1000024,1000037,
1417  &1000024,1000037,1000024,1000037,1000024,1000037,2*-1000005,
1418  &2*-2000005,1000002,1000004,1000012,1000014,2*1000016,-5,-6,21,11/
1419  data (kfdp(i,2),i= 525, 860)/-3,-4,-5,-6,-7,-8,-13,-15,-17,-4,-5,
1420  &-11,-13,-15,-5,-3,12,14,16,-24,-52,-24,-52,-1,-2,-3,-4,-5,-6,-7,
1421  &-8,-11,-12,-13,-14,-15,-16,-17,-18,23,51,23,51,2,4,6,8,2,4,6,8,2,
1422  &4,6,8,2,4,6,8,12,14,16,18,2*51,-1,-2,-3,-4,-5,-6,-7,-8,-11,-12,
1423  &-13,-14,-15,-16,-17,-18,-82,12,14,-1,-3,11,13,15,1,4,3,4,1,3,22,
1424  &11,14,12,-211,2*22,-13,-11,111,22,111,-321,211,-311,-211,2*211,
1425  &111,-211,211,111,211,-321,130,310,22,2*111,3*211,-311,22,211,111,
1426  &211,-211,211,2*111,2*211,111,-211,211,111,-323,323,111,211,111,
1427  &2*211,-313,323,211,3*-211,113,4*111,-211,213,2*211,-211,111,-211,
1428  &111,211,22,111,-211,211,11,-211,13,22,-211,111,-211,22,111,-11,
1429  &-211,111,2*-211,-321,130,310,221,111,-211,111,113,223,22,111,
1430  &-321,310,211,111,2*-211,221,22,-11,-13,-211,-321,130,310,221,
1431  &-211,111,-211,111,-321,130,310,-211,111,211,-211,111,-321,130,
1432  &310,2*-211,-311,311,-321,321,111,-211,211,113,-211,111,2*211,111,
1433  &-211,211,111,211,-321,2*-311,-321,22,-311,311,-321,321,-211,111,
1434  &321,130,310,-211,211,111,-323,323,-321,321,-211,211,2*111,-211,
1435  &11,-11,13,-13,-211,111,22,-211,111,2*0,14,111,211,111,-11,-13,
1436  &-211,111,22,211,111,22,-213,113,-211,111,223,213,113,211,111,223,
1437  &-211,111,-211,111,-211,211,-213,113,223,221,211,111,211,111,
1438  &2*211,213,113,223,221,22,-211,111,211,111,-211,111,211,111,-211/
1439  data (kfdp(i,2),i= 861,1266)/111,-211,111,-213,113,211,111,211,
1440  &111,213,113,-211,111,-213,113,-211,111,211,111,213,113,211,111,
1441  &-213,113,-211,111,-211,111,221,-211,111,213,113,211,111,211,111,
1442  &221,211,111,-211,111,-211,111,225,333,223,211,111,211,111,225,
1443  &333,223,-211,111,-211,111,225,333,223,211,111,211,111,225,333,
1444  &223,11*12,11*14,2*211,2*213,211,20213,2*321,2*323,211,213,211,
1445  &213,211,213,211,213,211,213,211,213,3*211,213,211,2*321,8*211,
1446  &2*113,3*211,111,22,211,111,211,111,4*211,8*12,8*14,2*211,2*213,
1447  &2*111,221,2*113,223,333,20213,211,2*321,323,2*311,313,-211,111,
1448  &113,2*211,321,2*211,311,321,310,211,-211,4*211,321,4*211,113,
1449  &2*211,-321,111,22,-211,111,-211,111,-211,211,-211,211,16,5*12,
1450  &5*14,3*211,3*213,211,2*111,2*113,2*-311,2*-313,-2112,3*321,323,
1451  &2*-1,22,111,321,311,321,311,-82,-11,-13,-82,22,-82,211,111,211,
1452  &111,-211,111,-211,111,311,2*321,311,22,2*-82,211,111,-211,111,
1453  &321,311,22,-82,-11,-13,-82,211,111,221,111,4*22,6*-11,6*-13,
1454  &2*-15,211,213,20213,211,213,20213,431,433,431,433,311,313,311,
1455  &313,311,313,-1,-4,-3,-4,-1,-3,22,-211,111,-211,111,-211,211,-211,
1456  &211,6*-11,6*-13,2*-15,211,213,20213,211,213,20213,431,433,431,
1457  &433,321,323,321,323,321,323,-1,-4,-3,-4,-1,-3,22,211,111,211,111,
1458  &4*211,6*-11,6*-13,2*-15,211,213,20213,211,213,20213,431,433,431/
1459  data (kfdp(i,2),i=1267,1572)/433,221,331,333,221,331,333,221,331,
1460  &333,-1,-4,-3,-4,-1,-3,22,-321,-311,-321,-311,-15,-3,-1,2*-11,
1461  &2*-13,2*-15,-1,-4,-3,-4,-3,-4,-1,-4,2*12,2*14,2,3,2,3,2*12,2*14,
1462  &2,1,22,411,421,411,421,21,-11,-13,-15,-1,-2,-3,-4,2*21,22,21,
1463  &-211,111,-211,111,211,111,211,111,-321,-311,-321,-311,411,421,
1464  &411,421,22,2*21,-211,111,211,111,-321,-311,411,421,22,21,-11,-13,
1465  &-15,-1,-2,-3,-4,2*21,211,111,3*22,10220,111,-211,211,111,-211,
1466  &113,-213,10220,111,211,-211,111,211,113,213,10220,111,-211,211,
1467  &111,-211,113,-213,221,10220,111,211,-211,111,211,113,213,221,
1468  &10220,111,-211,211,111,-211,113,-213,221,10220,111,211,-211,111,
1469  &211,113,213,221,10220,111,-211,211,111,-211,113,-213,221,311,
1470  &10220,111,211,-211,111,211,113,213,221,321,10220,111,-211,211,
1471  &111,-211,113,-213,221,311,10220,111,211,-211,111,211,113,213,221,
1472  &321,10220,111,-211,211,111,-211,113,-213,10220,111,211,-211,111,
1473  &211,113,213,10220,111,-211,211,111,-211,113,-213,311,10220,111,
1474  &211,-211,111,211,113,213,321,10220,111,-211,211,111,-211,113,
1475  &-213,221,2*311,321,10220,111,211,-211,111,211,113,213,221,2*321,
1476  &311,10220,111,-211,211,111,-211,113,-213,221,2*311,321,223,10220,
1477  &111,211,-211,111,211,113,213,221,2*321,311,223,10220,111,-211,
1478  &211,111,-211,113,-213,221,2*311,321,223,331,10220,111,211,-211/
1479  data (kfdp(i,2),i=1573,1870)/111,211,113,213,221,2*321,311,223,
1480  &331,2*-211,111,22,111,211,22,211,2*-211,111,-211,-213,111,-211,
1481  &111,-211,111,-211,211,113,-213,111,211,111,211,111,211,-211,113,
1482  &213,2*211,111,211,213,2*-211,111,-211,-213,111,-211,111,-211,111,
1483  &-211,211,113,-213,111,211,111,211,111,211,-211,113,213,2*211,111,
1484  &211,213,-211,111,-211,-213,311,111,-211,111,-211,211,113,-213,
1485  &311,321,111,211,111,211,-211,113,213,321,311,211,111,211,213,321,
1486  &2*-211,111,-211,-213,111,-211,111,-211,111,-211,211,113,-213,111,
1487  &211,111,211,111,211,-211,113,213,2*211,111,211,213,-211,111,-211,
1488  &-213,311,111,-211,111,-211,211,113,-213,311,321,111,211,111,211,
1489  &-211,113,213,321,311,211,111,211,213,321,2*-211,111,-211,-213,
1490  &111,-211,111,-211,111,-211,211,113,-213,111,211,111,211,111,211,
1491  &-211,113,213,2*211,111,211,213,2*-211,111,-211,311,111,-211,111,
1492  &-211,111,-211,211,311,321,111,211,111,211,111,211,-211,321,311,
1493  &2*211,111,211,321,-211,-213,111,-211,113,-213,111,211,113,213,
1494  &211,213,2*-211,111,-211,-213,311,111,-211,111,-211,111,-211,211,
1495  &113,-213,311,321,111,211,111,211,111,211,-211,113,213,321,311,
1496  &2*211,111,211,213,321,-211,111,22,11,111,-211,211,-311,-321,111,
1497  &-211,211,111,-211,211,-311,-321,111,-211,211,111,-211,211,-311,
1498  &-321,111,-211,211,221,111,-211,211,-311,-321,111,-211,211,221/
1499  data (kfdp(i,2),i=1871,2156)/111,-211,211,-311,-321,111,-211,211,
1500  &111,-211,211,-313,-323,-311,-321,111,-211,211,111,-211,211,-313,
1501  &-323,-311,-321,111,-211,211,221,111,-211,211,-313,-323,-311,-321,
1502  &111,-211,211,221,111,-211,211,-311,-321,111,-211,211,111,-211,
1503  &211,223,-313,-323,-311,-321,111,-211,211,221,311,321,111,-211,
1504  &211,223,-313,-323,-311,-321,111,-211,211,111,-211,211,223,-313,
1505  &-323,-211,11,22,111,211,22,2*-211,2*111,-211,2*211,111,211,-321,
1506  &2*-211,111,-311,-321,311,-311,-321,111,-211,211,-311,-321,311,
1507  &321,-311,211,111,211,-321,-311,321,-321,2*-211,111,3*-211,-311,
1508  &-321,111,-211,211,-211,211,2*111,-311,211,111,211,111,2*211,-321,
1509  &2*-211,111,221,-211,-311,-321,-311,-321,111,-211,211,221,111,
1510  &-311,-321,-311,211,111,211,221,211,-321,-311,-321,2*-211,111,
1511  &2*-211,-311,-321,111,-211,211,-211,211,111,-311,211,111,211,111,
1512  &211,-321,2*-211,111,221,311,-211,-311,-321,311,-311,-321,111,
1513  &-211,211,221,311,321,-211,211,-311,-321,311,321,-311,211,111,211,
1514  &221,321,111,-321,-311,321,-321,2*-211,111,221,2*-211,-311,-321,
1515  &311,-311,-321,111,-211,211,221,-211,211,111,-311,-321,311,321,
1516  &-311,211,111,211,221,111,211,-321,-311,321,-321,2*-211,111,311,
1517  &2*-211,-311,-321,311,-311,-321,111,-211,211,311,321,-211,211,111,
1518  &-311,-321,311,321,-311,211,111,211,321,111,211,-321,-311,321/
1519  data (kfdp(i,2),i=2157,2583)/-211,22,11,13,11,111,2*22,-211,
1520  &2*111,211,111,-211,2*-321,-311,111,-211,111,211,2*-311,-321,111,
1521  &211,111,-211,2*-321,-311,111,-211,111,211,2*-311,-321,111,211,
1522  &111,-211,2*-321,-311,111,-211,111,211,2*-311,-321,111,211,2*-321,
1523  &3*-311,2*-321,-211,111,11,2*-211,7*12,7*14,-321,-323,-311,-313,
1524  &-311,-313,211,213,211,213,211,213,111,221,331,113,223,111,221,
1525  &113,223,321,323,321,-211,-213,111,221,331,113,223,333,10221,111,
1526  &221,331,113,223,211,213,211,213,321,323,321,323,321,323,311,313,
1527  &311,313,2*-1,-3,-1,2203,3201,3203,2203,2101,2103,12,14,-1,-3,
1528  &2*111,2*211,12,14,-1,-3,22,111,2*22,111,22,12,14,-1,-3,22,12,14,
1529  &-1,-3,12,14,-1,-3,12,14,-1,-3,12,14,-1,-3,12,14,-1,-3,12,14,-1,
1530  &-3,12,14,-1,-3,2*-211,11,13,15,-211,-213,-20213,-431,-433,3*3122,
1531  &1,4,3,4,1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,1,3,2*111,
1532  &2*211,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,1,3,4*22,11,13,15,1,
1533  &4,3,4,1,3,22,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,1,3,11,13,15,
1534  &1,4,3,4,1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,
1535  &4,3,4,1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,
1536  &3,4,1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,
1537  &4,1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,
1538  &1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,1/
1539  data (kfdp(i,2),i=2584,3059)/3,1,2*2,4*1,2*-24,2*-37,1,2,2*1,4*2,
1540  &2*24,2*37,2,3,2*4,4*3,2*-24,2*-37,3,4,2*3,4*4,2*24,2*37,4,5,2*6,
1541  &4*5,2*-24,2*-37,5,6,2*5,4*6,2*24,2*37,6,4,11,2*12,4*11,2*-24,
1542  &2*-37,12,2*11,4*12,2*24,2*37,13,2*14,4*13,2*-24,2*-37,14,2*13,
1543  &4*14,2*24,2*37,15,2*16,4*15,2*-24,2*-37,16,2*15,4*16,2*24,2*37,
1544  &21,-1,1,-1,1,-2,2,-2,2,-3,3,-3,3,-4,4,-4,4,-5,5,-5,5,-6,6,-6,6,1,
1545  &3,5,2,4,6,1,3,5,2,4,6,1,3,5,2,4,6,1,3,5,2,4,6,1,-1,3,-3,5,-5,1,
1546  &-1,3,-3,5,-5,22,23,25,35,36,22,23,25,35,36,22,23,11,13,15,12,14,
1547  &16,1,3,5,2,4,25,35,36,-24,24,11,-11,13,-13,15,-15,1,-1,3,-3,-24,
1548  &24,11,-11,13,-13,15,-15,1,-1,3,-3,-37,37,-37,37,-1,1,-1,1,-2,2,
1549  &-2,2,-3,3,-3,3,-4,4,-4,4,-5,5,-5,5,-6,6,-6,6,-11,11,-11,11,-12,
1550  &12,-12,12,-13,13,-13,13,-14,14,-14,14,-15,15,-15,15,-16,16,-16,
1551  &16,1,3,5,2,4,24,37,24,-11,-13,-15,-1,-3,24,-11,-13,-15,-1,-3,24,
1552  &-11,-13,-15,-1,-3,24,-11,-13,-15,-1,-3,4*37,2*-1,2*2,2*-3,2*4,
1553  &2*-5,2*6,2*-11,2*12,2*-13,2*14,2*-15,2*16,-1,-3,22,23,25,35,36,
1554  &22,23,11,13,15,12,14,16,1,3,5,2,4,25,35,36,22,23,11,13,15,12,14,
1555  &16,1,3,5,2,4,25,35,36,-24,24,11,-11,13,-13,15,-15,1,-1,3,-3,-24,
1556  &24,11,-11,13,-13,15,-15,1,-1,3,-3,-37,37,-37,37,-1,1,-1,1,-2,2,
1557  &-2,2,-3,3,-3,3,-4,4,-4,4,-5,5,-5,5,-6,6,-6,6,-11,11,-11,11,-12,
1558  &12,-12,12,-13,13,-13,13,-14,14,-14,14,-15,15,-15,15,-16,16,-16/
1559  data (kfdp(i,2),i=3060,4000)/16,1,3,5,2,4,22,23,25,35,36,22,23,
1560  &11,13,15,12,14,16,1,3,5,2,4,25,35,36,22,23,11,13,15,12,14,16,1,3,
1561  &5,2,4,25,35,36,22,23,11,13,15,12,14,16,1,3,5,2,4,25,35,36,-24,24,
1562  &11,-11,13,-13,15,-15,1,-1,3,-3,-24,24,11,-11,13,-13,15,-15,1,-1,
1563  &3,-3,-37,37,-37,37,-1,1,-1,1,-2,2,-2,2,-3,3,-3,3,-4,4,-4,4,-5,5,
1564  &-5,5,-6,6,-6,6,-11,11,-11,11,-12,12,-12,12,-13,13,-13,13,-14,14,
1565  &-14,14,-15,15,-15,15,-16,16,-16,16,1,3,5,2,4,24,37,23,11,13,15,
1566  &12,14,16,1,3,5,2,4,25,35,36,24,-11,-13,-15,-1,-3,24,-11,-13,-15,
1567  &-1,-3,24,-11,-13,-15,-1,-3,24,-11,-13,-15,-1,-3,4*37,2*-1,2*2,
1568  &2*-3,2*4,2*-5,2*6,2*-11,2*12,2*-13,2*14,2*-15,2*16,-1,-3,1,2*2,
1569  &4*1,23,25,35,36,2*-24,2*-37,1,2,2*1,4*2,23,25,35,36,2*24,2*37,2,
1570  &3,2*4,4*3,23,25,35,36,2*-24,2*-37,3,4,2*3,4*4,23,25,35,36,2*24,
1571  &2*37,4,5,2*6,4*5,23,25,35,36,2*-24,2*-37,5,6,2*5,4*6,23,25,35,36,
1572  &2*24,2*37,6,11,2*12,4*11,23,25,35,36,2*-24,2*-37,13,2*14,4*13,23,
1573  &25,35,36,2*-24,2*-37,15,2*16,4*15,23,25,35,36,2*-24,2*-37,3*1,
1574  &4*2,1,2*11,2*12,11,576*0/
1575  data (kfdp(i,3),i= 1,1276)/79*0,14,6*0,2*16,2*0,6*111,310,130,
1576  &2*0,3*111,310,130,321,113,211,223,221,2*113,2*211,2*223,2*221,
1577  &2*113,221,2*113,2*213,-213,113,2*111,310,130,310,130,2*310,130,
1578  &470*0,4*3,4*4,1,4,3,2*2,0,-11,20*0,-211,8*0,111,26*0,111,211,4*0,
1579  &111,211,3*0,211,111,4*0,2*111,-211,-11,11,-13,22,111,3*0,22,3*0,
1580  &111,211,4*0,111,2*221,3*0,111,5*0,111,24*0,211,111,11*0,-211,
1581  &2*111,-211,3*111,-211,111,211,20*0,2*111,211,-211,211,-211,9*0,
1582  &-211,3*111,20*0,111,-211,8*0,111,-211,83*0,-311,-313,-311,-321,
1583  &-313,-323,111,221,331,113,223,-311,-313,-311,-321,-313,-323,111,
1584  &221,331,113,223,22*0,111,113,2*211,-211,-311,211,111,3*211,-211,
1585  &7*211,7*0,111,-211,111,-211,-321,-323,-311,-321,-313,-323,-211,
1586  &-213,-321,-323,-311,-321,-313,-323,-211,-213,22*0,111,113,-311,
1587  &2*-211,211,-211,310,-211,2*111,211,2*-211,-321,-211,2*211,-211,
1588  &111,-211,2*211,6*0,111,-211,111,-211,0,221,331,333,321,311,221,
1589  &331,333,321,311,20*0,3,39*0,-211,111,6*0,-411,-413,-10413,-10411,
1590  &-20413,-415,-411,-413,-10413,-10411,-20413,-415,-411,-413,16*0,
1591  &-4,-1,-4,-3,2*-2,5*0,111,-211,111,-211,-421,-423,-10423,-10421,
1592  &-20423,-425,-421,-423,-10423,-10421,-20423,-425,-421,-423,16*0,
1593  &-4,-1,-4,-3,2*-2,5*0,111,-211,111,-211,-431,-433,-10433,-10431,
1594  &-20433,-435,-431,-433,-10433,-10431,-20433,-435,-431,-433,19*0/
1595  data (kfdp(i,3),i=1277,3098)/-4,-1,-4,-3,2*-2,8*0,441,443,441,
1596  &443,441,443,-4,-1,-4,-3,-4,-3,-4,-1,531,533,531,533,3,2,3,2,511,
1597  &513,511,513,1,2,13*0,2*21,38*0,2*21,-211,111,460*0,2212,114*0,
1598  &2112,207*0,2*3122,3212,58*0,3322,2*0,3122,3212,3214,2112,2114,
1599  &2212,2112,3122,3212,3214,2112,2114,2212,2112,52*0,3*3,1,6*0,4*3,
1600  &4*0,4*3,6*0,4*3,0,28*3,2*0,3*4122,8*0,4,1,4,3,2*2,4*4,1,4,3,2*2,
1601  &4*4,1,4,3,2*2,4*0,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*0,4*4,1,4,3,2*2,
1602  &0,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,
1603  &4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,
1604  &4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,
1605  &3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,
1606  &164*0,-1,-3,-5,-2,-4,-6,-1,-3,-5,-2,-4,-6,-1,-3,-5,-2,-4,-6,-1,
1607  &-3,-5,-2,-4,-6,-2,2,-4,4,-6,6,-2,2,-4,4,-6,6,12*0,-11,-13,-15,
1608  &-12,-14,-16,-1,-3,-5,-2,-4,5*0,-12,12,-14,14,-16,16,-2,2,-4,4,
1609  &2*0,-12,12,-14,14,-16,16,-2,2,-4,4,52*0,-1,-3,-5,-2,-4,3*0,12,14,
1610  &16,2,4,0,12,14,16,2,4,0,12,14,16,2,4,0,12,14,16,2,4,28*0,2,4,7*0,
1611  &-11,-13,-15,-12,-14,-16,-1,-3,-5,-2,-4,5*0,-11,-13,-15,-12,-14,
1612  &-16,-1,-3,-5,-2,-4,5*0,-12,12,-14,14,-16,16,-2,2,-4,4,2*0,-12,12,
1613  &-14,14,-16,16,-2,2,-4,4,52*0,-1,-3,-5,-2,-4,7*0,-11,-13,-15,-12,
1614  &-14,-16,-1,-3,-5,-2,-4,5*0,-11,-13,-15,-12,-14,-16,-1,-3,-5,-2/
1615  data (kfdp(i,3),i=3099,4000)/-4,5*0,-11,-13,-15,-12,-14,-16,-1,
1616  &-3,-5,-2,-4,5*0,-12,12,-14,14,-16,16,-2,2,-4,4,2*0,-12,12,-14,14,
1617  &-16,16,-2,2,-4,4,52*0,-1,-3,-5,-2,-4,3*0,-11,-13,-15,-12,-14,-16,
1618  &-1,-3,-5,-2,-4,4*0,12,14,16,2,4,0,12,14,16,2,4,0,12,14,16,2,4,0,
1619  &12,14,16,2,4,28*0,2,4,730*0/
1620  data (kfdp(i,4),i= 1,4000)/92*0,4*111,6*0,111,2*0,-211,0,-211,
1621  &3*0,111,2*-211,0,111,0,2*111,113,221,2*111,-213,-211,211,113,
1622  &6*111,310,2*130,470*0,13*81,58*0,111,-211,4*0,2*111,13*0,-11,
1623  &10*0,111,-211,4*0,111,36*0,-211,111,177*0,111,211,111,211,7*0,
1624  &111,211,111,211,35*0,2*-211,2*111,211,111,-211,2*211,2*-211,13*0,
1625  &-211,111,-211,111,4*0,-211,111,-211,111,34*0,111,-211,3*111,
1626  &3*-211,2*111,3*-211,14*0,-321,-311,3*0,-321,-311,20*0,-3,77*0,
1627  &6*1,39*0,6*2,42*0,6*3,14*0,8*4,4*0,4*-5,4*0,2*-5,909*0,-211,111,
1628  &5*0,-211,111,52*0,2101,2103,2*2101,6*0,4*81,4*0,4*81,6*0,4*81,0,
1629  &28*81,13*0,6*2101,18*81,4*0,18*81,4*0,9*81,0,162*81,1416*0/
1630  data (kfdp(i,5),i= 1,4000)/94*0,2*111,17*0,111,7*0,2*111,0,
1631  &3*111,0,111,856*0,-211,2*111,-211,111,-211,111,65*0,111,-211,
1632  &3*111,-211,111,2936*0/
1633 
1634 !...PYDAT4, with particle names (character strings).
1635  data (chaf(i,1),i= 1, 184)/'d','u','s','c','b','t','b''','t''',
1636  &2*' ','e-','nu_e','mu-','nu_mu','tau-','nu_tau','tau''-',
1637  &'nu''_tau',2*' ','g','gamma','Z0','W+','h0',2*' ','reggeon',
1638  &'pomeron',2*' ','Z''0','Z"0','W''+','H0','A0','H+','eta_tech0',
1639  &'LQ_ue','R0',10*' ','pi_tech0','pi_tech+','pi''_tech0',
1640  &'rho_tech0','rho_tech+','omega_tech',24*' ','specflav',
1641  &'rndmflav','phasespa','c-hadron','b-hadron',5*' ','cluster',
1642  &'string','indep.','CMshower','SPHEaxis','THRUaxis','CLUSjet',
1643  &'CELLjet','table',' ','pi_diffr+','rho_diff0','phi_diff',
1644  &'omega_di','n_diffr0','p_diffr+','dd_1','ud_0','ud_1','uu_1',
1645  &'sd_0','sd_1','su_0','su_1','ss_1','cd_0','cd_1','cu_0','cu_1',
1646  &'cc_1','cs_0','cs_1','bd_0','bd_1','bu_0','bu_1','bs_0','bs_1',
1647  &'bc_0','bc_1','bb_1','pi0','pi+','rho0','rho+','a_00','a_0+',
1648  &'a_10','a_1+','a_20','a_2+','b_10','b_1+','pi(1300)0',
1649  &'pi(1300)+','pi_2(1670)0','pi_2(1670)+','rho(1465)0',
1650  &'rho(1465)+','rho(1700)0','rho(1700)+','sigma-m','eta','omega',
1651  &'f_2','eta''','phi','f''_2','f_0','h_1','f''_0','h''_1',
1652  &'eta(1295)','f_1','f''_1','f_0(1300)','omega(1420)','phi(1680)',
1653  &'f_1(1510)','omega(1600)','K_L0','K_S0','K0','K+','K*0','K*+',
1654  &'K_10','K_1+','K*_20','K*_2+','K*_00','K*_0+','K*_10','K*_1+'/
1655  data (chaf(i,1),i= 185, 305)/'K*(1410)0','K*(1410)+','K*(1680)0',
1656  &'K*(1680)+','K_3(1780)0','K_3(1780)+','K_2(1770)0','K_2(1770)+',
1657  &'K_2(1820)0','K_2(1820)+','D+','D*+','D*_2+','D0','D*0','D*_20',
1658  &'D_s+','D*_s+','D*_2s+','J/psi_di','eta_c','J/psi','chi_2c',
1659  &'D*_0+','D_1+','D*_00','D_10','D*_0s+','D_1s+','chi_0c','h_1c',
1660  &'D*_1+','D*_10','D*_1s+','chi_1c','psi''','B0','B*0','B*_20',
1661  &'B+','B*+','B*_2+','B_s0','B*_s0','B*_2s0','B_c+','B*_c+',
1662  &'B*_2c+','eta_b','Upsilon','chi_2b','B*_00','B_10','B*_0+',
1663  &'B_1+','B*_0s0','B_1s0','B*_0c+','B_1c+','chi_0b','h_1b','B*_10',
1664  &'B*_1+','B*_1s0','B*_1c+','chi_1b','Upsilon''','n0','p+',
1665  &'N(1440)0','N(1440)+','N(1520)0','N(1520)+','N(1535)0',
1666  &'N(1535)+','N(1650)0','N(1650)+','N(1675)0','N(1675)+',
1667  &'N(1680)0','N(1680)+','N(1700)0','N(1700)+','N(1710)0',
1668  &'N(1710)+','N(1720)0','N(1720)+','N(1990)0','N(1990)+','Delta-',
1669  &'Delta0','Delta+','Delta++','D(1600)-','D(1600)0','D(1600)+',
1670  &'D(1600)++','D(1620)-','D(1620)0','D(1620)+','D(1620)++',
1671  &'D(1700)-','D(1700)0','D(1700)+','D(1700)++','D(1900)-',
1672  &'D(1900)0','D(1900)+','D(1900)++','D(1905)-','D(1905)0',
1673  &'D(1905)+','D(1905)++','D(1910)-','D(1910)0','D(1910)+',
1674  &'D(1910)++','D(1920)-','D(1920)0','D(1920)+','D(1920)++'/
1675  data (chaf(i,1),i= 306, 413)/'D(1930)-','D(1930)0','D(1930)+',
1676  &'D(1930)++','D(1950)-','D(1950)0','D(1950)+','D(1950)++',
1677  &'Lambda0','L(1405)0','L(1520)0','L(1600)0','L(1670)0','L(1690)0',
1678  &'L(1800)0','L(1810)0','L(1820)0','L(1830)0','L(1890)0',
1679  &'L(2100)0','L(2110)0','Sigma-','Sigma0','Sigma+','Sigma*-',
1680  &'Sigma*0','Sigma*+','S(1660)-','S(1660)0','S(1660)+','S(1670)-',
1681  &'S(1670)0','S(1670)+','S(1750)-','S(1750)0','S(1750)+',
1682  &'S(1775)-','S(1775)0','S(1775)+','S(1915)-','S(1915)0',
1683  &'S(1915)+','S(1940)-','S(1940)0','S(1940)+','S(2030)-',
1684  &'S(2030)0','S(2030)+','Xi-','Xi0','Xi*-','Xi*0','X(1690)-',
1685  &'X(1690)0','X(1820)-','X(1820)0','X(1950)-','X(1950)0',
1686  &'X(2030)-','X(2030)0','Omega-','Sigma_c0','Sigma*_c0',
1687  &'Lambda_c+','Xi_c0','Sigma_c+','Sigma*_c+','Sigma_c++',
1688  &'Sigma*_c++','Xi_c+','Xi''_c0','Xi*_c0','Xi''_c+','Xi*_c+',
1689  &'Omega_c0','Omega*_c0','Xi_cc+','Xi*_cc+','Xi_cc++','Xi*_cc++',
1690  &'Omega_cc+','Omega*_cc+','Omega*_ccc++','Sigma_b-','Sigma*_b-',
1691  &'Lambda_b0','Xi_b-','Xi_bc0','Sigma_b0','Sigma*_b0','Sigma_b+',
1692  &'Sigma*_b+','Xi_b0','Xi_bc+','Xi''_b-','Xi*_b-','Xi''_b0',
1693  &'Xi*_b0','Omega_b-','Omega*_b-','Omega_bc0','Xi''_bc0','Xi*_bc0',
1694  &'Xi''_bc+','Xi*_bc+','Omega''_bc0','Omega*_bc0','Omega_bcc+'/
1695  data (chaf(i,1),i= 414, 500)/'Omega*_bcc+','Xi_bb-','Xi*_bb-',
1696  &'Xi_bb0','Xi*_bb0','Omega_bb-','Omega*_bb-','Omega_bbc0',
1697  &'Omega*_bbc0','Omega*_bbb-','~d_L','~u_L','~s_L','~c_L','~b_1',
1698  &'~t_1','~e_L-','~nu_eL','~mu_L-','~nu_muL','~tau_1-','~nu_tauL',
1699  &'~g','~chi_10','~chi_20','~chi_1+','~chi_30','~chi_40','~chi_2+',
1700  &'~gravitino','~d_R','~u_R','~s_R','~c_R','~b_2','~t_2','~e_R-',
1701  &'~nu_eR','~mu_R-','~nu_muR','~tau_2-','~nu_tauR','d*','u*','e*-',
1702  &'nu*_e0',41*' '/
1703  data (chaf(i,2),i= 1, 227)/'dbar','ubar','sbar','cbar','bbar',
1704  &'tbar','b''bar','t''bar',2*' ','e+','nu_ebar','mu+','nu_mubar',
1705  &'tau+','nu_taubar','tau''+','nu''_taubar',5*' ','W-',9*' ',
1706  &'W''-',2*' ','H-',' ','LQ_uebar','Rbar0',11*' ','pi_tech-',2*' ',
1707  &'rho_tech-',26*' ','rndmflavbar',' ','c-hadronbar','b-hadronbar',
1708  &15*' ','pi_diffr-',3*' ','n_diffrbar0','p_diffrbar-','dd_1bar',
1709  &'ud_0bar','ud_1bar','uu_1bar','sd_0bar','sd_1bar','su_0bar',
1710  &'su_1bar','ss_1bar','cd_0bar','cd_1bar','cu_0bar','cu_1bar',
1711  &'cc_1bar','cs_0bar','cs_1bar','bd_0bar','bd_1bar','bu_0bar',
1712  &'bu_1bar','bs_0bar','bs_1bar','bc_0bar','bc_1bar','bb_1bar',' ',
1713  &'pi-',' ','rho-',' ','a_0-',' ','a_1-',' ','a_2-',' ','b_1-',' ',
1714  &'pi(1300)-',' ','pi_2(1670)-',' ','rho(1465)-',' ','rho(1700)-',
1715  &21*' ','Kbar0','K-','K*bar0','K*-','K_1bar0','K_1-','K*_2bar0',
1716  &'K*_2-','K*_0bar0','K*_0-','K*_1bar0','K*_1-','K*(1410)bar0',
1717  &'K*(1410)-','K*(1680)bar0','K*(1680)-','K_3(1780)bar0',
1718  &'K_3(1780)-','K_2(1770)bar0','K_2(1770)-','K_2(1820)bar0',
1719  &'K_2(1820)-','D-','D*-','D*_2-','Dbar0','D*bar0','D*_2bar0',
1720  &'D_s-','D*_s-','D*_2s-',4*' ','D*_0-','D_1-','D*_0bar0',
1721  &'D_1bar0','D*_0s-','D_1s-',2*' ','D*_1-','D*_1bar0','D*_1s-',
1722  &2*' ','Bbar0','B*bar0','B*_2bar0','B-','B*-','B*_2-','B_sbar0'/
1723  data (chaf(i,2),i= 228, 319)/'B*_sbar0','B*_2sbar0','B_c-',
1724  &'B*_c-','B*_2c-',3*' ','B*_0bar0','B_1bar0','B*_0-','B_1-',
1725  &'B*_0sbar0','B_1sbar0','B*_0c-','B_1c-',2*' ','B*_1bar0','B*_1-',
1726  &'B*_1sbar0','B*_1c-',2*' ','nbar0','pbar-','N(1440)bar0',
1727  &'N(1440)bar-','N(1520)bar0','N(1520)bar-','N(1535)bar0',
1728  &'N(1535)bar-','N(1650)bar0','N(1650)bar-','N(1675)bar0',
1729  &'N(1675)bar-','N(1680)bar0','N(1680)bar-','N(1700)bar0',
1730  &'N(1700)bar-','N(1710)bar0','N(1710)bar-','N(1720)bar0',
1731  &'N(1720)bar-','N(1990)bar0','N(1990)bar-','Deltabar+',
1732  &'Deltabar0','Deltabar-','Deltabar--','D(1600)bar+','D(1600)bar0',
1733  &'D(1600)bar-','D(1600)bar--','D(1620)bar+','D(1620)bar0',
1734  &'D(1620)bar-','D(1620)bar--','D(1700)bar+','D(1700)bar0',
1735  &'D(1700)bar-','D(1700)bar--','D(1900)bar+','D(1900)bar0',
1736  &'D(1900)bar-','D(1900)bar--','D(1905)bar+','D(1905)bar0',
1737  &'D(1905)bar-','D(1905)bar--','D(1910)bar+','D(1910)bar0',
1738  &'D(1910)bar-','D(1910)bar--','D(1920)bar+','D(1920)bar0',
1739  &'D(1920)bar-','D(1920)bar--','D(1930)bar+','D(1930)bar0',
1740  &'D(1930)bar-','D(1930)bar--','D(1950)bar+','D(1950)bar0',
1741  &'D(1950)bar-','D(1950)bar--','Lambdabar0','L(1405)bar0',
1742  &'L(1520)bar0','L(1600)bar0','L(1670)bar0','L(1690)bar0'/
1743  data (chaf(i,2),i= 320, 401)/'L(1800)bar0','L(1810)bar0',
1744  &'L(1820)bar0','L(1830)bar0','L(1890)bar0','L(2100)bar0',
1745  &'L(2110)bar0','Sigmabar+','Sigmabar0','Sigmabar-','Sigma*bar+',
1746  &'Sigma*bar0','Sigma*bar-','S(1660)bar+','S(1660)bar0',
1747  &'S(1660)bar-','S(1670)bar+','S(1670)bar0','S(1670)bar-',
1748  &'S(1750)bar+','S(1750)bar0','S(1750)bar-','S(1775)bar+',
1749  &'S(1775)bar0','S(1775)bar-','S(1915)bar+','S(1915)bar0',
1750  &'S(1915)bar-','S(1940)bar+','S(1940)bar0','S(1940)bar-',
1751  &'S(2030)bar+','S(2030)bar0','S(2030)bar-','Xibar+','Xibar0',
1752  &'Xi*bar+','Xi*bar0','X(1690)bar+','X(1690)bar0','X(1820)bar+',
1753  &'X(1820)bar0','X(1950)bar+','X(1950)bar0','X(2030)bar+',
1754  &'X(2030)bar0','Omegabar+','Sigma_cbar0','Sigma*_cbar0',
1755  &'Lambda_cbar-','Xi_cbar0','Sigma_cbar-','Sigma*_cbar-',
1756  &'Sigma_cbar--','Sigma*_cbar--','Xi_cbar-','Xi''_cbar0',
1757  &'Xi*_cbar0','Xi''_cbar-','Xi*_cbar-','Omega_cbar0',
1758  &'Omega*_cbar0','Xi_ccbar-','Xi*_ccbar-','Xi_ccbar--',
1759  &'Xi*_ccbar--','Omega_ccbar-','Omega*_ccbar-','Omega*_cccbar-',
1760  &'Sigma_bbar+','Sigma*_bbar+','Lambda_bbar0','Xi_bbar+',
1761  &'Xi_bcbar0','Sigma_bbar0','Sigma*_bbar0','Sigma_bbar-',
1762  &'Sigma*_bbar-','Xi_bbar0','Xi_bcbar-','Xi''_bbar+','Xi*_bbar+'/
1763  data (chaf(i,2),i= 402, 500)/'Xi''_bbar0','Xi*_bbar0',
1764  &'Omega_bbar+','Omega*_bbar+','Omega_bcbar0','Xi''_bcbar0',
1765  &'Xi*_bcbar0','Xi''_bcbar-','Xi*_bcbar-','Omega''_bcba',
1766  &'Omega*_bcbar0','Omega_bccbar-','Omega*_bccbar-','Xi_bbbar+',
1767  &'Xi*_bbbar+','Xi_bbbar0','Xi*_bbbar0','Omega_bbbar+',
1768  &'Omega*_bbbar+','Omega_bbcbar0','Omega*_bbcbar0',
1769  &'Omega*_bbbbar+','~d_Lbar','~u_Lbar','~s_Lbar','~c_Lbar',
1770  &'~b_1bar','~t_1bar','~e_L+','~nu_eLbar','~mu_L+','~nu_muLbar',
1771  &'~tau_1+','~nu_tauLbar',3*' ','~chi_1-',2*' ','~chi_2-',' ',
1772  &'~d_Rbar','~u_Rbar','~s_Rbar','~c_Rbar','~b_2bar','~t_2bar',
1773  &'~e_R+','~nu_eRbar','~mu_R+','~nu_muRbar','~tau_2+',
1774  &'~nu_tauRbar','d*bar','u*bar','e*bar+','nu*_ebar0',41*' '/
1775 
1776 !...Treatment of resonances.
1777  data (mwid(i) ,i= 1, 500)/5*0,3*1,8*0,1,5*0,3*1,6*0,1,0,7*1,
1778  &10*0,2*1,0,3*1,367*0,19*2,0,7*2,0,2,0,2,0,4*1,41*0/
1779  end
z
nodes z
Definition: cstdatmos0.fNew2.h:93
i
nodes i
Definition: cstdatmos0.fNew2.h:93
parp
subroutine mydummy common pjpars parp(200)
j
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code kiron data j
Definition: cblkHeavy.h:36
d0
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
Definition: cblkEvhnp.h:5
o
struct ob o[NpMax]
Definition: Zprivate.h:34
msti
subroutine mydummy common pjpars msti(200)
d
dE dx *! Nuc Int sampling table d
Definition: cblkMuInt.h:130
a
nodes a
Definition: cstdatmos0.fNew2.h:93
mstp
subroutine mydummy common pjpars mstp(200)
h
dE dx *! Nuc Int sampling table h
Definition: cblkMuInt.h:130