pythiaeRHIC/pythia/pydata.f

0001 C*********************************************************************
0002 C*********************************************************************
0003 C*                                                                  **
0004 C*                                                 September 2007   **
0005 C*                                                                  **
0006 C*                       The Lund Monte Carlo                       **
0007 C*                                                                  **
0008 C*                        PYTHIA version 6.4                        **
0009 C*                                                                  **
0010 C*                        Torbjorn Sjostrand                        **
0011 C*               CERN/PH, CH-1211 Geneva, Switzerland               **
0012 C*                    phone +41 - 22 - 767 82 27                    **
0013 C*                               and                                **
0014 C*                 Department of Theoretical Physics                **
0015 C*                         Lund University                          **
0016 C*               Solvegatan 14A, S-223 62 Lund, Sweden              **
0017 C*                    E-mail torbjorn@thep.lu.se                    **
0018 C*                                                                  **
0019 C*                  SUSY and Technicolor parts by                   **
0020 C*                         Stephen Mrenna                           **
0021 C*                       Computing Division                         **
0022 C*            Generators and Detector Simulation Group              **
0023 C*              Fermi National Accelerator Laboratory               **
0024 C*                 MS 234, Batavia, IL  60510, USA                  **
0025 C*                   phone + 1 - 630 - 840 - 2556                   **
0026 C*                      E-mail mrenna@fnal.gov                      **
0027 C*                                                                  **
0028 C*         New multiple interactions and more SUSY parts by         **
0029 C*                          Peter Skands                            **
0030 C*                  Theoretical Physics Department                  **
0031 C*              Fermi National Accelerator Laboratory               **
0032 C*                 MS 106, Batavia, IL  60510, USA                  **
0033 C*                               and                                **
0034 C*               CERN/PH, CH-1211 Geneva, Switzerland               **
0035 C*                    phone +41 - 22 - 767 24 59                    **
0036 C*                      E-mail skands@fnal.gov                      **
0037 C*                                                                  **
0038 C*         Several parts are written by Hans-Uno Bengtsson          **
0039 C*          PYSHOW is written together with Mats Bengtsson          **
0040 C*               PYMAEL is written by Emanuel Norrbin               **
0041 C*     advanced popcorn baryon production written by Patrik Eden    **
0042 C*    code for virtual photons mainly written by Christer Friberg   **
0043 C*    code for low-mass strings mainly written by Emanuel Norrbin   **
0044 C*        Bose-Einstein code mainly written by Leif Lonnblad        **
0045 C*      CTEQ  parton distributions are by the CTEQ collaboration    **
0046 C*      GRV 94 parton distributions are by Glueck, Reya and Vogt    **
0047 C*   SaS photon parton distributions together with Gerhard Schuler  **
0048 C*     g + g and q + qbar -> t + tbar + H code by Zoltan Kunszt     **
0049 C*         MSSM Higgs mass calculation code by M. Carena,           **
0050 C*           J.R. Espinosa, M. Quiros and C.E.M. Wagner             **
0051 C*         PYGAUS adapted from CERN library (K.S. Kolbig)           **
0052 C*        NRQCD/colour octet production of onium by S. Wolf         **
0053 C*                                                                  **
0054 C*   The latest program version and documentation is found on WWW   **
0055 C*            http://www.thep.lu.se/~torbjorn/Pythia.html           **
0056 C*                                                                  **
0057 C*        Copyright Torbjorn Sjostrand, Lund (and CERN) 2007        **
0058 C*                                                                  **
0059 C*********************************************************************
0060 C*********************************************************************
0061 C                                                                    *
0062 C  List of subprograms in order of appearance, with main purpose     *
0063 C  (S = subroutine, F = function, B = block data)                    *
0064 C                                                                    *
0065 C  B   PYDATA   to contain all default values                        *
0066 C  S   PYCKBD   to check that BLOCK DATA has been correctly loaded   *
0067 C  S   PYTEST   to test the proper functioning of the package        *
0068 C  S   PYHEPC   to convert between /PYJETS/ and /HEPEVT/ records     *
0069 C                                                                    *
0070 C  S   PYINIT   to administer the initialization procedure           *
0071 C  S   PYEVNT   to administer the generation of an event             *
0072 C  S   PYEVNW   ditto, for new multiple interactions scenario        *
0073 C  S   PYSTAT   to print cross-section and other information         *
0074 C  S   PYUPEV   to administer the generation of an LHA hard process  *
0075 C  S   PYUPIN   to provide initialization needed for LHA input       *
0076 C  S   PYLHEF   to produce a Les Houches Event File from run         *
0077 C  S   PYINRE   to initialize treatment of resonances                *
0078 C  S   PYINBM   to read in beam, target and frame choices            *
0079 C  S   PYINKI   to initialize kinematics of incoming particles       *
0080 C  S   PYINPR   to set up the selection of included processes        *
0081 C  S   PYXTOT   to give total, elastic and diffractive cross-sect.   *
0082 C  S   PYMAXI   to find differential cross-section maxima            *
0083 C  S   PYPILE   to select multiplicity of pileup events              *
0084 C  S   PYSAVE   to save alternatives for gamma-p and gamma-gamma     *
0085 C  S   PYGAGA   to handle lepton -> lepton + gamma branchings        *
0086 C  S   PYRAND   to select subprocess and kinematics for event        *
0087 C  S   PYSCAT   to set up kinematics and colour flow of event        *
0088 C  S   PYEVOL   handler for pT-ordered ISR and multiple interactions *
0089 C  S   PYSSPA   to simulate initial state spacelike showers          *
0090 C  S   PYPTIS   to do pT-ordered initial state spacelike showers     *
0091 C  S   PYMEMX   auxiliary to PYSSPA/PYPTIS for ME correction maximum *
0092 C  S   PYMEWT   auxiliary to PYSSPA/.. for matrix element correction *
0093 C  S   PYPTMI   to do pT-ordered multiple interactions               *
0094 C  F   PYFCMP   to give companion quark x*f distribution             *
0095 C  F   PYPCMP   to calculate momentum integral for companion quarks  *
0096 C  S   PYUPRE   to rearranges contents of the HEPEUP commonblock     *
0097 C  S   PYADSH   to administrate sequential final-state showers       *
0098 C  S   PYVETO   to allow the generation of an event to be aborted    *
0099 C  S   PYRESD   to perform resonance decays                          *
0100 C  S   PYMULT   to generate multiple interactions - old scheme       *
0101 C  S   PYREMN   to add on target remnants - old scheme               *
0102 C  S   PYMIGN   to generate multiple interactions - new scheme       *
0103 C  S   PYMIHK   to connect colours in mult. int. - new scheme        *
0104 C  S   PYCTTR   to translate PYTHIA colour information to LHA1 tags  *
0105 C  S   PYMIHG   to collapse two pairs of LHA1 colour tags.           *
0106 C  S   PYMIRM   to add on target remnants in mult. int.- new scheme  *
0107 C  S   PYFSCR   to perform final state colour reconnections - -"-    *
0108 C  S   PYDIFF   to set up kinematics for diffractive events          *
0109 C  S   PYDISG   to set up kinematics, remnant and showers for DIS    *
0110 C  S   PYDOCU   to compute cross-sections and handle documentation   *
0111 C  S   PYFRAM   to perform boosts between different frames           *
0112 C  S   PYWIDT   to calculate full and partial widths of resonances   *
0113 C  S   PYOFSH   to calculate partial width into off-shell channels   *
0114 C  S   PYRECO   to handle colour reconnection in W+W- events         *
0115 C  S   PYKLIM   to calculate borders of allowed kinematical region   *
0116 C  S   PYKMAP   to construct value of kinematical variable           *
0117 C  S   PYSIGH   to calculate differential cross-sections             *
0118 C  S   PYSGQC   auxiliary to PYSIGH for QCD processes                *
0119 C  S   PYSGHF   auxiliary to PYSIGH for heavy flavour processes      *
0120 C  S   PYSGWZ   auxiliary to PYSIGH for W and Z processes            *
0121 C  S   PYSGHG   auxiliary to PYSIGH for Higgs processes              *
0122 C  S   PYSGSU   auxiliary to PYSIGH for supersymmetry processes      *
0123 C  S   PYSGTC   auxiliary to PYSIGH for technicolor processes        *
0124 C  S   PYSGEX   auxiliary to PYSIGH for various exotic processes     *
0125 C  S   PYPDFU   to evaluate parton distributions                     *
0126 C  S   PYPDFL   to evaluate parton distributions at low x and Q^2    *
0127 C  S   PYPDEL   to evaluate electron parton distributions            *
0128 C  S   PYPDGA   to evaluate photon parton distributions (generic)    *
0129 C  S   PYGGAM   to evaluate photon parton distributions (SaS sets)   *
0130 C  S   PYGVMD   to evaluate VMD part of photon parton distributions  *
0131 C  S   PYGANO   to evaluate anomalous part of photon PDFs            *
0132 C  S   PYGBEH   to evaluate Bethe-Heitler part of photon PDFs        *
0133 C  S   PYGDIR   to evaluate direct contribution to photon PDFs       *
0134 C  S   PYPDPI   to evaluate pion parton distributions                *
0135 C  S   PYPDPR   to evaluate proton parton distributions              *
0136 C  F   PYCTEQ   to evaluate the CTEQ 3 proton parton distributions   *
0137 C  S   PYGRVL   to evaluate the GRV 94L proton parton distributions  *
0138 C  S   PYGRVM   to evaluate the GRV 94M proton parton distributions  *
0139 C  S   PYGRVD   to evaluate the GRV 94D proton parton distributions  *
0140 C  F   PYGRVV   auxiliary to the PYGRV* routines                     *
0141 C  F   PYGRVW   auxiliary to the PYGRV* routines                     *
0142 C  F   PYGRVS   auxiliary to the PYGRV* routines                     *
0143 C  F   PYCT5L   to evaluate the CTEQ 5L proton parton distributions  *
0144 C  F   PYCT5M   to evaluate the CTEQ 5M1 proton parton distributions *
0145 C  S   PYPDPO   to evaluate old proton parton distributions          *
0146 C  F   PYHFTH   to evaluate threshold factor for heavy flavour       *
0147 C  S   PYSPLI   to find flavours left in hadron when one removed     *
0148 C  F   PYGAMM   to evaluate ordinary Gamma function Gamma(x)         *
0149 C  S   PYWAUX   to evaluate auxiliary functions W1(s) and W2(s)      *
0150 C  S   PYI3AU   to evaluate auxiliary function I3(s,t,u,v)           *
0151 C  F   PYSPEN   to evaluate Spence (dilogarithm) function Sp(x)      *
0152 C  S   PYQQBH   to evaluate matrix element for g + g -> Q + Qbar + H *
0153 C  S   PYSTBH   to evaluate matrix element for t + b + H processes   *
0154 C  S   PYTBHB   auxiliary to PYSTBH                                  *
0155 C  S   PYTBHG   auxiliary to PYSTBH                                  *
0156 C  S   PYTBHQ   auxiliary to PYSTBH                                  *
0157 C  F   PYTBHS   auxiliary to PYSTBH                                  *
0158 C                                                                    *
0159 C  S   PYMSIN   to initialize the supersymmetry simulation           *
0160 C  S   PYSLHA   to interface to SUSY spectrum and decay calculators  *
0161 C  S   PYAPPS   to determine MSSM parameters from SUGRA input        *
0162 C  S   PYSUGI   to determine MSSM parameters using ISASUSY           *
0163 C  S   PYFEYN   to determine MSSM Higgs parameters using FEYNHIGGS   *
0164 C  F   PYRNMQ   to determine running squark masses                   *
0165 C  S   PYTHRG   to calculate sfermion third-gen. mass eigenstates    *
0166 C  S   PYINOM   to calculate neutralino/chargino mass eigenstates    *
0167 C  F   PYRNM3   to determine running M3, gluino mass                 *
0168 C  S   PYEIG4   to calculate eigenvalues and -vectors in 4*4 matrix  *
0169 C  S   PYHGGM   to determine Higgs mass spectrum                     *
0170 C  S   PYSUBH   to determine Higgs masses in the MSSM                *
0171 C  S   PYPOLE   to determine Higgs masses in the MSSM                *
0172 C  S   PYRGHM   auxiliary to PYPOLE                                  *
0173 C  S   PYGFXX   auxiliary to PYRGHM                                  *
0174 C  F   PYFINT   auxiliary to PYPOLE                                  *
0175 C  F   PYFISB   auxiliary to PYFINT                                  *
0176 C  S   PYSFDC   to calculate sfermion decay partial widths           *
0177 C  S   PYGLUI   to calculate gluino decay partial widths             *
0178 C  S   PYTBBN   to calculate 3-body decay of gluino to neutralino    *
0179 C  S   PYTBBC   to calculate 3-body decay of gluino to chargino      *
0180 C  S   PYNJDC   to calculate neutralino decay partial widths         *
0181 C  S   PYCJDC   to calculate chargino decay partial widths           *
0182 C  F   PYXXZ6   auxiliary for ino 3-body decays                      *
0183 C  F   PYXXGA   auxiliary for ino -> ino + gamma decay               *
0184 C  F   PYX2XG   auxiliary for ino -> ino + gauge boson decay         *
0185 C  F   PYX2XH   auxiliary for ino -> ino + Higgs decay               *
0186 C  S   PYHEXT   to calculate non-SM Higgs decay partial widths       *
0187 C  F   PYH2XX   auxiliary for H -> ino + ino decay                   *
0188 C  F   PYGAUS   to perform Gaussian integration                      *
0189 C  F   PYGAU2   copy of PYGAUS to allow two-dimensional integration  *
0190 C  F   PYSIMP   to perform Simpson integration                       *
0191 C  F   PYLAMF   to evaluate the lambda kinematics function           *
0192 C  S   PYTBDY   to perform 3-body decay of gauginos                  *
0193 C  S   PYTECM   to calculate techni_rho/omega masses                 *
0194 C  S   PYEICG   to calculate eigenvalues of a 4*4 complex matrix     *
0195 C  S   PYCMQR   auxiliary to PYEICG                                  *
0196 C  S   PYCMQ2   auxiliary to PYEICG                                  *
0197 C  S   PYCDIV   auxiliary to PYCMQR                                  *
0198 C  S   PYCSRT   auxiliary to PYCMQR                                  *
0199 C  S   PYTHAG   auxiliary to PYCMQR                                  *
0200 C  S   PYCBAL   auxiliary to PYEICG                                  *
0201 C  S   PYCBA2   auxiliary to PYEICG                                  *
0202 C  S   PYCRTH   auxiliary to PYEICG                                  *
0203 C  S   PYLDCM   auxiliary to PYSIGH, for technicolor in QCD 2 -> 2   *
0204 C  S   PYBKSB   auxiliary to PYSIGH, for technicolor in QCD 2 -> 2   *
0205 C  S   PYWIDX   to calculate decay widths from within PYWIDT         *
0206 C  S   PYRVSF   to calculate R-violating sfermion decay widths       *
0207 C  S   PYRVNE   to calculate R-violating neutralino decay widths     *
0208 C  S   PYRVCH   to calculate R-violating chargino decay widths       *
0209 C  S   PYRVGL   to calculate R-violating gluino decay widths         *
0210 C  F   PYRVSB   auxiliary to PYRVSF                                  *
0211 C  S   PYRVGW   to calculate R-Violating 3-body widths               *
0212 C  F   PYRVI1   auxiliary to PYRVGW, to do PS integration for res.   *
0213 C  F   PYRVI2   auxiliary to PYRVGW, to do PS integration for LR-int.*
0214 C  F   PYRVI3   auxiliary to PYRVGW, to do PS X integral for int.    *
0215 C  F   PYRVG1   auxiliary to PYRVI1, general matrix element, res.    *
0216 C  F   PYRVG2   auxiliary to PYRVI2, general matrix element, LR-int. *
0217 C  F   PYRVG3   auxiliary to PYRVI3, to do PS Y integral for int.    *
0218 C  F   PYRVG4   auxiliary to PYRVG3, general matrix element, int.    *
0219 C  F   PYRVR    auxiliary to PYRVG1, Breit-Wigner                    *
0220 C  F   PYRVS    auxiliary to PYRVG2 & PYRVG4                         *
0221 C                                                                    *
0222 C  S   PY1ENT   to fill one entry (= parton or particle)             *
0223 C  S   PY2ENT   to fill two entries                                  *
0224 C  S   PY3ENT   to fill three entries                                *
0225 C  S   PY4ENT   to fill four entries                                 *
0226 C  S   PY2FRM   to interface to generic two-fermion generator        *
0227 C  S   PY4FRM   to interface to generic four-fermion generator       *
0228 C  S   PY6FRM   to interface to generic six-fermion generator        *
0229 C  S   PY4JET   to generate a shower from a given 4-parton config    *
0230 C  S   PY4JTW   to evaluate the weight od a shower history for above *
0231 C  S   PY4JTS   to set up the parton configuration for above         *
0232 C  S   PYJOIN   to connect entries with colour flow information      *
0233 C  S   PYGIVE   to fill (or query) commonblock variables             *
0234 C  S   PYONOF   to allow easy control of particle decay modes        *
0235 C  S   PYTUNE   to select a predefined 'tune' for min-bias and UE    *
0236 C  S   PYEXEC   to administrate fragmentation and decay chain        *
0237 C  S   PYPREP   to rearrange showered partons along strings          *
0238 C  S   PYSTRF   to do string fragmentation of jet system             *
0239 C  S   PYJURF   to find boost to string junction rest frame          *
0240 C  S   PYINDF   to do independent fragmentation of one or many jets  *
0241 C  S   PYDECY   to do the decay of a particle                        *
0242 C  S   PYDCYK   to select parton and hadron flavours in decays       *
0243 C  S   PYKFDI   to select parton and hadron flavours in fragm        *
0244 C  S   PYNMES   to select number of popcorn mesons                   *
0245 C  S   PYKFIN   to calculate falvour prod. ratios from input params. *
0246 C  S   PYPTDI   to select transverse momenta in fragm                *
0247 C  S   PYZDIS   to select longitudinal scaling variable in fragm     *
0248 C  S   PYSHOW   to do m-ordered timelike parton shower evolution     *
0249 C  S   PYPTFS   to do pT-ordered timelike parton shower evolution    *
0250 C  F   PYMAEL   auxiliary to PYSHOW & PYPTFS: gluon emission ME's    *
0251 C  S   PYBOEI   to include Bose-Einstein effects (crudely)           *
0252 C  S   PYBESQ   auxiliary to PYBOEI                                  *
0253 C  F   PYMASS   to give the mass of a particle or parton             *
0254 C  F   PYMRUN   to give the running MSbar mass of a quark            *
0255 C  S   PYNAME   to give the name of a particle or parton             *
0256 C  F   PYCHGE   to give three times the electric charge              *
0257 C  F   PYCOMP   to compress standard KF flavour code to internal KC  *
0258 C  S   PYERRM   to write error messages and abort faulty run         *
0259 C  F   PYALEM   to give the alpha_electromagnetic value              *
0260 C  F   PYALPS   to give the alpha_strong value                       *
0261 C  F   PYANGL   to give the angle from known x and y components      *
0262 C  F   PYR      to provide a random number generator                 *
0263 C  S   PYRGET   to save the state of the random number generator     *
0264 C  S   PYRSET   to set the state of the random number generator      *
0265 C  S   PYROBO   to rotate and/or boost an event                      *
0266 C  S   PYEDIT   to remove unwanted entries from record               *
0267 C  S   PYLIST   to list event record or particle data                *
0268 C  S   PYLOGO   to write a logo                                      *
0269 C  S   PYUPDA   to update particle data                              *
0270 C  F   PYK      to provide integer-valued event information          *
0271 C  F   PYP      to provide real-valued event information             *
0272 C  S   PYSPHE   to perform sphericity analysis                       *
0273 C  S   PYTHRU   to perform thrust analysis                           *
0274 C  S   PYCLUS   to perform three-dimensional cluster analysis        *
0275 C  S   PYCELL   to perform cluster analysis in (eta, phi, E_T)       *
0276 C  S   PYJMAS   to give high and low jet mass of event               *
0277 C  S   PYFOWO   to give Fox-Wolfram moments                          *
0278 C  S   PYTABU   to analyze events, with tabular output               *
0279 C                                                                    *
0280 C  S   PYEEVT   to administrate the generation of an e+e- event      *
0281 C  S   PYXTEE   to give the total cross-section at given CM energy   *
0282 C  S   PYRADK   to generate initial state photon radiation           *
0283 C  S   PYXKFL   to select flavour of primary qqbar pair              *
0284 C  S   PYXJET   to select (matrix element) jet multiplicity          *
0285 C  S   PYX3JT   to select kinematics of three-jet event              *
0286 C  S   PYX4JT   to select kinematics of four-jet event               *
0287 C  S   PYXDIF   to select angular orientation of event               *
0288 C  S   PYONIA   to perform generation of onium decay to gluons       *
0289 C                                                                    *
0290 C  S   PYBOOK   to book a histogram                                  *
0291 C  S   PYFILL   to fill an entry in a histogram                      *
0292 C  S   PYFACT   to multiply histogram contents by a factor           *
0293 C  S   PYOPER   to perform operations between histograms             *
0294 C  S   PYHIST   to print and reset all histograms                    *
0295 C  S   PYPLOT   to print a single histogram                          *
0296 C  S   PYNULL   to reset contents of a single histogram              *
0297 C  S   PYDUMP   to dump histogram contents onto a file               *
0298 C                                                                    *
0299 C  S   PYSTOP   routine to handle Fortran STOP condition             *
0300 C                                                                    *
0301 C  S   PYKCUT   dummy routine for user kinematical cuts              *
0302 C  S   PYEVWT   dummy routine for weighting events                   *
0303 C  S   UPINIT   dummy routine to initialize user processes           *
0304 C  S   UPEVNT   dummy routine to generate a user process event       *
0305 C  S   UPVETO   dummy routine to abort event at parton level         *
0306 C  S   PDFSET   dummy routine to be removed when using PDFLIB        *
0307 C  S   STRUCTM  dummy routine to be removed when using PDFLIB        *
0308 C  S   STRUCTP  dummy routine to be removed when using PDFLIB        *
0309 C  S   SUGRA    dummy routine to be removed when linking with ISAJET *
0310 C  F   VISAJE   dummy functn. to be removed when linking with ISAJET *
0311 C  S   SSMSSM   dummy routine to be removed when linking with ISAJET *
0312 C  S   FHSETFLAGS  dummy routine          -"-              FEYNHIGGS *
0313 C  S   FHSETPARA   dummy routine          -"-              FEYNHIGGS *
0314 C  S   FHHIGGSCORR dummy routine          -"-              FEYNHIGGS *
0315 C  S   PYTAUD   dummy routine for interface to tau decay libraries   *
0316 C  S   PYTIME   dummy routine for giving date and time               *
0317 C                                                                    *
0318 C*********************************************************************
0319
0320 C...PYDATA
0321 C...Default values for switches and parameters,
0322 C...and particle, decay and process data.
0323
0324       BLOCK DATA PYDATA
0325
0326 C...Double precision and integer declarations.
0327       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
0328       IMPLICIT INTEGER(I-N)
0329       INTEGER PYK,PYCHGE,PYCOMP
0330 C...Commonblocks.
0331       COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
0332       COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
0333       COMMON/PYDAT3/MDCY(500,3),MDME(8000,2),BRAT(8000),KFDP(8000,5)
0334       COMMON/PYDAT4/CHAF(500,2)
0335       CHARACTER CHAF*16
0336       COMMON/PYDATR/MRPY(6),RRPY(100)
0337       COMMON/PYSUBS/MSEL,MSELPD,MSUB(500),KFIN(2,-40:40),CKIN(200)
0338       COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
0339       COMMON/PYINT1/MINT(400),VINT(400)
0340       COMMON/PYINT2/ISET(500),KFPR(500,2),COEF(500,20),ICOL(40,4,2)
0341       COMMON/PYINT3/XSFX(2,-40:40),ISIG(1000,3),SIGH(1000)
0342       COMMON/PYINT4/MWID(500),WIDS(500,5)
0343       COMMON/PYINT5/NGENPD,NGEN(0:500,3),XSEC(0:500,3)
0344       COMMON/PYINT6/PROC(0:500)
0345       CHARACTER PROC*28
0346       COMMON/PYINT7/SIGT(0:6,0:6,0:5)
0347       COMMON/PYMSSM/IMSS(0:99),RMSS(0:99)
0348       COMMON/PYSSMT/ZMIX(4,4),UMIX(2,2),VMIX(2,2),SMZ(4),SMW(2),
0349      &SFMIX(16,4),ZMIXI(4,4),UMIXI(2,2),VMIXI(2,2)
0350       COMMON/PYMSRV/RVLAM(3,3,3), RVLAMP(3,3,3), RVLAMB(3,3,3)
0351       COMMON/PYTCSM/ITCM(0:99),RTCM(0:99)
0352       COMMON/PYBINS/IHIST(4),INDX(1000),BIN(20000)
0353       COMMON/PYLH3P/MODSEL(200),PARMIN(100),PAREXT(200),RMSOFT(0:100),
0354      &     AU(3,3),AD(3,3),AE(3,3)
0355       COMMON/PYLH3C/CPRO(2),CVER(2)
0356       CHARACTER CPRO*12,CVER*12
0357       SAVE /PYDAT1/,/PYDAT2/,/PYDAT3/,/PYDAT4/,/PYDATR/,/PYSUBS/,
0358      &/PYPARS/,/PYINT1/,/PYINT2/,/PYINT3/,/PYINT4/,/PYINT5/,
0359      &/PYINT6/,/PYINT7/,/PYMSSM/,/PYSSMT/,/PYMSRV/,/PYTCSM/,
0360      &/PYBINS/,/PYLH3P/,/PYLH3C/
0361
0362 C...PYDAT1, containing status codes and most parameters.
0363       DATA MSTU/
0364      &   0,    0,    0, 4000,10000,  500, 8000,    0,    0,    2,
0365      1   6,    0,    1,    0,    0,    1,    0,    0,    0,    0,
0366      2   2,   10,    0,    0,    1,   10,    0,    0,    0,    0,
0367      3   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
0368      4   2,    2,    1,    4,    2,    1,    1,    0,    0,    0,
0369      5  25,   24,    0,    1,    0,    0,    0,    0,    0,    0,
0370      6   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
0371      7  30*0,
0372      1   1,    0,    0,    0,    0,    0,    0,    0,    0,    0,
0373      2   1,    5,    3,    5,    0,    0,    0,    0,    0,    0,
0374      &  80*0/
0375       DATA (PARU(I),I=1,100)/
0376      &  3.141592653589793D0, 6.283185307179586D0,
0377      &  0.197327D0, 5.06773D0, 0.389380D0, 2.56819D0,  4*0D0,
0378      1  0.001D0, 0.09D0, 0.01D0, 2D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0,
0379      2  0D0,   0D0,   0D0,   0D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0,
0380      3  0D0,   0D0,   0D0,   0D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0,
0381      4  2.0D0,  1.0D0, 0.25D0,  2.5D0, 0.05D0,
0382      4  0D0,   0D0, 0.0001D0, 0D0,   0D0,
0383      5  2.5D0,1.5D0,7.0D0,1.0D0,0.5D0,2.0D0,3.2D0, 0D0, 0D0, 0D0,
0384      6  40*0D0/
0385       DATA (PARU(I),I=101,200)/
0386      &  0.00729735D0, 0.232D0, 0.007764D0, 1.0D0, 1.16639D-5,
0387      &  0D0, 0D0, 0D0, 0D0,  0D0,
0388      1  0.20D0, 0.25D0, 1.0D0, 4.0D0, 10D0, 0D0, 0D0,  0D0, 0D0, 0D0,
0389      2 -0.693D0, -1.0D0, 0.387D0, 1.0D0, -0.08D0,
0390      2 -1.0D0,  1.0D0,  1.0D0,  1.0D0,  0D0,
0391      3  1.0D0,-1.0D0, 1.0D0,-1.0D0, 1.0D0,  0D0,  0D0, 0D0, 0D0, 0D0,
0392      4  5.0D0, 1.0D0, 1.0D0,  0D0, 1.0D0, 1.0D0,  0D0, 0D0, 0D0, 0D0,
0393      5  1.0D0,   0D0,   0D0,   0D0,   0D0,   0D0, 0D0, 0D0, 0D0, 0D0,
0394      6  1.0D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0,  0D0,  0D0, 0D0, 0D0, 0D0,
0395      7  1.0D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0, 0D0,0D0,0D0,
0396      8  1.0D0, 1.0D0, 1.0D0, 0.0D0, 0.0D0, 1.0D0, 1.0D0, 0D0,0D0,0D0,
0397      9  0D0,  0D0,  0D0,  0D0, 1.0D0,  0D0,  0D0, 0D0, 0D0, 0D0/
0398       DATA MSTJ/
0399      &  1,    3,    0,    0,    0,    0,    0,    0,    0,    0,
0400      1  4,    2,    0,    1,    0,    2,    2,   20,    0,    0,
0401      2  2,    1,    1,    2,    1,    2,    2,    0,    0,    0,
0402      3  0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
0403      4  2,    2,    4,    2,    5,    3,    3,    0,    0,    3,
0404      5  0,    3,    0,    2,    0,    0,    1,    0,    0,    0,
0405      6  40*0,
0406      &  5,    2,    7,    5,    1,    1,    0,    2,    0,    2,
0407      1  0,    0,    0,    0,    1,    1,    0,    0,    0,    0,
0408      2  80*0/
0409       DATA PARJ/
0410      &  0.10D0, 0.30D0, 0.40D0, 0.05D0, 0.50D0,
0411      &  0.50D0, 0.50D0,   0.6D0,   1.2D0,   0.6D0,
0412      1  0.50D0,0.60D0,0.75D0, 0D0, 0D0, 0D0, 0D0, 1.0D0, 1.0D0, 0D0,
0413      2  0.36D0, 1.0D0,0.01D0, 2.0D0,1.0D0,0.4D0, 0D0, 0D0, 0D0, 0D0,
0414      3  0.10D0, 1.0D0, 0.8D0, 1.5D0,0D0,2.0D0,0.2D0, 0D0,0.08D0,1D0,
0415      4  0.3D0, 0.58D0, 0.5D0, 0.9D0,0.5D0,1.0D0,1.0D0,1.5D0,1D0,10D0,
0416      5  0.77D0, 0.77D0, 0.77D0, -0.05D0, -0.005D0,
0417      5  0D0, 0D0, 0D0, 1.0D0, 0D0,
0418      6  4.5D0, 0.7D0, 0D0,0.003D0, 0.5D0, 0.5D0, 0D0, 0D0, 0D0, 0D0,
0419      7  10D0, 1000D0, 100D0, 1000D0, 0D0, 0.7D0,10D0, 0D0,0D0,0.5D0,
0420      8  0.29D0, 1.0D0, 1.0D0,  0D0,  10D0, 10D0, 0D0, 0D0, 0D0,1D-4,
0421      9  0.02D0, 1.0D0, 0.2D0,  0D0,  0D0,  0D0,  0D0, 0D0, 0D0, 0D0,
0422      &  0D0,  0D0,  0D0,  0D0,   0D0,   0D0,  0D0,  0D0,  0D0,  0D0,
0423      1  0D0,  0D0,  0D0,  0D0,   0D0,   0D0,  0D0,  0D0,  0D0,  0D0,
0424      2  1.0D0, 0.25D0,91.187D0,2.489D0, 0.01D0,
0425      2  2.0D0,  1.0D0, 0.25D0,0.002D0,   0D0,
0426      3  0D0, 0D0, 0D0, 0D0, 0.01D0, 0.99D0, 0D0, 0D0,  0.2D0,   0D0,
0427      4  10*0D0,
0428      5  10*0D0,
0429      6  10*0D0,
0430      7  0D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0, -0.693D0,
0431      8 -1.0D0, 0.387D0, 1.0D0, -0.08D0, -1.0D0,
0432      8  1.0D0,  1.0D0, -0.693D0, -1.0D0, 0.387D0,
0433      9  1.0D0, -0.08D0, -1.0D0,   1.0D0, 1.0D0,
0434      9  5*0D0/
0435
0436 C...PYDAT2, with particle data and flavour treatment parameters.
0437       DATA (KCHG(I,1),I=   1, 500)/-1,2,-1,2,-1,2,-1,2,2*0,-3,0,-3,0,
0438      &-3,0,-3,6*0,3,9*0,3,2*0,3,4*0,-1,41*0,2,-1,20*0,3*3,7*0,3*3,3*0,
0439      &3*3,3*0,3*3,6*0,3*3,3*0,3*3,4*0,-2,-3,2*1,2*0,4,2*3,6,2*-2,2*-3,
0440      &0,2*1,2*0,2*3,-2,2*-3,2*0,-3,2*1,2*0,3,0,2*4,2*3,2*6,3,2*1,2*0,
0441      &2*3,2*0,4,2*3,2*6,2*3,6,2*-2,2*-3,0,-3,0,2*1,2*0,2*3,0,3,2*-2,
0442      &2*-3,2*0,2*-3,0,2*1,2*0,2*3,2*0,2*3,-2,2*-3,2*0,2*-3,2*0,-3,2*0,
0443      &2*3,4*0,2*3,2*0,2*3,2*0,2*3,4*0,2*3,2*0,2*3,3*0,3,2*0,3,0,3,0,3,
0444      &2*0,3,0,3,3*0,-1,2,-1,2,-1,2,-3,0,-3,0,-3,4*0,3,2*0,3,0,-1,2,-1,
0445      &2,-1,2,-3,0,-3,0,-3,2*0,3,3*0,3,8*0,-1,2,-3,6*0,3,2*6,0,3,4*0,3,
0446      &139*0/
0447       DATA (KCHG(I,2),I=   1, 500)/8*1,12*0,2,20*0,1,107*0,-1,0,2*-1,
0448      &2*0,-1,3*0,2*-1,3*0,2*-1,4*0,-1,5*0,2*-1,4*0,2*-1,5*0,2*-1,6*0,
0449      &-1,7*0,2*-1,5*0,2*-1,6*0,2*-1,7*0,2*-1,8*0,-1,56*0,6*1,6*0,2,7*0,
0450      &6*1,9*0,2,3*0,2,0,5*2,2*1,17*0,6*2,133*0/
0451       DATA (KCHG(I,3),I=   1, 500)/8*1,2*0,8*1,5*0,1,9*0,1,2*0,1,3*0,
0452      &2*1,39*0,1,0,2*1,20*0,3*1,4*0,6*1,3*0,9*1,3*0,12*1,4*0,100*1,2*0,
0453      &2*1,2*0,4*1,2*0,6*1,2*0,8*1,3*0,1,0,2*1,0,3*1,0,4*1,3*0,12*1,3*0,
0454      &1,2*0,1,0,12*1,0,1,3*0,1,8*0,4*1,5*0,3*1,0,1,3*0,2*1,139*0/
0455       DATA (KCHG(I,4),I=   1, 290)/1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,
0456      &16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,
0457      &37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,
0458      &58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,
0459      &79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,
0460      &100,110,111,113,115,130,211,213,215,221,223,225,310,311,313,315,
0461      &321,323,325,331,333,335,411,413,415,421,423,425,431,433,435,441,
0462      &443,445,511,513,515,521,523,525,531,533,535,541,543,545,551,553,
0463      &555,990,1103,1114,2101,2103,2112,2114,2203,2212,2214,2224,3101,
0464      &3103,3112,3114,3122,3201,3203,3212,3214,3222,3224,3303,3312,3314,
0465      &3322,3324,3334,4101,4103,4112,4114,4122,4132,4201,4203,4212,4214,
0466      &4222,4224,4232,4301,4303,4312,4314,4322,4324,4332,4334,4403,4412,
0467      &4414,4422,4424,4432,4434,4444,5101,5103,5112,5114,5122,5132,5142,
0468      &5201,5203,5212,5214,5222,5224,5232,5242,5301,5303,5312,5314,5322,
0469      &5324,5332,5334,5342,5401,5403,5412,5414,5422,5424,5432,5434,5442,
0470      &5444,5503,5512,5514,5522,5524,5532,5534,5542,5544,5554,10111,
0471      &10113,10211,10213,10221,10223,10311,10313,10321,10323,10331,
0472      &10333,10411,10413,10421,10423,10431,10433,10441,10443,10511,
0473      &10513,10521,10523,10531,10533,10541,10543,10551,10553,20113,
0474      &20213,20223,20313,20323,20333,20413,20423,20433,20443,20513/
0475       DATA (KCHG(I,4),I= 291, 500)/20523,20533,20543,20553,100443,
0476      &100553,1000001,1000002,1000003,1000004,1000005,1000006,1000011,
0477      &1000012,1000013,1000014,1000015,1000016,1000021,1000022,1000023,
0478      &1000024,1000025,1000035,1000037,1000039,2000001,2000002,2000003,
0479      &2000004,2000005,2000006,2000011,2000012,2000013,2000014,2000015,
0480      &2000016,3000111,3000211,3000221,3000331,3000113,3000213,3000223,
0481      &3100021,3100111,3200111,3100113,3200113,3300113,3400113,4000001,
0482      &4000002,4000011,4000012,5000039,9900012,9900014,9900016,9900023,
0483      &9900024,9900041,9900042,9900110,9900210,9900220,9900330,9900440,
0484      &9902110,9902210,9900443,9900441,9910441,9900553,9900551,9910551,
0485      &133*0/
0486       DATA (PMAS(I,1),I=   1, 217)/2*0.33D0,0.5D0,1.5D0,4.8D0,175D0,
0487      &2*400D0,2*0D0,0.00051D0,0D0,0.10566D0,0D0,1.777D0,0D0,400D0,
0488      &5*0D0,91.188D0,80.45D0,115D0,6*0D0,500D0,900D0,500D0,3*300D0,
0489      &3*0D0,5000D0,200D0,40*0D0,1D0,2D0,5D0,16*0D0,0.13498D0,0.7685D0,
0490      &1.318D0,0.49767D0,0.13957D0,0.7669D0,1.318D0,0.54745D0,0.78194D0,
0491      &1.275D0,2*0.49767D0,0.8961D0,1.432D0,0.4936D0,0.8916D0,1.425D0,
0492      &0.95777D0,1.0194D0,1.525D0,1.8693D0,2.01D0,2.46D0,1.8645D0,
0493      &2.0067D0,2.46D0,1.9685D0,2.1124D0,2.5735D0,2.9798D0,3.09688D0,
0494      &3.5562D0,5.2792D0,5.3248D0,5.83D0,5.2789D0,5.3248D0,5.83D0,
0495      &5.3693D0,5.4163D0,6.07D0,6.594D0,6.602D0,7.35D0,9.4D0,9.4603D0,
0496      &9.9132D0,0D0,0.77133D0,1.234D0,0.57933D0,0.77133D0,0.93957D0,
0497      &1.233D0,0.77133D0,0.93827D0,1.232D0,1.231D0,0.80473D0,0.92953D0,
0498      &1.19744D0,1.3872D0,1.11568D0,0.80473D0,0.92953D0,1.19255D0,
0499      &1.3837D0,1.18937D0,1.3828D0,1.09361D0,1.3213D0,1.535D0,1.3149D0,
0500      &1.5318D0,1.67245D0,1.96908D0,2.00808D0,2.4521D0,2.5D0,2.2849D0,
0501      &2.4703D0,1.96908D0,2.00808D0,2.4535D0,2.5D0,2.4529D0,2.5D0,
0502      &2.4656D0,2.15432D0,2.17967D0,2.55D0,2.63D0,2.55D0,2.63D0,2.704D0,
0503      &2.8D0,3.27531D0,3.59798D0,3.65648D0,3.59798D0,3.65648D0,
0504      &3.78663D0,3.82466D0,4.91594D0,5.38897D0,5.40145D0,5.8D0,5.81D0,
0505      &5.641D0,5.84D0,7.00575D0,5.38897D0,5.40145D0,5.8D0,5.81D0,5.8D0/
0506       DATA (PMAS(I,1),I= 218, 500)/5.81D0,5.84D0,7.00575D0,5.56725D0,
0507      &5.57536D0,5.96D0,5.97D0,5.96D0,5.97D0,6.12D0,6.13D0,7.19099D0,
0508      &6.67143D0,6.67397D0,7.03724D0,7.0485D0,7.03724D0,7.0485D0,
0509      &7.21101D0,7.219D0,8.30945D0,8.31325D0,10.07354D0,10.42272D0,
0510      &10.44144D0,10.42272D0,10.44144D0,10.60209D0,10.61426D0,
0511      &11.70767D0,11.71147D0,15.11061D0,0.9835D0,1.231D0,0.9835D0,
0512      &1.231D0,1D0,1.17D0,1.429D0,1.29D0,1.429D0,1.29D0,2*1.4D0,2.272D0,
0513      &2.424D0,2.272D0,2.424D0,2.5D0,2.536D0,3.4151D0,3.46D0,5.68D0,
0514      &5.73D0,5.68D0,5.73D0,5.92D0,5.97D0,7.25D0,7.3D0,9.8598D0,9.875D0,
0515      &2*1.23D0,1.282D0,2*1.402D0,1.427D0,2*2.372D0,2.56D0,3.5106D0,
0516      &2*5.78D0,6.02D0,7.3D0,9.8919D0,3.686D0,10.0233D0,32*500D0,
0517      &3*110D0,350D0,3*210D0,500D0,125D0,250D0,400D0,2*350D0,300D0,
0518      &4*400D0,1000D0,3*500D0,1200D0,750D0,2*200D0,7*0D0,3*3.1D0,
0519      &3*9.5D0,133*0D0/
0520       DATA (PMAS(I,2),I=   1, 500)/5*0D0,1.39816D0,16*0D0,2.47813D0,
0521      &2.07115D0,0.00367D0,6*0D0,14.54029D0,0D0,16.66099D0,8.38842D0,
0522      &3.3752D0,4.17669D0,3*0D0,417.29147D0,0.39162D0,60*0D0,0.151D0,
0523      &0.107D0,2*0D0,0.149D0,0.107D0,0D0,0.00843D0,0.185D0,2*0D0,
0524      &0.0505D0,0.109D0,0D0,0.0498D0,0.098D0,0.0002D0,0.00443D0,0.076D0,
0525      &2*0D0,0.023D0,2*0D0,0.023D0,2*0D0,0.015D0,0.0013D0,0D0,0.002D0,
0526      &2*0D0,0.02D0,2*0D0,0.02D0,2*0D0,0.02D0,2*0D0,0.02D0,5*0D0,0.12D0,
0527      &3*0D0,0.12D0,2*0D0,2*0.12D0,3*0D0,0.0394D0,4*0D0,0.036D0,0D0,
0528      &0.0358D0,2*0D0,0.0099D0,0D0,0.0091D0,74*0D0,0.06D0,0.142D0,
0529      &0.06D0,0.142D0,0D0,0.36D0,0.287D0,0.09D0,0.287D0,0.09D0,0.25D0,
0530      &0.08D0,0.05D0,0.02D0,0.05D0,0.02D0,0.05D0,0D0,0.014D0,0.01D0,
0531      &8*0.05D0,0D0,0.01D0,2*0.4D0,0.025D0,2*0.174D0,0.053D0,3*0.05D0,
0532      &0.0009D0,4*0.05D0,3*0D0,19*1D0,0D0,7*1D0,0D0,1D0,0D0,1D0,0D0,
0533      &0.02911D0,0.01741D0,0.04536D0,0.09511D0,0.8686D0,0.62395D0,
0534      &0.19192D0,123.27638D0,0.02296D0,0.18886D0,23.26819D0,2.86306D0,
0535      &0D0,3.45903D0,2.59359D0,2.59687D0,0.42896D0,0.41912D0,0.14153D0,
0536      &2*0.00098D0,0.00097D0,26.7245D0,21.74916D0,0.88159D0,0.88001D0,
0537      &7*0D0,6*0.01D0,133*0D0/
0538       DATA (PMAS(I,3),I=   1, 500)/5*0D0,13.98156D0,16*0D0,24.78129D0,
0539      &20.71149D0,0.03669D0,6*0D0,145.40294D0,0D0,166.60993D0,
0540      &83.88423D0,33.75195D0,41.76694D0,3*0D0,4172.91467D0,3.91621D0,
0541      &60*0D0,0.4D0,0.25D0,2*0D0,0.4D0,0.25D0,0D0,0.1D0,0.17D0,2*0D0,
0542      &0.2D0,0.12D0,0D0,0.2D0,0.12D0,0.002D0,0.015D0,0.2D0,2*0D0,0.12D0,
0543      &2*0D0,0.12D0,2*0D0,0.05D0,0.005D0,0D0,0.01D0,2*0D0,0.05D0,2*0D0,
0544      &0.05D0,2*0D0,0.05D0,2*0D0,0.05D0,5*0D0,0.14D0,3*0D0,0.14D0,2*0D0,
0545      &2*0.14D0,3*0D0,0.04D0,4*0D0,0.035D0,0D0,0.035D0,2*0D0,0.05D0,0D0,
0546      &0.05D0,74*0D0,0.05D0,0.25D0,0.05D0,0.25D0,0D0,0.2D0,0.4D0,
0547      &0.005D0,0.4D0,0.01D0,0.35D0,0.001D0,0.1D0,0.08D0,0.1D0,0.08D0,
0548      &0.1D0,0D0,0.05D0,0.02D0,6*0.1D0,0.05D0,0.1D0,0D0,0.02D0,2*0.3D0,
0549      &0.05D0,2*0.3D0,0.02D0,2*0.1D0,0.03D0,0.001D0,4*0.1D0,3*0D0,
0550      &19*10D0,0.00001D0,7*10D0,0.00001D0,10D0,0.00001D0,10D0,0.00001D0,
0551      &0.29108D0,0.17412D0,0.45362D0,0.95114D0,8.68604D0,6.23946D0,
0552      &1.91923D0,450D0,0.22959D0,1.88863D0,232.68185D0,28.63059D0,0D0,
0553      &34.59032D0,25.93594D0,25.96873D0,4.28961D0,4.19124D0,1.41528D0,
0554      &0.00977D0,0.00976D0,0.00973D0,267.24501D0,217.49162D0,8.81592D0,
0555      &8.80013D0,13*0D0,133*0D0/
0556       DATA (PMAS(I,4),I=   1, 500)/12*0D0,658654D0,0D0,0.0872D0,68*0D0,
0557      &0.1D0,0.387D0,16*0D0,0.00003D0,2*0D0,15500D0,7804.5D0,5*0D0,
0558      &26.762D0,3*0D0,3709D0,5*0D0,0.317D0,2*0D0,0.1244D0,2*0D0,0.14D0,
0559      &5*0D0,0.468D0,2*0D0,0.462D0,2*0D0,0.483D0,2*0D0,0.15D0,18*0D0,
0560      &44.34D0,0D0,78.88D0,4*0D0,23.96D0,2*0D0,49.1D0,0D0,87.1D0,0D0,
0561      &24.6D0,4*0D0,0.0618D0,0.029D0,6*0D0,0.106D0,6*0D0,0.019D0,2*0D0,
0562      &7*0.1D0,4*0D0,0.342D0,2*0.387D0,6*0D0,2*0.387D0,6*0D0,0.387D0,
0563      &0D0,0.387D0,2*0D0,8*0.387D0,0D0,9*0.387D0,118*0D0,133*0D0/
0564       DATA PARF/
0565      &  0.5D0,0.25D0, 0.5D0,0.25D0, 1D0, 0.5D0,  0D0,  0D0,  0D0, 0D0,
0566      1  0.5D0,  0D0, 0.5D0,  0D0,  1D0,  1D0,  0D0,  0D0,  0D0, 0D0,
0567      2  0.5D0,  0D0, 0.5D0,  0D0,  1D0,  1D0,  0D0,  0D0,  0D0, 0D0,
0568      3  0.5D0,  0D0, 0.5D0,  0D0,  1D0,  1D0,  0D0,  0D0,  0D0, 0D0,
0569      4  0.5D0,  0D0, 0.5D0,  0D0,  1D0,  1D0,  0D0,  0D0,  0D0, 0D0,
0570      5  0.5D0,  0D0, 0.5D0,  0D0,  1D0,  1D0,  0D0,  0D0,  0D0, 0D0,
0571      6  0.75D0, 0.5D0, 0D0,0.1667D0,0.0833D0,0.1667D0,0D0,0D0,0D0, 0D0,
0572      7  0D0,  0D0,  1D0,0.3333D0,0.6667D0,0.3333D0,0D0,0D0,0D0, 0D0,
0573      8  0D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0, 0D0,
0574      9  0.0099D0, 0.0056D0, 0.199D0, 1.23D0, 4.17D0, 165D0,  4*0D0,
0575      & 0.325D0,0.325D0,0.5D0,1.6D0, 5.0D0,  0D0,  0D0,  0D0,  0D0, 0D0,
0576      1 0D0,0.11D0,0.16D0,0.048D0,0.50D0,0.45D0,0.55D0,0.60D0,0D0,0D0,
0577      2 0.2D0, 0.1D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0, 0D0,
0578      3 60*0D0,
0579      4 0.2D0,  0.5D0,  8*0D0,
0580      5 1800*0D0/
0581       DATA ((VCKM(I,J),J=1,4),I=1,4)/
0582      &  0.95113D0,  0.04884D0,  0.00003D0,  0.00000D0,
0583      &  0.04884D0,  0.94940D0,  0.00176D0,  0.00000D0,
0584      &  0.00003D0,  0.00176D0,  0.99821D0,  0.00000D0,
0585      &  0.00000D0,  0.00000D0,  0.00000D0,  1.00000D0/
0586
0587 C...PYDAT3, with particle decay parameters and data.
0588
0589       DATA (MDCY(I,1),I=   1, 500)/5*0,3*1,6*0,1,0,1,5*0,3*1,6*0,1,0,
0590      &4*1,3*0,2*1,40*0,3*1,16*0,3*1,2*0,9*1,0,32*1,2*0,1,3*0,1,2*0,2*1,
0591      &2*0,3*1,2*0,4*1,0,5*1,2*0,4*1,2*0,5*1,2*0,6*1,0,7*1,2*0,5*1,2*0,
0592      &6*1,2*0,7*1,2*0,8*1,0,75*1,0,7*1,0,1,0,1,0,26*1,7*0,6*1,133*0/
0593       DATA (MDCY(I,2),I=   1, 351)/1,9,17,25,33,41,56,66,2*0,76,80,82,
0594      &87,89,143,145,150,2*0,153,162,174,190,210,6*0,289,0,311,334,420,
0595      &503,3*0,530,539,40*0,540,541,545,16*0,554,556,561,570,579,581,
0596      &583,590,598,604,613,615,617,620,630,636,639,650,656,667,673,736,
0597      &739,747,808,810,818,851,853,857,858,861,863,899,900,908,944,945,
0598      &953,992,993,997,1028,1029,1033,1034,1043,2*0,1045,3*0,1046,2*0,
0599      &1049,1052,2*0,1053,1055,1058,2*0,1062,1063,1066,1069,0,1072,1077,
0600      &1079,1082,1084,2*0,1088,1089,1090,1166,2*0,1170,1171,1172,1173,
0601      &1174,2*0,1178,1179,1181,1182,1184,1188,0,1189,1193,1197,1201,
0602      &1205,1209,1213,2*0,1217,1218,1219,1236,1245,2*0,1254,1255,1256,
0603      &1257,1258,1267,2*0,1276,1277,1278,1279,1280,1289,1290,2*0,1299,
0604      &1308,1317,1326,1335,1344,1353,1362,0,1371,1380,1389,1398,1407,
0605      &1416,1425,1434,1443,1452,1453,1454,1455,1456,1461,1464,1466,1471,
0606      &1473,1478,1485,1489,1491,1493,1495,1497,1499,1501,1503,1504,1506,
0607      &1508,1510,1512,1514,1516,1518,1520,1522,1523,1525,1527,1541,1543,
0608      &1545,1549,1551,1553,1555,1557,1559,1561,1563,1565,1567,1578,1592,
0609      &1637,1661,1706,1730,1775,1802,1833,1859,1891,1917,1949,1975,2162,
0610      &2331,2595,2826,3106,3402,0,3657,3706,3734,3783,3811,3860,3888,0,
0611      &3924,0,3960,0,3996,4004,4012,4020,4023,4047,4073,4097,4103,4110,
0612      &4117,4124,4130,4136,4145,4149,4153,4156,4158,4178,4200,4222,4244/
0613       DATA (MDCY(I,2),I= 352, 500)/4259,4271,4278,7*0,4285,4286,4287,
0614      &4288,4289,4290,133*0/
0615       DATA (MDCY(I,3),I=   1, 500)/5*8,15,2*10,2*0,4,2,5,2,54,2,5,3,
0616      &2*0,9,12,16,20,79,6*0,22,0,23,86,83,27,3*0,9,1,40*0,1,4,9,16*0,2,
0617      &5,2*9,2*2,7,8,6,9,2*2,3,10,6,3,11,6,11,6,63,3,8,61,2,8,33,2,4,1,
0618      &3,2,36,1,8,36,1,8,39,1,4,31,1,4,1,9,2,2*0,1,3*0,3,2*0,3,1,2*0,2,
0619      &3,4,2*0,1,3*3,0,5,2,3,2,4,2*0,2*1,76,4,2*0,4*1,4,2*0,1,2,1,2,4,1,
0620      &0,7*4,2*0,2*1,17,2*9,2*0,4*1,2*9,2*0,4*1,9,1,9,2*0,8*9,0,9*9,4*1,
0621      &5,3,2,5,2,5,7,4,7*2,1,9*2,1,2*2,14,2*2,4,9*2,11,14,45,24,45,24,
0622      &45,27,31,26,32,26,32,26,187,169,264,231,280,296,255,0,49,28,49,
0623      &28,49,28,36,0,36,0,36,0,3*8,3,24,26,24,6,3*7,2*6,9,2*4,3,2,20,
0624      &3*22,15,12,2*7,7*0,6*1,133*0/
0625       DATA (MDME(I,1),I=   1,8000)/6*1,-1,7*1,-1,7*1,-1,7*1,-1,7*1,-1,
0626      &7*1,-1,1,7*-1,8*1,2*-1,8*1,2*-1,73*1,-1,2*1,-1,5*1,0,2*-1,6*1,0,
0627      &2*-1,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,
0628      &2*-1,6*1,8*-1,3*1,-1,3*1,-1,3*1,5*-1,3*1,4*-1,200*1,2*-1,2*1,-1,
0629      &1249*1,2*-1,377*1,2*-1,1868*1,2*-1,6*1,2*-1,9*1,-1,3*1,-1,3*1,
0630      &5*-1,3*1,-1,14*1,2*-1,6*1,2*-1,67*1,2*-1,6*1,2*-1,117*1,3710*0/
0631       DATA (MDME(I,2),I=   1,8000)/43*102,4*0,102,0,6*53,3*102,4*0,102,
0632      &2*0,3*102,4*0,102,2*0,6*102,42,6*102,2*42,2*0,8*41,2*0,36*41,
0633      &8*102,0,102,0,102,2*0,21*102,8*32,8*0,16*32,4*0,8*32,9*0,62*53,
0634      &8*32,14*0,16*32,7*0,8*32,16*0,62*53,8*32,13*0,62*53,4*32,5*0,
0635      &18*53,6*32,4*0,12,2*42,2*11,9*42,0,2,3,15*0,4*42,5*0,3,12*0,2,
0636      &3*0,1,0,3,16*0,2*3,15*0,2*42,2*3,18*0,2*3,3*0,1,11*0,22*42,41*0,
0637      &2*3,9*0,16*42,45*0,3,10*0,10*42,20*0,2*13,6*0,12,2*0,12,0,12,
0638      &14*42,16*0,48,3*13,2*42,9*0,14*42,16*0,48,3*13,2*42,9*0,14*42,
0639      &19*0,48,3*13,2*42,6*0,2*11,28*42,5*0,32,3*0,4*32,2*4,0,32,45*0,
0640      &14*42,52*0,10*13,2*42,2*11,4*0,2*42,2*11,6*0,2*42,2*11,0,2*42,
0641      &2*11,2*42,2*11,2*42,2*11,2*42,2*11,2*42,2*11,2*42,2*11,2*42,2*11,
0642      &2*0,3*42,8*0,48,3*13,20*42,4*0,18*42,4*0,9*42,0,162*42,50*0,2*12,
0643      &17*0,2*32,33*0,12,9*0,32,2*0,12,11*0,4*32,2*4,5*0,2404*53,4*32,
0644      &3*0,6*32,3*0,4*32,3*0,4*32,8*0,8*32,14*0,16*32,12*0,8*32,8*0,
0645      &46*32,3*53,12*0,8*32,12*0,66*51,6*32,9*0,9*32,17*0,6*51,3710*0/
0646       DATA (BRAT(I)  ,I=   1, 346)/43*0D0,0.00003D0,0.001765D0,
0647      &0.998205D0,35*0D0,1D0,6*0D0,0.1783D0,0.1735D0,0.1131D0,0.2494D0,
0648      &0.003D0,0.09D0,0.0027D0,0.01D0,0.0014D0,0.0012D0,2*0.00025D0,
0649      &0.0071D0,0.012D0,0.0004D0,0.00075D0,0.00006D0,2*0.00078D0,
0650      &0.0034D0,0.08D0,0.011D0,0.0191D0,0.00006D0,0.005D0,0.0133D0,
0651      &0.0067D0,0.0005D0,0.0035D0,0.0006D0,0.0015D0,0.00021D0,0.0002D0,
0652      &0.00075D0,0.0001D0,0.0002D0,0.0011D0,3*0.0002D0,0.00022D0,
0653      &0.0004D0,0.0001D0,2*0.00205D0,2*0.00069D0,0.00025D0,0.00051D0,
0654      &0.00025D0,35*0D0,0.153995D0,0.11942D0,0.153984D0,0.119259D0,
0655      &0.152272D0,3*0D0,0.033576D0,0.066806D0,0.033576D0,0.066806D0,
0656      &0.0335D0,0.066806D0,2*0D0,0.321369D0,0.016494D0,2*0D0,0.016502D0,
0657      &0.320615D0,2*0D0,0.00001D0,0.000591D0,6*0D0,2*0.108166D0,
0658      &0.108087D0,0D0,0.000001D0,0D0,0.000349D0,0.048707D0,0.768308D0,
0659      &4*0D0,0.000227D0,0.064048D0,0D0,0.040621D0,0.002043D0,0.000615D0,
0660      &0.006981D0,0.068099D0,62*0D0,0.145835D0,0.113276D0,0.145835D0,
0661      &0.113271D0,0.145781D0,0.049002D0,2*0D0,0.032025D0,0.063642D0,
0662      &0.032025D0,0.063642D0,0.032022D0,0.063642D0,8*0D0,0.251225D0,
0663      &0.0129D0,0.000006D0,0D0,0.0129D0,0.250764D0,0.00038D0,0D0,
0664      &0.000008D0,0.000465D0,0.215418D0,5*0D0,2*0.085312D0,0.08531D0,
0665      &7*0D0,0.000049D0,0.000774D0,5*0D0,0.000074D0,0D0,0.000417D0/
0666       DATA (BRAT(I)  ,I= 347, 651)/0.000015D0,0.000061D0,0.30671D0,
0667      &0.689011D0,0D0,0.002889D0,69*0D0,0.000001D0,0.000121D0,
0668      &0.001924D0,4*0D0,0.000001D0,0.000184D0,0D0,0.003106D0,0.000015D0,
0669      &0.000003D0,2*0D0,0.994646D0,66*0D0,0.000021D0,0.090135D0,2*0D0,
0670      &0.000013D0,0.003714D0,0D0,0.906117D0,18*0D0,3*0.215119D0,
0671      &0.214724D0,2*0D0,0.06996D0,0.069959D0,0D0,2*1D0,2*0.08D0,0.76D0,
0672      &0.08D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,
0673      &0.005D0,0.988D0,0.012D0,0.998739D0,0.00079D0,0.00038D0,
0674      &0.000046D0,0.000045D0,2*0.34725D0,0.144D0,0.104D0,0.0245D0,
0675      &2*0.01225D0,0.0028D0,0.0057D0,0.2112D0,0.1256D0,2*0.1939D0,
0676      &2*0.1359D0,0.002D0,0.001D0,0.0006D0,0.999877D0,0.000123D0,
0677      &0.99955D0,0.00045D0,2*0.34725D0,0.144D0,0.104D0,0.049D0,0.0028D0,
0678      &0.0057D0,0.3923D0,0.321D0,0.2317D0,0.0478D0,0.0049D0,0.0013D0,
0679      &0.0003D0,0.0007D0,0.89D0,0.08693D0,0.0221D0,0.00083D0,
0680      &2*0.00007D0,0.564D0,0.282D0,0.072D0,0.028D0,0.023D0,2*0.0115D0,
0681      &0.005D0,0.003D0,0.6861D0,0.3139D0,2*0.5D0,0.665D0,0.333D0,
0682      &0.002D0,0.333D0,0.166D0,0.168D0,0.084D0,0.087D0,0.043D0,0.059D0,
0683      &2*0.029D0,0.002D0,0.6352D0,0.2116D0,0.0559D0,0.0173D0,0.0482D0,
0684      &0.0318D0,0.666D0,0.333D0,0.001D0,0.332D0,0.166D0,0.168D0,0.084D0,
0685      &0.086D0,0.043D0,0.059D0,2*0.029D0,2*0.002D0,0.437D0,0.208D0/
0686       DATA (BRAT(I)  ,I= 652, 823)/0.302D0,0.0302D0,0.0212D0,0.0016D0,
0687      &0.48947D0,0.34D0,3*0.043D0,0.027D0,0.0126D0,0.0013D0,0.0003D0,
0688      &0.00025D0,0.00008D0,0.444D0,2*0.222D0,0.104D0,2*0.004D0,0.07D0,
0689      &0.065D0,2*0.005D0,2*0.011D0,5*0.001D0,0.07D0,0.065D0,2*0.005D0,
0690      &2*0.011D0,5*0.001D0,0.026D0,0.019D0,0.066D0,0.041D0,0.045D0,
0691      &0.076D0,0.0073D0,2*0.0047D0,0.026D0,0.001D0,0.0006D0,0.0066D0,
0692      &0.005D0,2*0.003D0,2*0.0006D0,2*0.001D0,0.006D0,0.005D0,0.012D0,
0693      &0.0057D0,0.067D0,0.008D0,0.0022D0,0.027D0,0.004D0,0.019D0,
0694      &0.012D0,0.002D0,0.009D0,0.0218D0,0.001D0,0.022D0,0.087D0,0.001D0,
0695      &0.0019D0,0.0015D0,0.0028D0,0.683D0,0.306D0,0.011D0,0.3D0,0.15D0,
0696      &0.16D0,0.08D0,0.13D0,0.06D0,0.08D0,0.04D0,0.034D0,0.027D0,
0697      &2*0.002D0,2*0.004D0,2*0.002D0,0.034D0,0.027D0,2*0.002D0,
0698      &2*0.004D0,2*0.002D0,0.0365D0,0.045D0,0.073D0,0.062D0,3*0.021D0,
0699      &0.0061D0,0.015D0,0.025D0,0.0088D0,0.074D0,0.0109D0,0.0041D0,
0700      &0.002D0,0.0035D0,0.0011D0,0.001D0,0.0027D0,2*0.0016D0,0.0018D0,
0701      &0.011D0,0.0063D0,0.0052D0,0.018D0,0.016D0,0.0034D0,0.0036D0,
0702      &0.0009D0,0.0006D0,0.015D0,0.0923D0,0.018D0,0.022D0,0.0077D0,
0703      &0.009D0,0.0075D0,0.024D0,0.0085D0,0.067D0,0.0511D0,0.017D0,
0704      &0.0004D0,0.0028D0,0.619D0,0.381D0,0.3D0,0.15D0,0.16D0,0.08D0,
0705      &0.13D0,0.06D0,0.08D0,0.04D0,0.01D0,2*0.02D0,0.03D0,2*0.005D0/
0706       DATA (BRAT(I)  ,I= 824, 991)/2*0.02D0,0.03D0,2*0.005D0,0.015D0,
0707      &0.037D0,0.028D0,0.079D0,0.095D0,0.052D0,0.0078D0,4*0.001D0,
0708      &0.028D0,0.033D0,0.026D0,0.05D0,0.01D0,4*0.005D0,0.25D0,0.0952D0,
0709      &0.94D0,0.06D0,2*0.4D0,2*0.1D0,1D0,0.0602D0,0.0601D0,0.8797D0,
0710      &0.135D0,0.865D0,0.02D0,0.055D0,2*0.005D0,0.008D0,0.012D0,0.02D0,
0711      &0.055D0,2*0.005D0,0.008D0,0.012D0,0.01D0,0.03D0,0.0035D0,0.011D0,
0712      &0.0055D0,0.0042D0,0.009D0,0.018D0,0.015D0,0.0185D0,0.0135D0,
0713      &0.025D0,0.0004D0,0.0007D0,0.0008D0,0.0014D0,0.0019D0,0.0025D0,
0714      &0.4291D0,0.08D0,0.07D0,0.02D0,0.015D0,0.005D0,1D0,0.3D0,0.15D0,
0715      &0.16D0,0.08D0,0.13D0,0.06D0,0.08D0,0.04D0,0.02D0,0.055D0,
0716      &2*0.005D0,0.008D0,0.012D0,0.02D0,0.055D0,2*0.005D0,0.008D0,
0717      &0.012D0,0.01D0,0.03D0,0.0035D0,0.011D0,0.0055D0,0.0042D0,0.009D0,
0718      &0.018D0,0.015D0,0.0185D0,0.0135D0,0.025D0,0.0004D0,0.0007D0,
0719      &0.0008D0,0.0014D0,0.0019D0,0.0025D0,0.4291D0,0.08D0,0.07D0,
0720      &0.02D0,0.015D0,0.005D0,1D0,0.3D0,0.15D0,0.16D0,0.08D0,0.13D0,
0721      &0.06D0,0.08D0,0.04D0,0.02D0,0.055D0,2*0.005D0,0.008D0,0.012D0,
0722      &0.02D0,0.055D0,2*0.005D0,0.008D0,0.012D0,0.01D0,0.03D0,0.0035D0,
0723      &0.011D0,0.0055D0,0.0042D0,0.009D0,0.018D0,0.015D0,0.0185D0,
0724      &0.0135D0,0.025D0,2*0.0002D0,0.0007D0,2*0.0004D0,0.0014D0,0.001D0,
0725      &0.0009D0,0.0025D0,0.4291D0,0.08D0,0.07D0,0.02D0,0.015D0,0.005D0/
0726       DATA (BRAT(I)  ,I= 992,1183)/1D0,2*0.3D0,2*0.2D0,0.047D0,0.122D0,
0727      &0.006D0,0.012D0,0.035D0,0.012D0,0.035D0,0.003D0,0.007D0,0.15D0,
0728      &0.037D0,0.008D0,0.002D0,0.05D0,0.015D0,0.003D0,0.001D0,0.014D0,
0729      &0.042D0,0.014D0,0.042D0,0.24D0,0.065D0,0.012D0,0.003D0,0.001D0,
0730      &0.002D0,0.001D0,0.002D0,0.014D0,0.003D0,1D0,2*0.3D0,2*0.2D0,1D0,
0731      &0.0252D0,0.0248D0,0.0267D0,0.015D0,0.045D0,0.015D0,0.045D0,
0732      &0.7743D0,0.029D0,0.22D0,0.78D0,1D0,0.331D0,0.663D0,0.006D0,
0733      &0.663D0,0.331D0,0.006D0,1D0,0.999D0,0.001D0,0.88D0,2*0.06D0,
0734      &0.639D0,0.358D0,0.002D0,0.001D0,1D0,0.88D0,2*0.06D0,0.516D0,
0735      &0.483D0,0.001D0,0.88D0,2*0.06D0,0.9988D0,0.0001D0,0.0006D0,
0736      &0.0004D0,0.0001D0,0.667D0,0.333D0,0.9954D0,0.0011D0,0.0035D0,
0737      &0.333D0,0.667D0,0.676D0,0.234D0,0.085D0,0.005D0,2*1D0,0.018D0,
0738      &2*0.005D0,0.003D0,0.002D0,2*0.006D0,0.018D0,2*0.005D0,0.003D0,
0739      &0.002D0,2*0.006D0,0.0066D0,0.025D0,0.016D0,0.0088D0,2*0.005D0,
0740      &0.0058D0,0.005D0,0.0055D0,4*0.004D0,2*0.002D0,2*0.004D0,0.003D0,
0741      &0.002D0,2*0.003D0,3*0.002D0,2*0.001D0,0.002D0,2*0.001D0,
0742      &2*0.002D0,0.0013D0,0.0018D0,5*0.001D0,4*0.003D0,2*0.005D0,
0743      &2*0.002D0,2*0.001D0,2*0.002D0,2*0.001D0,0.2432D0,0.057D0,
0744      &2*0.035D0,0.15D0,2*0.075D0,0.03D0,2*0.015D0,2*0.08D0,0.76D0,
0745      &0.08D0,4*1D0,2*0.08D0,0.76D0,0.08D0,1D0,2*0.5D0,1D0,2*0.5D0/
0746       DATA (BRAT(I)  ,I=1184,1377)/2*0.08D0,0.76D0,0.08D0,1D0,2*0.08D0,
0747      &0.76D0,3*0.08D0,0.76D0,3*0.08D0,0.76D0,3*0.08D0,0.76D0,3*0.08D0,
0748      &0.76D0,3*0.08D0,0.76D0,3*0.08D0,0.76D0,0.08D0,2*1D0,2*0.105D0,
0749      &0.04D0,0.0077D0,0.02D0,0.0235D0,0.0285D0,0.0435D0,0.0011D0,
0750      &0.0022D0,0.0044D0,0.4291D0,0.08D0,0.07D0,0.02D0,0.015D0,0.005D0,
0751      &2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,
0752      &2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,
0753      &4*1D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,
0754      &0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,
0755      &0.005D0,4*1D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
0756      &0.015D0,0.005D0,1D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
0757      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
0758      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
0759      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
0760      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
0761      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
0762      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
0763      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
0764      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
0765      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0/
0766       DATA (BRAT(I)  ,I=1378,1580)/0.015D0,0.005D0,2*0.105D0,0.04D0,
0767      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,0.04D0,
0768      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,0.04D0,
0769      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,0.04D0,
0770      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,0.04D0,
0771      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,0.04D0,
0772      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,0.04D0,
0773      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,0.04D0,
0774      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,4*1D0,0.52D0,0.26D0,
0775      &0.11D0,2*0.055D0,0.333D0,0.334D0,0.333D0,0.667D0,0.333D0,0.28D0,
0776      &0.14D0,0.313D0,0.157D0,0.11D0,0.667D0,0.333D0,0.28D0,0.14D0,
0777      &0.313D0,0.157D0,0.11D0,0.36D0,0.18D0,0.03D0,2*0.015D0,2*0.2D0,
0778      &4*0.25D0,0.667D0,0.333D0,0.667D0,0.333D0,0.667D0,0.333D0,0.667D0,
0779      &0.333D0,4*0.5D0,0.007D0,0.993D0,1D0,0.667D0,0.333D0,0.667D0,
0780      &0.333D0,0.667D0,0.333D0,0.667D0,0.333D0,8*0.5D0,0.02D0,0.98D0,
0781      &1D0,4*0.5D0,3*0.146D0,3*0.05D0,0.15D0,2*0.05D0,4*0.024D0,0.066D0,
0782      &0.667D0,0.333D0,0.667D0,0.333D0,4*0.25D0,0.667D0,0.333D0,0.667D0,
0783      &0.333D0,2*0.5D0,0.273D0,0.727D0,0.667D0,0.333D0,0.667D0,0.333D0,
0784      &4*0.5D0,0.35D0,0.65D0,2*0.0083D0,0.1866D0,0.324D0,0.184D0,
0785      &0.027D0,0.001D0,0.093D0,0.087D0,0.078D0,0.0028D0,3*0.014D0/
0786       DATA (BRAT(I)  ,I=1581,4149)/0.008D0,0.024D0,0.008D0,0.024D0,
0787      &0.425D0,0.02D0,0.185D0,0.088D0,0.043D0,0.067D0,0.066D0,2404*0D0,
0788      &0.017431D0,0.054048D0,0.857694D0,2*0D0,0.00025D0,0.070578D0,0D0,
0789      &0.022748D0,0.026576D0,0.359486D0,0.561581D0,2*0D0,0.000104D0,
0790      &0.029504D0,0.011185D0,0.034681D0,0.550354D0,2*0D0,0.00016D0,
0791      &0.045287D0,0.358333D0,0.445781D0,0D0,0.554219D0,0.144051D0,
0792      &2*0.351902D0,0D0,0.082107D0,0.029566D0,0.001511D0,0.000726D0,
0793      &0.004518D0,0.006522D0,0.004518D0,0.006522D0,0.004513D0,3*0D0,
0794      &0.002908D0,0.000973D0,0.002908D0,0.000973D0,0.002908D0,
0795      &0.000973D0,2*0D0,0.143982D0,0.489888D0,0.1951D0,0D0,0.114302D0,
0796      &0.008426D0,0.014868D0,0.000763D0,2*0D0,0.000763D0,0.01484D0,
0797      &0.000003D0,2*0D0,0.000027D0,0.001945D0,5*0D0,3*0.00503D0,0D0,
0798      &0.133776D0,0.003284D0,0.37169D0,0.006838D0,2*0.030954D0,
0799      &0.00163D0,0D0,0.047224D0,0.073737D0,0.047224D0,0.073732D0,
0800      &0.047179D0,3*0D0,0.034761D0,0.009166D0,0.034761D0,0.009166D0,
0801      &0.034759D0,0.009166D0,2*0D0,4*0.009069D0,0.510147D0,0.453576D0,
0802      &6*0D0,1D0,6*0D0,1D0,4*0.001128D0,0.571047D0,0.382288D0,
0803      &0.042153D0,4*0.016597D0,0.93361D0,0D0,4*0.016597D0,0.93361D0,0D0,
0804      &4*0.05515D0,0.34469D0,0D0,0.228998D0,0.164208D0,0.041503D0,
0805      &0.850973D0,0.005411D0,0.045025D0,0.098591D0,0.849898D0/
0806       DATA (BRAT(I)  ,I=4150,4280)/0.021617D0,0.030018D0,0.098466D0,
0807      &0.294448D0,0.10945D0,0.596102D0,0.389906D0,0.610094D0,3*0.0633D0,
0808      &0.063299D0,0.063295D0,0.056281D0,2*0D0,6*0.020495D0,2*0D0,
0809      &0.327919D0,0.04099D0,0.045236D0,0.090112D0,0.19874D0,0.010204D0,
0810      &0.000003D0,0.010205D0,0.198356D0,0.000151D0,0.000006D0,
0811      &0.000367D0,0.081967D0,0.19874D0,0.010204D0,0.000003D0,0.010205D0,
0812      &0.198356D0,0.000151D0,0.000006D0,0.000367D0,0.081967D0,4*0D0,
0813      &0.198776D0,0.010206D0,0.000003D0,0.010207D0,0.19839D0,0.000151D0,
0814      &0.000006D0,0.000367D0,0.081893D0,0.198776D0,0.010206D0,
0815      &0.000003D0,0.010207D0,0.19839D0,0.000151D0,0.000006D0,0.000367D0,
0816      &0.081893D0,4*0D0,0.199344D0,0.010234D0,0.000003D0,0.010236D0,
0817      &0.198928D0,0.000149D0,0.000006D0,0.000368D0,0.080733D0,
0818      &0.199344D0,0.010234D0,0.000003D0,0.010236D0,0.198928D0,
0819      &0.000149D0,0.000006D0,0.000368D0,0.080733D0,4*0D0,0.184738D0,
0820      &0.104588D0,0.184738D0,0.104587D0,0.184731D0,0.09582D0,0.022902D0,
0821      &0.008429D0,0.015602D0,0.022902D0,0.008429D0,0.015602D0,
0822      &0.022902D0,0.008429D0,0.015602D0,0.28959D0,0.01487D0,0.000008D0,
0823      &0.01487D0,0.289061D0,0.000492D0,0.000009D0,0.000536D0,0.27911D0,
0824      &2*0.037151D0,0.03715D0,0.090266D0,2*0.001805D0,0.090266D0,
0825      &0.001805D0,0.812263D0,0.00179D0,0.090428D0,0.001809D0,0.001808D0/
0826       DATA (BRAT(I)  ,I=4281,8000)/0.090428D0,0.001808D0,0.81372D0,0D0,
0827      &6*1D0,3710*0D0/
0828       DATA (KFDP(I,1),I=   1, 377)/21,22,23,4*-24,25,21,22,23,4*24,25,
0829      &21,22,23,4*-24,25,21,22,23,4*24,25,21,22,23,4*-24,25,21,22,23,
0830      &4*24,25,37,1000022,1000023,1000025,1000035,1000021,1000039,21,22,
0831      &23,4*-24,25,2*-37,21,22,23,4*24,25,2*37,22,23,-24,25,23,24,-12,
0832      &22,23,-24,25,23,24,-12,-14,48*16,22,23,-24,25,23,24,22,23,-24,25,
0833      &-37,23,24,37,1,2,3,4,5,6,7,8,21,1,2,3,4,5,6,7,8,11,13,15,17,1,2,
0834      &3,4,5,6,7,8,11,12,13,14,15,16,17,18,4*-1,4*-3,4*-5,4*-7,-11,-13,
0835      &-15,-17,1,2,3,4,5,6,7,8,11,13,15,17,21,2*22,23,24,1000022,
0836      &2*1000023,3*1000025,4*1000035,2*1000024,2*1000037,1000001,
0837      &2000001,1000001,-1000001,1000002,2000002,1000002,-1000002,
0838      &1000003,2000003,1000003,-1000003,1000004,2000004,1000004,
0839      &-1000004,1000005,2000005,1000005,-1000005,1000006,2000006,
0840      &1000006,-1000006,1000011,2000011,1000011,-1000011,1000012,
0841      &2000012,1000012,-1000012,1000013,2000013,1000013,-1000013,
0842      &1000014,2000014,1000014,-1000014,1000015,2000015,1000015,
0843      &-1000015,1000016,2000016,1000016,-1000016,1,2,3,4,5,6,7,8,11,12,
0844      &13,14,15,16,17,18,24,37,2*23,25,35,4*-1,4*-3,4*-5,4*-7,-11,-13,
0845      &-15,-17,3*24,1,2,3,4,5,6,7,8,11,13,15,17,21,2*22,23,24,23,25,24,
0846      &37,23,25,36,1000022,2*1000023,3*1000025,4*1000035,2*1000024,
0847      &2*1000037,1000001,2000001,1000001,-1000001,1000002,2000002/
0848       DATA (KFDP(I,1),I= 378, 580)/1000002,-1000002,1000003,2000003,
0849      &1000003,-1000003,1000004,2000004,1000004,-1000004,1000005,
0850      &2000005,1000005,-1000005,1000006,2000006,1000006,-1000006,
0851      &1000011,2000011,1000011,-1000011,1000012,2000012,1000012,
0852      &-1000012,1000013,2000013,1000013,-1000013,1000014,2000014,
0853      &1000014,-1000014,1000015,2000015,1000015,-1000015,1000016,
0854      &2000016,1000016,-1000016,1,2,3,4,5,6,7,8,11,13,15,17,21,2*22,23,
0855      &24,23,25,24,37,1000022,2*1000023,3*1000025,4*1000035,2*1000024,
0856      &2*1000037,1000001,2000001,1000001,-1000001,1000002,2000002,
0857      &1000002,-1000002,1000003,2000003,1000003,-1000003,1000004,
0858      &2000004,1000004,-1000004,1000005,2000005,1000005,-1000005,
0859      &1000006,2000006,1000006,-1000006,1000011,2000011,1000011,
0860      &-1000011,1000012,2000012,1000012,-1000012,1000013,2000013,
0861      &1000013,-1000013,1000014,2000014,1000014,-1000014,1000015,
0862      &2000015,1000015,-1000015,1000016,2000016,1000016,-1000016,-1,-3,
0863      &-5,-7,-11,-13,-15,-17,24,2*1000022,2*1000023,2*1000025,2*1000035,
0864      &1000006,2000006,1000006,2000006,-1000001,-1000003,-1000011,
0865      &-1000013,-1000015,-2000015,1,2,3,4,5,6,11,13,15,2,82,-11,-13,2*2,
0866      &-12,-14,-16,2*-2,2*-4,-2,-4,2*22,211,111,221,13,11,213,-213,221,
0867      &223,321,130,310,111,331,111,211,-12,12,-14,14,211,111,22,-13,-11/
0868       DATA (KFDP(I,1),I= 581, 992)/2*211,213,113,221,223,321,211,331,
0869      &22,111,211,2*22,211,22,111,211,22,211,221,111,11,211,111,2*211,
0870      &321,130,310,221,111,211,111,130,310,321,2*311,321,311,323,313,
0871      &323,313,321,3*311,-13,3*211,12,14,311,2*321,311,321,313,323,313,
0872      &323,311,4*321,211,111,3*22,111,321,130,-213,113,213,211,22,111,
0873      &11,13,211,321,130,310,221,211,111,11*-11,11*-13,-311,-313,-311,
0874      &-313,-20313,2*-311,-313,-311,-313,2*111,2*221,2*331,2*113,2*223,
0875      &2*333,-311,-313,2*-321,211,-311,-321,333,-311,-313,-321,211,
0876      &2*-321,2*-311,-321,211,113,421,2*411,421,411,423,413,423,413,421,
0877      &411,8*-11,8*-13,-321,-323,-321,-323,-311,2*-313,-311,-313,2*-311,
0878      &-321,-10323,-321,-323,-321,-311,2*-313,211,111,333,3*-321,-311,
0879      &-313,-321,-313,310,333,211,2*-321,-311,-313,-311,211,-321,3*-311,
0880      &211,113,321,2*421,411,421,413,423,413,423,411,421,-15,5*-11,
0881      &5*-13,221,331,333,221,331,333,10221,211,213,211,213,321,323,321,
0882      &323,2212,221,331,333,221,2*2,2*431,421,411,423,413,82,11,13,82,
0883      &443,82,6*12,6*14,2*16,3*-411,3*-413,2*-411,2*-413,2*441,2*443,
0884      &2*20443,2*2,2*4,2,4,511,521,511,523,513,523,513,521,511,6*12,
0885      &6*14,2*16,3*-421,3*-423,2*-421,2*-423,2*441,2*443,2*20443,2*2,
0886      &2*4,2,4,521,511,521,513,523,513,523,511,521,6*12,6*14,2*16,
0887      &3*-431,3*-433,2*-431,2*-433,3*441,3*443,3*20443,2*2,2*4,2,4,531/
0888       DATA (KFDP(I,1),I= 993,1402)/521,511,523,513,16,2*4,2*12,2*14,
0889      &2*16,4*2,4*4,2*-11,2*-13,2*-1,2*-3,2*-11,2*-13,2*-1,541,511,521,
0890      &513,523,21,11,13,15,1,2,3,4,21,22,553,21,2112,2212,2*2112,2212,
0891      &2112,2*2212,2112,-12,3122,3212,3112,2212,2*2112,-12,2*3122,3222,
0892      &3112,2212,2112,2212,3122,3222,3212,3122,3112,-12,-14,-12,3322,
0893      &3312,2*3122,3212,3322,3312,3122,3322,3312,-12,2*4122,7*-11,7*-13,
0894      &2*2224,2*2212,2*2214,2*3122,2*3212,2*3214,5*3222,4*3224,2*3322,
0895      &3324,2*2224,7*2212,5*2214,2*2112,2*2114,2*3122,2*3212,2*3214,
0896      &2*3222,2*3224,4*2,3,2*2,1,2*2,-11,-13,2*2,4*4122,-11,-13,2*2,
0897      &3*4132,3*4232,-11,-13,2*2,4332,-11,-13,2*2,-11,-13,2*2,-11,-13,
0898      &2*2,-11,-13,2*2,-11,-13,2*2,-11,-13,2*2,-11,-13,2*2,2*5122,-12,
0899      &-14,-16,5*4122,441,443,20443,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,
0900      &2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,4*5122,-12,-14,-16,2*-2,
0901      &2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,2*5132,2*5232,-12,-14,-16,
0902      &2*-2,2*-4,-2,-4,5332,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,
0903      &2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,
0904      &2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,
0905      &-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,
0906      &-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,
0907      &2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2/
0908       DATA (KFDP(I,1),I=1403,1713)/2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,
0909      &-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,
0910      &-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,221,223,221,
0911      &223,211,111,321,130,310,213,113,-213,321,311,321,311,323,313,
0912      &2*311,321,311,321,313,323,321,211,111,321,130,310,2*211,313,-313,
0913      &323,-323,421,411,423,413,411,421,413,423,411,421,423,413,443,
0914      &2*82,521,511,523,513,511,521,513,523,521,511,523,513,511,521,513,
0915      &523,553,2*21,213,-213,113,213,10211,10111,-10211,2*221,213,2*113,
0916      &-213,2*321,2*311,113,323,2*313,323,313,-313,323,-323,423,2*413,
0917      &2*423,413,443,82,523,2*513,2*523,2*513,523,553,21,11,13,82,4*443,
0918      &10441,20443,445,441,11,13,15,1,2,3,4,21,22,2*553,10551,20553,555,
0919      &1000039,-1000024,-1000037,1000022,1000023,1000025,1000035,
0920      &1000002,2000002,1000002,2000002,1000021,3*-12,3*-14,3*-16,12,11,
0921      &12,11,12,11,14,13,14,13,14,13,16,15,16,15,16,15,2*-2,2*-4,2*-6,
0922      &1000039,1000024,1000037,1000022,1000023,1000025,1000035,1000001,
0923      &2000001,1000001,2000001,1000021,3*-11,3*-13,3*-15,2*-1,-3,
0924      &1000039,-1000024,-1000037,1000022,1000023,1000025,1000035,
0925      &1000004,2000004,1000004,2000004,1000021,3*-12,3*-14,3*-16,12,11,
0926      &12,11,12,11,14,13,14,13,14,13,16,15,16,15,16,15,2*-2,2*-4,2*-6,
0927      &1000039,1000024,1000037,1000022,1000023,1000025,1000035,1000003/
0928       DATA (KFDP(I,1),I=1714,1984)/2000003,1000003,2000003,1000021,
0929      &3*-11,3*-13,3*-15,2*-1,-3,1000039,-1000024,-1000037,1000022,
0930      &1000023,1000025,1000035,1000006,2000006,1000006,2000006,1000021,
0931      &3*-12,3*-14,3*-16,12,11,12,11,12,11,14,13,14,13,14,13,16,15,16,
0932      &15,16,15,2*-2,2*-4,2*-6,1000039,1000024,1000037,1000022,1000023,
0933      &1000025,1000035,1000005,2000005,1000005,2000005,1000021,1000022,
0934      &1000016,-1000015,3*-11,3*-13,3*-15,2*-1,-3,1000039,-1000024,
0935      &-1000037,1000022,1000023,1000025,1000035,1000012,2000012,1000012,
0936      &2*12,2*14,2*16,3*-14,3*-16,3*-2,3*-4,3*-6,1000039,1000024,
0937      &1000037,1000022,1000023,1000025,1000035,1000011,2000011,1000011,
0938      &2000011,3*-13,3*-15,3*-1,3*-3,3*-5,1000039,-1000024,-1000037,
0939      &1000022,1000023,1000025,1000035,1000014,2000014,1000014,2000014,
0940      &2*12,2*14,2*16,3*-12,3*-16,3*-2,3*-4,3*-6,1000039,1000024,
0941      &1000037,1000022,1000023,1000025,1000035,1000013,2000013,1000013,
0942      &2000013,3*-11,3*-15,3*-1,3*-3,3*-5,1000039,-1000024,-1000037,
0943      &1000022,1000023,1000025,1000035,1000016,2000016,1000016,2000016,
0944      &2*12,2*14,2*16,3*-12,3*-14,3*-2,3*-4,3*-6,1000039,1000024,
0945      &1000037,1000022,1000023,1000025,1000035,1000015,2000015,1000015,
0946      &2000015,3*-11,3*-13,3*-1,3*-3,3*-5,1000039,1000001,-1000001,
0947      &2000001,-2000001,1000002,-1000002,2000002,-2000002,1000003/
0948       DATA (KFDP(I,1),I=1985,2321)/-1000003,2000003,-2000003,1000004,
0949      &-1000004,2000004,-2000004,1000005,-1000005,2000005,-2000005,
0950      &1000006,-1000006,2000006,-2000006,6*1000022,6*1000023,6*1000025,
0951      &6*1000035,1000024,-1000024,1000024,-1000024,1000024,-1000024,
0952      &1000037,-1000037,1000037,-1000037,1000037,-1000037,-12,12,-11,11,
0953      &-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,
0954      &-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-14,14,-13,13,
0955      &-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,
0956      &-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-16,16,-15,15,
0957      &-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,
0958      &-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-2,2,-2,2,-2,2,
0959      &-4,4,-4,4,-4,4,-6,6,-6,6,-6,6,5*1000039,4,1,-12,12,-12,12,-12,12,
0960      &-12,12,-12,12,-12,12,-14,14,-14,14,-14,14,-14,14,-14,14,-14,14,
0961      &-16,16,-16,16,-16,16,-16,16,-16,16,-16,16,-12,12,-11,11,-12,12,
0962      &-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,
0963      &-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-14,14,-13,13,-14,14,
0964      &-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,
0965      &-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-16,16,-15,15,-16,16,
0966      &-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,
0967      &-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-2,2,-2,2,-2,2,-4,4,-4/
0968       DATA (KFDP(I,1),I=2322,2573)/4,-4,4,-6,6,-6,6,-6,6,5*1000039,
0969      &16*1000022,1000024,-1000024,1000024,-1000024,1000024,-1000024,
0970      &1000024,-1000024,1000024,-1000024,1000024,-1000024,1000037,
0971      &-1000037,1000037,-1000037,1000037,-1000037,1000037,-1000037,
0972      &1000037,-1000037,1000037,-1000037,1000024,-1000024,1000037,
0973      &-1000037,1000001,-1000001,2000001,-2000001,1000002,-1000002,
0974      &2000002,-2000002,1000003,-1000003,2000003,-2000003,1000004,
0975      &-1000004,2000004,-2000004,1000005,-1000005,2000005,-2000005,
0976      &1000006,-1000006,2000006,-2000006,1000011,-1000011,2000011,
0977      &-2000011,1000012,-1000012,2000012,-2000012,1000013,-1000013,
0978      &2000013,-2000013,1000014,-1000014,2000014,-2000014,1000015,
0979      &-1000015,2000015,-2000015,1000016,-1000016,2000016,-2000016,
0980      &5*1000021,-12,12,-12,12,-12,12,-12,12,-12,12,-12,12,-14,14,-14,
0981      &14,-14,14,-14,14,-14,14,-14,14,-16,16,-16,16,-16,16,-16,16,-16,
0982      &16,-16,16,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,
0983      &11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,
0984      &12,-11,11,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,
0985      &13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,
0986      &14,-13,13,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,
0987      &15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16/
0988       DATA (KFDP(I,1),I=2574,2892)/16,-15,15,-2,2,-2,2,-2,2,-4,4,-4,4,
0989      &-4,4,-6,6,-6,6,-6,6,2*1000039,6*1000022,6*1000023,6*1000025,
0990      &6*1000035,1000022,1000023,1000025,1000035,1000002,2000002,
0991      &-1000001,-2000001,1000004,2000004,-1000003,-2000003,1000006,
0992      &2000006,-1000005,-2000005,1000012,2000012,-1000011,-2000011,
0993      &1000014,2000014,-1000013,-2000013,1000016,2000016,-1000015,
0994      &-2000015,2*1000021,-12,12,-11,-12,12,-11,-12,12,-11,-12,12,-11,
0995      &-12,12,-11,-12,12,-11,-14,-13,-14,-13,-14,-13,-14,14,-13,-14,14,
0996      &-13,-14,14,-13,-16,-15,-16,-15,-16,-15,-16,-15,-16,-15,-16,-15,
0997      &-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,
0998      &-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,-14,2*-13,14,
0999      &-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,
1000      &-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-16,2*-15,16,-16,2*-15,16,
1001      &-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,
1002      &-16,2*-15,16,-16,2*-15,16,2,-1,2,-1,2*2,-1,2,-1,3*2,-1,2*4,-3,
1003      &3*4,-3,2*6,5*1000039,16*1000022,16*1000023,1000024,-1000024,
1004      &1000024,-1000024,1000024,-1000024,1000024,-1000024,1000024,
1005      &-1000024,1000024,-1000024,1000037,-1000037,1000037,-1000037,
1006      &1000037,-1000037,1000037,-1000037,1000037,-1000037,1000037,
1007      &-1000037,1000024,-1000024,1000037,-1000037,1000001,-1000001/
1008       DATA (KFDP(I,1),I=2893,3182)/2000001,-2000001,1000002,-1000002,
1009      &2000002,-2000002,1000003,-1000003,2000003,-2000003,1000004,
1010      &-1000004,2000004,-2000004,1000005,-1000005,2000005,-2000005,
1011      &1000006,-1000006,2000006,-2000006,1000011,-1000011,2000011,
1012      &-2000011,1000012,-1000012,2000012,-2000012,1000013,-1000013,
1013      &2000013,-2000013,1000014,-1000014,2000014,-2000014,1000015,
1014      &-1000015,2000015,-2000015,1000016,-1000016,2000016,-2000016,
1015      &5*1000021,-12,12,-12,12,-12,12,-12,12,-12,12,-12,12,-14,14,-14,
1016      &14,-14,14,-14,14,-14,14,-14,14,-16,16,-16,16,-16,16,-16,16,-16,
1017      &16,-16,16,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,
1018      &11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,
1019      &12,-11,11,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,
1020      &13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,
1021      &14,-13,13,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,
1022      &15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,
1023      &16,-15,15,-2,2,-2,2,-2,2,-4,4,-4,4,-4,4,-6,6,-6,6,-6,6,5*1000039,
1024      &16*1000022,16*1000023,16*1000025,1000024,-1000024,1000024,
1025      &-1000024,1000024,-1000024,1000024,-1000024,1000024,-1000024,
1026      &1000024,-1000024,1000037,-1000037,1000037,-1000037,1000037,
1027      &-1000037,1000037,-1000037,1000037,-1000037,1000037,-1000037/
1028       DATA (KFDP(I,1),I=3183,3459)/1000024,-1000024,1000037,-1000037,
1029      &1000001,-1000001,2000001,-2000001,1000002,-1000002,2000002,
1030      &-2000002,1000003,-1000003,2000003,-2000003,1000004,-1000004,
1031      &2000004,-2000004,1000005,-1000005,2000005,-2000005,1000006,
1032      &-1000006,2000006,-2000006,1000011,-1000011,2000011,-2000011,
1033      &1000012,-1000012,2000012,-2000012,1000013,-1000013,2000013,
1034      &-2000013,1000014,-1000014,2000014,-2000014,1000015,-1000015,
1035      &2000015,-2000015,1000016,-1000016,2000016,-2000016,5*1000021,-12,
1036      &12,-12,12,-12,12,-12,12,-12,12,-12,12,-14,14,-14,14,-14,14,-14,
1037      &14,-14,14,-14,14,-16,16,-16,16,-16,16,-16,16,-16,16,-16,16,-12,
1038      &12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,
1039      &11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-14,
1040      &14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,
1041      &13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-16,
1042      &16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,
1043      &15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-2,2,
1044      &-2,2,-2,2,-4,4,-4,4,-4,4,-6,6,-6,6,-6,6,2*1000039,15*1000024,
1045      &6*1000022,6*1000023,6*1000025,6*1000035,1000022,1000023,1000025,
1046      &1000035,1000002,2000002,-1000001,-2000001,1000004,2000004,
1047      &-1000003,-2000003,1000006,2000006,-1000005,-2000005,1000012/
1048       DATA (KFDP(I,1),I=3460,3782)/2000012,-1000011,-2000011,1000014,
1049      &2000014,-1000013,-2000013,1000016,2000016,-1000015,-2000015,
1050      &2*1000021,-12,12,-11,-12,12,-11,-12,12,-11,-12,12,-11,-12,12,-11,
1051      &-12,12,-11,-14,14,-13,-14,14,-13,-14,14,-13,-14,14,-13,-14,14,
1052      &-13,-14,14,-13,-16,16,-15,-16,16,-15,-16,16,-15,-16,16,-15,-16,
1053      &16,-15,-16,16,-15,-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,-12,
1054      &2*-11,12,-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,-12,
1055      &2*-11,12,-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-14,
1056      &2*-13,14,-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-16,
1057      &2*-15,16,-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,-16,
1058      &2*-15,16,-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,2,-1,2,-1,2*2,-1,
1059      &2,-1,3*2,-1,2*4,-3,3*4,-3,2*6,1000039,-1000024,-1000037,1000022,
1060      &1000023,1000025,1000035,4*1000001,1000002,2000002,1000002,
1061      &2000002,1000021,3*-12,3*-14,3*-16,12,11,12,11,12,11,14,13,14,13,
1062      &14,13,16,15,16,15,16,15,2*-2,2*-4,2*-6,1000039,1000024,1000037,
1063      &1000022,1000023,1000025,1000035,4*1000002,1000001,2000001,
1064      &1000001,2000001,1000021,3*-11,3*-13,3*-15,2*-1,-3,1000039,
1065      &-1000024,-1000037,1000022,1000023,1000025,1000035,4*1000003,
1066      &1000004,2000004,1000004,2000004,1000021,3*-12,3*-14,3*-16,12,11,
1067      &12,11,12,11,14,13,14,13,14,13,16,15,16,15,16,15,2*-2,2*-4,2*-6/
1068       DATA (KFDP(I,1),I=3783,4127)/1000039,1000024,1000037,1000022,
1069      &1000023,1000025,1000035,4*1000004,1000003,2000003,1000003,
1070      &2000003,1000021,3*-11,3*-13,3*-15,2*-1,-3,1000039,-1000024,
1071      &-1000037,1000022,1000023,1000025,1000035,4*1000005,1000006,
1072      &2000006,1000006,2000006,1000021,3*-12,3*-14,3*-16,12,11,12,11,12,
1073      &11,14,13,14,13,14,13,16,15,16,15,16,15,2*-2,2*-4,2*-6,1000039,
1074      &1000024,1000037,1000022,1000023,1000025,1000035,4*1000006,
1075      &1000005,2000005,1000005,2000005,1000021,3*-11,3*-13,3*-15,2*-1,
1076      &-3,1000039,-1000024,-1000037,1000022,1000023,1000025,1000035,
1077      &4*1000011,1000012,2000012,1000012,2000012,2*12,2*14,2*16,3*-14,
1078      &3*-16,3*-2,3*-4,3*-6,1000039,-1000024,-1000037,1000022,1000023,
1079      &1000025,1000035,4*1000013,1000014,2000014,1000014,2000014,2*12,
1080      &2*14,2*16,3*-12,3*-16,3*-2,3*-4,3*-6,1000039,-1000024,-1000037,
1081      &1000022,1000023,1000025,1000035,4*1000015,1000016,2000016,
1082      &1000016,2000016,2*12,2*14,2*16,3*-12,3*-14,3*-2,3*-4,3*-6,3,4,5,
1083      &6,11,13,15,21,2*4,2,4,24,-11,-13,-15,3,4,5,6,11,13,15,21,5,6,21,
1084      &2*24,2*3000211,2*22,2*23,1,2,3,4,5,6,7,8,11,12,13,14,15,16,17,18,
1085      &2*24,3*3000211,24,4*-1,4*-3,4*-5,4*-7,-11,-13,-15,-17,22,23,22,
1086      &23,24,3000211,24,3000211,1,2,3,4,5,6,7,8,11,12,13,14,15,16,17,18,
1087      &1,2,3,4,5,6,1,2,3,4,5,6,21,1,2,3,4,5,6,21,1,2,3,4,5,6,21,1,2,3,4/
1088       DATA (KFDP(I,1),I=4128,8000)/5,6,1,2,3,4,5,6,1,2,3,4,5,6,21,
1089      &3100111,3200111,21,22,23,-24,21,22,23,24,22,23,-24,23,24,1,2,3,4,
1090      &5,6,7,8,11,12,13,14,15,16,17,18,21,22,23,24,9*11,9*-11,2*11,
1091      &2*-11,9*13,9*-13,2*13,2*-13,9*15,9*-15,2*15,2*-15,1,2,3,4,5,6,11,
1092      &12,9900012,13,14,9900014,15,16,9900016,3*-1,3*-3,3*-5,-11,-13,
1093      &-15,3*-11,2*-13,-15,24,3*-11,2*-13,-15,9900024,3*443,3*553,
1094      &3710*0/
1095       DATA (KFDP(I,2),I=   1, 339)/3*1,2,4,6,8,1,3*2,1,3,5,7,2,3*3,2,4,
1096      &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,6*1000006,3*7,
1097      &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,
1098      &13,11,13,-211,-213,-211,-213,-211,-213,-211,-213,2*-211,-321,
1099      &-323,-321,2*-323,3*-321,4*-211,-213,-211,-213,-211,-213,-211,
1100      &-213,-211,-213,3*-211,-213,4*-211,-323,-321,2*-211,2*-321,3*-211,
1101      &2*15,16,15,16,15,2*17,18,17,2*18,2*17,-1,-2,-3,-4,-5,-6,-7,-8,21,
1102      &-1,-2,-3,-4,-5,-6,-7,-8,-11,-13,-15,-17,-1,-2,-3,-4,-5,-6,-7,-8,
1103      &-11,-12,-13,-14,-15,-16,-17,-18,2,4,6,8,2,4,6,8,2,4,6,8,2,4,6,8,
1104      &12,14,16,18,-1,-2,-3,-4,-5,-6,-7,-8,-11,-13,-15,-17,21,22,2*23,
1105      &-24,2*1000022,1000023,1000022,1000023,1000025,1000022,1000023,
1106      &1000025,1000035,-1000024,-1000037,-1000024,-1000037,-1000001,
1107      &2*-2000001,2000001,-1000002,2*-2000002,2000002,-1000003,
1108      &2*-2000003,2000003,-1000004,2*-2000004,2000004,-1000005,
1109      &2*-2000005,2000005,-1000006,2*-2000006,2000006,-1000011,
1110      &2*-2000011,2000011,-1000012,2*-2000012,2000012,-1000013,
1111      &2*-2000013,2000013,-1000014,2*-2000014,2000014,-1000015,
1112      &2*-2000015,2000015,-1000016,2*-2000016,2000016,-1,-2,-3,-4,-5,-6,
1113      &-7,-8,-11,-12,-13,-14,-15,-16,-17,-18,-24,-37,22,25,2*36,2,4,6,8,
1114      &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/
1115       DATA (KFDP(I,2),I= 340, 533)/-7,-8,-11,-13,-15,-17,21,22,2*23,
1116      &-24,2*25,-37,-24,3*36,2*1000022,1000023,1000022,1000023,1000025,
1117      &1000022,1000023,1000025,1000035,-1000024,-1000037,-1000024,
1118      &-1000037,-1000001,2*-2000001,2000001,-1000002,2*-2000002,2000002,
1119      &-1000003,2*-2000003,2000003,-1000004,2*-2000004,2000004,-1000005,
1120      &2*-2000005,2000005,-1000006,2*-2000006,2000006,-1000011,
1121      &2*-2000011,2000011,-1000012,2*-2000012,2000012,-1000013,
1122      &2*-2000013,2000013,-1000014,2*-2000014,2000014,-1000015,
1123      &2*-2000015,2000015,-1000016,2*-2000016,2000016,-1,-2,-3,-4,-5,-6,
1124      &-7,-8,-11,-13,-15,-17,21,22,2*23,-24,2*25,-37,-24,2*1000022,
1125      &1000023,1000022,1000023,1000025,1000022,1000023,1000025,1000035,
1126      &-1000024,-1000037,-1000024,-1000037,-1000001,2*-2000001,2000001,
1127      &-1000002,2*-2000002,2000002,-1000003,2*-2000003,2000003,-1000004,
1128      &2*-2000004,2000004,-1000005,2*-2000005,2000005,-1000006,
1129      &2*-2000006,2000006,-1000011,2*-2000011,2000011,-1000012,
1130      &2*-2000012,2000012,-1000013,2*-2000013,2000013,-1000014,
1131      &2*-2000014,2000014,-1000015,2*-2000015,2000015,-1000016,
1132      &2*-2000016,2000016,2,4,6,8,12,14,16,18,25,1000024,1000037,
1133      &1000024,1000037,1000024,1000037,1000024,1000037,2*-1000005,
1134      &2*-2000005,1000002,1000004,1000012,1000014,2*1000016,-3,-4,-5,-6/
1135       DATA (KFDP(I,2),I= 534, 938)/-7,-8,-13,-15,-17,11,-82,12,14,-1,
1136      &-3,11,13,15,1,4,3,4,1,3,22,11,-211,2*22,-13,-11,-211,211,111,211,
1137      &-321,130,310,22,2*111,-211,11,-11,13,-13,-211,111,22,14,12,111,
1138      &22,111,3*211,-311,22,211,22,111,-211,211,11,-211,13,22,-211,111,
1139      &-211,22,111,-11,-211,111,2*-211,-321,130,310,221,111,-211,111,
1140      &2*0,-211,111,22,-211,111,-211,111,-211,211,-213,113,223,221,14,
1141      &111,211,111,-11,-13,211,111,22,211,111,211,111,2*211,213,113,223,
1142      &221,22,-211,111,113,223,22,111,-321,310,211,111,2*-211,221,22,
1143      &-11,-13,-211,-321,130,310,221,-211,111,11*12,11*14,2*211,2*213,
1144      &211,20213,2*321,2*323,211,213,211,213,211,213,211,213,211,213,
1145      &211,213,3*211,213,211,2*321,8*211,2*113,3*211,111,22,211,111,211,
1146      &111,4*211,8*12,8*14,2*211,2*213,2*111,221,2*113,223,333,20213,
1147      &211,2*321,323,2*311,313,-211,111,113,2*211,321,2*211,311,321,310,
1148      &211,-211,4*211,321,4*211,113,2*211,-321,111,22,-211,111,-211,111,
1149      &-211,211,-211,211,16,5*12,5*14,3*211,3*213,211,2*111,2*113,
1150      &2*-311,2*-313,-2112,3*321,323,2*-1,22,111,321,311,321,311,-82,
1151      &-11,-13,-82,22,-82,6*-11,6*-13,2*-15,211,213,20213,211,213,20213,
1152      &431,433,431,433,311,313,311,313,311,313,-1,-4,-3,-4,-1,-3,22,
1153      &-211,111,-211,111,-211,211,-211,211,6*-11,6*-13,2*-15,211,213,
1154      &20213,211,213,20213,431,433,431,433,321,323,321,323,321,323,-1/
1155       DATA (KFDP(I,2),I= 939,1352)/-4,-3,-4,-1,-3,22,211,111,211,111,
1156      &4*211,6*-11,6*-13,2*-15,211,213,20213,211,213,20213,431,433,431,
1157      &433,221,331,333,221,331,333,221,331,333,-1,-4,-3,-4,-1,-3,22,
1158      &-321,-311,-321,-311,-15,-3,-1,2*-11,2*-13,2*-15,-1,-4,-3,-4,-3,
1159      &-4,-1,-4,2*12,2*14,2,3,2,3,2*12,2*14,2,1,22,411,421,411,421,21,
1160      &-11,-13,-15,-1,-2,-3,-4,2*21,22,21,2*-211,111,22,111,211,22,211,
1161      &-211,11,2*-211,111,-211,111,22,11,22,111,-211,211,111,211,22,211,
1162      &111,211,-211,22,11,13,11,-211,2*111,2*22,111,211,-321,-211,111,
1163      &11,2*-211,7*12,7*14,-321,-323,-311,-313,-311,-313,211,213,211,
1164      &213,211,213,111,221,331,113,223,111,221,113,223,321,323,321,-211,
1165      &-213,111,221,331,113,223,333,10221,111,221,331,113,223,211,213,
1166      &211,213,321,323,321,323,321,323,311,313,311,313,2*-1,-3,-1,2203,
1167      &3201,3203,2203,2101,2103,12,14,-1,-3,2*111,2*211,12,14,-1,-3,22,
1168      &111,2*22,111,22,12,14,-1,-3,22,12,14,-1,-3,12,14,-1,-3,12,14,-1,
1169      &-3,12,14,-1,-3,12,14,-1,-3,12,14,-1,-3,12,14,-1,-3,2*-211,11,13,
1170      &15,-211,-213,-20213,-431,-433,3*3122,1,4,3,4,1,3,11,13,15,1,4,3,
1171      &4,1,3,11,13,15,1,4,3,4,1,3,2*111,2*211,11,13,15,1,4,3,4,1,3,11,
1172      &13,15,1,4,3,4,1,3,4*22,11,13,15,1,4,3,4,1,3,22,11,13,15,1,4,3,4,
1173      &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,
1174      &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,3/
1175       DATA (KFDP(I,2),I=1353,1815)/11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,
1176      &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,
1177      &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,
1178      &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,3,
1179      &2*111,2*211,-211,111,-321,130,310,-211,111,211,-211,111,-213,113,
1180      &-211,111,223,211,111,213,113,211,111,223,-211,111,-321,130,310,
1181      &2*-211,-311,311,-321,321,211,111,211,111,-211,111,-211,111,311,
1182      &2*321,311,22,2*-82,-211,111,-211,111,211,111,211,111,-321,-311,
1183      &-321,-311,411,421,411,421,22,2*21,-211,2*211,111,-211,111,2*211,
1184      &111,-211,211,111,211,-321,2*-311,-321,22,-211,111,211,111,-311,
1185      &311,-321,321,211,111,-211,111,321,311,22,-82,-211,111,211,111,
1186      &-321,-311,411,421,22,21,-11,-13,-82,211,111,221,111,4*22,-11,-13,
1187      &-15,-1,-2,-3,-4,2*21,211,111,3*22,1,2*2,4*1,2*-24,2*-37,2*1,3,5,
1188      &1,3,5,1,3,5,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6,-3,-5,-3,-5,-3,
1189      &-5,2,2*1,4*2,2*24,2*37,2,1,3,5,1,3,5,1,3,5,-3,2*-5,3,2*4,4*3,
1190      &2*-24,2*-37,3,1,3,5,1,3,5,1,3,5,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,
1191      &5,6,-1,-5,-1,-5,-1,-5,4,2*3,4*4,2*24,2*37,4,1,3,5,1,3,5,1,3,5,-3,
1192      &2*-5,5,2*6,4*5,2*-24,2*-37,5,1,3,5,1,3,5,1,3,5,1,2,3,4,5,6,1,2,3,
1193      &4,5,6,1,2,3,4,5,6,-1,-3,-1,-3,-1,-3,6,2*5,4*6,2*24,2*37,6,4,-15,
1194      &16,1,3,5,1,3,5,1,3,5,-3,2*-5,11,2*12,4*11,2*-24,-37,13,15,11,15/
1195       DATA (KFDP(I,2),I=1816,2317)/11,13,11,13,15,11,13,15,1,3,5,1,3,5,
1196      &1,3,5,12,2*11,4*12,2*24,2*37,11,13,15,11,13,15,1,3,5,1,3,5,1,3,5,
1197      &13,2*14,4*13,2*-24,2*-37,13,15,11,15,11,13,11,13,15,11,13,15,1,3,
1198      &5,1,3,5,1,3,5,14,2*13,4*14,2*24,2*37,11,13,15,11,13,15,1,3,5,1,3,
1199      &5,1,3,5,15,2*16,4*15,2*-24,2*-37,13,15,11,15,11,13,11,13,15,11,
1200      &13,15,1,3,5,1,3,5,1,3,5,16,2*15,4*16,2*24,2*37,11,13,15,11,13,15,
1201      &1,3,5,1,3,5,1,3,5,21,-1,1,-1,1,-2,2,-2,2,-3,3,-3,3,-4,4,-4,4,-5,
1202      &5,-5,5,-6,6,-6,6,1,3,5,2,4,6,1,3,5,2,4,6,1,3,5,2,4,6,1,3,5,2,4,6,
1203      &1,-1,3,-3,5,-5,1,-1,3,-3,5,-5,-1,1,-2,2,-1,1,-2,2,-1,1,-2,2,-3,3,
1204      &-4,4,-3,3,-4,4,-3,3,-4,4,-5,5,-6,6,-5,5,-6,6,-5,5,-6,6,-1,1,-2,2,
1205      &-1,1,-2,2,-1,1,-2,2,-3,3,-4,4,-3,3,-4,4,-3,3,-4,4,-5,5,-6,6,-5,5,
1206      &-6,6,-5,5,-6,6,-1,1,-2,2,-1,1,-2,2,-1,1,-2,2,-3,3,-4,4,-3,3,-4,4,
1207      &-3,3,-4,4,-5,5,-6,6,-5,5,-6,6,-5,5,-6,6,-1,1,-1,1,-3,3,-1,1,-1,1,
1208      &-3,3,-1,1,-1,1,-3,3,22,23,25,35,36,-1,-3,-13,13,-13,13,-13,13,
1209      &-15,15,-15,15,-15,15,-11,11,-11,11,-11,11,-15,15,-15,15,-15,15,
1210      &-11,11,-11,11,-11,11,-13,13,-13,13,-13,13,-1,1,-2,2,-1,1,-2,2,-1,
1211      &1,-2,2,-3,3,-4,4,-3,3,-4,4,-3,3,-4,4,-5,5,-6,6,-5,5,-6,6,-5,5,-6,
1212      &6,-1,1,-2,2,-1,1,-2,2,-1,1,-2,2,-3,3,-4,4,-3,3,-4,4,-3,3,-4,4,-5,
1213      &5,-6,6,-5,5,-6,6,-5,5,-6,6,-1,1,-2,2,-1,1,-2,2,-1,1,-2,2,-3,3,-4,
1214      &4,-3,3,-4,4,-3,3,-4,4,-5,5,-6,6,-5,5,-6,6,-5,5,-6,6,-1,1,-1,1,-3/
1215       DATA (KFDP(I,2),I=2318,2770)/3,-1,1,-1,1,-3,3,-1,1,-1,1,-3,3,22,
1216      &23,25,35,36,22,23,11,13,15,12,14,16,1,3,5,2,4,25,35,36,-24,24,11,
1217      &-11,13,-13,15,-15,1,-1,3,-3,-24,24,11,-11,13,-13,15,-15,1,-1,3,
1218      &-3,-37,37,-37,37,-1,1,-1,1,-2,2,-2,2,-3,3,-3,3,-4,4,-4,4,-5,5,-5,
1219      &5,-6,6,-6,6,-11,11,-11,11,-12,12,-12,12,-13,13,-13,13,-14,14,-14,
1220      &14,-15,15,-15,15,-16,16,-16,16,1,3,5,2,4,-13,13,-13,13,-13,13,
1221      &-15,15,-15,15,-15,15,-11,11,-11,11,-11,11,-15,15,-15,15,-15,15,
1222      &-11,11,-11,11,-11,11,-13,13,-13,13,-13,13,-1,1,-2,2,-1,1,-2,2,-1,
1223      &1,-2,2,-3,3,-4,4,-3,3,-4,4,-3,3,-4,4,-5,5,-6,6,-5,5,-6,6,-5,5,-6,
1224      &6,-1,1,-2,2,-1,1,-2,2,-1,1,-2,2,-3,3,-4,4,-3,3,-4,4,-3,3,-4,4,-5,
1225      &5,-6,6,-5,5,-6,6,-5,5,-6,6,-1,1,-2,2,-1,1,-2,2,-1,1,-2,2,-3,3,-4,
1226      &4,-3,3,-4,4,-3,3,-4,4,-5,5,-6,6,-5,5,-6,6,-5,5,-6,6,-1,1,-1,1,-3,
1227      &3,-1,1,-1,1,-3,3,-1,1,-1,1,-3,3,24,37,24,-11,-13,-15,-1,-3,24,
1228      &-11,-13,-15,-1,-3,24,-11,-13,-15,-1,-3,24,-11,-13,-15,-1,-3,4*37,
1229      &2*-1,2*2,2*-3,2*4,2*-5,2*6,2*-11,2*12,2*-13,2*14,2*-15,2*16,-1,
1230      &-3,-13,14,2*-13,14,2*-13,14,-13,-15,16,2*-15,16,2*-15,16,-15,
1231      &6*-11,-15,16,2*-15,16,2*-15,16,-15,6*-11,6*-13,-1,-2,-1,2,-1,-2,
1232      &-1,2,-1,-2,-1,2,-3,-4,-3,4,-3,-4,-3,4,-3,-4,-3,4,-5,-6,-5,6,-5,
1233      &-6,-5,6,-5,-6,-5,6,-1,-2,-1,2,-1,-2,-1,2,-1,-2,-1,2,-3,-4,-3,4,
1234      &-3,-4,-3,4,-3,-4,-3,4,-5,-6,-5,6,-5,-6,-5,6,-5,-6,-5,6,-1,-2,-1/
1235       DATA (KFDP(I,2),I=2771,3221)/2,-1,-2,-1,2,-1,-2,-1,2,-3,-4,-3,4,
1236      &-3,-4,-3,4,-3,-4,-3,4,-5,-6,-5,6,-5,-6,-5,6,-5,-6,-5,6,2,-1,2,-1,
1237      &2*4,-3,4,-3,3*6,-5,2*4,-3,3*6,-5,2*6,22,23,25,35,36,22,23,11,13,
1238      &15,12,14,16,1,3,5,2,4,25,35,36,22,23,11,13,15,12,14,16,1,3,5,2,4,
1239      &25,35,36,-24,24,11,-11,13,-13,15,-15,1,-1,3,-3,-24,24,11,-11,13,
1240      &-13,15,-15,1,-1,3,-3,-37,37,-37,37,-1,1,-1,1,-2,2,-2,2,-3,3,-3,3,
1241      &-4,4,-4,4,-5,5,-5,5,-6,6,-6,6,-11,11,-11,11,-12,12,-12,12,-13,13,
1242      &-13,13,-14,14,-14,14,-15,15,-15,15,-16,16,-16,16,1,3,5,2,4,-13,
1243      &13,-13,13,-13,13,-15,15,-15,15,-15,15,-11,11,-11,11,-11,11,-15,
1244      &15,-15,15,-15,15,-11,11,-11,11,-11,11,-13,13,-13,13,-13,13,-1,1,
1245      &-2,2,-1,1,-2,2,-1,1,-2,2,-3,3,-4,4,-3,3,-4,4,-3,3,-4,4,-5,5,-6,6,
1246      &-5,5,-6,6,-5,5,-6,6,-1,1,-2,2,-1,1,-2,2,-1,1,-2,2,-3,3,-4,4,-3,3,
1247      &-4,4,-3,3,-4,4,-5,5,-6,6,-5,5,-6,6,-5,5,-6,6,-1,1,-2,2,-1,1,-2,2,
1248      &-1,1,-2,2,-3,3,-4,4,-3,3,-4,4,-3,3,-4,4,-5,5,-6,6,-5,5,-6,6,-5,5,
1249      &-6,6,-1,1,-1,1,-3,3,-1,1,-1,1,-3,3,-1,1,-1,1,-3,3,22,23,25,35,36,
1250      &22,23,11,13,15,12,14,16,1,3,5,2,4,25,35,36,22,23,11,13,15,12,14,
1251      &16,1,3,5,2,4,25,35,36,22,23,11,13,15,12,14,16,1,3,5,2,4,25,35,36,
1252      &-24,24,11,-11,13,-13,15,-15,1,-1,3,-3,-24,24,11,-11,13,-13,15,
1253      &-15,1,-1,3,-3,-37,37,-37,37,-1,1,-1,1,-2,2,-2,2,-3,3,-3,3,-4,4,
1254      &-4,4,-5,5,-5,5,-6,6,-6,6,-11,11,-11,11,-12,12,-12,12,-13,13,-13/
1255       DATA (KFDP(I,2),I=3222,3669)/13,-14,14,-14,14,-15,15,-15,15,-16,
1256      &16,-16,16,1,3,5,2,4,-13,13,-13,13,-13,13,-15,15,-15,15,-15,15,
1257      &-11,11,-11,11,-11,11,-15,15,-15,15,-15,15,-11,11,-11,11,-11,11,
1258      &-13,13,-13,13,-13,13,-1,1,-2,2,-1,1,-2,2,-1,1,-2,2,-3,3,-4,4,-3,
1259      &3,-4,4,-3,3,-4,4,-5,5,-6,6,-5,5,-6,6,-5,5,-6,6,-1,1,-2,2,-1,1,-2,
1260      &2,-1,1,-2,2,-3,3,-4,4,-3,3,-4,4,-3,3,-4,4,-5,5,-6,6,-5,5,-6,6,-5,
1261      &5,-6,6,-1,1,-2,2,-1,1,-2,2,-1,1,-2,2,-3,3,-4,4,-3,3,-4,4,-3,3,-4,
1262      &4,-5,5,-6,6,-5,5,-6,6,-5,5,-6,6,-1,1,-1,1,-3,3,-1,1,-1,1,-3,3,-1,
1263      &1,-1,1,-3,3,24,37,23,11,13,15,12,14,16,1,3,5,2,4,25,35,36,24,-11,
1264      &-13,-15,-1,-3,24,-11,-13,-15,-1,-3,24,-11,-13,-15,-1,-3,24,-11,
1265      &-13,-15,-1,-3,4*37,2*-1,2*2,2*-3,2*4,2*-5,2*6,2*-11,2*12,2*-13,
1266      &2*14,2*-15,2*16,-1,-3,-13,14,2*-13,14,2*-13,14,-13,-15,16,2*-15,
1267      &16,2*-15,16,-15,-11,12,2*-11,12,2*-11,12,-11,-15,16,2*-15,16,
1268      &2*-15,16,-15,-11,12,2*-11,12,2*-11,12,-11,-13,14,2*-13,14,2*-13,
1269      &14,-13,-1,-2,-1,2,-1,-2,-1,2,-1,-2,-1,2,-3,-4,-3,4,-3,-4,-3,4,-3,
1270      &-4,-3,4,-5,-6,-5,6,-5,-6,-5,6,-5,-6,-5,6,-1,-2,-1,2,-1,-2,-1,2,
1271      &-1,-2,-1,2,-3,-4,-3,4,-3,-4,-3,4,-3,-4,-3,4,-5,-6,-5,6,-5,-6,-5,
1272      &6,-5,-6,-5,6,-1,-2,-1,2,-1,-2,-1,2,-1,-2,-1,2,-3,-4,-3,4,-3,-4,
1273      &-3,4,-3,-4,-3,4,-5,-6,-5,6,-5,-6,-5,6,-5,-6,-5,6,2,-1,2,-1,2*4,
1274      &-3,4,-3,3*6,-5,2*4,-3,3*6,-5,2*6,1,2*2,4*1,23,25,35,36,2*-24/
1275       DATA (KFDP(I,2),I=3670,4136)/2*-37,2*1,3,5,1,3,5,1,3,5,1,2,3,4,5,
1276      &6,1,2,3,4,5,6,1,2,3,4,5,6,-3,-5,-3,-5,-3,-5,2,2*1,4*2,23,25,35,
1277      &36,2*24,2*37,2,1,3,5,1,3,5,1,3,5,-3,2*-5,3,2*4,4*3,23,25,35,36,
1278      &2*-24,2*-37,3,1,3,5,1,3,5,1,3,5,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,
1279      &5,6,-1,-5,-1,-5,-1,-5,4,2*3,4*4,23,25,35,36,2*24,2*37,4,1,3,5,1,
1280      &3,5,1,3,5,-3,2*-5,5,2*6,4*5,23,25,35,36,2*-24,2*-37,5,1,3,5,1,3,
1281      &5,1,3,5,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6,-1,-3,-1,-3,-1,-3,6,
1282      &2*5,4*6,23,25,35,36,2*24,2*37,6,1,3,5,1,3,5,1,3,5,-3,2*-5,11,
1283      &2*12,4*11,23,25,35,36,2*-24,2*-37,13,15,11,15,11,13,11,13,15,11,
1284      &13,15,1,3,5,1,3,5,1,3,5,13,2*14,4*13,23,25,35,36,2*-24,2*-37,13,
1285      &15,11,15,11,13,11,13,15,11,13,15,1,3,5,1,3,5,1,3,5,15,2*16,4*15,
1286      &23,25,35,36,2*-24,2*-37,13,15,11,15,11,13,11,13,15,11,13,15,1,3,
1287      &5,1,3,5,1,3,5,-3,-4,-5,-6,-11,-13,-15,21,-1,-3,2*-5,5,12,14,16,
1288      &-3,-4,-5,-6,-11,-13,-15,21,-5,-6,21,-24,-3000211,-24,-3000211,
1289      &3000111,3000221,3000111,3000221,-1,-2,-3,-4,-5,-6,-7,-8,-11,-12,
1290      &-13,-14,-15,-16,-17,-18,23,3000111,23,3000111,22,3000221,2,4,6,8,
1291      &2,4,6,8,2,4,6,8,2,4,6,8,12,14,16,18,2*3000111,2*3000221,-3000211,
1292      &2*-24,-3000211,-1,-2,-3,-4,-5,-6,-7,-8,-11,-12,-13,-14,-15,-16,
1293      &-17,-18,-1,-2,-3,-4,-5,-6,-1,-2,-3,-4,-5,-6,21,-1,-2,-3,-4,-5,-6,
1294      &21,-1,-2,-3,-4,-5,-6,21,-1,-2,-3,-4,-5,-6,-1,-2,-3,-4,-5,-6,-1/
1295       DATA (KFDP(I,2),I=4137,8000)/-2,-3,-4,-5,-6,3*21,3*1,4*2,1,2*11,
1296      &2*12,11,-1,-2,-3,-4,-5,-6,-7,-8,-11,-12,-13,-14,-15,-16,-17,-18,
1297      &21,22,23,-24,3*-1,3*-3,3*-5,3*1,3*3,3*5,2*-13,2*15,3*-1,3*-3,
1298      &3*-5,3*1,3*3,3*5,2*-11,2*15,3*-1,3*-3,3*-5,3*1,3*3,3*5,2*-11,
1299      &2*13,-1,-2,-3,-4,-5,-6,-11,-12,9900012,-13,-14,9900014,-15,-16,
1300      &9900016,2,4,6,2,4,6,2,4,6,9900012,9900014,9900016,-11,-13,-15,
1301      &-13,2*-15,24,-11,-13,-15,-13,2*-15,9900024,6*21,3710*0/
1302       DATA (KFDP(I,3),I=   1,1021)/81*0,14,6*0,2*16,2*0,6*111,310,130,
1303      &2*0,3*111,310,130,321,113,211,223,221,2*113,2*211,2*223,2*221,
1304      &2*113,221,2*113,2*213,-213,113,2*111,310,130,310,130,2*310,130,
1305      &402*0,4*3,4*4,1,4,3,2*2,0,-11,8*0,-211,5*0,2*111,211,-211,211,
1306      &-211,10*0,111,4*0,2*111,-211,-11,11,-13,22,111,3*0,22,3*0,111,
1307      &211,4*0,111,11*0,111,-211,6*0,-211,3*111,7*0,111,-211,5*0,2*221,
1308      &3*0,111,5*0,111,11*0,-311,-313,-311,-321,-313,-323,111,221,331,
1309      &113,223,-311,-313,-311,-321,-313,-323,111,221,331,113,223,22*0,
1310      &111,113,2*211,-211,-311,211,111,3*211,-211,7*211,7*0,111,-211,
1311      &111,-211,-321,-323,-311,-321,-313,-323,-211,-213,-321,-323,-311,
1312      &-321,-313,-323,-211,-213,22*0,111,113,-311,2*-211,211,-211,310,
1313      &-211,2*111,211,2*-211,-321,-211,2*211,-211,111,-211,2*211,6*0,
1314      &111,-211,111,-211,0,221,331,333,321,311,221,331,333,321,311,20*0,
1315      &3,13*0,-411,-413,-10413,-10411,-20413,-415,-411,-413,-10413,
1316      &-10411,-20413,-415,-411,-413,16*0,-4,-1,-4,-3,2*-2,5*0,111,-211,
1317      &111,-211,-421,-423,-10423,-10421,-20423,-425,-421,-423,-10423,
1318      &-10421,-20423,-425,-421,-423,16*0,-4,-1,-4,-3,2*-2,5*0,111,-211,
1319      &111,-211,-431,-433,-10433,-10431,-20433,-435,-431,-433,-10433,
1320      &-10431,-20433,-435,-431,-433,19*0,-4,-1,-4,-3,2*-2,8*0,441,443,
1321      &441,443,441,443,-4,-1,-4,-3,-4,-3,-4,-1,531,533,531,533,3,2,3,2/
1322       DATA (KFDP(I,3),I=1022,2223)/511,513,511,513,1,2,13*0,2*21,11*0,
1323      &2112,6*0,2212,12*0,2*3122,3212,10*0,3322,2*0,3122,3212,3214,2112,
1324      &2114,2212,2112,3122,3212,3214,2112,2114,2212,2112,52*0,3*3,1,6*0,
1325      &4*3,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,
1326      &2*2,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,
1327      &2*2,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,
1328      &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,
1329      &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,
1330      &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,
1331      &3,2*2,31*0,211,111,45*0,-211,2*111,-211,3*111,-211,111,211,30*0,
1332      &-211,111,13*0,2*21,-211,111,199*0,2*5,210*0,-1,-3,-5,-2,-4,-6,-1,
1333      &-3,-5,-2,-4,-6,-1,-3,-5,-2,-4,-6,-1,-3,-5,-2,-4,-6,-2,2,-4,4,-6,
1334      &6,-2,2,-4,4,-6,6,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,
1335      &-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,
1336      &-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,
1337      &-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,
1338      &-1,1,-1,3,-3,3,-3,5,-5,5,-5,-3,3,-5,5,-5,5,-3,3,-5,5,-5,5,-3,3,
1339      &-5,5,-5,5,5*0,11,12,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,
1340      &-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,
1341      &-11,13,-13,15,-15,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3/
1342       DATA (KFDP(I,3),I=2224,2783)/-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,
1343      &-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,
1344      &-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,
1345      &-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,-3,3,
1346      &-5,5,-5,5,-3,3,-5,5,-5,5,-3,3,-5,5,-5,5,7*0,-11,-13,-15,-12,-14,
1347      &-16,-1,-3,-5,-2,-4,5*0,-12,12,-14,14,-16,16,-2,2,-4,4,2*0,-12,12,
1348      &-14,14,-16,16,-2,2,-4,4,52*0,-1,-3,-5,-2,-4,11,-11,13,-13,15,-15,
1349      &11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,
1350      &11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,1,-1,1,-1,3,-3,3,-3,5,
1351      &-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,
1352      &-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,
1353      &-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,
1354      &-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,-3,3,-5,5,
1355      &-5,5,-3,3,-5,5,-5,5,-3,3,-5,5,-5,5,3*0,12,14,16,2,4,0,12,14,16,2,
1356      &4,0,12,14,16,2,4,0,12,14,16,2,4,28*0,2,4,12,-11,11,14,-13,13,16,
1357      &-15,15,12,-11,11,14,-13,13,16,-15,15,12,11,14,13,16,15,12,-11,11,
1358      &14,-13,13,16,-15,15,12,11,14,13,16,15,12,11,14,13,16,15,2*2,1,-1,
1359      &2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3,
1360      &2*6,5,-5,2*2,1,-1,2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3,2*6,5,-5,
1361      &2*2,1,-1,2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3,2*6,5,-5,2*2,1,-1/
1362       DATA (KFDP(I,3),I=2784,3354)/2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3,
1363      &2*6,5,-5,3,-3,5,-5,1,3,-3,5,-5,1,3,5,-5,1,5,-5,1,3,5,-5,1,3,7*0,
1364      &-11,-13,-15,-12,-14,-16,-1,-3,-5,-2,-4,5*0,-11,-13,-15,-12,-14,
1365      &-16,-1,-3,-5,-2,-4,5*0,-12,12,-14,14,-16,16,-2,2,-4,4,2*0,-12,12,
1366      &-14,14,-16,16,-2,2,-4,4,52*0,-1,-3,-5,-2,-4,11,-11,13,-13,15,-15,
1367      &11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,
1368      &11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,1,-1,1,-1,3,-3,3,-3,5,
1369      &-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,
1370      &-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,
1371      &-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,
1372      &-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,-3,3,-5,5,
1373      &-5,5,-3,3,-5,5,-5,5,-3,3,-5,5,-5,5,7*0,-11,-13,-15,-12,-14,-16,
1374      &-1,-3,-5,-2,-4,5*0,-11,-13,-15,-12,-14,-16,-1,-3,-5,-2,-4,5*0,
1375      &-11,-13,-15,-12,-14,-16,-1,-3,-5,-2,-4,5*0,-12,12,-14,14,-16,16,
1376      &-2,2,-4,4,2*0,-12,12,-14,14,-16,16,-2,2,-4,4,52*0,-1,-3,-5,-2,-4,
1377      &11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,
1378      &11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,1,
1379      &-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,
1380      &-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,
1381      &-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3/
1382       DATA (KFDP(I,3),I=3355,8000)/-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,
1383      &-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,-3,3,-5,5,-5,5,-3,3,-5,5,
1384      &-5,5,-3,3,-5,5,-5,5,3*0,-11,-13,-15,-12,-14,-16,-1,-3,-5,-2,-4,
1385      &4*0,12,14,16,2,4,0,12,14,16,2,4,0,12,14,16,2,4,0,12,14,16,2,4,
1386      &28*0,2,4,12,-11,11,14,-13,13,16,-15,15,12,-11,11,14,-13,13,16,
1387      &-15,15,12,-11,11,14,-13,13,16,-15,15,12,-11,11,14,-13,13,16,-15,
1388      &15,12,-11,11,14,-13,13,16,-15,15,12,-11,11,14,-13,13,16,-15,15,
1389      &2*2,1,-1,2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3,2*6,5,-5,2*2,1,-1,
1390      &2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3,
1391      &2*6,5,-5,2*2,1,-1,2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3,2*6,5,-5,
1392      &2*2,1,-1,2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3,2*6,5,-5,3,-3,5,-5,
1393      &1,3,-3,5,-5,1,3,5,-5,1,5,-5,1,3,5,-5,1,3,351*0,-5,169*0,2,4,6,2,
1394      &4,6,2,4,6,-2,-4,-6,-2,-4,-6,-2,-4,-6,2*9900014,2*9900016,2,4,6,2,
1395      &4,6,2,4,6,-2,-4,-6,-2,-4,-6,-2,-4,-6,2*9900012,2*9900016,2,4,6,2,
1396      &4,6,2,4,6,-2,-4,-6,-2,-4,-6,-2,-4,-6,2*9900012,2*9900014,3757*0/
1397       DATA (KFDP(I,4),I=   1,8000)/94*0,4*111,6*0,111,2*0,-211,0,-211,
1398      &3*0,111,2*-211,0,111,0,2*111,113,221,2*111,-213,-211,211,113,
1399      &6*111,310,2*130,402*0,13*81,41*0,-11,10*0,111,-211,4*0,111,62*0,
1400      &111,211,111,211,7*0,111,211,111,211,35*0,2*-211,2*111,211,111,
1401      &-211,2*211,2*-211,13*0,-211,111,-211,111,4*0,-211,111,-211,111,
1402      &34*0,111,-211,3*111,3*-211,2*111,3*-211,14*0,-321,-311,3*0,-321,
1403      &-311,20*0,-3,43*0,6*1,39*0,6*2,42*0,6*3,14*0,8*4,4*0,4*-5,4*0,
1404      &2*-5,67*0,-211,111,5*0,-211,111,52*0,2101,2103,2*2101,6*0,4*81,
1405      &4*0,4*81,6*0,4*81,0,28*81,13*0,6*2101,18*81,4*0,18*81,4*0,9*81,0,
1406      &162*81,31*0,-211,111,6516*0/
1407       DATA (KFDP(I,5),I=   1,8000)/96*0,2*111,17*0,111,7*0,2*111,0,
1408      &3*111,0,111,597*0,-211,2*111,-211,111,-211,111,65*0,111,-211,
1409      &3*111,-211,111,7193*0/
1410
1411 C...PYDAT4, with particle names (character strings).
1412
1413       DATA (CHAF(I,1),I=   1, 202)/'d','u','s','c','b','t','b''','t''',
1414      &2*' ','e-','nu_e','mu-','nu_mu','tau-','nu_tau','tau''-',
1415      &'nu''_tau',2*' ','g','gamma','Z0','W+','h0',6*' ','Z''0','Z"0',
1416      &'W''+','H0','A0','H+',' ','Graviton',' ','R0','LQ_ue',38*' ',
1417      &'specflav','rndmflav','phasespa','c-hadron','b-hadron',2*' ',
1418      &'junction',' ','system','cluster','string','indep.','CMshower',
1419      &'SPHEaxis','THRUaxis','CLUSjet','CELLjet','table',' ','reggeon',
1420      &'pi0','rho0','a_20','K_L0','pi+','rho+','a_2+','eta','omega',
1421      &'f_2','K_S0','K0','K*0','K*_20','K+','K*+','K*_2+','eta''','phi',
1422      &'f''_2','D+','D*+','D*_2+','D0','D*0','D*_20','D_s+','D*_s+',
1423      &'D*_2s+','eta_c','J/psi','chi_2c','B0','B*0','B*_20','B+','B*+',
1424      &'B*_2+','B_s0','B*_s0','B*_2s0','B_c+','B*_c+','B*_2c+','eta_b',
1425      &'Upsilon','chi_2b','pomeron','dd_1','Delta-','ud_0','ud_1','n0',
1426      &'Delta0','uu_1','p+','Delta+','Delta++','sd_0','sd_1','Sigma-',
1427      &'Sigma*-','Lambda0','su_0','su_1','Sigma0','Sigma*0','Sigma+',
1428      &'Sigma*+','ss_1','Xi-','Xi*-','Xi0','Xi*0','Omega-','cd_0',
1429      &'cd_1','Sigma_c0','Sigma*_c0','Lambda_c+','Xi_c0','cu_0','cu_1',
1430      &'Sigma_c+','Sigma*_c+','Sigma_c++','Sigma*_c++','Xi_c+','cs_0',
1431      &'cs_1','Xi''_c0','Xi*_c0','Xi''_c+','Xi*_c+','Omega_c0',
1432      &'Omega*_c0','cc_1','Xi_cc+','Xi*_cc+','Xi_cc++','Xi*_cc++'/
1433       DATA (CHAF(I,1),I= 203, 332)/'Omega_cc+','Omega*_cc+',
1434      &'Omega*_ccc++','bd_0','bd_1','Sigma_b-','Sigma*_b-','Lambda_b0',
1435      &'Xi_b-','Xi_bc0','bu_0','bu_1','Sigma_b0','Sigma*_b0','Sigma_b+',
1436      &'Sigma*_b+','Xi_b0','Xi_bc+','bs_0','bs_1','Xi''_b-','Xi*_b-',
1437      &'Xi''_b0','Xi*_b0','Omega_b-','Omega*_b-','Omega_bc0','bc_0',
1438      &'bc_1','Xi''_bc0','Xi*_bc0','Xi''_bc+','Xi*_bc+','Omega''_bc0',
1439      &'Omega*_bc0','Omega_bcc+','Omega*_bcc+','bb_1','Xi_bb-',
1440      &'Xi*_bb-','Xi_bb0','Xi*_bb0','Omega_bb-','Omega*_bb-',
1441      &'Omega_bbc0','Omega*_bbc0','Omega*_bbb-','a_00','b_10','a_0+',
1442      &'b_1+','f_0','h_1','K*_00','K_10','K*_0+','K_1+','f''_0','h''_1',
1443      &'D*_0+','D_1+','D*_00','D_10','D*_0s+','D_1s+','chi_0c','h_1c',
1444      &'B*_00','B_10','B*_0+','B_1+','B*_0s0','B_1s0','B*_0c+','B_1c+',
1445      &'chi_0b','h_1b','a_10','a_1+','f_1','K*_10','K*_1+','f''_1',
1446      &'D*_1+','D*_10','D*_1s+','chi_1c','B*_10','B*_1+','B*_1s0',
1447      &'B*_1c+','chi_1b','psi''','Upsilon''','~d_L','~u_L','~s_L',
1448      &'~c_L','~b_1','~t_1','~e_L-','~nu_eL','~mu_L-','~nu_muL',
1449      &'~tau_1-','~nu_tauL','~g','~chi_10','~chi_20','~chi_1+',
1450      &'~chi_30','~chi_40','~chi_2+','~Gravitino','~d_R','~u_R','~s_R',
1451      &'~c_R','~b_2','~t_2','~e_R-','~nu_eR','~mu_R-','~nu_muR',
1452      &'~tau_2-','~nu_tauR','pi_tc0','pi_tc+','pi''_tc0','eta_tc0'/
1453       DATA (CHAF(I,1),I= 333, 500)/'rho_tc0','rho_tc+','omega_tc',
1454      &'V8_tc','pi_22_1_tc','pi_22_8_tc','rho_11_tc','rho_12_tc',
1455      &'rho_21_tc','rho_22_tc','d*','u*','e*-','nu*_e0','Graviton*',
1456      &'nu_Re','nu_Rmu','nu_Rtau','Z_R0','W_R+','H_L++','H_R++',
1457      &'rho_diff0','pi_diffr+','omega_di','phi_diff','J/psi_di',
1458      &'n_diffr0','p_diffr+','cc~[3S18]','cc~[1S08]','cc~[3P08]',
1459      &'bb~[3S18]','bb~[1S08]','bb~[3P08]',133*' '/
1460       DATA (CHAF(I,2),I=   1, 205)/'dbar','ubar','sbar','cbar','bbar',
1461      &'tbar','b''bar','t''bar',2*' ','e+','nu_ebar','mu+','nu_mubar',
1462      &'tau+','nu_taubar','tau''+','nu''_taubar',5*' ','W-',9*' ',
1463      &'W''-',2*' ','H-',3*' ','Rbar0','LQ_uebar',39*' ','rndmflavbar',
1464      &' ','c-hadronbar','b-hadronbar',20*' ','pi-','rho-','a_2-',4*' ',
1465      &'Kbar0','K*bar0','K*_2bar0','K-','K*-','K*_2-',3*' ','D-','D*-',
1466      &'D*_2-','Dbar0','D*bar0','D*_2bar0','D_s-','D*_s-','D*_2s-',
1467      &3*' ','Bbar0','B*bar0','B*_2bar0','B-','B*-','B*_2-','B_sbar0',
1468      &'B*_sbar0','B*_2sbar0','B_c-','B*_c-','B*_2c-',4*' ','dd_1bar',
1469      &'Deltabar+','ud_0bar','ud_1bar','nbar0','Deltabar0','uu_1bar',
1470      &'pbar-','Deltabar-','Deltabar--','sd_0bar','sd_1bar','Sigmabar+',
1471      &'Sigma*bar+','Lambdabar0','su_0bar','su_1bar','Sigmabar0',
1472      &'Sigma*bar0','Sigmabar-','Sigma*bar-','ss_1bar','Xibar+',
1473      &'Xi*bar+','Xibar0','Xi*bar0','Omegabar+','cd_0bar','cd_1bar',
1474      &'Sigma_cbar0','Sigma*_cbar0','Lambda_cbar-','Xi_cbar0','cu_0bar',
1475      &'cu_1bar','Sigma_cbar-','Sigma*_cbar-','Sigma_cbar--',
1476      &'Sigma*_cbar--','Xi_cbar-','cs_0bar','cs_1bar','Xi''_cbar0',
1477      &'Xi*_cbar0','Xi''_cbar-','Xi*_cbar-','Omega_cbar0',
1478      &'Omega*_cbar0','cc_1bar','Xi_ccbar-','Xi*_ccbar-','Xi_ccbar--',
1479      &'Xi*_ccbar--','Omega_ccbar-','Omega*_ccbar-','Omega*_cccbar-'/
1480       DATA (CHAF(I,2),I= 206, 325)/'bd_0bar','bd_1bar','Sigma_bbar+',
1481      &'Sigma*_bbar+','Lambda_bbar0','Xi_bbar+','Xi_bcbar0','bu_0bar',
1482      &'bu_1bar','Sigma_bbar0','Sigma*_bbar0','Sigma_bbar-',
1483      &'Sigma*_bbar-','Xi_bbar0','Xi_bcbar-','bs_0bar','bs_1bar',
1484      &'Xi''_bbar+','Xi*_bbar+','Xi''_bbar0','Xi*_bbar0','Omega_bbar+',
1485      &'Omega*_bbar+','Omega_bcbar0','bc_0bar','bc_1bar','Xi''_bcbar0',
1486      &'Xi*_bcbar0','Xi''_bcbar-','Xi*_bcbar-','Omega''_bcba',
1487      &'Omega*_bcbar0','Omega_bccbar-','Omega*_bccbar-','bb_1bar',
1488      &'Xi_bbbar+','Xi*_bbbar+','Xi_bbbar0','Xi*_bbbar0','Omega_bbbar+',
1489      &'Omega*_bbbar+','Omega_bbcbar0','Omega*_bbcbar0',
1490      &'Omega*_bbbbar+',2*' ','a_0-','b_1-',2*' ','K*_0bar0','K_1bar0',
1491      &'K*_0-','K_1-',2*' ','D*_0-','D_1-','D*_0bar0','D_1bar0',
1492      &'D*_0s-','D_1s-',2*' ','B*_0bar0','B_1bar0','B*_0-','B_1-',
1493      &'B*_0sbar0','B_1sbar0','B*_0c-','B_1c-',3*' ','a_1-',' ',
1494      &'K*_1bar0','K*_1-',' ','D*_1-','D*_1bar0','D*_1s-',' ',
1495      &'B*_1bar0','B*_1-','B*_1sbar0','B*_1c-',3*' ','~d_Lbar',
1496      &'~u_Lbar','~s_Lbar','~c_Lbar','~b_1bar','~t_1bar','~e_L+',
1497      &'~nu_eLbar','~mu_L+','~nu_muLbar','~tau_1+','~nu_tauLbar',3*' ',
1498      &'~chi_1-',2*' ','~chi_2-',' ','~d_Rbar','~u_Rbar','~s_Rbar',
1499      &'~c_Rbar','~b_2bar','~t_2bar','~e_R+','~nu_eRbar','~mu_R+'/
1500       DATA (CHAF(I,2),I= 326, 500)/'~nu_muRbar','~tau_2+',
1501      &'~nu_tauRbar',' ','pi_tc-',3*' ','rho_tc-',8*' ','d*bar','u*bar',
1502      &'e*bar+','nu*_ebar0',5*' ','W_R-','H_L--','H_R--',' ',
1503      &'pi_diffr-',3*' ','n_diffrbar0','p_diffrbar-',139*' '/
1504
1505 C...PYDATR, with initial values for the random number generator.
1506       DATA MRPY/19780503,0,0,97,33,0/
1507
1508 C...Default values for allowed processes and kinematics constraints.
1509       DATA MSEL/1/
1510       DATA MSUB/500*0/
1511       DATA ((KFIN(I,J),J=-40,40),I=1,2)/16*0,4*1,4*0,6*1,5*0,5*1,0,
1512      &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,
1513      &6*1,4*0,4*1,16*0/
1514       DATA CKIN/
1515      &  2.0D0, -1.0D0,  0.0D0, -1.0D0,  1.0D0,
1516      &  1.0D0,  -10D0,   10D0,  -40D0,   40D0,
1517      1  -40D0,   40D0,  -40D0,   40D0,  -40D0,
1518      1   40D0, -1.0D0,  1.0D0, -1.0D0,  1.0D0,
1519      2  0.0D0,  1.0D0,  0.0D0,  1.0D0, -1.0D0,
1520      2  1.0D0, -1.0D0,  1.0D0,    0D0,    0D0,
1521      3  2.0D0, -1.0D0,    0D0,    0D0,  0.0D0,
1522      3 -1.0D0,  0.0D0, -1.0D0,  4.0D0, -1.0D0,
1523      4 12.0D0, -1.0D0, 12.0D0, -1.0D0, 12.0D0,
1524      4 -1.0D0, 12.0D0, -1.0D0,    0D0,    0D0,
1525      5  0.0D0, -1.0D0,  0.0D0, -1.0D0,  0.0D0,
1526      5 -1.0D0,    0D0,    0D0,    0D0,    0D0,
1527      6 0.0001D0, 0.99D0, 0.0001D0, 0.99D0,    0D0,
1528      6   -1D0,    0D0,   -1D0,    0D0,   -1D0,
1529      7    0D0,   -1D0, 0.0001D0, 0.99D0, 0.0001D0,
1530      7 0.99D0,    2D0,   -1D0,    0D0,    0D0,
1531      8  120*0D0/
1532
1533 C...Default values for main switches and parameters. Reset information.
1534       DATA (MSTP(I),I=1,100)/
1535      &  3,    1,    2,    0,    0,    0,    0,    0,    0,    0,
1536      1  1,    0,    1,   30,    0,    1,    4,    3,    4,    3,
1537      2  1,    0,    1,    0,    0,    0,    0,    0,    0,    1,
1538      3  1,    8,    0,    1,    0,    2,    1,    5,    2,    0,
1539      4  2,    1,    3,    7,    3,    1,    1,    0,    1,    0,
1540      5  7,    1,    3,    1,    5,    1,    1,    5,    1,    7,
1541      6  2,    3,    2,    2,    1,    5,    2,    3,    0,    0,
1542      7  1,    1,    0,    0,    0,    0,    0,    0,    0,    0,
1543      8  1,    4,  100,    1,    1,    2,    4,    1,    1,    0,
1544      9  1,    3,    1,    3,    1,    0,    0,    0,    0,    0/
1545       DATA (MSTP(I),I=101,200)/
1546      &  3,    1,    0,    0,    0,    0,    0,    0,    0,    0,
1547      1  1,    1,    1,    0,    0,    0,    0,    0,    0,    0,
1548      2  0,    1,    2,    1,    1,  100,    0,    0,   10,    0,
1549      3  0,    4,    0,    1,    0,    0,    0,    0,    0,    0,
1550      4  0,    0,    0,    0,    0,    1,    0,    0,    0,    0,
1551      5  0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1552      6  0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1553      7  0,    2,    0,    0,    0,    0,    0,    0,    0,    0,
1554      8  6,  413, 2007,   09,   12,    0,    0,    0,    0,    0,
1555      9  0,    0,    0,    0,    0,    0,    0,    0,    0,    0/
1556       DATA (PARP(I),I=1,100)/
1557      &  0.25D0,  10D0, 8*0D0,
1558      1  0D0, 0D0, 1.0D0, 0.01D0, 0.5D0, 1.0D0, 1.0D0, 0.4D0, 2*0D0,
1559      2  10*0D0,
1560      3  1.5D0,2.0D0,0.075D0,1.0D0,0.2D0,0D0,1.0D0,0.70D0,0.006D0,0D0,
1561      4  0.02D0,2.0D0,0.10D0,1000D0,2054D0,123D0,246D0,50D0,0D0,0.054D0,
1562      5  10*0D0,
1563      6  0.25D0, 1.0D0,0.25D0, 1.0D0, 2.0D0,1D-3, 4.0D0,1D-3,2*0D0,
1564      7  4.0D0, 0.25D0, 5*0D0, 0.025D0, 2.0D0, 0.1D0,
1565      8  1.90D0, 2.0D0, 0.5D0, 0.4D0, 0.90D0,
1566      8  0.95D0, 0.7D0, 0.5D0, 1800D0, 0.16D0,
1567      9  2.0D0,0.40D0,5.0D0,1.0D0,0.0D0,3.0D0,1.0D0,0.75D0,1.0D0,5.0D0/
1568       DATA (PARP(I),I=101,200)/
1569      &  0.5D0, 0.28D0,  1.0D0, 0.8D0, 0D0, 0D0, 0D0, 0D0, 0D0, 1D0,
1570      1  2.0D0, 3*0D0, 1.5D0, 0.5D0, 0.6D0, 2.5D0, 2.0D0, 1.0D0,
1571      2  1.0D0,  0.4D0, 8*0D0,
1572      3  0.01D0, 9*0D0,
1573      4  1.16D0, 0.0119D0, 0.01D0, 0.01D0, 0.05D0,
1574      4  9.28D0, 0.15D0, 0.02D0, 0.48D0, 0.09D0,
1575      5  0D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0,
1576      6  2.20D0, 23.6D0, 18.4D0, 11.5D0, 0.5D0, 0D0, 0D0, 0D0, 2*0D0,
1577      7  0D0,   0D0,   0D0,  1.0D0, 6*0D0,
1578      8  0.1D0, 0.01D0, 0.01D0, 0.01D0, 0.1D0, 0.01D0, 0.01D0, 0.01D0,
1579      8  0.3D0, 0.64D0,
1580      9  0.64D0, 5.0D0, 1.0D4, 1.0D4, 6*0D0/
1581       DATA MSTI/200*0/
1582       DATA PARI/200*0D0/
1583       DATA MINT/400*0/
1584       DATA VINT/400*0D0/
1585
1586 C...Constants for the generation of the various processes.
1587       DATA (ISET(I),I=1,100)/
1588      &  1,    1,    1,   -1,    3,   -1,   -1,    3,   -2,    2,
1589      1  2,    2,    2,    2,    2,    2,   -1,    2,    2,    2,
1590      2 -1,    2,    2,    2,    2,    2,   -1,    2,    2,    2,
1591      3  2,    2,    2,    2,    2,    2,   -1,   -1,   -1,   -1,
1592      4 -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,
1593      5 -1,   -1,    2,    2,   -1,   -1,   -1,    2,   -1,   -1,
1594      6 -1,   -1,   -1,   -1,   -1,   -1,   -1,    2,    2,    2,
1595      7  4,    4,    4,   -1,   -1,    4,    4,   -1,   -1,    2,
1596      8  2,    2,    2,    2,    2,    2,    2,    2,    2,   -2,
1597      9  0,    0,    0,    0,    0,    9,   -2,   -2,    8,   -2/
1598       DATA (ISET(I),I=101,200)/
1599      & -1,    1,    1,    1,    1,    2,    2,    2,   -2,    2,
1600      1  2,    2,    2,    2,    2,   -1,   -1,   -1,   -2,   -2,
1601      2  5,    5,    5,    5,   -2,   -2,   -2,   -2,   -2,   -2,
1602      3  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1603      4  1,    1,    1,    1,    1,    1,    1,    1,    1,   -2,
1604      5  1,    1,    1,   -2,   -2,    1,    1,    1,   -2,   -2,
1605      6  2,    2,    2,    2,    2,    2,    2,    2,    2,   -2,
1606      7  2,    2,    5,    5,   -2,    2,    2,    5,    5,   -2,
1607      8  5,    5,    2,    2,    2,    5,    5,    2,    2,    2,
1608      9  1,    1,    1,    2,    2,   -2,   -2,   -2,   -2,   -2/
1609       DATA (ISET(I),I=201,300)/
1610      &  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1611      1  2,    2,    2,    2,   -2,    2,    2,    2,    2,    2,
1612      2  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1613      3  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1614      4  2,    2,    2,    2,   -1,    2,    2,    2,    2,    2,
1615      5  2,    2,    2,    2,   -1,    2,   -1,    2,    2,   -2,
1616      6  2,    2,    2,    2,    2,   -1,   -1,   -1,   -1,   -1,
1617      7  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1618      8  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1619      9  2,    2,    2,    2,    2,    2,    2,    2,    2,    2/
1620       DATA (ISET(I),I=301,500)/
1621      &  2,   39*-2,
1622      4  1,    1,    2,    2,    2,    2,    2,    2,    2,    2,
1623      5  5,    5,    1,    1,   -1,   -1,   -1,   -1,   -1,   -1,
1624      6  2,    2,    2,    2,    2,    2,    2,    2,   -1,    2,
1625      7  2,    2,    2,    2,    2,    2,    2,   -1,   -1,   -1,
1626      8  2,    2,    2,    2,    2,    2,    2,    2,   -2,   -2,
1627      9  1,    1,    2,    2,    2, 5*-2,
1628      &  5,    5, 18*-2,
1629      2  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1630      3  2,    2,    2,    2,    2,    2,    2,    2,    2, 21*-2,
1631      6  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1632      7  2,    2,    2,    2,    2,    2,    2,    2,    2, 21*-2/
1633       DATA ((KFPR(I,J),J=1,2),I=1,50)/
1634      &  23,    0,   24,    0,   25,    0,   24,    0,   25,    0,
1635      &  24,    0,   23,    0,   25,    0,    0,    0,    0,    0,
1636      1   0,    0,    0,    0,   21,   21,   21,   22,   21,   23,
1637      1  21,   24,   21,   25,   22,   22,   22,   23,   22,   24,
1638      2  22,   25,   23,   23,   23,   24,   23,   25,   24,   24,
1639      2  24,   25,   25,   25,    0,   21,    0,   22,    0,   23,
1640      3   0,   24,    0,   25,    0,   21,    0,   22,    0,   23,
1641      3   0,   24,    0,   25,    0,   21,    0,   22,    0,   23,
1642      4   0,   24,    0,   25,    0,   21,    0,   22,    0,   23,
1643      4   0,   24,    0,   25,    0,   21,    0,   22,    0,   23/
1644       DATA ((KFPR(I,J),J=1,2),I=51,100)/
1645      5   0,   24,    0,   25,    0,    0,    0,    0,    0,    0,
1646      5   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1647      6   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1648      6   0,    0,    0,    0,   21,   21,   24,   24,   23,   24,
1649      7  23,   23,   24,   24,   23,   24,   23,   25,   22,   22,
1650      7  23,   23,   24,   24,   24,   25,   25,   25,    0,  211,
1651      8   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1652      8 443,   21,10441,   21,20443,   21,  445,   21,    0,    0,
1653      9   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1654      9   0,    0,    0,    0,    0,    0,    0,    0,    0,    0/
1655       DATA ((KFPR(I,J),J=1,2),I=101,150)/
1656      &  23,    0,   25,    0,   25,    0,10441,    0,  445,    0,
1657      & 443,   22,  443,   21,  443,   22,    0,    0,   22,   25,
1658      1  21,   25,    0,   25,   21,   25,   22,   22,   21,   22,
1659      1  22,   23,   23,   23,   24,   24,    0,    0,    0,    0,
1660      2  25,    6,   25,    6,   25,    0,   25,    0,    0,    0,
1661      2   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1662      3   0,   21,    0,   21,    0,   22,    0,   22,    0,    0,
1663      3   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1664      4  32,    0,   34,    0,   37,    0,   41,    0,   42,    0,
1665      4 4000011, 0, 4000001, 0, 4000002, 0, 3000331, 0,   0,    0/
1666       DATA ((KFPR(I,J),J=1,2),I=151,200)/
1667      5  35,    0,   35,    0,   35,    0,    0,    0,    0,    0,
1668      5  36,    0,   36,    0,   36,    0,    0,    0,    0,    0,
1669      6   6,   37,   42,    0,   42,   42,   42,   42,   11,    0,
1670      6  11,    0, 0, 4000001, 0, 4000002, 0, 4000011,    0,    0,
1671      7  23,   35,   24,   35,   35,    0,   35,    0,    0,    0,
1672      7  23,   36,   24,   36,   36,    0,   36,    0,    0,    0,
1673      8  35,    6,   35,    6,   21,   35,    0,   35,   21,   35,
1674      8  36,    6,   36,    6,   21,   36,    0,   36,   21,   36,
1675      9  3000113, 0, 3000213, 0, 3000223, 0, 11,    0,   11,    0,
1676      9   0,    0,    0,    0,    0,    0,    0,    0,    0,    0/
1677       DATA ((KFPR(I,J),J=1,2),I=201,240)/
1678      &  1000011,   1000011,   2000011,   2000011,   1000011,
1679      &  2000011,   1000013,   1000013,   2000013,   2000013,
1680      &  1000013,   2000013,   1000015,   1000015,   2000015,
1681      &  2000015,   1000015,   2000015,   1000011,   1000012,
1682      1  1000015,   1000016,   2000015,   1000016,   1000012,
1683      1  1000012,   1000016,   1000016,         0,         0,
1684      1  1000022,   1000022,   1000023,   1000023,   1000025,
1685      1  1000025,   1000035,   1000035,   1000022,   1000023,
1686      2  1000022,   1000025,   1000022,   1000035,   1000023,
1687      2  1000025,   1000023,   1000035,   1000025,   1000035,
1688      2  1000024,   1000024,   1000037,   1000037,   1000024,
1689      2  1000037,   1000022,   1000024,   1000023,   1000024,
1690      3  1000025,   1000024,   1000035,   1000024,   1000022,
1691      3  1000037,   1000023,   1000037,   1000025,   1000037,
1692      3  1000035,   1000037,   1000021,   1000022,   1000021,
1693      3  1000023,   1000021,   1000025,   1000021,   1000035/
1694       DATA ((KFPR(I,J),J=1,2),I=241,280)/
1695      4  1000021,   1000024,   1000021,   1000037,   1000021,
1696      4  1000021,   1000021,   1000021,         0,         0,
1697      4  1000002,   1000022,   2000002,   1000022,   1000002,
1698      4  1000023,   2000002,   1000023,   1000002,   1000025,
1699      5  2000002,   1000025,   1000002,   1000035,   2000002,
1700      5  1000035,   1000001,   1000024,   2000005,   1000024,
1701      5  1000001,   1000037,   2000005,   1000037,   1000002,
1702      5  1000021,   2000002,   1000021,         0,         0,
1703      6  1000006,   1000006,   2000006,   2000006,   1000006,
1704      6  2000006,   1000006,   1000006,   2000006,   2000006,
1705      6        0,         0,         0,         0,         0,
1706      6        0,         0,         0,         0,         0,
1707      7  1000002,   1000002,   2000002,   2000002,   1000002,
1708      7  2000002,   1000002,   1000002,   2000002,   2000002,
1709      7  1000002,   2000002,   1000002,   1000002,   2000002,
1710      7  2000002,   1000002,   1000002,   2000002,   2000002/
1711       DATA ((KFPR(I,J),J=1,2),I=281,350)/
1712      8  1000005,   1000002,   2000005,   2000002,   1000005,
1713      8  2000002,   1000005,   1000002,   2000005,   2000002,
1714      8  1000005,   2000002,   1000005,   1000005,   2000005,
1715      8  2000005,   1000005,   1000005,   2000005,   2000005,
1716      9  1000005,   1000005,   2000005,   2000005,   1000005,
1717      9  2000005,   1000005,   1000021,   2000005,   1000021,
1718      9  1000005,   2000005,        37,        25,        37,
1719      9       35,        36,        25,        36,        35,
1720      &       37,        37,      78*0,
1721      4  9900041,         0,   9900042,         0,   9900041,
1722      4       11,   9900042,        11,   9900041,        13,
1723      4  9900042,        13,   9900041,        15,   9900042,
1724      4       15,   9900041,   9900041,   9900042,   9900042/
1725       DATA ((KFPR(I,J),J=1,2),I=351,400)/
1726      5  9900041,         0,   9900042,         0,   9900023,
1727      5        0,   9900024,         0,         0,         0,
1728      5        0,         0,         0,         0,         0,
1729      5        0,         0,         0,         0,         0,
1730      6       24,        24,        24,   3000211,   3000211,
1731      6  3000211,        22,   3000111,        22,   3000221,
1732      6       23,   3000111,        23,   3000221,        24,
1733      6  3000211,         0,         0,        24,        23,
1734      7       24,   3000111,   3000211,        23,   3000211,
1735      7  3000111,        22,   3000211,        23,   3000211,
1736      7       24,   3000111,        24,   3000221,         0,
1737      7        0,         0,         0,         0,         0,
1738      8   0,    0,    0,    0,   21,   21,    0,   21,    0,    0,
1739      8  21,   21,    0,    0,    0,    0,    0,    0,    0,    0,
1740      9  5000039,         0,   5000039,         0,        21,
1741      9  5000039,         0,   5000039,        21,   5000039,
1742      9     10*0/
1743       DATA ((KFPR(I,J),J=1,2),I=401,500)/
1744      &  37,    6,   37,    6,    36*0,
1745      2      443,        21,   9900443,        21,   9900441,
1746      2       21,   9910441,        21,         0,   9900443,
1747      2        0,   9900441,         0,   9910441,        21,
1748      2  9900443,        21,   9900441,        21,   9910441,
1749      3 10441, 21, 20443,  21,  445,   21,    0, 10441,   0, 20443,
1750      3   0,  445,   21, 10441,  21, 20443,  21,  445,  42*0,
1751      6      553,        21,   9900553,        21,   9900551,
1752      6       21,   9910551,        21,         0,   9900553,
1753      6        0,   9900551,         0,   9910551,        21,
1754      6  9900553,        21,   9900551,        21,   9910551,
1755      7 10551, 21, 20553,  21,  555,   21,    0, 10551,   0, 20553,
1756      7   0,  555,   21, 10551,  21, 20553,  21,  555, 42*0/
1757       DATA COEF/10000*0D0/
1758       DATA (((ICOL(I,J,K),K=1,2),J=1,4),I=1,40)/
1759      &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,
1760      &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,
1761      &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,
1762      &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,
1763      &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,
1764      &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,
1765      &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,
1766      &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,
1767      &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,
1768      &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/
1769
1770 C...Treatment of resonances.
1771       DATA (MWID(I)  ,I=   1, 500)/5*0,3*1,8*0,1,5*0,3*1,6*0,1,0,4*1,
1772      &3*0,2*1,254*0,19*2,0,7*2,0,2,0,2,0,26*1,7*0,6*2,133*0/
1773
1774 C...Character constants: name of processes.
1775       DATA PROC(0)/                    'All included subprocesses   '/
1776       DATA (PROC(I),I=1,20)/
1777      &'f + fbar -> gamma*/Z0       ',  'f + fbar'' -> W+/-           ',
1778      &'f + fbar -> h0              ',  'gamma + W+/- -> W+/-        ',
1779      &'Z0 + Z0 -> h0               ',  'Z0 + W+/- -> W+/-           ',
1780      &'                            ',  'W+ + W- -> h0               ',
1781      &'                            ',  'f + f'' -> f + f'' (QFD)      ',
1782      1'f + f'' -> f + f'' (QCD)      ','f + fbar -> f'' + fbar''      ',
1783      1'f + fbar -> g + g           ',  'f + fbar -> g + gamma       ',
1784      1'f + fbar -> g + Z0          ',  'f + fbar'' -> g + W+/-       ',
1785      1'f + fbar -> g + h0          ',  'f + fbar -> gamma + gamma   ',
1786      1'f + fbar -> gamma + Z0      ',  'f + fbar'' -> gamma + W+/-   '/
1787       DATA (PROC(I),I=21,40)/
1788      2'f + fbar -> gamma + h0      ',  'f + fbar -> Z0 + Z0         ',
1789      2'f + fbar'' -> Z0 + W+/-      ', 'f + fbar -> Z0 + h0         ',
1790      2'f + fbar -> W+ + W-         ',  'f + fbar'' -> W+/- + h0      ',
1791      2'f + fbar -> h0 + h0         ',  'f + g -> f + g              ',
1792      2'f + g -> f + gamma          ',  'f + g -> f + Z0             ',
1793      3'f + g -> f'' + W+/-          ', 'f + g -> f + h0             ',
1794      3'f + gamma -> f + g          ',  'f + gamma -> f + gamma      ',
1795      3'f + gamma -> f + Z0         ',  'f + gamma -> f'' + W+/-      ',
1796      3'f + gamma -> f + h0         ',  'f + Z0 -> f + g             ',
1797      3'f + Z0 -> f + gamma         ',  'f + Z0 -> f + Z0            '/
1798       DATA (PROC(I),I=41,60)/
1799      4'f + Z0 -> f'' + W+/-         ', 'f + Z0 -> f + h0            ',
1800      4'f + W+/- -> f'' + g          ', 'f + W+/- -> f'' + gamma      ',
1801      4'f + W+/- -> f'' + Z0         ', 'f + W+/- -> f'' + W+/-       ',
1802      4'f + W+/- -> f'' + h0         ', 'f + h0 -> f + g             ',
1803      4'f + h0 -> f + gamma         ',  'f + h0 -> f + Z0            ',
1804      5'f + h0 -> f'' + W+/-         ', 'f + h0 -> f + h0            ',
1805      5'g + g -> f + fbar           ',  'g + gamma -> f + fbar       ',
1806      5'g + Z0 -> f + fbar          ',  'g + W+/- -> f + fbar''       ',
1807      5'g + h0 -> f + fbar          ',  'gamma + gamma -> f + fbar   ',
1808      5'gamma + Z0 -> f + fbar      ',  'gamma + W+/- -> f + fbar''   '/
1809       DATA (PROC(I),I=61,80)/
1810      6'gamma + h0 -> f + fbar      ',  'Z0 + Z0 -> f + fbar         ',
1811      6'Z0 + W+/- -> f + fbar''      ', 'Z0 + h0 -> f + fbar         ',
1812      6'W+ + W- -> f + fbar         ',  'W+/- + h0 -> f + fbar''      ',
1813      6'h0 + h0 -> f + fbar         ',  'g + g -> g + g              ',
1814      6'gamma + gamma -> W+ + W-    ',  'gamma + W+/- -> Z0 + W+/-   ',
1815      7'Z0 + Z0 -> Z0 + Z0          ',  'Z0 + Z0 -> W+ + W-          ',
1816      7'Z0 + W+/- -> Z0 + W+/-      ',  'Z0 + Z0 -> Z0 + h0          ',
1817      7'W+ + W- -> gamma + gamma    ',  'W+ + W- -> Z0 + Z0          ',
1818      7'W+/- + W+/- -> W+/- + W+/-  ',  'W+/- + h0 -> W+/- + h0      ',
1819      7'h0 + h0 -> h0 + h0          ',  'q + gamma -> q'' + pi+/-     '/
1820       DATA (PROC(I),I=81,100)/
1821      8'q + qbar -> Q + Qbar, mass  ',  'g + g -> Q + Qbar, massive  ',
1822      8'f + q -> f'' + Q, massive    ', 'g + gamma -> Q + Qbar, mass ',
1823      8'gamma + gamma -> F + Fbar, m',  'g + g -> J/Psi + g          ',
1824      8'g + g -> chi_0c + g         ',  'g + g -> chi_1c + g         ',
1825      8'g + g -> chi_2c + g         ',  '                            ',
1826      9'Elastic scattering          ',  'Single diffractive (XB)     ',
1827      9'Single diffractive (AX)     ',  'Double  diffractive         ',
1828      9'Low-pT scattering           ',  'Semihard QCD 2 -> 2         ',
1829      9'                            ',  '                            ',
1830      9'q + gamma* -> q             ',  '                            '/
1831       DATA (PROC(I),I=101,120)/
1832      &'g + g -> gamma*/Z0          ',  'g + g -> h0                 ',
1833      &'gamma + gamma -> h0         ',  'g + g -> chi_0c             ',
1834      &'g + g -> chi_2c             ',  'g + g -> J/Psi + gamma      ',
1835      &'gamma + g -> J/Psi + g      ',  'gamma+gamma -> J/Psi + gamma',
1836      &'                            ',  'f + fbar -> gamma + h0      ',
1837      1'q + qbar -> g + h0          ',  'q + g -> q + h0             ',
1838      1'g + g -> g + h0             ',  'g + g -> gamma + gamma      ',
1839      1'g + g -> g + gamma          ',  'g + g -> gamma + Z0         ',
1840      1'g + g -> Z0 + Z0            ',  'g + g -> W+ + W-            ',
1841      1'                            ',  '                            '/
1842       DATA (PROC(I),I=121,140)/
1843      2'g + g -> Q + Qbar + h0      ',  'q + qbar -> Q + Qbar + h0   ',
1844      2'f + f'' -> f + f'' + h0       ',
1845      2'f + f'' -> f" + f"'' + h0     ',
1846      2'                            ',  '                            ',
1847      2'                            ',  '                            ',
1848      2'                            ',  '                            ',
1849      3'f + gamma*_T -> f + g       ',  'f + gamma*_L -> f + g       ',
1850      3'f + gamma*_T -> f + gamma   ',  'f + gamma*_L -> f + gamma   ',
1851      3'g + gamma*_T -> f + fbar    ',  'g + gamma*_L -> f + fbar    ',
1852      3'gamma*_T+gamma*_T -> f+fbar ',  'gamma*_T+gamma*_L -> f+fbar ',
1853      3'gamma*_L+gamma*_T -> f+fbar ',  'gamma*_L+gamma*_L -> f+fbar '/
1854       DATA (PROC(I),I=141,160)/
1855      4'f + fbar -> gamma*/Z0/Z''0   ', 'f + fbar'' -> W''+/-          ',
1856      4'f + fbar'' -> H+/-           ', 'f + fbar'' -> R              ',
1857      4'q + l -> LQ                 ',  'e + gamma -> e*             ',
1858      4'd + g -> d*                 ',  'u + g -> u*                 ',
1859      4'g + g -> eta_tc             ',  '                            ',
1860      5'f + fbar -> H0              ',  'g + g -> H0                 ',
1861      5'gamma + gamma -> H0         ',  '                            ',
1862      5'                            ',  'f + fbar -> A0              ',
1863      5'g + g -> A0                 ',  'gamma + gamma -> A0         ',
1864      5'                            ',  '                            '/
1865       DATA (PROC(I),I=161,180)/
1866      6'f + g -> f'' + H+/-          ', 'q + g -> LQ + lbar          ',
1867      6'g + g -> LQ + LQbar         ',  'q + qbar -> LQ + LQbar      ',
1868      6'f + fbar -> f'' + fbar'' (g/Z)',
1869      6'f +fbar'' -> f" + fbar"'' (W) ',
1870      6'q + q'' -> q" + d*           ',  'q + q'' -> q" + u*           ',
1871      6'q + qbar -> e + e*          ',  '                            ',
1872      7'f + fbar -> Z0 + H0         ', 'f + fbar'' -> W+/- + H0      ',
1873      7'f + f'' -> f + f'' + H0       ',
1874      7'f + f'' -> f" + f"'' + H0     ',
1875      7'                            ',  'f + fbar -> Z0 + A0         ',
1876      7'f + fbar'' -> W+/- + A0      ',
1877      7'f + f'' -> f + f'' + A0       ',
1878      7'f + f'' -> f" + f"'' + A0     ',
1879      7'                            '/
1880       DATA (PROC(I),I=181,200)/
1881      8'g + g -> Q + Qbar + H0      ',  'q + qbar -> Q + Qbar + H0   ',
1882      8'q + qbar -> g + H0          ',  'q + g -> q + H0             ',
1883      8'g + g -> g + H0             ',  'g + g -> Q + Qbar + A0      ',
1884      8'q + qbar -> Q + Qbar + A0   ',  'q + qbar -> g + A0          ',
1885      8'q + g -> q + A0             ',  'g + g -> g + A0             ',
1886      9'f + fbar -> rho_tc0         ',  'f + f'' -> rho_tc+/-         ',
1887      9'f + fbar -> omega_tc0      ',  'f+fbar -> f''+fbar'' (ETC)  ',
1888      9'f+fbar'' -> f"+fbar"'' (ETC)','                          ',
1889      9'                            ',  '                            ',
1890      9'                            ',  '                            '/
1891       DATA (PROC(I),I=201,220)/
1892      &'f + fbar -> ~e_L + ~e_Lbar  ',  'f + fbar -> ~e_R + ~e_Rbar  ',
1893      &'f + fbar -> ~e_R + ~e_Lbar  ',  'f + fbar -> ~mu_L + ~mu_Lbar',
1894      &'f + fbar -> ~mu_R + ~mu_Rbar',  'f + fbar -> ~mu_L + ~mu_Rbar',
1895      &'f+fbar -> ~tau_1 + ~tau_1bar',  'f+fbar -> ~tau_2 + ~tau_2bar',
1896      &'f+fbar -> ~tau_1 + ~tau_2bar',  'q + qbar'' -> ~l_L + ~nulbar ',
1897      1'q+qbar''-> ~tau_1 + ~nutaubar', 'q+qbar''-> ~tau_2 + ~nutaubar',
1898      1'f + fbar -> ~nul + ~nulbar  ',  'f+fbar -> ~nutau + ~nutaubar',
1899      1'                            ',  'f + fbar -> ~chi1 + ~chi1   ',
1900      1'f + fbar -> ~chi2 + ~chi2   ',  'f + fbar -> ~chi3 + ~chi3   ',
1901      1'f + fbar -> ~chi4 + ~chi4   ',  'f + fbar -> ~chi1 + ~chi2   '/
1902       DATA (PROC(I),I=221,240)/
1903      2'f + fbar -> ~chi1 + ~chi3   ',  'f + fbar -> ~chi1 + ~chi4   ',
1904      2'f + fbar -> ~chi2 + ~chi3   ',  'f + fbar -> ~chi2 + ~chi4   ',
1905      2'f + fbar -> ~chi3 + ~chi4   ',  'f+fbar -> ~chi+-1 + ~chi-+1 ',
1906      2'f+fbar -> ~chi+-2 + ~chi-+2 ',  'f+fbar -> ~chi+-1 + ~chi-+2 ',
1907      2'q + qbar'' -> ~chi1 + ~chi+-1', 'q + qbar'' -> ~chi2 + ~chi+-1',
1908      3'q + qbar'' -> ~chi3 + ~chi+-1', 'q + qbar'' -> ~chi4 + ~chi+-1',
1909      3'q + qbar'' -> ~chi1 + ~chi+-2', 'q + qbar'' -> ~chi2 + ~chi+-2',
1910      3'q + qbar'' -> ~chi3 + ~chi+-2', 'q + qbar'' -> ~chi4 + ~chi+-2',
1911      3'q + qbar -> ~chi1 + ~g      ',  'q + qbar -> ~chi2 + ~g      ',
1912      3'q + qbar -> ~chi3 + ~g      ',  'q + qbar -> ~chi4 + ~g      '/
1913       DATA (PROC(I),I=241,260)/
1914      4'q + qbar'' -> ~chi+-1 + ~g   ', 'q + qbar'' -> ~chi+-2 + ~g  ',
1915      4'q + qbar -> ~g + ~g         ',  'g + g -> ~g + ~g            ',
1916      4'                            ',  'qj + g -> ~qj_L + ~chi1     ',
1917      4'qj + g -> ~qj_R + ~chi1     ',  'qj + g -> ~qj_L + ~chi2     ',
1918      4'qj + g -> ~qj_R + ~chi2     ',  'qj + g -> ~qj_L + ~chi3     ',
1919      5'qj + g -> ~qj_R + ~chi3     ',  'qj + g -> ~qj_L + ~chi4     ',
1920      5'qj + g -> ~qj_R + ~chi4     ',  'qj + g -> ~qk_L + ~chi+-1   ',
1921      5'qj + g -> ~qk_R + ~chi+-1   ',  'qj + g -> ~qk_L + ~chi+-2   ',
1922      5'qj + g -> ~qk_R + ~chi+-2   ',  'qj + g -> ~qj_L + ~g        ',
1923      5'qj + g -> ~qj_R + ~g        ',  '                            '/
1924       DATA (PROC(I),I=261,300)/
1925      6'f + fbar -> ~t_1 + ~t_1bar  ',  'f + fbar -> ~t_2 + ~t_2bar  ',
1926      6'f + fbar -> ~t_1 + ~t_2bar  ',  'g + g -> ~t_1 + ~t_1bar     ',
1927      6'g + g -> ~t_2 + ~t_2bar     ',  '                            ',
1928      6'                            ',  '                            ',
1929      6'                            ',  '                            ',
1930      7'qi + qj -> ~qi_L + ~qj_L    ',  'qi + qj -> ~qi_R + ~qj_R    ',
1931      7'qi + qj -> ~qi_L + ~qj_R    ',  'qi+qjbar -> ~qi_L + ~qj_Lbar',
1932      7'qi+qjbar -> ~qi_R + ~qj_Rbar',  'qi+qjbar -> ~qi_L + ~qj_Rbar',
1933      7'f + fbar -> ~qi_L + ~qi_Lbar',  'f + fbar -> ~qi_R + ~qi_Rbar',
1934      7'g + g -> ~qi_L + ~qi_Lbar   ',  'g + g -> ~qi_R + ~qi_Rbar   ',
1935      8'b + qj -> ~b_1 + ~qj_L      ',  'b + qj -> ~b_2 + ~qj_R      ',
1936      8'b + qj -> ~b_1 + ~qj_R      ',  'b + qjbar -> ~b_1 + ~qj_Lbar',
1937      8'b + qjbar -> ~b_2 + ~qj_Rbar',  'b + qjbar -> ~b_1 + ~qj_Rbar',
1938      8'f + fbar -> ~b_1 + ~b_1bar  ',  'f + fbar -> ~b_2 + ~b_2bar  ',
1939      8'g + g -> ~b_1 + ~b_1bar     ',  'g + g -> ~b_2 + ~b_2bar     ',
1940      9'b + b -> ~b_1 + ~b_1        ',  'b + b -> ~b_2 + ~b_2        ',
1941      9'b + b -> ~b_1 + ~b_2        ',  'b + g -> ~b_1 + ~g          ',
1942      9'b + g -> ~b_2 + ~g          ',  'b + bbar -> ~b_1 + ~b_2bar  ',
1943      9'f + fbar'' -> H+/- + h0     ',  'f + fbar -> H+/- + H0       ',
1944      9'f + fbar -> A0 + h0         ',  'f + fbar -> A0 + H0         '/
1945       DATA (PROC(I),I=301,340)/
1946      &'f + fbar -> H+ + H-         ', 39*'                          '/
1947       DATA (PROC(I),I=341,380)/
1948      4'l + l -> H_L++/--           ',  'l + l -> H_R++/--           ',
1949      4'l + gamma -> H_L++/-- e-/+  ',  'l + gamma -> H_R++/-- e-/+  ',
1950      4'l + gamma -> H_L++/-- mu-/+ ',  'l + gamma -> H_R++/-- mu-/+ ',
1951      4'l + gamma -> H_L++/-- tau-/+',  'l + gamma -> H_R++/-- tau-/+',
1952      4'f + fbar -> H_L++ + H_L--   ',  'f + fbar -> H_R++ + H_R--   ',
1953      5'f + f -> f'' + f'' + H_L++/-- ',
1954      5'f + f -> f'' + f'' + H_R++/-- ','f + fbar -> Z_R0            ',
1955      5'f + fbar'' -> W_R+/-         ',5*'                            ',
1956      6'                            ',  'f + fbar -> W_L+ W_L-       ',
1957      6'f + fbar -> W_L+/- pi_T-/+  ',  'f + fbar -> pi_T+ pi_T-     ',
1958      6'f + fbar -> gamma pi_T0     ',  'f + fbar -> gamma pi_T0''    ',
1959      6'f + fbar -> Z0 pi_T0        ',  'f + fbar -> Z0 pi_T0''       ',
1960      6'f + fbar -> W+/- pi_T-/+    ',  '                            ',
1961      7'f + fbar'' -> W_L+/- Z_L0    ', 'f + fbar'' -> W_L+/- pi_T0   ',
1962      7'f + fbar'' -> pi_T+/- Z_L0   ', 'f + fbar'' -> pi_T+/- pi_T0  ',
1963      7'f + fbar'' -> gamma pi_T+/-  ', 'f + fbar'' -> Z0 pi_T+/-     ',
1964      7'f + fbar'' -> W+/- pi_T0     ',
1965      7'f + fbar'' -> W+/- pi_T0''    ',
1966      7'                            ',  '                            ',
1967      7'                            '/
1968       DATA (PROC(I),I=381,420)/
1969      8'f + f'' -> f + f'' (ETC)      ','f + fbar -> f'' + fbar'' (ETC)',
1970      8'f + fbar -> g + g (ETC)     ',  'f + g -> f + g (ETC)        ',
1971      8'g + g -> f + fbar (ETC)     ',  'g + g -> g + g (ETC)        ',
1972      8'q + qbar -> Q + Qbar (ETC)  ',  'g + g -> Q + Qbar (ETC)     ',
1973      8'                            ',  '                            ',
1974      9'f + fbar -> G*              ',  'g + g -> G*                 ',
1975      9'q + qbar -> g + G*          ',  'q + g -> q + G*             ',
1976      9'g + g -> g + G*             ',  '                            ',
1977      9 4*'                         ',
1978      &'g + g -> t + b + H+/-       ',  'q + qbar -> t + b + H+/-    ',
1979      & 18*'                            '/
1980       DATA (PROC(I),I=421,460)/
1981      2'g + g  -> cc~[3S1(1)] + g   ',  'g + g  -> cc~[3S1(8)] + g   ',
1982      2'g + g  -> cc~[1S0(8)] + g   ',  'g + g  -> cc~[3PJ(8)] + g   ',
1983      2'g + q  -> q + cc~[3S1(8)]   ',  'g + q  -> q + cc~[1S0(8)]   ',
1984      2'g + q  -> q + cc~[3PJ(8)]   ',  'q + q~ -> g + cc~[3S1(8)]   ',
1985      2'q + q~ -> g + cc~[1S0(8)]   ',  'q + q~ -> g + cc~[3PJ(8)]   ',
1986      3'g + g  -> cc~[3P0(1)] + g   ',  'g + g  -> cc~[3P1(1)] + g   ',
1987      3'g + g  -> cc~[3P2(1)] + g   ',  'q + g  -> q + cc~[3P0(1)]   ',
1988      3'q + g  -> q + cc~[3P1(1)]   ',  'q + g  -> q + cc~[3P2(1)]   ',
1989      3'q + q~ -> g + cc~[3P0(1)]   ',  'q + q~ -> g + cc~[3P1(1)]   ',
1990      3'q + q~ -> g + cc~[3P2(1)]   ',
1991      3     21 *'                            '/
1992       DATA (PROC(I),I=461,500)/
1993      6'g + g  -> bb~[3S1(1)] + g   ',  'g + g  -> bb~[3S1(8)] + g   ',
1994      6'g + g  -> bb~[1S0(8)] + g   ',  'g + g  -> bb~[3PJ(8)] + g   ',
1995      6'g + q  -> q + bb~[3S1(8)]   ',  'g + q  -> q + bb~[1S0(8)]   ',
1996      6'g + q  -> q + bb~[3PJ(8)]   ',  'q + q~ -> g + bb~[3S1(8)]   ',
1997      6'q + q~ -> g + bb~[1S0(8)]   ',  'q + q~ -> g + bb~[3PJ(8)]   ',
1998      7'g + g  -> bb~[3P0(1)] + g   ',  'g + g  -> bb~[3P1(1)] + g   ',
1999      7'g + g  -> bb~[3P2(1)] + g   ',  'q + g  -> q + bb~[3P0(1)]   ',
2000      7'q + g  -> q + bb~[3P1(1)]   ',  'q + g  -> q + bb~[3P2(1)]   ',
2001      7'q + q~ -> g + bb~[3P0(1)]   ',  'q + q~ -> g + bb~[3P1(1)]   ',
2002      7'q + q~ -> g + bb~[3P2(1)]   ',
2003      7     21 *'                            '/
2004
2005 C...Cross sections and slope offsets.
2006       DATA SIGT/294*0D0/
2007
2008 C...Supersymmetry switches and parameters.
2009       DATA IMSS/0,
2010      &  0,  0,  0,  1,  0,  0,  0,  0,  0,  0,
2011      1  89*0/
2012       DATA RMSS/0D0,
2013      &  80D0,160D0,500D0,800D0,2D0,250D0,200D0,800D0,700D0,800D0,
2014      1  700D0,500D0,250D0,200D0,800D0,400D0,0D0,0.1D0,850D0,0.041D0,
2015      2   1D0,800D0,1D4,1D4,1D4,0D0,0D0,0D0,24D17,0D0,
2016      3  10*0D0,
2017      4  0D0,1D0,8*0D0,
2018      5  49*0D0/
2019 C...Initial values for R-violating SUSY couplings.
2020 C...Should not be changed here. See PYMSIN.
2021       DATA RVLAM/27*0D0/
2022       DATA RVLAMP/27*0D0/
2023       DATA RVLAMB/27*0D0/
2024
2025 C...Technicolor switches and parameters
2026       DATA ITCM/0,
2027      &  4,  0,  0,  0,  0,  0,  0,  0,  0,  0,
2028      1  89*0/
2029       DATA RTCM/0D0,
2030      &  82D0,1.333D0,.333D0,0.408D0,1D0,1D0,.0182D0,1D0,0D0,1.333D0,
2031      1  .05D0,200D0,200D0,0D0,0D0,0D0,0D0,0D0,0D0,0D0,
2032      2  .283D0,.707D0,0D0,0D0,0D0,1.667D0,250D0,250D0,.707D0,0D0,
2033      3  .707D0,0D0,1D0,0D0,0D0,0D0,0D0,0D0,0D0,0D0,
2034      4  1000D0, 1D0, 1D0, 1D0, 1D0, 0D0, 4*0D0,
2035      4  49*0D0/
2036
2037 C...Data for histogramming routines.
2038       DATA IHIST/1000,20000,55,1/
2039       DATA INDX/1000*0/
2040
2041 C...Data for SUSY Les Houches Accord.
2042       DATA CPRO/'PYTHIA      ','PYTHIA      '/
2043       DATA CVER/'6.4         ','6.4         '/
2044       DATA MODSEL/200*0/
2045       DATA PARMIN/100*0D0/
2046       DATA RMSOFT/101*0D0/
2047       DATA AU/9*0D0/
2048       DATA AD/9*0D0/
2049       DATA AE/9*0D0/
2050
2051       END