Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 13:29:39

0001 C*********************************************************************
0002 C*********************************************************************
0003 C*                                                                  **
0004 C*                                                      June 2010   **
0005 C*                                                                  **
0006 C*                       The Lund Monte Carlo                       **
0007 C*                                                                  **
0008 C*                        PYTHIA version 6.4                        **
0009 C*                                                                  **
0010 C*                        Torbjorn Sjostrand                        **
0011 C*                 Department of Theoretical Physics                **
0012 C*                         Lund University                          **
0013 C*               Solvegatan 14A, S-223 62 Lund, Sweden              **
0014 C*                    E-mail torbjorn@thep.lu.se                    **
0015 C*                                                                  **
0016 C*                  SUSY and Technicolor parts by                   **
0017 C*                         Stephen Mrenna                           **
0018 C*                       Computing Division                         ** 
0019 C*            Generators and Detector Simulation Group              **
0020 C*              Fermi National Accelerator Laboratory               **
0021 C*                 MS 234, Batavia, IL  60510, USA                  **
0022 C*                   phone + 1 - 630 - 840 - 2556                   **
0023 C*                      E-mail mrenna@fnal.gov                      **
0024 C*                                                                  **
0025 C*         New multiple interactions and more SUSY parts by         **
0026 C*                          Peter Skands                            **
0027 C*               CERN/PH, CH-1211 Geneva, Switzerland               **
0028 C*                    phone +41 - 22 - 767 2447                     **
0029 C*                   E-mail peter.skands@cern.ch                    **
0030 C*                                                                  **
0031 C*         Several parts are written by Hans-Uno Bengtsson          **
0032 C*          PYSHOW is written together with Mats Bengtsson          **
0033 C*               PYMAEL is written by Emanuel Norrbin               **
0034 C*     advanced popcorn baryon production written by Patrik Eden    **
0035 C*    code for virtual photons mainly written by Christer Friberg   **
0036 C*    code for low-mass strings mainly written by Emanuel Norrbin   **
0037 C*        Bose-Einstein code mainly written by Leif Lonnblad        **
0038 C*      CTEQ  parton distributions are by the CTEQ collaboration    **
0039 C*      GRV 94 parton distributions are by Glueck, Reya and Vogt    **
0040 C*   SaS photon parton distributions together with Gerhard Schuler  **
0041 C*     g + g and q + qbar -> t + tbar + H code by Zoltan Kunszt     **
0042 C*         MSSM Higgs mass calculation code by M. Carena,           **
0043 C*           J.R. Espinosa, M. Quiros and C.E.M. Wagner             **
0044 C*  UED implementation by M. Elkacimi, D. Goujdami, H. Przysiezniak **
0045 C*         PYGAUS adapted from CERN library (K.S. Kolbig)           **
0046 C*        NRQCD/colour octet production of onium by S. Wolf         **
0047 C*                                                                  **
0048 C*   The latest program version and documentation is found on WWW   **
0049 C*            http://www.thep.lu.se/~torbjorn/Pythia.html           **
0050 C*                                                                  **
0051 C*              Copyright Torbjorn Sjostrand, Lund 2010             **
0052 C*                                                                  **
0053 C*********************************************************************
0054 C*********************************************************************
0055 C                                                                    *
0056 C  List of subprograms in order of appearance, with main purpose     *
0057 C  (S = subroutine, F = function, B = block data)                    *
0058 C                                                                    *
0059 C  B   PYDATA   to contain all default values                        *
0060 C  S   PYCKBD   to check that BLOCK DATA has been correctly loaded   *
0061 C  S   PYTEST   to test the proper functioning of the package        *
0062 C  S   PYHEPC   to convert between /PYJETS/ and /HEPEVT/ records     *
0063 C                                                                    *
0064 C  S   PYINIT   to administer the initialization procedure           *
0065 C  S   PYEVNT   to administer the generation of an event             *
0066 C  S   PYEVNW   ditto, for new multiple interactions scenario        *
0067 C  S   PYSTAT   to print cross-section and other information         *
0068 C  S   PYUPEV   to administer the generation of an LHA hard process  *
0069 C  S   PYUPIN   to provide initialization needed for LHA input       *
0070 C  S   PYLHEF   to produce a Les Houches Event File from run         *
0071 C  S   PYINRE   to initialize treatment of resonances                *
0072 C  S   PYINBM   to read in beam, target and frame choices            *
0073 C  S   PYINKI   to initialize kinematics of incoming particles       *
0074 C  S   PYINPR   to set up the selection of included processes        *
0075 C  S   PYXTOT   to give total, elastic and diffractive cross-sect.   *
0076 C  S   PYMAXI   to find differential cross-section maxima            *
0077 C  S   PYPILE   to select multiplicity of pileup events              *
0078 C  S   PYSAVE   to save alternatives for gamma-p and gamma-gamma     *
0079 C  S   PYGAGA   to handle lepton -> lepton + gamma branchings        *
0080 C  S   PYRAND   to select subprocess and kinematics for event        *
0081 C  S   PYSCAT   to set up kinematics and colour flow of event        *
0082 C  S   PYEVOL   handler for pT-ordered ISR and multiple interactions *
0083 C  S   PYSSPA   to simulate initial state spacelike showers          *
0084 C  S   PYPTIS   to do pT-ordered initial state spacelike showers     *
0085 C  S   PYMEMX   auxiliary to PYSSPA/PYPTIS for ME correction maximum *
0086 C  S   PYMEWT   auxiliary to PYSSPA/.. for matrix element correction *
0087 C  S   PYPTMI   to do pT-ordered multiple interactions               *
0088 C  F   PYFCMP   to give companion quark x*f distribution             *
0089 C  F   PYPCMP   to calculate momentum integral for companion quarks  *
0090 C  S   PYUPRE   to rearranges contents of the HEPEUP commonblock     *
0091 C  S   PYADSH   to administrate sequential final-state showers       *
0092 C  S   PYVETO   to allow the generation of an event to be aborted    *
0093 C  S   PYRESD   to perform resonance decays                          *
0094 C  S   PYMULT   to generate multiple interactions - old scheme       *
0095 C  S   PYREMN   to add on target remnants - old scheme               *
0096 C  S   PYMIGN   to generate multiple interactions - new scheme       *
0097 C  S   PYMIHK   to connect colours in mult. int. - new scheme        *
0098 C  S   PYCTTR   to translate PYTHIA colour information to LHA1 tags  *
0099 C  S   PYMIHG   to collapse two pairs of LHA1 colour tags.           *
0100 C  S   PYMIRM   to add on target remnants in mult. int.- new scheme  *
0101 C  S   PYFSCR   to perform final state colour reconnections - -"-    *
0102 C  S   PYDIFF   to set up kinematics for diffractive events          *
0103 C  S   PYDISG   to set up kinematics, remnant and showers for DIS    *
0104 C  S   PYDOCU   to compute cross-sections and handle documentation   *
0105 C  S   PYFRAM   to perform boosts between different frames           *
0106 C  S   PYWIDT   to calculate full and partial widths of resonances   *
0107 C  S   PYOFSH   to calculate partial width into off-shell channels   *
0108 C  S   PYRECO   to handle colour reconnection in W+W- events         *
0109 C  S   PYKLIM   to calculate borders of allowed kinematical region   *
0110 C  S   PYKMAP   to construct value of kinematical variable           *
0111 C  S   PYSIGH   to calculate differential cross-sections             *
0112 C  S   PYSGQC   auxiliary to PYSIGH for QCD processes                *
0113 C  S   PYSGHF   auxiliary to PYSIGH for heavy flavour processes      *
0114 C  S   PYSGWZ   auxiliary to PYSIGH for W and Z processes            *
0115 C  S   PYSGHG   auxiliary to PYSIGH for Higgs processes              *
0116 C  S   PYSGSU   auxiliary to PYSIGH for supersymmetry processes      *
0117 C  S   PYSGTC   auxiliary to PYSIGH for technicolor processes        *
0118 C  S   PYSGEX   auxiliary to PYSIGH for various exotic processes     *
0119 C  S   PYPDFU   to evaluate parton distributions                     *
0120 C  S   PYPDFL   to evaluate parton distributions at low x and Q^2    *
0121 C  S   PYPDEL   to evaluate electron parton distributions            *
0122 C  S   PYPDGA   to evaluate photon parton distributions (generic)    *
0123 C  S   PYGGAM   to evaluate photon parton distributions (SaS sets)   *
0124 C  S   PYGVMD   to evaluate VMD part of photon parton distributions  *
0125 C  S   PYGANO   to evaluate anomalous part of photon PDFs            *
0126 C  S   PYGBEH   to evaluate Bethe-Heitler part of photon PDFs        *
0127 C  S   PYGDIR   to evaluate direct contribution to photon PDFs       *
0128 C  S   PYPDPI   to evaluate pion parton distributions                *
0129 C  S   PYPDPR   to evaluate proton parton distributions              *
0130 C  F   PYCTEQ   to evaluate the CTEQ 3 proton parton distributions   *
0131 C  S   PYGRVL   to evaluate the GRV 94L proton parton distributions  *
0132 C  S   PYGRVM   to evaluate the GRV 94M proton parton distributions  *
0133 C  S   PYGRVD   to evaluate the GRV 94D proton parton distributions  *
0134 C  F   PYGRVV   auxiliary to the PYGRV* routines                     *
0135 C  F   PYGRVW   auxiliary to the PYGRV* routines                     *
0136 C  F   PYGRVS   auxiliary to the PYGRV* routines                     *
0137 C  F   PYCT5L   to evaluate the CTEQ 5L proton parton distributions  *
0138 C  F   PYCT5M   to evaluate the CTEQ 5M1 proton parton distributions *
0139 C  S   PYPDPO   to evaluate old proton parton distributions          *
0140 C  F   PYHFTH   to evaluate threshold factor for heavy flavour       *
0141 C  S   PYSPLI   to find flavours left in hadron when one removed     *
0142 C  F   PYGAMM   to evaluate ordinary Gamma function Gamma(x)         *
0143 C  S   PYWAUX   to evaluate auxiliary functions W1(s) and W2(s)      *
0144 C  S   PYI3AU   to evaluate auxiliary function I3(s,t,u,v)           *
0145 C  F   PYSPEN   to evaluate Spence (dilogarithm) function Sp(x)      *
0146 C  S   PYQQBH   to evaluate matrix element for g + g -> Q + Qbar + H *
0147 C  S   PYSTBH   to evaluate matrix element for t + b + H processes   *
0148 C  S   PYTBHB   auxiliary to PYSTBH                                  *
0149 C  S   PYTBHG   auxiliary to PYSTBH                                  *
0150 C  S   PYTBHQ   auxiliary to PYSTBH                                  *
0151 C  F   PYTBHS   auxiliary to PYSTBH                                  *
0152 C                                                                    *
0153 C  S   PYMSIN   to initialize the supersymmetry simulation           *
0154 C  S   PYSLHA   to interface to SUSY spectrum and decay calculators  *
0155 C  S   PYAPPS   to determine MSSM parameters from SUGRA input        *
0156 C  S   PYSUGI   to determine MSSM parameters using ISASUSY           *
0157 C  S   PYFEYN   to determine MSSM Higgs parameters using FEYNHIGGS   *
0158 C  F   PYRNMQ   to determine running squark masses                   *
0159 C  S   PYTHRG   to calculate sfermion third-gen. mass eigenstates    *
0160 C  S   PYINOM   to calculate neutralino/chargino mass eigenstates    *
0161 C  F   PYRNM3   to determine running M3, gluino mass                 *
0162 C  S   PYEIG4   to calculate eigenvalues and -vectors in 4*4 matrix  *
0163 C  S   PYHGGM   to determine Higgs mass spectrum                     *
0164 C  S   PYSUBH   to determine Higgs masses in the MSSM                *
0165 C  S   PYPOLE   to determine Higgs masses in the MSSM                *
0166 C  S   PYRGHM   auxiliary to PYPOLE                                  *
0167 C  S   PYGFXX   auxiliary to PYRGHM                                  *
0168 C  F   PYFINT   auxiliary to PYPOLE                                  *
0169 C  F   PYFISB   auxiliary to PYFINT                                  *
0170 C  S   PYSFDC   to calculate sfermion decay partial widths           *
0171 C  S   PYGLUI   to calculate gluino decay partial widths             *
0172 C  S   PYTBBN   to calculate 3-body decay of gluino to neutralino    *
0173 C  S   PYTBBC   to calculate 3-body decay of gluino to chargino      *
0174 C  S   PYNJDC   to calculate neutralino decay partial widths         *
0175 C  S   PYCJDC   to calculate chargino decay partial widths           *
0176 C  F   PYXXZ6   auxiliary for ino 3-body decays                      *
0177 C  F   PYXXGA   auxiliary for ino -> ino + gamma decay               *
0178 C  F   PYX2XG   auxiliary for ino -> ino + gauge boson decay         *
0179 C  F   PYX2XH   auxiliary for ino -> ino + Higgs decay               *
0180 C  S   PYHEXT   to calculate non-SM Higgs decay partial widths       *
0181 C  F   PYH2XX   auxiliary for H -> ino + ino decay                   *
0182 C  F   PYGAUS   to perform Gaussian integration                      *
0183 C  F   PYGAU2   copy of PYGAUS to allow two-dimensional integration  *
0184 C  F   PYSIMP   to perform Simpson integration                       *
0185 C  F   PYLAMF   to evaluate the lambda kinematics function           *
0186 C  S   PYTBDY   to perform 3-body decay of gauginos                  *
0187 C  S   PYTECM   to calculate techni_rho/omega masses                 *
0188 C  S   PYXDIN   to initialize Universal Extra Dimensions             *
0189 C  S   PYUEDC   to compute UED mass radiative corrections            *
0190 C  S   PYXUED   to compute UED cross sections                        *
0191 C  S   PYGRAM   to generate UED G* (excited graviton) mass spectrum  *
0192 C  F   PYGRAW   to compute UED partial widths to G*                  *
0193 C  F   PYWDKK   to compute UED differential partial widths to G*     *
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/PYPUED/IUED(0:99),RUED(0:99)
0353       COMMON/PYBINS/IHIST(4),INDX(1000),BIN(20000)
0354       COMMON/PYLH3P/MODSEL(200),PARMIN(100),PAREXT(200),RMSOFT(0:100),
0355      &     AU(3,3),AD(3,3),AE(3,3)
0356       COMMON/PYLH3C/CPRO(2),CVER(2)
0357       CHARACTER CPRO*12,CVER*12
0358       SAVE /PYDAT1/,/PYDAT2/,/PYDAT3/,/PYDAT4/,/PYDATR/,/PYSUBS/,
0359      &/PYPARS/,/PYINT1/,/PYINT2/,/PYINT3/,/PYINT4/,/PYINT5/,
0360      &/PYINT6/,/PYINT7/,/PYMSSM/,/PYSSMT/,/PYMSRV/,/PYTCSM/,/PYPUED/,
0361      &/PYBINS/,/PYLH3P/,/PYLH3C/
0362  
0363 C...PYDAT1, containing status codes and most parameters.
0364       DATA MSTU/
0365      &   0,    0,    0, 4000,10000,  500, 8000,    0,    0,    2,
0366      1   6,    0,    1,    0,    0,    1,    0,    0,    0,    0,
0367      2   2,   10,    0,    0,    1,   10,    0,    0,    0,    0,
0368      3   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
0369      4   2,    2,    1,    4,    2,    1,    1,    0,    0,    0,
0370      5  25,   24,    0,    1,    0,    0,    0,    0,    0,    0,
0371      6   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
0372      7  30*0,
0373      1   1,    0,    0,    0,    0,    0,    0,    0,    0,    0,
0374      2   1,    5,    3,    5,    0,    0,    0,    0,    0,    0,
0375      &  80*0/
0376       DATA (PARU(I),I=1,100)/
0377      &  3.141592653589793D0, 6.283185307179586D0,
0378      &  0.197327D0, 5.06773D0, 0.389380D0, 2.56819D0,  4*0D0,
0379      1  0.001D0, 0.09D0, 0.01D0, 2D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0,
0380      2  0D0,   0D0,   0D0,   0D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0,
0381      3  0D0,   0D0,   0D0,   0D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0,
0382      4  2.0D0,  1.0D0, 0.25D0,  2.5D0, 0.05D0,
0383      4  0D0,   0D0, 0.0001D0, 0D0,   0D0,
0384      5  2.5D0,1.5D0,7.0D0,1.0D0,0.5D0,2.0D0,3.2D0, 0D0, 0D0, 0D0,
0385      6  40*0D0/
0386       DATA (PARU(I),I=101,200)/
0387      &  0.00729735D0, 0.232D0, 0.007764D0, 1.0D0, 1.16639D-5,
0388      &  0D0, 0D0, 0D0, 0D0,  0D0,
0389      1  0.20D0, 0.25D0, 1.0D0, 4.0D0, 10D0, 0D0, 0D0,  0D0, 0D0, 0D0,
0390      2 -0.693D0, -1.0D0, 0.387D0, 1.0D0, -0.08D0,
0391      2 -1.0D0,  1.0D0,  1.0D0,  1.0D0,  0D0,
0392      3  1.0D0,-1.0D0, 1.0D0,-1.0D0, 1.0D0,  0D0,  0D0, 0D0, 0D0, 0D0,
0393      4  5.0D0, 1.0D0, 1.0D0,  0D0, 1.0D0, 1.0D0,  0D0, 0D0, 0D0, 0D0,
0394      5  1.0D0,   0D0,   0D0,   0D0,   0D0,   0D0, 0D0, 0D0, 0D0, 0D0,
0395      6  1.0D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0,  0D0,  0D0, 0D0, 0D0, 0D0,
0396      7  1.0D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0, 0D0,0D0,0D0,
0397      8  1.0D0, 1.0D0, 1.0D0, 0.0D0, 0.0D0, 1.0D0, 1.0D0, 0D0,0D0,0D0,
0398      9  0D0,  0D0,  0D0,  0D0, 1.0D0,  0D0,  0D0, 0D0, 0D0, 0D0/
0399       DATA MSTJ/
0400      &  1,    3,    0,    0,    0,    0,    0,    0,    0,    0,
0401      1  4,    2,    0,    1,    0,    2,    2,   20,    0,    0,
0402      2  2,    1,    1,    2,    1,    2,    2,    0,    0,    0,
0403      3  0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
0404      4  2,    2,    4,    2,    5,    3,    3,    0,    0,    3,
0405      5  0,    3,    0,    2,    0,    0,    1,    0,    0,    0,
0406      6  40*0,
0407      &  5,    2,    7,    5,    1,    1,    0,    2,    0,    2,
0408      1  0,    0,    0,    0,    1,    1,    0,    0,    0,    0,
0409      2  80*0/
0410       DATA PARJ/
0411      &  0.10D0, 0.30D0, 0.40D0, 0.05D0, 0.50D0,
0412      &  0.50D0, 0.50D0,   0.6D0,   1.2D0,   0.6D0,
0413      1  0.50D0,0.60D0,0.75D0, 0D0, 0D0, 0D0, 0D0, 1.0D0, 1.0D0, 0D0,
0414      2  0.36D0, 1.0D0,0.01D0, 2.0D0,1.0D0,0.4D0, 0D0, 0D0, 0D0, 0D0,
0415      3  0.10D0, 1.0D0, 0.8D0, 1.5D0,0D0,2.0D0,0.2D0, 0D0,0.08D0,1D0,
0416      4  0.3D0, 0.58D0, 0.5D0, 0.9D0,0.5D0,1.0D0,1.0D0,1.5D0,1D0,10D0,
0417      5  0.77D0, 0.77D0, 0.77D0, -0.05D0, -0.005D0,
0418      5  0D0, 0D0, 0D0, 1.0D0, 0D0,
0419      6  4.5D0, 0.7D0, 0D0,0.003D0, 0.5D0, 0.5D0, 0D0, 0D0, 0D0, 0D0,
0420      7  10D0, 1000D0, 100D0, 1000D0, 0D0, 0.7D0,10D0, 0D0,0D0,0.5D0,
0421      8  0.29D0, 1.0D0, 1.0D0,  0D0,  10D0, 10D0, 0D0, 0D0, 0D0,1D-4,
0422      9  0.02D0, 1.0D0, 0.2D0,  0D0,  0D0,  0D0,  0D0, 0D0, 0D0, 0D0,
0423      &  0D0,  0D0,  0D0,  0D0,   0D0,   0D0,  0D0,  0D0,  0D0,  0D0,
0424      1  0D0,  0D0,  0D0,  0D0,   0D0,   0D0,  0D0,  0D0,  0D0,  0D0,
0425      2  1.0D0, 0.25D0,91.187D0,2.489D0, 0.01D0,
0426      2  2.0D0,  1.0D0, 0.25D0,0.002D0,   0D0,
0427      3  0D0, 0D0, 0D0, 0D0, 0.01D0, 0.99D0, 0D0, 0D0,  0.2D0,   0D0,
0428      4  10*0D0,
0429      5  10*0D0,
0430      6  10*0D0,
0431      7  0D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0, -0.693D0,
0432      8 -1.0D0, 0.387D0, 1.0D0, -0.08D0, -1.0D0,
0433      8  1.0D0,  1.0D0, -0.693D0, -1.0D0, 0.387D0,
0434      9  1.0D0, -0.08D0, -1.0D0,   1.0D0, 1.0D0,
0435      9  5*0D0/
0436  
0437 C...PYDAT2, with particle data and flavour treatment parameters.
0438       DATA (KCHG(I,1),I=   1, 500)/-1,2,-1,2,-1,2,-1,2,2*0,-3,0,-3,0,   
0439      &-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,  
0440      &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,  
0441      &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,   
0442      &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,    
0443      &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,  
0444      &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,  
0445      &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,  
0446      &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,  
0447      &7*0,3,
0448 C...UED singlet and doublet quarks, leptons, and KK g, gamma, Z, and W
0449      &81*0,-1,2,-1,2,-1,2,-1,2,-1,2,-1,2, 
0450      &3*-3,0,-3,0,-3,0,-3,
0451      &3*0,3, 
0452      &25*0/
0453       DATA (KCHG(I,2),I=   1, 500)/8*1,12*0,2,20*0,1,107*0,-1,0,2*-1,   
0454      &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,   
0455      &-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, 
0456      &6*1,9*0,2,3*0,2,0,5*2,2*1,17*0,6*2,
0457      &83*0,12*1,9*0,2,3*0,25*0/
0458       DATA (KCHG(I,3),I=   1, 500)/8*1,2*0,8*1,5*0,1,9*0,1,2*0,1,3*0,   
0459      &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, 
0460      &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, 
0461      &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,7*0,1,
0462      &81*0,21*1,3*0,1,25*0/
0463       DATA (KCHG(I,4),I=   1, 290)/1,2,3,4,5,6,7,8,9,10,11,12,13,14,15, 
0464      &16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,   
0465      &37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,   
0466      &58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,   
0467      &79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,   
0468      &100,110,111,113,115,130,211,213,215,221,223,225,310,311,313,315,  
0469      &321,323,325,331,333,335,411,413,415,421,423,425,431,433,435,441,  
0470      &443,445,511,513,515,521,523,525,531,533,535,541,543,545,551,553,  
0471      &555,990,1103,1114,2101,2103,2112,2114,2203,2212,2214,2224,3101,   
0472      &3103,3112,3114,3122,3201,3203,3212,3214,3222,3224,3303,3312,3314, 
0473      &3322,3324,3334,4101,4103,4112,4114,4122,4132,4201,4203,4212,4214, 
0474      &4222,4224,4232,4301,4303,4312,4314,4322,4324,4332,4334,4403,4412, 
0475      &4414,4422,4424,4432,4434,4444,5101,5103,5112,5114,5122,5132,5142, 
0476      &5201,5203,5212,5214,5222,5224,5232,5242,5301,5303,5312,5314,5322, 
0477      &5324,5332,5334,5342,5401,5403,5412,5414,5422,5424,5432,5434,5442, 
0478      &5444,5503,5512,5514,5522,5524,5532,5534,5542,5544,5554,10111,     
0479      &10113,10211,10213,10221,10223,10311,10313,10321,10323,10331,      
0480      &10333,10411,10413,10421,10423,10431,10433,10441,10443,10511,      
0481      &10513,10521,10523,10531,10533,10541,10543,10551,10553,20113,      
0482      &20213,20223,20313,20323,20333,20413,20423,20433,20443,20513/      
0483       DATA (KCHG(I,4),I= 291, 500)/20523,20533,20543,20553,100443,      
0484      &100553,1000001,1000002,1000003,1000004,1000005,1000006,1000011,   
0485      &1000012,1000013,1000014,1000015,1000016,1000021,1000022,1000023,  
0486      &1000024,1000025,1000035,1000037,1000039,2000001,2000002,2000003,  
0487      &2000004,2000005,2000006,2000011,2000012,2000013,2000014,2000015,  
0488      &2000016,3000111,3000211,3000221,3000331,3000113,3000213,3000223,  
0489      &3100021,3100111,3200111,3100113,3200113,3300113,3400113,4000001,  
0490      &4000002,4000011,4000012,5000039,9900012,9900014,9900016,9900023,  
0491      &9900024,9900041,9900042,9900110,9900210,9900220,9900330,9900440,  
0492      &9902110,9902210,9900443,9900441,9910441,9900553,9900551,9910551,  
0493      &3000115,3000215,
0494      &81*0,
0495 C...UED singlet and doublet quarks and leptons, and KK g, gamma, Z, and W.
0496      &6100001,6100002,6100003,6100004,6100005,6100006, 
0497      &5100001,5100002,5100003,5100004,5100005,5100006, 
0498      &6100011,6100013,6100015,
0499      &5100012,5100011,5100014,5100013,5100016,5100015, 
0500      &5100021,5100022,5100023,5100024,
0501      &25*0/ 
0502       DATA (PMAS(I,1),I=   1, 217)/2*0.33D0,0.5D0,1.5D0,4.8D0,175D0,    
0503      &2*400D0,2*0D0,0.00051D0,0D0,0.10566D0,0D0,1.777D0,0D0,400D0,      
0504      &5*0D0,91.188D0,80.45D0,115D0,6*0D0,500D0,900D0,500D0,3*300D0,     
0505      &3*0D0,5000D0,200D0,40*0D0,1D0,2D0,5D0,16*0D0,0.13498D0,0.7685D0,  
0506      &1.318D0,0.49767D0,0.13957D0,0.7669D0,1.318D0,0.54745D0,0.78194D0, 
0507      &1.275D0,2*0.49767D0,0.8961D0,1.432D0,0.4936D0,0.8916D0,1.425D0,   
0508      &0.95777D0,1.0194D0,1.525D0,1.8693D0,2.01D0,2.46D0,1.8645D0,       
0509      &2.0067D0,2.46D0,1.9685D0,2.1124D0,2.5735D0,2.9798D0,3.09688D0,    
0510      &3.5562D0,5.2792D0,5.3248D0,5.83D0,5.2789D0,5.3248D0,5.83D0,       
0511      &5.3693D0,5.4163D0,6.07D0,6.594D0,6.602D0,7.35D0,9.4D0,9.4603D0,   
0512      &9.9132D0,0D0,0.77133D0,1.234D0,0.57933D0,0.77133D0,0.93957D0,     
0513      &1.233D0,0.77133D0,0.93827D0,1.232D0,1.231D0,0.80473D0,0.92953D0,  
0514      &1.19744D0,1.3872D0,1.11568D0,0.80473D0,0.92953D0,1.19255D0,       
0515      &1.3837D0,1.18937D0,1.3828D0,1.09361D0,1.3213D0,1.535D0,1.3149D0,  
0516      &1.5318D0,1.67245D0,1.96908D0,2.00808D0,2.4521D0,2.5D0,2.2849D0,   
0517      &2.4703D0,1.96908D0,2.00808D0,2.4535D0,2.5D0,2.4529D0,2.5D0,       
0518      &2.4656D0,2.15432D0,2.17967D0,2.55D0,2.63D0,2.55D0,2.63D0,2.704D0, 
0519      &2.8D0,3.27531D0,3.59798D0,3.65648D0,3.59798D0,3.65648D0,          
0520      &3.78663D0,3.82466D0,4.91594D0,5.38897D0,5.40145D0,5.8D0,5.81D0,   
0521      &5.641D0,5.84D0,7.00575D0,5.38897D0,5.40145D0,5.8D0,5.81D0,5.8D0/  
0522       DATA (PMAS(I,1),I= 218, 500)/5.81D0,5.84D0,7.00575D0,5.56725D0,   
0523      &5.57536D0,5.96D0,5.97D0,5.96D0,5.97D0,6.12D0,6.13D0,7.19099D0,    
0524      &6.67143D0,6.67397D0,7.03724D0,7.0485D0,7.03724D0,7.0485D0,        
0525      &7.21101D0,7.219D0,8.30945D0,8.31325D0,10.07354D0,10.42272D0,      
0526      &10.44144D0,10.42272D0,10.44144D0,10.60209D0,10.61426D0,           
0527      &11.70767D0,11.71147D0,15.11061D0,0.9835D0,1.231D0,0.9835D0,       
0528      &1.231D0,1D0,1.17D0,1.429D0,1.29D0,1.429D0,1.29D0,2*1.4D0,2.272D0, 
0529      &2.424D0,2.272D0,2.424D0,2.5D0,2.536D0,3.4151D0,3.46D0,5.68D0,     
0530      &5.73D0,5.68D0,5.73D0,5.92D0,5.97D0,7.25D0,7.3D0,9.8598D0,9.875D0, 
0531      &2*1.23D0,1.282D0,2*1.402D0,1.427D0,2*2.372D0,2.56D0,3.5106D0,     
0532      &2*5.78D0,6.02D0,7.3D0,9.8919D0,3.686D0,10.0233D0,32*500D0,        
0533      &3*110D0,350D0,3*210D0,500D0,125D0,250D0,400D0,2*350D0,300D0,      
0534      &4*400D0,1000D0,3*500D0,1200D0,750D0,2*200D0,7*0D0,3*3.1D0,        
0535      &3*9.5D0,2*250D0,
0536      &81*0,
0537 C...UED
0538      &586.,588.,586.,588.,586.,586.,6*598.,
0539      &3*505.,6*516.,640.,501.,536.,536.,25*0.D0/
0540       DATA (PMAS(I,2),I=   1, 500)/5*0D0,1.39816D0,16*0D0,2.47813D0,    
0541      &2.07115D0,0.00367D0,6*0D0,14.54029D0,0D0,16.66099D0,8.38842D0,    
0542      &3.3752D0,4.17669D0,3*0D0,417.29147D0,0.39162D0,60*0D0,0.151D0,   
0543      &0.107D0,2*0D0,0.149D0,0.107D0,0D0,0.00843D0,0.185D0,2*0D0,        
0544      &0.0505D0,0.109D0,0D0,0.0498D0,0.098D0,0.0002D0,0.00443D0,0.076D0, 
0545      &2*0D0,0.023D0,2*0D0,0.023D0,2*0D0,0.015D0,0.0013D0,0D0,0.002D0,   
0546      &2*0D0,0.02D0,2*0D0,0.02D0,2*0D0,0.02D0,2*0D0,0.02D0,5*0D0,0.12D0, 
0547      &3*0D0,0.12D0,2*0D0,2*0.12D0,3*0D0,0.0394D0,4*0D0,0.036D0,0D0,     
0548      &0.0358D0,2*0D0,0.0099D0,0D0,0.0091D0,74*0D0,0.06D0,0.142D0,       
0549      &0.06D0,0.142D0,0D0,0.36D0,0.287D0,0.09D0,0.287D0,0.09D0,0.25D0,   
0550      &0.08D0,0.05D0,0.02D0,0.05D0,0.02D0,0.05D0,0D0,0.014D0,0.01D0,     
0551      &8*0.05D0,0D0,0.01D0,2*0.4D0,0.025D0,2*0.174D0,0.053D0,3*0.05D0,   
0552      &0.0009D0,4*0.05D0,3*0D0,19*1D0,0D0,7*1D0,0D0,1D0,0D0,1D0,0D0,     
0553      &0.0208D0,0.01195D0,0.03705D0,0.09511D0,1.89978D0,1.60746D0,       
0554      &0.13396D0,200.47294D0,0.02296D0,0.18886D0,94.66794D0,6.08718D0,   
0555      &0D0,2.17482D0,2.59359D0,2.59687D0,0.42896D0,0.41912D0,0.14153D0,  
0556      &2*0.00098D0,0.00097D0,26.7245D0,21.74916D0,0.88159D0,0.88001D0,   
0557      &7*0D0,6*0.01D0,0.25499D0,0.28446D0,131*0D0/                       
0558       DATA (PMAS(I,3),I=   1, 500)/5*0D0,13.98156D0,16*0D0,24.78129D0,  
0559      &20.71149D0,0.03669D0,6*0D0,145.40294D0,0D0,166.60993D0,           
0560      &83.88423D0,33.75195D0,41.76694D0,3*0D0,4172.91467D0,3.91621D0,    
0561      &60*0D0,0.4D0,0.25D0,2*0D0,0.4D0,0.25D0,0D0,0.1D0,0.17D0,2*0D0,    
0562      &0.2D0,0.12D0,0D0,0.2D0,0.12D0,0.002D0,0.015D0,0.2D0,2*0D0,0.12D0, 
0563      &2*0D0,0.12D0,2*0D0,0.05D0,0.005D0,0D0,0.01D0,2*0D0,0.05D0,2*0D0,  
0564      &0.05D0,2*0D0,0.05D0,2*0D0,0.05D0,5*0D0,0.14D0,3*0D0,0.14D0,2*0D0, 
0565      &2*0.14D0,3*0D0,0.04D0,4*0D0,0.035D0,0D0,0.035D0,2*0D0,0.05D0,0D0, 
0566      &0.05D0,74*0D0,0.05D0,0.25D0,0.05D0,0.25D0,0D0,0.2D0,0.4D0,        
0567      &0.005D0,0.4D0,0.01D0,0.35D0,0.001D0,0.1D0,0.08D0,0.1D0,0.08D0,    
0568      &0.1D0,0D0,0.05D0,0.02D0,6*0.1D0,0.05D0,0.1D0,0D0,0.02D0,2*0.3D0,  
0569      &0.05D0,2*0.3D0,0.02D0,2*0.1D0,0.03D0,0.001D0,4*0.1D0,3*0D0,       
0570      &19*10D0,0.00001D0,7*10D0,0.00001D0,10D0,0.00001D0,10D0,0.00001D0, 
0571      &0.20797D0,0.11949D0,0.37048D0,0.95114D0,18.99785D0,16.07463D0,    
0572      &1.33964D0,450D0,0.22959D0,1.88863D0,360D0,60.8718D0,0D0,          
0573      &21.74824D0,25.93594D0,25.96873D0,4.28961D0,4.19124D0,1.41528D0,   
0574      &0.00977D0,0.00976D0,0.00973D0,267.24501D0,217.49162D0,8.81592D0,  
0575      &8.80013D0,13*0D0,2.54987D0,2.84456D0,
0576      &81*0,
0577 C...UED
0578      &12*0.2D0,9*0.1D0,0.2,10.,0.07,0.3,25*0.D0/
0579       DATA (PMAS(I,4),I=   1, 500)/12*0D0,658654D0,0D0,0.0872D0,68*0D0, 
0580      &0.1D0,0.387D0,16*0D0,0.00003D0,2*0D0,15500D0,7804.5D0,5*0D0,      
0581      &26.762D0,3*0D0,3709D0,5*0D0,0.317D0,2*0D0,0.1244D0,2*0D0,0.14D0,  
0582      &5*0D0,0.468D0,2*0D0,0.462D0,2*0D0,0.483D0,2*0D0,0.15D0,18*0D0,    
0583      &44.34D0,0D0,78.88D0,4*0D0,23.96D0,2*0D0,49.1D0,0D0,87.1D0,0D0,    
0584      &24.6D0,4*0D0,0.0618D0,0.029D0,6*0D0,0.106D0,6*0D0,0.019D0,2*0D0,  
0585      &7*0.1D0,4*0D0,0.342D0,2*0.387D0,6*0D0,2*0.387D0,6*0D0,0.387D0,    
0586      &0D0,0.387D0,2*0D0,8*0.387D0,0D0,9*0.387D0,120*0D0,131*0D0/        
0587 
0588       DATA PARF/
0589      &  0.5D0,0.25D0, 0.5D0,0.25D0, 1D0, 0.5D0,  0D0,  0D0,  0D0, 0D0,
0590      1  0.5D0,  0D0, 0.5D0,  0D0,  1D0,  1D0,  0D0,  0D0,  0D0, 0D0,
0591      2  0.5D0,  0D0, 0.5D0,  0D0,  1D0,  1D0,  0D0,  0D0,  0D0, 0D0,
0592      3  0.5D0,  0D0, 0.5D0,  0D0,  1D0,  1D0,  0D0,  0D0,  0D0, 0D0,
0593      4  0.5D0,  0D0, 0.5D0,  0D0,  1D0,  1D0,  0D0,  0D0,  0D0, 0D0,
0594      5  0.5D0,  0D0, 0.5D0,  0D0,  1D0,  1D0,  0D0,  0D0,  0D0, 0D0,
0595      6  0.75D0, 0.5D0, 0D0,0.1667D0,0.0833D0,0.1667D0,0D0,0D0,0D0, 0D0,
0596      7  0D0,  0D0,  1D0,0.3333D0,0.6667D0,0.3333D0,0D0,0D0,0D0, 0D0,
0597      8  0D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0, 0D0,
0598      9  0.0099D0, 0.0056D0, 0.199D0, 1.23D0, 4.17D0, 165D0,  4*0D0,
0599      & 0.325D0,0.325D0,0.5D0,1.6D0, 5.0D0,  0D0,  0D0,  0D0,  0D0, 0D0,
0600      1 0D0,0.11D0,0.16D0,0.048D0,0.50D0,0.45D0,0.55D0,0.60D0,0D0,0D0,
0601      2 0.2D0, 0.1D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0, 0D0,
0602      3 60*0D0,
0603      4 0.2D0,  0.5D0,  8*0D0,
0604      5 1800*0D0/
0605       DATA ((VCKM(I,J),J=1,4),I=1,4)/
0606      &  0.95113D0,  0.04884D0,  0.00003D0,  0.00000D0,
0607      &  0.04884D0,  0.94940D0,  0.00176D0,  0.00000D0,
0608      &  0.00003D0,  0.00176D0,  0.99821D0,  0.00000D0,
0609      &  0.00000D0,  0.00000D0,  0.00000D0,  1.00000D0/
0610  
0611 C...PYDAT3, with particle decay parameters and data.
0612       DATA (MDCY(I,1),I=   1, 500)/5*0,3*1,6*0,1,0,1,5*0,3*1,6*0,1,0,   
0613      &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, 
0614      &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,  
0615      &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,8*1,
0616      &81*0,
0617 C...UED
0618      &5*1,0,5*1,0,13*1,25*0/
0619       DATA (MDCY(I,2),I=   1, 351)/1,9,17,25,33,41,56,66,2*0,76,80,82,  
0620      &87,89,143,145,150,2*0,153,162,174,190,210,6*0,289,0,311,334,420,  
0621      &503,3*0,530,539,40*0,540,541,545,16*0,554,556,561,570,579,581,    
0622      &583,590,598,604,613,615,617,620,630,636,639,650,656,667,673,736,  
0623      &739,747,808,810,818,851,853,857,858,861,863,899,900,908,944,945,  
0624      &953,992,993,997,1028,1029,1033,1034,1043,2*0,1045,3*0,1046,2*0,   
0625      &1049,1052,2*0,1053,1055,1058,2*0,1062,1063,1066,1069,0,1072,1077, 
0626      &1079,1082,1084,2*0,1088,1089,1090,1166,2*0,1170,1171,1172,1173,   
0627      &1174,2*0,1178,1179,1181,1182,1184,1188,0,1189,1193,1197,1201,     
0628      &1205,1209,1213,2*0,1217,1218,1219,1236,1245,2*0,1254,1255,1256,   
0629      &1257,1258,1267,2*0,1276,1277,1278,1279,1280,1289,1290,2*0,1299,   
0630      &1308,1317,1326,1335,1344,1353,1362,0,1371,1380,1389,1398,1407,    
0631      &1416,1425,1434,1443,1452,1453,1454,1455,1456,1461,1464,1466,1471, 
0632      &1473,1478,1485,1489,1491,1493,1495,1497,1499,1501,1503,1504,1506, 
0633      &1508,1510,1512,1514,1516,1518,1520,1522,1523,1525,1527,1541,1543, 
0634      &1545,1549,1551,1553,1555,1557,1559,1561,1563,1565,1567,1578,1592, 
0635      &1637,1661,1706,1730,1775,1802,1833,1859,1891,1917,1949,1975,2162, 
0636      &2331,2595,2826,3106,3402,0,3657,3706,3734,3783,3811,3860,3888,0,  
0637      &3924,0,3960,0,3996,4004,4012,4020,4217,4243,4270,4023,4029,4036,  
0638      &4043,4050,4056,4062,4071,4075,4079,4082,4084,4104,4126,4148,4170/ 
0639       DATA (MDCY(I,2),I= 352, 500)/4185,4197,4204,7*0,4211,4212,4213,   
0640      &4214,4215,4216,4296,4322,
0641      &81*0,
0642 C...UED
0643      %5001,5003,5005,5007,5009,5011,5013,5016,5019,5022,5025,5028,
0644      &5031,5032,5033,
0645      &5034,5035,5036,5037,5038,5039,5040,5064,5065,5083,
0646      &25*0/
0647       DATA (MDCY(I,3),I=   1, 500)/5*8,15,2*10,2*0,4,2,5,2,54,2,5,3,    
0648      &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, 
0649      &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,  
0650      &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,  
0651      &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, 
0652      &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, 
0653      &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,   
0654      &45,27,31,26,32,26,32,26,187,169,264,231,280,296,255,0,49,28,49,   
0655      &28,49,28,36,0,36,0,36,0,3*8,3,26,27,26,6,3*7,2*6,9,2*4,3,2,20,    
0656      &3*22,15,12,2*7,7*0,6*1,26,30,
0657      &81*0,
0658 C...UED
0659      &6*2,6*3,9*1,24,1,18,6,25*0/                                 
0660       DATA (MDME(I,1),I=   1,8000)/6*1,-1,7*1,-1,7*1,-1,7*1,-1,7*1,-1,  
0661      &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,  
0662      &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,6*1,  
0663      &2*-1,3*1,-1,5*1,62*1,6*1,2*-1,6*1,8*-1,3*1,-1,3*1,-1,3*1,5*-1,   
0664      &3*1,4*-1,6*1,2*-1,3*1,-1,12*1,62*1,6*1,2*-1,3*1,-1,9*1,62*1,    
0665      &3*1,-1,3*1,-1,1,18*1,4*1,2*-1,2*1,-1,1249*1,2*-1,377*1,2*-1,     
0666      &1921*1,2*-1,6*1,2*-1,133*1,2*-1,6*1,2*-1,10*1,-1,3*1,-1,3*1,5*-1, 
0667      &3*1,-1,16*1,2*-1,6*1,2*-1,16*1,2*-1,6*1,2*-1,13*1,-1,3*1,-1,3*1,  
0668      &5*-1,3*1,-1,
0669      &649*0,
0670 C...UED
0671      &10*1,2*0,15*1,3*0,9*1,5*1,0,5*1,0,5*1,0,5*1,0,
0672      &1,24*1,2912*0/
0673       DATA (MDME(I,2),I=   1,8000)/43*102,4*0,102,0,6*53,3*102,4*0,102, 
0674      &2*0,3*102,4*0,102,2*0,6*102,42,6*102,2*42,2*0,8*41,2*0,36*41,     
0675      &8*102,0,102,0,102,2*0,21*102,8*32,8*0,16*32,4*0,8*32,9*0,62*53,   
0676      &8*32,14*0,16*32,7*0,8*32,16*0,62*53,8*32,13*0,62*53,4*32,5*0,     
0677      &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,    
0678      &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,  
0679      &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,     
0680      &14*42,16*0,48,3*13,2*42,9*0,14*42,16*0,48,3*13,2*42,9*0,14*42,    
0681      &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,   
0682      &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,    
0683      &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, 
0684      &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, 
0685      &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,   
0686      &3*0,6*32,3*0,4*32,3*0,50*32,3*53,12*0,8*32,12*0,66*51,6*32,9*0,   
0687      &9*32,17*0,6*51,10*0,8*32,15*0,16*32,14*0,8*32,18*0,8*32,18*0,     
0688      &16*32,
0689 C...UED
0690      &653*0,30*0,9*0,12*0,37*0,2912*0/
0691       DATA (BRAT(I)  ,I=   1, 348)/43*0D0,0.00003D0,0.001765D0,         
0692      &0.998205D0,35*0D0,1D0,6*0D0,0.1783D0,0.1735D0,0.1131D0,0.2494D0,  
0693      &0.003D0,0.09D0,0.0027D0,0.01D0,0.0014D0,0.0012D0,2*0.00025D0,     
0694      &0.0071D0,0.012D0,0.0004D0,0.00075D0,0.00006D0,2*0.00078D0,        
0695      &0.0034D0,0.08D0,0.011D0,0.0191D0,0.00006D0,0.005D0,0.0133D0,      
0696      &0.0067D0,0.0005D0,0.0035D0,0.0006D0,0.0015D0,0.00021D0,0.0002D0,  
0697      &0.00075D0,0.0001D0,0.0002D0,0.0011D0,3*0.0002D0,0.00022D0,        
0698      &0.0004D0,0.0001D0,2*0.00205D0,2*0.00069D0,0.00025D0,0.00051D0,    
0699      &0.00025D0,35*0D0,0.153995D0,0.11942D0,0.153984D0,0.119259D0,      
0700      &0.152272D0,3*0D0,0.033576D0,0.066806D0,0.033576D0,0.066806D0,     
0701      &0.0335D0,0.066806D0,2*0D0,0.321369D0,0.016494D0,2*0D0,0.016502D0, 
0702      &0.320615D0,2*0D0,0.00001D0,0.000591D0,6*0D0,2*0.108166D0,         
0703      &0.108087D0,0D0,0.000001D0,0D0,0.000353D0,0.04359D0,0.795274D0,    
0704      &4*0D0,0.000339D0,0.095746D0,0D0,0.060724D0,0.003054D0,0.000919D0, 
0705      &64*0D0,0.145835D0,0.113276D0,0.145835D0,0.113271D0,0.145781D0,    
0706      &0.049002D0,2*0D0,0.032025D0,0.063642D0,0.032025D0,0.063642D0,     
0707      &0.032022D0,0.063642D0,8*0D0,0.251225D0,0.0129D0,0.000006D0,0D0,   
0708      &0.0129D0,0.250764D0,0.00038D0,0D0,0.000008D0,0.000465D0,          
0709      &0.215418D0,5*0D0,2*0.085312D0,0.08531D0,7*0D0,0.000029D0,         
0710      &0.000536D0,5*0D0,0.000074D0,0D0,0.000417D0,0.000015D0,0.000061D0/ 
0711       DATA (BRAT(I)  ,I= 349, 655)/0.306789D0,0.689189D0,0D0,0.00289D0, 
0712      &69*0D0,0.000001D0,0.000072D0,0.001333D0,4*0D0,0.000001D0,         
0713      &0.000184D0,0D0,0.003108D0,0.000015D0,0.000003D0,2*0D0,0.995284D0, 
0714      &66*0D0,0.000014D0,0.082234D0,2*0D0,0.000013D0,0.003746D0,0D0,     
0715      &0.913992D0,18*0D0,3*0.215119D0,0.214724D0,2*0D0,0.06996D0,        
0716      &0.069959D0,0D0,2*1D0,2*0.08D0,0.76D0,0.08D0,2*0.105D0,0.04D0,     
0717      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,0.988D0,0.012D0,       
0718      &0.998739D0,0.00079D0,0.00038D0,0.000046D0,0.000045D0,2*0.34725D0, 
0719      &0.144D0,0.104D0,0.0245D0,2*0.01225D0,0.0028D0,0.0057D0,0.2112D0,  
0720      &0.1256D0,2*0.1939D0,2*0.1359D0,0.002D0,0.001D0,0.0006D0,          
0721      &0.999877D0,0.000123D0,0.99955D0,0.00045D0,2*0.34725D0,0.144D0,    
0722      &0.104D0,0.049D0,0.0028D0,0.0057D0,0.3923D0,0.321D0,0.2317D0,      
0723      &0.0478D0,0.0049D0,0.0013D0,0.0003D0,0.0007D0,0.89D0,0.08693D0,    
0724      &0.0221D0,0.00083D0,2*0.00007D0,0.564D0,0.282D0,0.072D0,0.028D0,   
0725      &0.023D0,2*0.0115D0,0.005D0,0.003D0,0.6861D0,0.3139D0,2*0.5D0,     
0726      &0.665D0,0.333D0,0.002D0,0.333D0,0.166D0,0.168D0,0.084D0,0.087D0,  
0727      &0.043D0,0.059D0,2*0.029D0,0.002D0,0.6352D0,0.2116D0,0.0559D0,     
0728      &0.0173D0,0.0482D0,0.0318D0,0.666D0,0.333D0,0.001D0,0.332D0,       
0729      &0.166D0,0.168D0,0.084D0,0.086D0,0.043D0,0.059D0,2*0.029D0,        
0730      &2*0.002D0,0.437D0,0.208D0,0.302D0,0.0302D0,0.0212D0,0.0016D0/     
0731       DATA (BRAT(I)  ,I= 656, 831)/0.48947D0,0.34D0,3*0.043D0,0.027D0,  
0732      &0.0126D0,0.0013D0,0.0003D0,0.00025D0,0.00008D0,0.444D0,2*0.222D0, 
0733      &0.104D0,2*0.004D0,0.07D0,0.065D0,2*0.005D0,2*0.011D0,5*0.001D0,   
0734      &0.07D0,0.065D0,2*0.005D0,2*0.011D0,5*0.001D0,0.026D0,0.019D0,     
0735      &0.066D0,0.041D0,0.045D0,0.076D0,0.0073D0,2*0.0047D0,0.026D0,      
0736      &0.001D0,0.0006D0,0.0066D0,0.005D0,2*0.003D0,2*0.0006D0,2*0.001D0, 
0737      &0.006D0,0.005D0,0.012D0,0.0057D0,0.067D0,0.008D0,0.0022D0,        
0738      &0.027D0,0.004D0,0.019D0,0.012D0,0.002D0,0.009D0,0.0218D0,0.001D0, 
0739      &0.022D0,0.087D0,0.001D0,0.0019D0,0.0015D0,0.0028D0,0.683D0,       
0740      &0.306D0,0.011D0,0.3D0,0.15D0,0.16D0,0.08D0,0.13D0,0.06D0,0.08D0,  
0741      &0.04D0,0.034D0,0.027D0,2*0.002D0,2*0.004D0,2*0.002D0,0.034D0,     
0742      &0.027D0,2*0.002D0,2*0.004D0,2*0.002D0,0.0365D0,0.045D0,0.073D0,   
0743      &0.062D0,3*0.021D0,0.0061D0,0.015D0,0.025D0,0.0088D0,0.074D0,      
0744      &0.0109D0,0.0041D0,0.002D0,0.0035D0,0.0011D0,0.001D0,0.0027D0,     
0745      &2*0.0016D0,0.0018D0,0.011D0,0.0063D0,0.0052D0,0.018D0,0.016D0,    
0746      &0.0034D0,0.0036D0,0.0009D0,0.0006D0,0.015D0,0.0923D0,0.018D0,     
0747      &0.022D0,0.0077D0,0.009D0,0.0075D0,0.024D0,0.0085D0,0.067D0,       
0748      &0.0511D0,0.017D0,0.0004D0,0.0028D0,0.619D0,0.381D0,0.3D0,0.15D0,  
0749      &0.16D0,0.08D0,0.13D0,0.06D0,0.08D0,0.04D0,0.01D0,2*0.02D0,0.03D0, 
0750      &2*0.005D0,2*0.02D0,0.03D0,2*0.005D0,0.015D0,0.037D0,0.028D0/      
0751       DATA (BRAT(I)  ,I= 832, 997)/0.079D0,0.095D0,0.052D0,0.0078D0,    
0752      &4*0.001D0,0.028D0,0.033D0,0.026D0,0.05D0,0.01D0,4*0.005D0,0.25D0, 
0753      &0.0952D0,0.94D0,0.06D0,2*0.4D0,2*0.1D0,1D0,0.0602D0,0.0601D0,     
0754      &0.8797D0,0.135D0,0.865D0,0.02D0,0.055D0,2*0.005D0,0.008D0,        
0755      &0.012D0,0.02D0,0.055D0,2*0.005D0,0.008D0,0.012D0,0.01D0,0.03D0,   
0756      &0.0035D0,0.011D0,0.0055D0,0.0042D0,0.009D0,0.018D0,0.015D0,       
0757      &0.0185D0,0.0135D0,0.025D0,0.0004D0,0.0007D0,0.0008D0,0.0014D0,    
0758      &0.0019D0,0.0025D0,0.4291D0,0.08D0,0.07D0,0.02D0,0.015D0,0.005D0,  
0759      &1D0,0.3D0,0.15D0,0.16D0,0.08D0,0.13D0,0.06D0,0.08D0,0.04D0,       
0760      &0.02D0,0.055D0,2*0.005D0,0.008D0,0.012D0,0.02D0,0.055D0,          
0761      &2*0.005D0,0.008D0,0.012D0,0.01D0,0.03D0,0.0035D0,0.011D0,         
0762      &0.0055D0,0.0042D0,0.009D0,0.018D0,0.015D0,0.0185D0,0.0135D0,      
0763      &0.025D0,0.0004D0,0.0007D0,0.0008D0,0.0014D0,0.0019D0,0.0025D0,    
0764      &0.4291D0,0.08D0,0.07D0,0.02D0,0.015D0,0.005D0,1D0,0.3D0,0.15D0,   
0765      &0.16D0,0.08D0,0.13D0,0.06D0,0.08D0,0.04D0,0.02D0,0.055D0,         
0766      &2*0.005D0,0.008D0,0.012D0,0.02D0,0.055D0,2*0.005D0,0.008D0,       
0767      &0.012D0,0.01D0,0.03D0,0.0035D0,0.011D0,0.0055D0,0.0042D0,0.009D0, 
0768      &0.018D0,0.015D0,0.0185D0,0.0135D0,0.025D0,2*0.0002D0,0.0007D0,    
0769      &2*0.0004D0,0.0014D0,0.001D0,0.0009D0,0.0025D0,0.4291D0,0.08D0,    
0770      &0.07D0,0.02D0,0.015D0,0.005D0,1D0,2*0.3D0,2*0.2D0,0.047D0/        
0771       DATA (BRAT(I)  ,I= 998,1188)/0.122D0,0.006D0,0.012D0,0.035D0,     
0772      &0.012D0,0.035D0,0.003D0,0.007D0,0.15D0,0.037D0,0.008D0,0.002D0,   
0773      &0.05D0,0.015D0,0.003D0,0.001D0,0.014D0,0.042D0,0.014D0,0.042D0,   
0774      &0.24D0,0.065D0,0.012D0,0.003D0,0.001D0,0.002D0,0.001D0,0.002D0,   
0775      &0.014D0,0.003D0,1D0,2*0.3D0,2*0.2D0,1D0,0.0252D0,0.0248D0,        
0776      &0.0267D0,0.015D0,0.045D0,0.015D0,0.045D0,0.7743D0,0.029D0,0.22D0, 
0777      &0.78D0,1D0,0.331D0,0.663D0,0.006D0,0.663D0,0.331D0,0.006D0,1D0,   
0778      &0.999D0,0.001D0,0.88D0,2*0.06D0,0.639D0,0.358D0,0.002D0,0.001D0,  
0779      &1D0,0.88D0,2*0.06D0,0.516D0,0.483D0,0.001D0,0.88D0,2*0.06D0,      
0780      &0.9988D0,0.0001D0,0.0006D0,0.0004D0,0.0001D0,0.667D0,0.333D0,     
0781      &0.9954D0,0.0011D0,0.0035D0,0.333D0,0.667D0,0.676D0,0.234D0,       
0782      &0.085D0,0.005D0,2*1D0,0.018D0,2*0.005D0,0.003D0,0.002D0,          
0783      &2*0.006D0,0.018D0,2*0.005D0,0.003D0,0.002D0,2*0.006D0,0.0066D0,   
0784      &0.025D0,0.016D0,0.0088D0,2*0.005D0,0.0058D0,0.005D0,0.0055D0,     
0785      &4*0.004D0,2*0.002D0,2*0.004D0,0.003D0,0.002D0,2*0.003D0,          
0786      &3*0.002D0,2*0.001D0,0.002D0,2*0.001D0,2*0.002D0,0.0013D0,         
0787      &0.0018D0,5*0.001D0,4*0.003D0,2*0.005D0,2*0.002D0,2*0.001D0,       
0788      &2*0.002D0,2*0.001D0,0.2432D0,0.057D0,2*0.035D0,0.15D0,2*0.075D0,  
0789      &0.03D0,2*0.015D0,2*0.08D0,0.76D0,0.08D0,4*1D0,2*0.08D0,0.76D0,    
0790      &0.08D0,1D0,2*0.5D0,1D0,2*0.5D0,2*0.08D0,0.76D0,0.08D0,1D0/        
0791       DATA (BRAT(I)  ,I=1189,1381)/2*0.08D0,0.76D0,3*0.08D0,0.76D0,     
0792      &3*0.08D0,0.76D0,3*0.08D0,0.76D0,3*0.08D0,0.76D0,3*0.08D0,0.76D0,  
0793      &3*0.08D0,0.76D0,0.08D0,2*1D0,2*0.105D0,0.04D0,0.0077D0,0.02D0,    
0794      &0.0235D0,0.0285D0,0.0435D0,0.0011D0,0.0022D0,0.0044D0,0.4291D0,   
0795      &0.08D0,0.07D0,0.02D0,0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,      
0796      &0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,      
0797      &0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,4*1D0,2*0.105D0,0.04D0,      
0798      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,0.04D0,      
0799      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,4*1D0,2*0.105D0,       
0800      &0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,1D0,2*0.105D0,  
0801      &0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,      
0802      &0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,      
0803      &0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,      
0804      &0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,      
0805      &0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,      
0806      &0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,      
0807      &0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,      
0808      &0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,      
0809      &0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,      
0810      &0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0/      
0811       DATA (BRAT(I)  ,I=1382,1582)/0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,   
0812      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,      
0813      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,      
0814      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,      
0815      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,      
0816      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,      
0817      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,      
0818      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,      
0819      &0.015D0,0.005D0,4*1D0,0.52D0,0.26D0,0.11D0,2*0.055D0,0.333D0,     
0820      &0.334D0,0.333D0,0.667D0,0.333D0,0.28D0,0.14D0,0.313D0,0.157D0,    
0821      &0.11D0,0.667D0,0.333D0,0.28D0,0.14D0,0.313D0,0.157D0,0.11D0,      
0822      &0.36D0,0.18D0,0.03D0,2*0.015D0,2*0.2D0,4*0.25D0,0.667D0,0.333D0,  
0823      &0.667D0,0.333D0,0.667D0,0.333D0,0.667D0,0.333D0,4*0.5D0,0.007D0,  
0824      &0.993D0,1D0,0.667D0,0.333D0,0.667D0,0.333D0,0.667D0,0.333D0,      
0825      &0.667D0,0.333D0,8*0.5D0,0.02D0,0.98D0,1D0,4*0.5D0,3*0.146D0,      
0826      &3*0.05D0,0.15D0,2*0.05D0,4*0.024D0,0.066D0,0.667D0,0.333D0,       
0827      &0.667D0,0.333D0,4*0.25D0,0.667D0,0.333D0,0.667D0,0.333D0,2*0.5D0, 
0828      &0.273D0,0.727D0,0.667D0,0.333D0,0.667D0,0.333D0,4*0.5D0,0.35D0,   
0829      &0.65D0,2*0.0083D0,0.1866D0,0.324D0,0.184D0,0.027D0,0.001D0,       
0830      &0.093D0,0.087D0,0.078D0,0.0028D0,3*0.014D0,0.008D0,0.024D0/       
0831       DATA (BRAT(I)  ,I=1583,4150)/0.008D0,0.024D0,0.425D0,0.02D0,      
0832      &0.185D0,0.088D0,0.043D0,0.067D0,0.066D0,2404*0D0,0.024396D0,      
0833      &0.045285D0,0.83119D0,2*0D0,0.000349D0,0.09878D0,0D0,0.019884D0,   
0834      &0.02341D0,0.362776D0,0.550787D0,2*0D0,0.000152D0,0.042991D0,      
0835      &0.013695D0,0.025421D0,0.466595D0,2*0D0,0.000196D0,0.055451D0,     
0836      &0.438642D0,0.445781D0,0D0,0.554219D0,4*0.00335D0,0.522257D0,      
0837      &0.464343D0,6*0D0,1D0,6*0D0,1D0,4*0.013853D0,0.562703D0,           
0838      &0.376702D0,0.00518D0,4*0.006254D0,0.974985D0,7*0D0,4*0.148299D0,  
0839      &0.015351D0,0D0,0.182109D0,0.167099D0,0.042247D0,0.850973D0,       
0840      &0.005411D0,0.045025D0,0.098591D0,0.849898D0,0.021617D0,           
0841      &0.030018D0,0.098466D0,0.294448D0,0.10945D0,0.596102D0,0.389906D0, 
0842      &0.610094D0,3*0.0633D0,0.063299D0,0.063295D0,0.056281D0,2*0D0,     
0843      &6*0.020495D0,2*0D0,0.327919D0,0.04099D0,0.045236D0,0.090112D0,    
0844      &0.19874D0,0.010204D0,0.000003D0,0.010205D0,0.198356D0,0.000151D0, 
0845      &0.000006D0,0.000367D0,0.081967D0,0.19874D0,0.010204D0,0.000003D0, 
0846      &0.010205D0,0.198356D0,0.000151D0,0.000006D0,0.000367D0,           
0847      &0.081967D0,4*0D0,0.198776D0,0.010206D0,0.000003D0,0.010207D0,     
0848      &0.19839D0,0.000151D0,0.000006D0,0.000367D0,0.081893D0,0.198776D0, 
0849      &0.010206D0,0.000003D0,0.010207D0,0.19839D0,0.000151D0,0.000006D0, 
0850      &0.000367D0,0.081893D0,4*0D0,0.199344D0,0.010234D0,0.000003D0/     
0851       DATA (BRAT(I)  ,I=4151,4281)/0.010236D0,0.198928D0,0.000149D0,    
0852      &0.000006D0,0.000368D0,0.080733D0,0.199344D0,0.010234D0,           
0853      &0.000003D0,0.010236D0,0.198928D0,0.000149D0,0.000006D0,           
0854      &0.000368D0,0.080733D0,4*0D0,0.184738D0,0.104588D0,0.184738D0,     
0855      &0.104587D0,0.184731D0,0.09582D0,0.022902D0,0.008429D0,0.015602D0, 
0856      &0.022902D0,0.008429D0,0.015602D0,0.022902D0,0.008429D0,           
0857      &0.015602D0,0.28959D0,0.01487D0,0.000008D0,0.01487D0,0.289061D0,   
0858      &0.000492D0,0.000009D0,0.000536D0,0.27911D0,2*0.037151D0,          
0859      &0.03715D0,0.090266D0,2*0.001805D0,0.090266D0,0.001805D0,          
0860      &0.812263D0,0.00179D0,0.090428D0,0.001809D0,0.001808D0,0.090428D0, 
0861      &0.001808D0,0.81372D0,0D0,6*1D0,0.095602D0,2*0.338272D0,           
0862      &0.156896D0,0.019193D0,0.017993D0,0.001168D0,0.001462D0,           
0863      &0.009608D0,0.003306D0,0.002132D0,0.003127D0,0.002132D0,           
0864      &0.003127D0,0.00213D0,3*0D0,0.001411D0,0.00045D0,0.001411D0,       
0865      &0.00045D0,0.001411D0,0.00045D0,2*0D0,0.097996D0,0.399787D0,       
0866      &0.262464D0,0.185427D0,0.022683D0,0.007648D0,0.004259D0,           
0867      &0.005925D0,0.000304D0,2*0D0,0.000304D0,0.005914D0,0.000002D0,     
0868      &2*0D0,0.000011D0,0.001258D0,5*0D0,3*0.002005D0,0D0,0.272178D0,    
0869      &0.022112D0,0.255165D0,0.015534D0,2*0.108965D0,0.031557D0,         
0870      &0.005562D0,0.044965D0,0.004674D0,0.007637D0,0.020597D0/           
0871       DATA (BRAT(I)  ,I=4282,8000)/0.007636D0,0.020595D0,0.007616D0,    
0872      &3*0D0,0.017298D0,0.004782D0,0.017298D0,0.004782D0,0.017297D0,     
0873      &0.004782D0,2*0D0,0.055332D0,2*0.319757D0,0.121576D0,2*0.001556D0, 
0874      &4*0D0,0.0277D0,0.021481D0,0.027699D0,0.021477D0,0.027658D0,3*0D0, 
0875      &0.006071D0,0.01208D0,0.006071D0,0.01208D0,0.006069D0,0.01208D0,   
0876      &2*0D0,0.035891D0,0.209476D0,0.129084D0,0.286631D0,0.10742D0,      
0877      &0.109486D0,4*0D0,0.035282D0,0.001812D0,2*0D0,0.001812D0,          
0878      &0.035215D0,0.000021D0,0D0,0.000001D0,0.000065D0,0.011965D0,5*0D0, 
0879      &2*0.011947D0,0.011946D0,0D0,
0880      &649*0.D0,
0881 C....UED
0882      &0.001D0,0.999D0,0.001D0,0.999D0,0.001D0,0.999D0,
0883      &0.001D0,0.999D0,0.001D0,0.999D0,0.001D0,0.999D0, 
0884      &0.33D0,0.66D0,0.01D0,0.33D0,0.66D0,0.01D0,0.33D0,0.66D0,0.01D0,
0885      &0.33D0,0.66D0,0.01D0,0.98D0,0.D0,0.02D0,0.33D0,0.66D0,0.01D0,
0886      &9*1.D0,              
0887      &24*0.0416667,        
0888      &1.,                  
0889      &3*0.D0,6*0.08333D0, 
0890      &3*0.D0,6*0.08333D0,
0891      &6*0.166667D0,        
0892      &2912*0.D0/
0893       DATA (KFDP(I,1),I=   1, 377)/21,22,23,4*-24,25,21,22,23,4*24,25,  
0894      &21,22,23,4*-24,25,21,22,23,4*24,25,21,22,23,4*-24,25,21,22,23,    
0895      &4*24,25,37,1000022,1000023,1000025,1000035,1000021,1000039,21,22, 
0896      &23,4*-24,25,2*-37,21,22,23,4*24,25,2*37,22,23,-24,25,23,24,-12,   
0897      &22,23,-24,25,23,24,-12,-14,48*16,22,23,-24,25,23,24,22,23,-24,25, 
0898      &-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,  
0899      &3,4,5,6,7,8,11,12,13,14,15,16,17,18,4*-1,4*-3,4*-5,4*-7,-11,-13,  
0900      &-15,-17,1,2,3,4,5,6,7,8,11,13,15,17,21,2*22,23,24,1000022,        
0901      &2*1000023,3*1000025,4*1000035,2*1000024,2*1000037,1000001,        
0902      &2000001,1000001,-1000001,1000002,2000002,1000002,-1000002,        
0903      &1000003,2000003,1000003,-1000003,1000004,2000004,1000004,         
0904      &-1000004,1000005,2000005,1000005,-1000005,1000006,2000006,        
0905      &1000006,-1000006,1000011,2000011,1000011,-1000011,1000012,        
0906      &2000012,1000012,-1000012,1000013,2000013,1000013,-1000013,        
0907      &1000014,2000014,1000014,-1000014,1000015,2000015,1000015,         
0908      &-1000015,1000016,2000016,1000016,-1000016,1,2,3,4,5,6,7,8,11,12,  
0909      &13,14,15,16,17,18,24,37,2*23,25,35,4*-1,4*-3,4*-5,4*-7,-11,-13,   
0910      &-15,-17,3*24,1,2,3,4,5,6,7,8,11,13,15,17,21,2*22,23,24,23,25,24,  
0911      &37,23,25,36,1000022,2*1000023,3*1000025,4*1000035,2*1000024,      
0912      &2*1000037,1000001,2000001,1000001,-1000001,1000002,2000002/       
0913       DATA (KFDP(I,1),I= 378, 580)/1000002,-1000002,1000003,2000003,    
0914      &1000003,-1000003,1000004,2000004,1000004,-1000004,1000005,        
0915      &2000005,1000005,-1000005,1000006,2000006,1000006,-1000006,        
0916      &1000011,2000011,1000011,-1000011,1000012,2000012,1000012,         
0917      &-1000012,1000013,2000013,1000013,-1000013,1000014,2000014,        
0918      &1000014,-1000014,1000015,2000015,1000015,-1000015,1000016,        
0919      &2000016,1000016,-1000016,1,2,3,4,5,6,7,8,11,13,15,17,21,2*22,23,  
0920      &24,23,25,24,37,1000022,2*1000023,3*1000025,4*1000035,2*1000024,   
0921      &2*1000037,1000001,2000001,1000001,-1000001,1000002,2000002,       
0922      &1000002,-1000002,1000003,2000003,1000003,-1000003,1000004,        
0923      &2000004,1000004,-1000004,1000005,2000005,1000005,-1000005,        
0924      &1000006,2000006,1000006,-1000006,1000011,2000011,1000011,         
0925      &-1000011,1000012,2000012,1000012,-1000012,1000013,2000013,        
0926      &1000013,-1000013,1000014,2000014,1000014,-1000014,1000015,        
0927      &2000015,1000015,-1000015,1000016,2000016,1000016,-1000016,-1,-3,  
0928      &-5,-7,-11,-13,-15,-17,24,2*1000022,2*1000023,2*1000025,2*1000035, 
0929      &1000006,2000006,1000006,2000006,-1000001,-1000003,-1000011,       
0930      &-1000013,-1000015,-2000015,1,2,3,4,5,6,11,13,15,2,82,-11,-13,2*2, 
0931      &-12,-14,-16,2*-2,2*-4,-2,-4,2*22,211,111,221,13,11,213,-213,221,  
0932      &223,321,130,310,111,331,111,211,-12,12,-14,14,211,111,22,-13,-11/ 
0933       DATA (KFDP(I,1),I= 581, 992)/2*211,213,113,221,223,321,211,331,   
0934      &22,111,211,2*22,211,22,111,211,22,211,221,111,11,211,111,2*211,   
0935      &321,130,310,221,111,211,111,130,310,321,2*311,321,311,323,313,    
0936      &323,313,321,3*311,-13,3*211,12,14,311,2*321,311,321,313,323,313,  
0937      &323,311,4*321,211,111,3*22,111,321,130,-213,113,213,211,22,111,   
0938      &11,13,211,321,130,310,221,211,111,11*-11,11*-13,-311,-313,-311,   
0939      &-313,-20313,2*-311,-313,-311,-313,2*111,2*221,2*331,2*113,2*223,  
0940      &2*333,-311,-313,2*-321,211,-311,-321,333,-311,-313,-321,211,      
0941      &2*-321,2*-311,-321,211,113,421,2*411,421,411,423,413,423,413,421, 
0942      &411,8*-11,8*-13,-321,-323,-321,-323,-311,2*-313,-311,-313,2*-311, 
0943      &-321,-10323,-321,-323,-321,-311,2*-313,211,111,333,3*-321,-311,   
0944      &-313,-321,-313,310,333,211,2*-321,-311,-313,-311,211,-321,3*-311, 
0945      &211,113,321,2*421,411,421,413,423,413,423,411,421,-15,5*-11,      
0946      &5*-13,221,331,333,221,331,333,10221,211,213,211,213,321,323,321,  
0947      &323,2212,221,331,333,221,2*2,2*431,421,411,423,413,82,11,13,82,   
0948      &443,82,6*12,6*14,2*16,3*-411,3*-413,2*-411,2*-413,2*441,2*443,    
0949      &2*20443,2*2,2*4,2,4,511,521,511,523,513,523,513,521,511,6*12,     
0950      &6*14,2*16,3*-421,3*-423,2*-421,2*-423,2*441,2*443,2*20443,2*2,    
0951      &2*4,2,4,521,511,521,513,523,513,523,511,521,6*12,6*14,2*16,       
0952      &3*-431,3*-433,2*-431,2*-433,3*441,3*443,3*20443,2*2,2*4,2,4,531/  
0953       DATA (KFDP(I,1),I= 993,1402)/521,511,523,513,16,2*4,2*12,2*14,    
0954      &2*16,4*2,4*4,2*-11,2*-13,2*-1,2*-3,2*-11,2*-13,2*-1,541,511,521,  
0955      &513,523,21,11,13,15,1,2,3,4,21,22,553,21,2112,2212,2*2112,2212,   
0956      &2112,2*2212,2112,-12,3122,3212,3112,2212,2*2112,-12,2*3122,3222,  
0957      &3112,2212,2112,2212,3122,3222,3212,3122,3112,-12,-14,-12,3322,    
0958      &3312,2*3122,3212,3322,3312,3122,3322,3312,-12,2*4122,7*-11,7*-13, 
0959      &2*2224,2*2212,2*2214,2*3122,2*3212,2*3214,5*3222,4*3224,2*3322,   
0960      &3324,2*2224,7*2212,5*2214,2*2112,2*2114,2*3122,2*3212,2*3214,     
0961      &2*3222,2*3224,4*2,3,2*2,1,2*2,-11,-13,2*2,4*4122,-11,-13,2*2,     
0962      &3*4132,3*4232,-11,-13,2*2,4332,-11,-13,2*2,-11,-13,2*2,-11,-13,   
0963      &2*2,-11,-13,2*2,-11,-13,2*2,-11,-13,2*2,-11,-13,2*2,2*5122,-12,   
0964      &-14,-16,5*4122,441,443,20443,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,    
0965      &2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,4*5122,-12,-14,-16,2*-2,   
0966      &2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,2*5132,2*5232,-12,-14,-16, 
0967      &2*-2,2*-4,-2,-4,5332,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,     
0968      &2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,     
0969      &2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,  
0970      &-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,   
0971      &-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,  
0972      &2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2/     
0973       DATA (KFDP(I,1),I=1403,1713)/2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2, 
0974      &-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,   
0975      &-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,221,223,221,  
0976      &223,211,111,321,130,310,213,113,-213,321,311,321,311,323,313,     
0977      &2*311,321,311,321,313,323,321,211,111,321,130,310,2*211,313,-313, 
0978      &323,-323,421,411,423,413,411,421,413,423,411,421,423,413,443,     
0979      &2*82,521,511,523,513,511,521,513,523,521,511,523,513,511,521,513, 
0980      &523,553,2*21,213,-213,113,213,10211,10111,-10211,2*221,213,2*113, 
0981      &-213,2*321,2*311,113,323,2*313,323,313,-313,323,-323,423,2*413,   
0982      &2*423,413,443,82,523,2*513,2*523,2*513,523,553,21,11,13,82,4*443, 
0983      &10441,20443,445,441,11,13,15,1,2,3,4,21,22,2*553,10551,20553,555, 
0984      &1000039,-1000024,-1000037,1000022,1000023,1000025,1000035,        
0985      &1000002,2000002,1000002,2000002,1000021,3*-12,3*-14,3*-16,12,11,  
0986      &12,11,12,11,14,13,14,13,14,13,16,15,16,15,16,15,2*-2,2*-4,2*-6,   
0987      &1000039,1000024,1000037,1000022,1000023,1000025,1000035,1000001,  
0988      &2000001,1000001,2000001,1000021,3*-11,3*-13,3*-15,2*-1,-3,        
0989      &1000039,-1000024,-1000037,1000022,1000023,1000025,1000035,        
0990      &1000004,2000004,1000004,2000004,1000021,3*-12,3*-14,3*-16,12,11,  
0991      &12,11,12,11,14,13,14,13,14,13,16,15,16,15,16,15,2*-2,2*-4,2*-6,   
0992      &1000039,1000024,1000037,1000022,1000023,1000025,1000035,1000003/  
0993       DATA (KFDP(I,1),I=1714,1984)/2000003,1000003,2000003,1000021,     
0994      &3*-11,3*-13,3*-15,2*-1,-3,1000039,-1000024,-1000037,1000022,      
0995      &1000023,1000025,1000035,1000006,2000006,1000006,2000006,1000021,  
0996      &3*-12,3*-14,3*-16,12,11,12,11,12,11,14,13,14,13,14,13,16,15,16,   
0997      &15,16,15,2*-2,2*-4,2*-6,1000039,1000024,1000037,1000022,1000023,  
0998      &1000025,1000035,1000005,2000005,1000005,2000005,1000021,1000022,  
0999      &1000016,-1000015,3*-11,3*-13,3*-15,2*-1,-3,1000039,-1000024,      
1000      &-1000037,1000022,1000023,1000025,1000035,1000012,2000012,1000012, 
1001      &2*12,2*14,2*16,3*-14,3*-16,3*-2,3*-4,3*-6,1000039,1000024,        
1002      &1000037,1000022,1000023,1000025,1000035,1000011,2000011,1000011,  
1003      &2000011,3*-13,3*-15,3*-1,3*-3,3*-5,1000039,-1000024,-1000037,     
1004      &1000022,1000023,1000025,1000035,1000014,2000014,1000014,2000014,  
1005      &2*12,2*14,2*16,3*-12,3*-16,3*-2,3*-4,3*-6,1000039,1000024,        
1006      &1000037,1000022,1000023,1000025,1000035,1000013,2000013,1000013,  
1007      &2000013,3*-11,3*-15,3*-1,3*-3,3*-5,1000039,-1000024,-1000037,     
1008      &1000022,1000023,1000025,1000035,1000016,2000016,1000016,2000016,  
1009      &2*12,2*14,2*16,3*-12,3*-14,3*-2,3*-4,3*-6,1000039,1000024,        
1010      &1000037,1000022,1000023,1000025,1000035,1000015,2000015,1000015,  
1011      &2000015,3*-11,3*-13,3*-1,3*-3,3*-5,1000039,1000001,-1000001,      
1012      &2000001,-2000001,1000002,-1000002,2000002,-2000002,1000003/       
1013       DATA (KFDP(I,1),I=1985,2321)/-1000003,2000003,-2000003,1000004,   
1014      &-1000004,2000004,-2000004,1000005,-1000005,2000005,-2000005,      
1015      &1000006,-1000006,2000006,-2000006,6*1000022,6*1000023,6*1000025,  
1016      &6*1000035,1000024,-1000024,1000024,-1000024,1000024,-1000024,     
1017      &1000037,-1000037,1000037,-1000037,1000037,-1000037,-12,12,-11,11, 
1018      &-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,   
1019      &-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-14,14,-13,13,   
1020      &-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,   
1021      &-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-16,16,-15,15,   
1022      &-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,   
1023      &-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-2,2,-2,2,-2,2,  
1024      &-4,4,-4,4,-4,4,-6,6,-6,6,-6,6,5*1000039,4,1,-12,12,-12,12,-12,12, 
1025      &-12,12,-12,12,-12,12,-14,14,-14,14,-14,14,-14,14,-14,14,-14,14,   
1026      &-16,16,-16,16,-16,16,-16,16,-16,16,-16,16,-12,12,-11,11,-12,12,   
1027      &-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,   
1028      &-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-14,14,-13,13,-14,14,   
1029      &-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,   
1030      &-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-16,16,-15,15,-16,16,   
1031      &-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,   
1032      &-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-2,2,-2,2,-2,2,-4,4,-4/ 
1033       DATA (KFDP(I,1),I=2322,2573)/4,-4,4,-6,6,-6,6,-6,6,5*1000039,     
1034      &16*1000022,1000024,-1000024,1000024,-1000024,1000024,-1000024,    
1035      &1000024,-1000024,1000024,-1000024,1000024,-1000024,1000037,       
1036      &-1000037,1000037,-1000037,1000037,-1000037,1000037,-1000037,      
1037      &1000037,-1000037,1000037,-1000037,1000024,-1000024,1000037,       
1038      &-1000037,1000001,-1000001,2000001,-2000001,1000002,-1000002,      
1039      &2000002,-2000002,1000003,-1000003,2000003,-2000003,1000004,       
1040      &-1000004,2000004,-2000004,1000005,-1000005,2000005,-2000005,      
1041      &1000006,-1000006,2000006,-2000006,1000011,-1000011,2000011,       
1042      &-2000011,1000012,-1000012,2000012,-2000012,1000013,-1000013,      
1043      &2000013,-2000013,1000014,-1000014,2000014,-2000014,1000015,       
1044      &-1000015,2000015,-2000015,1000016,-1000016,2000016,-2000016,      
1045      &5*1000021,-12,12,-12,12,-12,12,-12,12,-12,12,-12,12,-14,14,-14,   
1046      &14,-14,14,-14,14,-14,14,-14,14,-16,16,-16,16,-16,16,-16,16,-16,   
1047      &16,-16,16,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,   
1048      &11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,   
1049      &12,-11,11,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,   
1050      &13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,   
1051      &14,-13,13,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,   
1052      &15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16/   
1053       DATA (KFDP(I,1),I=2574,2892)/16,-15,15,-2,2,-2,2,-2,2,-4,4,-4,4,  
1054      &-4,4,-6,6,-6,6,-6,6,2*1000039,6*1000022,6*1000023,6*1000025,      
1055      &6*1000035,1000022,1000023,1000025,1000035,1000002,2000002,        
1056      &-1000001,-2000001,1000004,2000004,-1000003,-2000003,1000006,      
1057      &2000006,-1000005,-2000005,1000012,2000012,-1000011,-2000011,      
1058      &1000014,2000014,-1000013,-2000013,1000016,2000016,-1000015,       
1059      &-2000015,2*1000021,-12,12,-11,-12,12,-11,-12,12,-11,-12,12,-11,   
1060      &-12,12,-11,-12,12,-11,-14,-13,-14,-13,-14,-13,-14,14,-13,-14,14,  
1061      &-13,-14,14,-13,-16,-15,-16,-15,-16,-15,-16,-15,-16,-15,-16,-15,   
1062      &-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,-12,2*-11,12, 
1063      &-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,-14,2*-13,14, 
1064      &-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-14,2*-13,14, 
1065      &-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-16,2*-15,16,-16,2*-15,16, 
1066      &-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,-16,2*-15,16, 
1067      &-16,2*-15,16,-16,2*-15,16,2,-1,2,-1,2*2,-1,2,-1,3*2,-1,2*4,-3,    
1068      &3*4,-3,2*6,5*1000039,16*1000022,16*1000023,1000024,-1000024,      
1069      &1000024,-1000024,1000024,-1000024,1000024,-1000024,1000024,       
1070      &-1000024,1000024,-1000024,1000037,-1000037,1000037,-1000037,      
1071      &1000037,-1000037,1000037,-1000037,1000037,-1000037,1000037,       
1072      &-1000037,1000024,-1000024,1000037,-1000037,1000001,-1000001/      
1073       DATA (KFDP(I,1),I=2893,3182)/2000001,-2000001,1000002,-1000002,   
1074      &2000002,-2000002,1000003,-1000003,2000003,-2000003,1000004,       
1075      &-1000004,2000004,-2000004,1000005,-1000005,2000005,-2000005,      
1076      &1000006,-1000006,2000006,-2000006,1000011,-1000011,2000011,       
1077      &-2000011,1000012,-1000012,2000012,-2000012,1000013,-1000013,      
1078      &2000013,-2000013,1000014,-1000014,2000014,-2000014,1000015,       
1079      &-1000015,2000015,-2000015,1000016,-1000016,2000016,-2000016,      
1080      &5*1000021,-12,12,-12,12,-12,12,-12,12,-12,12,-12,12,-14,14,-14,   
1081      &14,-14,14,-14,14,-14,14,-14,14,-16,16,-16,16,-16,16,-16,16,-16,   
1082      &16,-16,16,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,   
1083      &11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,   
1084      &12,-11,11,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,   
1085      &13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,   
1086      &14,-13,13,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,   
1087      &15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,   
1088      &16,-15,15,-2,2,-2,2,-2,2,-4,4,-4,4,-4,4,-6,6,-6,6,-6,6,5*1000039, 
1089      &16*1000022,16*1000023,16*1000025,1000024,-1000024,1000024,        
1090      &-1000024,1000024,-1000024,1000024,-1000024,1000024,-1000024,      
1091      &1000024,-1000024,1000037,-1000037,1000037,-1000037,1000037,       
1092      &-1000037,1000037,-1000037,1000037,-1000037,1000037,-1000037/      
1093       DATA (KFDP(I,1),I=3183,3459)/1000024,-1000024,1000037,-1000037,   
1094      &1000001,-1000001,2000001,-2000001,1000002,-1000002,2000002,       
1095      &-2000002,1000003,-1000003,2000003,-2000003,1000004,-1000004,      
1096      &2000004,-2000004,1000005,-1000005,2000005,-2000005,1000006,       
1097      &-1000006,2000006,-2000006,1000011,-1000011,2000011,-2000011,      
1098      &1000012,-1000012,2000012,-2000012,1000013,-1000013,2000013,       
1099      &-2000013,1000014,-1000014,2000014,-2000014,1000015,-1000015,      
1100      &2000015,-2000015,1000016,-1000016,2000016,-2000016,5*1000021,-12, 
1101      &12,-12,12,-12,12,-12,12,-12,12,-12,12,-14,14,-14,14,-14,14,-14,   
1102      &14,-14,14,-14,14,-16,16,-16,16,-16,16,-16,16,-16,16,-16,16,-12,   
1103      &12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,   
1104      &11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-14,   
1105      &14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,   
1106      &13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-16,   
1107      &16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,   
1108      &15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-2,2,  
1109      &-2,2,-2,2,-4,4,-4,4,-4,4,-6,6,-6,6,-6,6,2*1000039,15*1000024,     
1110      &6*1000022,6*1000023,6*1000025,6*1000035,1000022,1000023,1000025,  
1111      &1000035,1000002,2000002,-1000001,-2000001,1000004,2000004,        
1112      &-1000003,-2000003,1000006,2000006,-1000005,-2000005,1000012/      
1113       DATA (KFDP(I,1),I=3460,3782)/2000012,-1000011,-2000011,1000014,   
1114      &2000014,-1000013,-2000013,1000016,2000016,-1000015,-2000015,      
1115      &2*1000021,-12,12,-11,-12,12,-11,-12,12,-11,-12,12,-11,-12,12,-11, 
1116      &-12,12,-11,-14,14,-13,-14,14,-13,-14,14,-13,-14,14,-13,-14,14,    
1117      &-13,-14,14,-13,-16,16,-15,-16,16,-15,-16,16,-15,-16,16,-15,-16,   
1118      &16,-15,-16,16,-15,-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,-12,     
1119      &2*-11,12,-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,-12, 
1120      &2*-11,12,-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-14, 
1121      &2*-13,14,-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-16, 
1122      &2*-15,16,-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,-16, 
1123      &2*-15,16,-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,2,-1,2,-1,2*2,-1, 
1124      &2,-1,3*2,-1,2*4,-3,3*4,-3,2*6,1000039,-1000024,-1000037,1000022,  
1125      &1000023,1000025,1000035,4*1000001,1000002,2000002,1000002,        
1126      &2000002,1000021,3*-12,3*-14,3*-16,12,11,12,11,12,11,14,13,14,13,  
1127      &14,13,16,15,16,15,16,15,2*-2,2*-4,2*-6,1000039,1000024,1000037,   
1128      &1000022,1000023,1000025,1000035,4*1000002,1000001,2000001,        
1129      &1000001,2000001,1000021,3*-11,3*-13,3*-15,2*-1,-3,1000039,        
1130      &-1000024,-1000037,1000022,1000023,1000025,1000035,4*1000003,      
1131      &1000004,2000004,1000004,2000004,1000021,3*-12,3*-14,3*-16,12,11,  
1132      &12,11,12,11,14,13,14,13,14,13,16,15,16,15,16,15,2*-2,2*-4,2*-6/   
1133       DATA (KFDP(I,1),I=3783,4156)/1000039,1000024,1000037,1000022,     
1134      &1000023,1000025,1000035,4*1000004,1000003,2000003,1000003,        
1135      &2000003,1000021,3*-11,3*-13,3*-15,2*-1,-3,1000039,-1000024,       
1136      &-1000037,1000022,1000023,1000025,1000035,4*1000005,1000006,       
1137      &2000006,1000006,2000006,1000021,3*-12,3*-14,3*-16,12,11,12,11,12, 
1138      &11,14,13,14,13,14,13,16,15,16,15,16,15,2*-2,2*-4,2*-6,1000039,    
1139      &1000024,1000037,1000022,1000023,1000025,1000035,4*1000006,        
1140      &1000005,2000005,1000005,2000005,1000021,3*-11,3*-13,3*-15,2*-1,   
1141      &-3,1000039,-1000024,-1000037,1000022,1000023,1000025,1000035,     
1142      &4*1000011,1000012,2000012,1000012,2000012,2*12,2*14,2*16,3*-14,   
1143      &3*-16,3*-2,3*-4,3*-6,1000039,-1000024,-1000037,1000022,1000023,   
1144      &1000025,1000035,4*1000013,1000014,2000014,1000014,2000014,2*12,   
1145      &2*14,2*16,3*-12,3*-16,3*-2,3*-4,3*-6,1000039,-1000024,-1000037,   
1146      &1000022,1000023,1000025,1000035,4*1000015,1000016,2000016,        
1147      &1000016,2000016,2*12,2*14,2*16,3*-12,3*-14,3*-2,3*-4,3*-6,3,4,5,  
1148      &6,11,13,15,21,2*4,2,4,24,-11,-13,-15,3,4,5,6,11,13,15,21,5,6,21,  
1149      &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, 
1150      &5,6,1,2,3,4,5,6,1,2,3,4,5,6,21,3100111,3200111,21,22,23,-24,21,   
1151      &22,23,24,22,23,-24,23,24,1,2,3,4,5,6,7,8,11,12,13,14,15,16,17,18, 
1152      &21,22,23,24,9*11,9*-11,11,-11,11,-11,9*13,9*-13,13,-13,13,-13,
1153      &9*15/     
1154       DATA (KFDP(I,1),I=4157,8000)/9*-15,15,-15,15,-15,1,2,3,4,5,6,11,
1155      &12,9900012,13,14,9900014,15,16,9900016,3*-1,3*-3,3*-5,-11,-13,-15,   
1156      &3*-11,2*-13,-15,24,3*-11,2*-13,-15,9900024,3*443,3*553,2*24,      
1157      &2*3000211,2*22,2*23,22,23,1,2,3,4,5,6,7,8,11,12,13,14,15,16,17,   
1158      &18,2*24,3*3000211,2*24,4*-1,4*-3,4*-5,4*-7,-11,-13,-15,-17,22,23, 
1159      &22,23,24,3000211,24,3000211,22,23,1,2,3,4,5,6,7,8,11,12,13,14,15, 
1160      &16,17,18,2*24,-24,23,2*22,24,-24,2*23,1,2,3,4,5,6,7,8,11,12,13,   
1161      &14,15,16,17,18,2*22,23,2*24,23,22,2*24,23,4*-1,4*-3,4*-5,4*-7,    
1162      &-11,-13,-15,-17,
1163      &649*0,
1164 C...UED
1165      &5100023,5100022,5100023,5100022,5100023,5100022,
1166      &5100023,5100022,5100023,5100022,5100023,5100022, 
1167      &5100023,-5100024,5100022,5100023,5100024,5100022,
1168      &5100023,-5100024,5100022,5100023,5100024,5100022,
1169      &5100023,-5100024,5100022,5100023,5100024,5100022, 
1170      &9*5100022, 
1171      &6100001,6100002,6100003,6100004,6100005,6100006,
1172      &5100001,5100002,5100003,5100004,5100005,5100006,
1173      &-6100001,-6100002,-6100003,-6100004,-6100005,-6100006,
1174      &-5100001,-5100002,-5100003,-5100004,-5100005,-5100006, 
1175      &39, 
1176      &6100011,6100013,6100015,
1177      &5100011,5100013,5100015,
1178      %5100012,5100014,5100016,
1179      &-6100011,-6100013,-6100015,
1180      &-5100011,-5100013,-5100015,
1181      %-5100012,-5100014,-5100016,
1182      &-5100011,-5100013,-5100015,
1183      &5100012,5100014,5100016,
1184      &2912*0/
1185       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, 
1186      &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,  
1187      &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, 
1188      &13,11,13,-211,-213,-211,-213,-211,-213,-211,-213,2*-211,-321,     
1189      &-323,-321,2*-323,3*-321,4*-211,-213,-211,-213,-211,-213,-211,     
1190      &-213,-211,-213,3*-211,-213,4*-211,-323,-321,2*-211,2*-321,3*-211, 
1191      &2*15,16,15,16,15,2*17,18,17,2*18,2*17,-1,-2,-3,-4,-5,-6,-7,-8,21, 
1192      &-1,-2,-3,-4,-5,-6,-7,-8,-11,-13,-15,-17,-1,-2,-3,-4,-5,-6,-7,-8,  
1193      &-11,-12,-13,-14,-15,-16,-17,-18,2,4,6,8,2,4,6,8,2,4,6,8,2,4,6,8,  
1194      &12,14,16,18,-1,-2,-3,-4,-5,-6,-7,-8,-11,-13,-15,-17,21,22,2*23,   
1195      &-24,2*1000022,1000023,1000022,1000023,1000025,1000022,1000023,    
1196      &1000025,1000035,-1000024,-1000037,-1000024,-1000037,-1000001,     
1197      &2*-2000001,2000001,-1000002,2*-2000002,2000002,-1000003,          
1198      &2*-2000003,2000003,-1000004,2*-2000004,2000004,-1000005,          
1199      &2*-2000005,2000005,-1000006,2*-2000006,2000006,-1000011,          
1200      &2*-2000011,2000011,-1000012,2*-2000012,2000012,-1000013,          
1201      &2*-2000013,2000013,-1000014,2*-2000014,2000014,-1000015,          
1202      &2*-2000015,2000015,-1000016,2*-2000016,2000016,-1,-2,-3,-4,-5,-6, 
1203      &-7,-8,-11,-12,-13,-14,-15,-16,-17,-18,-24,-37,22,25,2*36,2,4,6,8, 
1204      &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/   
1205       DATA (KFDP(I,2),I= 340, 533)/-7,-8,-11,-13,-15,-17,21,22,2*23,    
1206      &-24,2*25,-37,-24,3*36,2*1000022,1000023,1000022,1000023,1000025,  
1207      &1000022,1000023,1000025,1000035,-1000024,-1000037,-1000024,       
1208      &-1000037,-1000001,2*-2000001,2000001,-1000002,2*-2000002,2000002, 
1209      &-1000003,2*-2000003,2000003,-1000004,2*-2000004,2000004,-1000005, 
1210      &2*-2000005,2000005,-1000006,2*-2000006,2000006,-1000011,          
1211      &2*-2000011,2000011,-1000012,2*-2000012,2000012,-1000013,          
1212      &2*-2000013,2000013,-1000014,2*-2000014,2000014,-1000015,          
1213      &2*-2000015,2000015,-1000016,2*-2000016,2000016,-1,-2,-3,-4,-5,-6, 
1214      &-7,-8,-11,-13,-15,-17,21,22,2*23,-24,2*25,-37,-24,2*1000022,      
1215      &1000023,1000022,1000023,1000025,1000022,1000023,1000025,1000035,  
1216      &-1000024,-1000037,-1000024,-1000037,-1000001,2*-2000001,2000001,  
1217      &-1000002,2*-2000002,2000002,-1000003,2*-2000003,2000003,-1000004, 
1218      &2*-2000004,2000004,-1000005,2*-2000005,2000005,-1000006,          
1219      &2*-2000006,2000006,-1000011,2*-2000011,2000011,-1000012,          
1220      &2*-2000012,2000012,-1000013,2*-2000013,2000013,-1000014,          
1221      &2*-2000014,2000014,-1000015,2*-2000015,2000015,-1000016,          
1222      &2*-2000016,2000016,2,4,6,8,12,14,16,18,25,1000024,1000037,        
1223      &1000024,1000037,1000024,1000037,1000024,1000037,2*-1000005,       
1224      &2*-2000005,1000002,1000004,1000012,1000014,2*1000016,-3,-4,-5,-6/ 
1225       DATA (KFDP(I,2),I= 534, 938)/-7,-8,-13,-15,-17,11,-82,12,14,-1,   
1226      &-3,11,13,15,1,4,3,4,1,3,22,11,-211,2*22,-13,-11,-211,211,111,211, 
1227      &-321,130,310,22,2*111,-211,11,-11,13,-13,-211,111,22,14,12,111,   
1228      &22,111,3*211,-311,22,211,22,111,-211,211,11,-211,13,22,-211,111,  
1229      &-211,22,111,-11,-211,111,2*-211,-321,130,310,221,111,-211,111,    
1230      &2*0,-211,111,22,-211,111,-211,111,-211,211,-213,113,223,221,14,   
1231      &111,211,111,-11,-13,211,111,22,211,111,211,111,2*211,213,113,223, 
1232      &221,22,-211,111,113,223,22,111,-321,310,211,111,2*-211,221,22,    
1233      &-11,-13,-211,-321,130,310,221,-211,111,11*12,11*14,2*211,2*213,   
1234      &211,20213,2*321,2*323,211,213,211,213,211,213,211,213,211,213,    
1235      &211,213,3*211,213,211,2*321,8*211,2*113,3*211,111,22,211,111,211, 
1236      &111,4*211,8*12,8*14,2*211,2*213,2*111,221,2*113,223,333,20213,    
1237      &211,2*321,323,2*311,313,-211,111,113,2*211,321,2*211,311,321,310, 
1238      &211,-211,4*211,321,4*211,113,2*211,-321,111,22,-211,111,-211,111, 
1239      &-211,211,-211,211,16,5*12,5*14,3*211,3*213,211,2*111,2*113,       
1240      &2*-311,2*-313,-2112,3*321,323,2*-1,22,111,321,311,321,311,-82,    
1241      &-11,-13,-82,22,-82,6*-11,6*-13,2*-15,211,213,20213,211,213,20213, 
1242      &431,433,431,433,311,313,311,313,311,313,-1,-4,-3,-4,-1,-3,22,     
1243      &-211,111,-211,111,-211,211,-211,211,6*-11,6*-13,2*-15,211,213,    
1244      &20213,211,213,20213,431,433,431,433,321,323,321,323,321,323,-1/   
1245       DATA (KFDP(I,2),I= 939,1352)/-4,-3,-4,-1,-3,22,211,111,211,111,   
1246      &4*211,6*-11,6*-13,2*-15,211,213,20213,211,213,20213,431,433,431,  
1247      &433,221,331,333,221,331,333,221,331,333,-1,-4,-3,-4,-1,-3,22,     
1248      &-321,-311,-321,-311,-15,-3,-1,2*-11,2*-13,2*-15,-1,-4,-3,-4,-3,   
1249      &-4,-1,-4,2*12,2*14,2,3,2,3,2*12,2*14,2,1,22,411,421,411,421,21,   
1250      &-11,-13,-15,-1,-2,-3,-4,2*21,22,21,2*-211,111,22,111,211,22,211,  
1251      &-211,11,2*-211,111,-211,111,22,11,22,111,-211,211,111,211,22,211, 
1252      &111,211,-211,22,11,13,11,-211,2*111,2*22,111,211,-321,-211,111,   
1253      &11,2*-211,7*12,7*14,-321,-323,-311,-313,-311,-313,211,213,211,    
1254      &213,211,213,111,221,331,113,223,111,221,113,223,321,323,321,-211, 
1255      &-213,111,221,331,113,223,333,10221,111,221,331,113,223,211,213,   
1256      &211,213,321,323,321,323,321,323,311,313,311,313,2*-1,-3,-1,2203,  
1257      &3201,3203,2203,2101,2103,12,14,-1,-3,2*111,2*211,12,14,-1,-3,22,  
1258      &111,2*22,111,22,12,14,-1,-3,22,12,14,-1,-3,12,14,-1,-3,12,14,-1,  
1259      &-3,12,14,-1,-3,12,14,-1,-3,12,14,-1,-3,12,14,-1,-3,2*-211,11,13,  
1260      &15,-211,-213,-20213,-431,-433,3*3122,1,4,3,4,1,3,11,13,15,1,4,3,  
1261      &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,   
1262      &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,  
1263      &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, 
1264      &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/ 
1265       DATA (KFDP(I,2),I=1353,1815)/11,13,15,1,4,3,4,1,3,11,13,15,1,4,3, 
1266      &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, 
1267      &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, 
1268      &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, 
1269      &2*111,2*211,-211,111,-321,130,310,-211,111,211,-211,111,-213,113, 
1270      &-211,111,223,211,111,213,113,211,111,223,-211,111,-321,130,310,   
1271      &2*-211,-311,311,-321,321,211,111,211,111,-211,111,-211,111,311,   
1272      &2*321,311,22,2*-82,-211,111,-211,111,211,111,211,111,-321,-311,   
1273      &-321,-311,411,421,411,421,22,2*21,-211,2*211,111,-211,111,2*211,  
1274      &111,-211,211,111,211,-321,2*-311,-321,22,-211,111,211,111,-311,   
1275      &311,-321,321,211,111,-211,111,321,311,22,-82,-211,111,211,111,    
1276      &-321,-311,411,421,22,21,-11,-13,-82,211,111,221,111,4*22,-11,-13, 
1277      &-15,-1,-2,-3,-4,2*21,211,111,3*22,1,2*2,4*1,2*-24,2*-37,2*1,3,5,  
1278      &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,   
1279      &-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,     
1280      &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,  
1281      &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, 
1282      &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, 
1283      &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,  
1284      &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/  
1285       DATA (KFDP(I,2),I=1816,2317)/11,13,11,13,15,11,13,15,1,3,5,1,3,5, 
1286      &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, 
1287      &13,2*14,4*13,2*-24,2*-37,13,15,11,15,11,13,11,13,15,11,13,15,1,3, 
1288      &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, 
1289      &5,1,3,5,15,2*16,4*15,2*-24,2*-37,13,15,11,15,11,13,11,13,15,11,   
1290      &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, 
1291      &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,  
1292      &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, 
1293      &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, 
1294      &-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, 
1295      &-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, 
1296      &-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, 
1297      &-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, 
1298      &-3,3,-1,1,-1,1,-3,3,22,23,25,35,36,-1,-3,-13,13,-13,13,-13,13,    
1299      &-15,15,-15,15,-15,15,-11,11,-11,11,-11,11,-15,15,-15,15,-15,15,   
1300      &-11,11,-11,11,-11,11,-13,13,-13,13,-13,13,-1,1,-2,2,-1,1,-2,2,-1, 
1301      &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, 
1302      &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, 
1303      &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, 
1304      &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/ 
1305       DATA (KFDP(I,2),I=2318,2770)/3,-1,1,-1,1,-3,3,-1,1,-1,1,-3,3,22,  
1306      &23,25,35,36,22,23,11,13,15,12,14,16,1,3,5,2,4,25,35,36,-24,24,11, 
1307      &-11,13,-13,15,-15,1,-1,3,-3,-24,24,11,-11,13,-13,15,-15,1,-1,3,   
1308      &-3,-37,37,-37,37,-1,1,-1,1,-2,2,-2,2,-3,3,-3,3,-4,4,-4,4,-5,5,-5, 
1309      &5,-6,6,-6,6,-11,11,-11,11,-12,12,-12,12,-13,13,-13,13,-14,14,-14, 
1310      &14,-15,15,-15,15,-16,16,-16,16,1,3,5,2,4,-13,13,-13,13,-13,13,    
1311      &-15,15,-15,15,-15,15,-11,11,-11,11,-11,11,-15,15,-15,15,-15,15,   
1312      &-11,11,-11,11,-11,11,-13,13,-13,13,-13,13,-1,1,-2,2,-1,1,-2,2,-1, 
1313      &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, 
1314      &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, 
1315      &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, 
1316      &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, 
1317      &3,-1,1,-1,1,-3,3,-1,1,-1,1,-3,3,24,37,24,-11,-13,-15,-1,-3,24,    
1318      &-11,-13,-15,-1,-3,24,-11,-13,-15,-1,-3,24,-11,-13,-15,-1,-3,4*37, 
1319      &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,   
1320      &-3,-13,14,2*-13,14,2*-13,14,-13,-15,16,2*-15,16,2*-15,16,-15,     
1321      &6*-11,-15,16,2*-15,16,2*-15,16,-15,6*-11,6*-13,-1,-2,-1,2,-1,-2,  
1322      &-1,2,-1,-2,-1,2,-3,-4,-3,4,-3,-4,-3,4,-3,-4,-3,4,-5,-6,-5,6,-5,   
1323      &-6,-5,6,-5,-6,-5,6,-1,-2,-1,2,-1,-2,-1,2,-1,-2,-1,2,-3,-4,-3,4,   
1324      &-3,-4,-3,4,-3,-4,-3,4,-5,-6,-5,6,-5,-6,-5,6,-5,-6,-5,6,-1,-2,-1/  
1325       DATA (KFDP(I,2),I=2771,3221)/2,-1,-2,-1,2,-1,-2,-1,2,-3,-4,-3,4,  
1326      &-3,-4,-3,4,-3,-4,-3,4,-5,-6,-5,6,-5,-6,-5,6,-5,-6,-5,6,2,-1,2,-1, 
1327      &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,  
1328      &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, 
1329      &25,35,36,-24,24,11,-11,13,-13,15,-15,1,-1,3,-3,-24,24,11,-11,13,  
1330      &-13,15,-15,1,-1,3,-3,-37,37,-37,37,-1,1,-1,1,-2,2,-2,2,-3,3,-3,3, 
1331      &-4,4,-4,4,-5,5,-5,5,-6,6,-6,6,-11,11,-11,11,-12,12,-12,12,-13,13, 
1332      &-13,13,-14,14,-14,14,-15,15,-15,15,-16,16,-16,16,1,3,5,2,4,-13,   
1333      &13,-13,13,-13,13,-15,15,-15,15,-15,15,-11,11,-11,11,-11,11,-15,   
1334      &15,-15,15,-15,15,-11,11,-11,11,-11,11,-13,13,-13,13,-13,13,-1,1,  
1335      &-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, 
1336      &-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, 
1337      &-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, 
1338      &-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, 
1339      &-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, 
1340      &22,23,11,13,15,12,14,16,1,3,5,2,4,25,35,36,22,23,11,13,15,12,14,  
1341      &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, 
1342      &-24,24,11,-11,13,-13,15,-15,1,-1,3,-3,-24,24,11,-11,13,-13,15,    
1343      &-15,1,-1,3,-3,-37,37,-37,37,-1,1,-1,1,-2,2,-2,2,-3,3,-3,3,-4,4,   
1344      &-4,4,-5,5,-5,5,-6,6,-6,6,-11,11,-11,11,-12,12,-12,12,-13,13,-13/  
1345       DATA (KFDP(I,2),I=3222,3669)/13,-14,14,-14,14,-15,15,-15,15,-16,  
1346      &16,-16,16,1,3,5,2,4,-13,13,-13,13,-13,13,-15,15,-15,15,-15,15,    
1347      &-11,11,-11,11,-11,11,-15,15,-15,15,-15,15,-11,11,-11,11,-11,11,   
1348      &-13,13,-13,13,-13,13,-1,1,-2,2,-1,1,-2,2,-1,1,-2,2,-3,3,-4,4,-3,  
1349      &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, 
1350      &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, 
1351      &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, 
1352      &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, 
1353      &1,-1,1,-3,3,24,37,23,11,13,15,12,14,16,1,3,5,2,4,25,35,36,24,-11, 
1354      &-13,-15,-1,-3,24,-11,-13,-15,-1,-3,24,-11,-13,-15,-1,-3,24,-11,   
1355      &-13,-15,-1,-3,4*37,2*-1,2*2,2*-3,2*4,2*-5,2*6,2*-11,2*12,2*-13,   
1356      &2*14,2*-15,2*16,-1,-3,-13,14,2*-13,14,2*-13,14,-13,-15,16,2*-15,  
1357      &16,2*-15,16,-15,-11,12,2*-11,12,2*-11,12,-11,-15,16,2*-15,16,     
1358      &2*-15,16,-15,-11,12,2*-11,12,2*-11,12,-11,-13,14,2*-13,14,2*-13,  
1359      &14,-13,-1,-2,-1,2,-1,-2,-1,2,-1,-2,-1,2,-3,-4,-3,4,-3,-4,-3,4,-3, 
1360      &-4,-3,4,-5,-6,-5,6,-5,-6,-5,6,-5,-6,-5,6,-1,-2,-1,2,-1,-2,-1,2,   
1361      &-1,-2,-1,2,-3,-4,-3,4,-3,-4,-3,4,-3,-4,-3,4,-5,-6,-5,6,-5,-6,-5,  
1362      &6,-5,-6,-5,6,-1,-2,-1,2,-1,-2,-1,2,-1,-2,-1,2,-3,-4,-3,4,-3,-4,   
1363      &-3,4,-3,-4,-3,4,-5,-6,-5,6,-5,-6,-5,6,-5,-6,-5,6,2,-1,2,-1,2*4,   
1364      &-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/     
1365       DATA (KFDP(I,2),I=3670,4183)/2*-37,2*1,3,5,1,3,5,1,3,5,1,2,3,4,5, 
1366      &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,   
1367      &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,   
1368      &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,  
1369      &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,  
1370      &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,  
1371      &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,  
1372      &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,     
1373      &2*12,4*11,23,25,35,36,2*-24,2*-37,13,15,11,15,11,13,11,13,15,11,  
1374      &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,  
1375      &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,  
1376      &23,25,35,36,2*-24,2*-37,13,15,11,15,11,13,11,13,15,11,13,15,1,3,  
1377      &5,1,3,5,1,3,5,-3,-4,-5,-6,-11,-13,-15,21,-1,-3,2*-5,5,12,14,16,   
1378      &-3,-4,-5,-6,-11,-13,-15,21,-5,-6,21,-1,-2,-3,-4,-5,-6,-1,-2,-3,   
1379      &-4,-5,-6,21,-1,-2,-3,-4,-5,-6,21,-1,-2,-3,-4,-5,-6,21,-1,-2,-3,   
1380      &-4,-5,-6,-1,-2,-3,-4,-5,-6,-1,-2,-3,-4,-5,-6,3*21,3*1,4*2,1,2*11, 
1381      &2*12,11,-1,-2,-3,-4,-5,-6,-7,-8,-11,-12,-13,-14,-15,-16,-17,-18,  
1382      &21,22,23,-24,3*-1,3*-3,3*-5,3*1,3*3,3*5,-13,13,-15,15,3*-1,3*-3,     
1383      &3*-5,3*1,3*3,3*5,-11,11,-15,15,3*-1,3*-3,3*-5,3*1,3*3,3*5,-11,11,     
1384      &-13,13,-1,-2,-3,-4,-5,-6,-11,-12,9900012,-13,-14,9900014,-15,-16/   
1385       DATA (KFDP(I,2),I=4184,8000)/9900016,2,4,6,2,4,6,2,4,6,9900012,   
1386      &9900014,9900016,-11,-13,-15,-13,2*-15,24,-11,-13,-15,-13,2*-15,   
1387      &9900024,6*21,-24,-3000211,-24,-3000211,3000111,3000221,3000111,   
1388      &3000221,2*23,-1,-2,-3,-4,-5,-6,-7,-8,-11,-12,-13,-14,-15,-16,-17, 
1389      &-18,23,3000111,23,3000111,22,3000221,22,2,4,6,8,2,4,6,8,2,4,6,8,  
1390      &2,4,6,8,12,14,16,18,2*3000111,2*3000221,-3000211,2*-24,-3000211,  
1391      &2*23,-1,-2,-3,-4,-5,-6,-7,-8,-11,-12,-13,-14,-15,-16,-17,-18,-24, 
1392      &-3000211,3000211,3000221,3000113,3000223,-3000213,3000213,        
1393      &3000113,3000223,-1,-2,-3,-4,-5,-6,-7,-8,-11,-12,-13,-14,-15,-16,  
1394      &-17,-18,24,3000211,24,3000111,3000221,3000211,3000213,3000113,    
1395      &3000223,3000213,2,4,6,8,2,4,6,8,2,4,6,8,2,4,6,8,12,14,16,18,      
1396      &649*0,
1397 C...UED     
1398      &1,1,2,2,3,3,4,4,5,5,6,6, 
1399      &1,2,1,2,1,2,3,4,3,4,3,4,5,6,5,6,5,6,
1400      &11,13,15,12,11,14,13,16,15, 
1401      &-1,-2,-3,-4,-5,-6,-1,-2,-3,-4,-5,-6,
1402      &1,2,3,4,5,6,1,2,3,4,5,6, 
1403      &22, 
1404      &-11,-13,-15,-11,-13,-15,-12,-14,-16,
1405      &11,13,15,11,13,15,12,14,16,
1406      &12,14,16,-11,-13,-15, 
1407      &2912*0/
1408       DATA (KFDP(I,3),I=   1,1021)/81*0,14,6*0,2*16,2*0,6*111,310,130,  
1409      &2*0,3*111,310,130,321,113,211,223,221,2*113,2*211,2*223,2*221,    
1410      &2*113,221,2*113,2*213,-213,113,2*111,310,130,310,130,2*310,130,   
1411      &402*0,4*3,4*4,1,4,3,2*2,0,-11,8*0,-211,5*0,2*111,211,-211,211,    
1412      &-211,10*0,111,4*0,2*111,-211,-11,11,-13,22,111,3*0,22,3*0,111,    
1413      &211,4*0,111,11*0,111,-211,6*0,-211,3*111,7*0,111,-211,5*0,2*221,  
1414      &3*0,111,5*0,111,11*0,-311,-313,-311,-321,-313,-323,111,221,331,   
1415      &113,223,-311,-313,-311,-321,-313,-323,111,221,331,113,223,22*0,   
1416      &111,113,2*211,-211,-311,211,111,3*211,-211,7*211,7*0,111,-211,    
1417      &111,-211,-321,-323,-311,-321,-313,-323,-211,-213,-321,-323,-311,  
1418      &-321,-313,-323,-211,-213,22*0,111,113,-311,2*-211,211,-211,310,   
1419      &-211,2*111,211,2*-211,-321,-211,2*211,-211,111,-211,2*211,6*0,    
1420      &111,-211,111,-211,0,221,331,333,321,311,221,331,333,321,311,20*0, 
1421      &3,13*0,-411,-413,-10413,-10411,-20413,-415,-411,-413,-10413,      
1422      &-10411,-20413,-415,-411,-413,16*0,-4,-1,-4,-3,2*-2,5*0,111,-211,  
1423      &111,-211,-421,-423,-10423,-10421,-20423,-425,-421,-423,-10423,    
1424      &-10421,-20423,-425,-421,-423,16*0,-4,-1,-4,-3,2*-2,5*0,111,-211,  
1425      &111,-211,-431,-433,-10433,-10431,-20433,-435,-431,-433,-10433,    
1426      &-10431,-20433,-435,-431,-433,19*0,-4,-1,-4,-3,2*-2,8*0,441,443,   
1427      &441,443,441,443,-4,-1,-4,-3,-4,-3,-4,-1,531,533,531,533,3,2,3,2/  
1428       DATA (KFDP(I,3),I=1022,2223)/511,513,511,513,1,2,13*0,2*21,11*0,  
1429      &2112,6*0,2212,12*0,2*3122,3212,10*0,3322,2*0,3122,3212,3214,2112, 
1430      &2114,2212,2112,3122,3212,3214,2112,2114,2212,2112,52*0,3*3,1,6*0, 
1431      &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,  
1432      &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,  
1433      &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,    
1434      &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,  
1435      &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,    
1436      &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,  
1437      &3,2*2,31*0,211,111,45*0,-211,2*111,-211,3*111,-211,111,211,30*0,  
1438      &-211,111,13*0,2*21,-211,111,199*0,2*5,210*0,-1,-3,-5,-2,-4,-6,-1, 
1439      &-3,-5,-2,-4,-6,-1,-3,-5,-2,-4,-6,-1,-3,-5,-2,-4,-6,-2,2,-4,4,-6,  
1440      &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,  
1441      &-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, 
1442      &-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, 
1443      &-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, 
1444      &-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,   
1445      &-5,5,-5,5,5*0,11,12,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11, 
1446      &-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,   
1447      &-11,13,-13,15,-15,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3/ 
1448       DATA (KFDP(I,3),I=2224,2783)/-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,  
1449      &-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, 
1450      &-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, 
1451      &-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,   
1452      &-5,5,-5,5,-3,3,-5,5,-5,5,-3,3,-5,5,-5,5,7*0,-11,-13,-15,-12,-14,  
1453      &-16,-1,-3,-5,-2,-4,5*0,-12,12,-14,14,-16,16,-2,2,-4,4,2*0,-12,12, 
1454      &-14,14,-16,16,-2,2,-4,4,52*0,-1,-3,-5,-2,-4,11,-11,13,-13,15,-15, 
1455      &11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,   
1456      &11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,1,-1,1,-1,3,-3,3,-3,5,  
1457      &-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, 
1458      &-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, 
1459      &-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, 
1460      &-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,   
1461      &-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, 
1462      &4,0,12,14,16,2,4,0,12,14,16,2,4,28*0,2,4,12,-11,11,14,-13,13,16,  
1463      &-15,15,12,-11,11,14,-13,13,16,-15,15,12,11,14,13,16,15,12,-11,11, 
1464      &14,-13,13,16,-15,15,12,11,14,13,16,15,12,11,14,13,16,15,2*2,1,-1, 
1465      &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,   
1466      &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,   
1467      &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/   
1468       DATA (KFDP(I,3),I=2784,3354)/2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3, 
1469      &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,  
1470      &-11,-13,-15,-12,-14,-16,-1,-3,-5,-2,-4,5*0,-11,-13,-15,-12,-14,   
1471      &-16,-1,-3,-5,-2,-4,5*0,-12,12,-14,14,-16,16,-2,2,-4,4,2*0,-12,12, 
1472      &-14,14,-16,16,-2,2,-4,4,52*0,-1,-3,-5,-2,-4,11,-11,13,-13,15,-15, 
1473      &11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,   
1474      &11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,1,-1,1,-1,3,-3,3,-3,5,  
1475      &-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, 
1476      &-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, 
1477      &-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, 
1478      &-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,   
1479      &-5,5,-3,3,-5,5,-5,5,-3,3,-5,5,-5,5,7*0,-11,-13,-15,-12,-14,-16,   
1480      &-1,-3,-5,-2,-4,5*0,-11,-13,-15,-12,-14,-16,-1,-3,-5,-2,-4,5*0,    
1481      &-11,-13,-15,-12,-14,-16,-1,-3,-5,-2,-4,5*0,-12,12,-14,14,-16,16,  
1482      &-2,2,-4,4,2*0,-12,12,-14,14,-16,16,-2,2,-4,4,52*0,-1,-3,-5,-2,-4, 
1483      &11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,   
1484      &11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,1, 
1485      &-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, 
1486      &-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, 
1487      &-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/ 
1488       DATA (KFDP(I,3),I=3355,8000)/-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,  
1489      &-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,   
1490      &-5,5,-3,3,-5,5,-5,5,3*0,-11,-13,-15,-12,-14,-16,-1,-3,-5,-2,-4,   
1491      &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,    
1492      &28*0,2,4,12,-11,11,14,-13,13,16,-15,15,12,-11,11,14,-13,13,16,    
1493      &-15,15,12,-11,11,14,-13,13,16,-15,15,12,-11,11,14,-13,13,16,-15,  
1494      &15,12,-11,11,14,-13,13,16,-15,15,12,-11,11,14,-13,13,16,-15,15,   
1495      &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,   
1496      &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,   
1497      &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,   
1498      &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,  
1499      &1,3,-3,5,-5,1,3,5,-5,1,5,-5,1,3,5,-5,1,3,351*0,-5,95*0,2,4,6,2,4, 
1500      &6,2,4,6,-2,-4,-6,-2,-4,-6,-2,-4,-6,2*9900014,2*9900016,2,4,6,2,4, 
1501      &6,2,4,6,-2,-4,-6,-2,-4,-6,-2,-4,-6,2*9900012,2*9900016,2,4,6,2,4, 
1502      &6,2,4,6,-2,-4,-6,-2,-4,-6,-2,-4,-6,2*9900012,2*9900014,3831*0/    
1503       DATA (KFDP(I,4),I=   1,8000)/94*0,4*111,6*0,111,2*0,-211,0,-211,  
1504      &3*0,111,2*-211,0,111,0,2*111,113,221,2*111,-213,-211,211,113,     
1505      &6*111,310,2*130,402*0,13*81,41*0,-11,10*0,111,-211,4*0,111,62*0,  
1506      &111,211,111,211,7*0,111,211,111,211,35*0,2*-211,2*111,211,111,    
1507      &-211,2*211,2*-211,13*0,-211,111,-211,111,4*0,-211,111,-211,111,   
1508      &34*0,111,-211,3*111,3*-211,2*111,3*-211,14*0,-321,-311,3*0,-321,  
1509      &-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,    
1510      &2*-5,67*0,-211,111,5*0,-211,111,52*0,2101,2103,2*2101,6*0,4*81,   
1511      &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, 
1512      &162*81,31*0,-211,111,6516*0/                                      
1513       DATA (KFDP(I,5),I=   1,8000)/96*0,2*111,17*0,111,7*0,2*111,0,     
1514      &3*111,0,111,597*0,-211,2*111,-211,111,-211,111,65*0,111,-211,     
1515      &3*111,-211,111,7193*0/                                            
1516  
1517 C...PYDAT4, with particle names (character strings).
1518       DATA (CHAF(I,1),I=   1, 202)/'d','u','s','c','b','t','b''','t''', 
1519      &2*' ','e-','nu_e','mu-','nu_mu','tau-','nu_tau','tau''-',         
1520      &'nu''_tau',2*' ','g','gamma','Z0','W+','h0',6*' ','Z''0','Z"0',   
1521      &'W''+','H0','A0','H+',' ','Graviton',' ','R0','LQ_ue',38*' ',     
1522      &'specflav','rndmflav','phasespa','c-hadron','b-hadron',2*' ',     
1523      &'junction',' ','system','cluster','string','indep.','CMshower',   
1524      &'SPHEaxis','THRUaxis','CLUSjet','CELLjet','table',' ','reggeon',  
1525      &'pi0','rho0','a_20','K_L0','pi+','rho+','a_2+','eta','omega',     
1526      &'f_2','K_S0','K0','K*0','K*_20','K+','K*+','K*_2+','eta''','phi', 
1527      &'f''_2','D+','D*+','D*_2+','D0','D*0','D*_20','D_s+','D*_s+',     
1528      &'D*_2s+','eta_c','J/psi','chi_2c','B0','B*0','B*_20','B+','B*+',  
1529      &'B*_2+','B_s0','B*_s0','B*_2s0','B_c+','B*_c+','B*_2c+','eta_b',  
1530      &'Upsilon','chi_2b','pomeron','dd_1','Delta-','ud_0','ud_1','n0',  
1531      &'Delta0','uu_1','p+','Delta+','Delta++','sd_0','sd_1','Sigma-',   
1532      &'Sigma*-','Lambda0','su_0','su_1','Sigma0','Sigma*0','Sigma+',    
1533      &'Sigma*+','ss_1','Xi-','Xi*-','Xi0','Xi*0','Omega-','cd_0',       
1534      &'cd_1','Sigma_c0','Sigma*_c0','Lambda_c+','Xi_c0','cu_0','cu_1',  
1535      &'Sigma_c+','Sigma*_c+','Sigma_c++','Sigma*_c++','Xi_c+','cs_0',   
1536      &'cs_1','Xi''_c0','Xi*_c0','Xi''_c+','Xi*_c+','Omega_c0',          
1537      &'Omega*_c0','cc_1','Xi_cc+','Xi*_cc+','Xi_cc++','Xi*_cc++'/       
1538       DATA (CHAF(I,1),I= 203, 332)/'Omega_cc+','Omega*_cc+',            
1539      &'Omega*_ccc++','bd_0','bd_1','Sigma_b-','Sigma*_b-','Lambda_b0',  
1540      &'Xi_b-','Xi_bc0','bu_0','bu_1','Sigma_b0','Sigma*_b0','Sigma_b+', 
1541      &'Sigma*_b+','Xi_b0','Xi_bc+','bs_0','bs_1','Xi''_b-','Xi*_b-',    
1542      &'Xi''_b0','Xi*_b0','Omega_b-','Omega*_b-','Omega_bc0','bc_0',     
1543      &'bc_1','Xi''_bc0','Xi*_bc0','Xi''_bc+','Xi*_bc+','Omega''_bc0',   
1544      &'Omega*_bc0','Omega_bcc+','Omega*_bcc+','bb_1','Xi_bb-',          
1545      &'Xi*_bb-','Xi_bb0','Xi*_bb0','Omega_bb-','Omega*_bb-',            
1546      &'Omega_bbc0','Omega*_bbc0','Omega*_bbb-','a_00','b_10','a_0+',    
1547      &'b_1+','f_0','h_1','K*_00','K_10','K*_0+','K_1+','f''_0','h''_1', 
1548      &'D*_0+','D_1+','D*_00','D_10','D*_0s+','D_1s+','chi_0c','h_1c',   
1549      &'B*_00','B_10','B*_0+','B_1+','B*_0s0','B_1s0','B*_0c+','B_1c+',  
1550      &'chi_0b','h_1b','a_10','a_1+','f_1','K*_10','K*_1+','f''_1',      
1551      &'D*_1+','D*_10','D*_1s+','chi_1c','B*_10','B*_1+','B*_1s0',       
1552      &'B*_1c+','chi_1b','psi''','Upsilon''','~d_L','~u_L','~s_L',       
1553      &'~c_L','~b_1','~t_1','~e_L-','~nu_eL','~mu_L-','~nu_muL',         
1554      &'~tau_1-','~nu_tauL','~g','~chi_10','~chi_20','~chi_1+',          
1555      &'~chi_30','~chi_40','~chi_2+','~Gravitino','~d_R','~u_R','~s_R',  
1556      &'~c_R','~b_2','~t_2','~e_R-','~nu_eR','~mu_R-','~nu_muR',         
1557      &'~tau_2-','~nu_tauR','pi_tc0','pi_tc+','pi''_tc0','eta_tc0'/      
1558       DATA (CHAF(I,1),I= 333, 500)/'rho_tc0','rho_tc+','omega_tc',      
1559      &'V8_tc','pi_22_1_tc','pi_22_8_tc','rho_11_tc','rho_12_tc',        
1560      &'rho_21_tc','rho_22_tc','d*','u*','e*-','nu*_e0','Graviton*',     
1561      &'nu_Re','nu_Rmu','nu_Rtau','Z_R0','W_R+','H_L++','H_R++',         
1562      &'rho_diff0','pi_diffr+','omega_di','phi_diff','J/psi_di',         
1563      &'n_diffr0','p_diffr+','cc~[3S18]','cc~[1S08]','cc~[3P08]',        
1564      &'bb~[3S18]','bb~[1S08]','bb~[3P08]','a_tc0','a_tc+',
1565      &81*' ',
1566 C...UED    
1567      &'d*_S','u*_S','s*_S','c*_S','b*_S','t*_S',
1568      &'d*_D','u*_D','s*_D','c*_D','b*_D','t*_D',
1569      &'e*_S-','mu*_S-','tau*_S-',
1570      &'nu*_eD','e*_D-','nu*_muD','mu*_D-','nu*_tauD','tau*_D-',
1571      &'g*','gamma*','Z*0','W*+',25*' '/               
1572       DATA (CHAF(I,2),I=   1, 205)/'dbar','ubar','sbar','cbar','bbar',  
1573      &'tbar','b''bar','t''bar',2*' ','e+','nu_ebar','mu+','nu_mubar',   
1574      &'tau+','nu_taubar','tau''+','nu''_taubar',5*' ','W-',9*' ',       
1575      &'W''-',2*' ','H-',3*' ','Rbar0','LQ_uebar',39*' ','rndmflavbar',  
1576      &' ','c-hadronbar','b-hadronbar',20*' ','pi-','rho-','a_2-',4*' ', 
1577      &'Kbar0','K*bar0','K*_2bar0','K-','K*-','K*_2-',3*' ','D-','D*-',  
1578      &'D*_2-','Dbar0','D*bar0','D*_2bar0','D_s-','D*_s-','D*_2s-',      
1579      &3*' ','Bbar0','B*bar0','B*_2bar0','B-','B*-','B*_2-','B_sbar0',   
1580      &'B*_sbar0','B*_2sbar0','B_c-','B*_c-','B*_2c-',4*' ','dd_1bar',   
1581      &'Deltabar+','ud_0bar','ud_1bar','nbar0','Deltabar0','uu_1bar',    
1582      &'pbar-','Deltabar-','Deltabar--','sd_0bar','sd_1bar','Sigmabar+', 
1583      &'Sigma*bar+','Lambdabar0','su_0bar','su_1bar','Sigmabar0',        
1584      &'Sigma*bar0','Sigmabar-','Sigma*bar-','ss_1bar','Xibar+',         
1585      &'Xi*bar+','Xibar0','Xi*bar0','Omegabar+','cd_0bar','cd_1bar',     
1586      &'Sigma_cbar0','Sigma*_cbar0','Lambda_cbar-','Xi_cbar0','cu_0bar', 
1587      &'cu_1bar','Sigma_cbar-','Sigma*_cbar-','Sigma_cbar--',            
1588      &'Sigma*_cbar--','Xi_cbar-','cs_0bar','cs_1bar','Xi''_cbar0',      
1589      &'Xi*_cbar0','Xi''_cbar-','Xi*_cbar-','Omega_cbar0',               
1590      &'Omega*_cbar0','cc_1bar','Xi_ccbar-','Xi*_ccbar-','Xi_ccbar--',   
1591      &'Xi*_ccbar--','Omega_ccbar-','Omega*_ccbar-','Omega*_cccbar-'/    
1592       DATA (CHAF(I,2),I= 206, 325)/'bd_0bar','bd_1bar','Sigma_bbar+',   
1593      &'Sigma*_bbar+','Lambda_bbar0','Xi_bbar+','Xi_bcbar0','bu_0bar',   
1594      &'bu_1bar','Sigma_bbar0','Sigma*_bbar0','Sigma_bbar-',             
1595      &'Sigma*_bbar-','Xi_bbar0','Xi_bcbar-','bs_0bar','bs_1bar',        
1596      &'Xi''_bbar+','Xi*_bbar+','Xi''_bbar0','Xi*_bbar0','Omega_bbar+',  
1597      &'Omega*_bbar+','Omega_bcbar0','bc_0bar','bc_1bar','Xi''_bcbar0',  
1598      &'Xi*_bcbar0','Xi''_bcbar-','Xi*_bcbar-','Omega''_bcba',           
1599      &'Omega*_bcbar0','Omega_bccbar-','Omega*_bccbar-','bb_1bar',       
1600      &'Xi_bbbar+','Xi*_bbbar+','Xi_bbbar0','Xi*_bbbar0','Omega_bbbar+', 
1601      &'Omega*_bbbar+','Omega_bbcbar0','Omega*_bbcbar0',                 
1602      &'Omega*_bbbbar+',2*' ','a_0-','b_1-',2*' ','K*_0bar0','K_1bar0',  
1603      &'K*_0-','K_1-',2*' ','D*_0-','D_1-','D*_0bar0','D_1bar0',         
1604      &'D*_0s-','D_1s-',2*' ','B*_0bar0','B_1bar0','B*_0-','B_1-',       
1605      &'B*_0sbar0','B_1sbar0','B*_0c-','B_1c-',3*' ','a_1-',' ',         
1606      &'K*_1bar0','K*_1-',' ','D*_1-','D*_1bar0','D*_1s-',' ',           
1607      &'B*_1bar0','B*_1-','B*_1sbar0','B*_1c-',3*' ','~d_Lbar',          
1608      &'~u_Lbar','~s_Lbar','~c_Lbar','~b_1bar','~t_1bar','~e_L+',        
1609      &'~nu_eLbar','~mu_L+','~nu_muLbar','~tau_1+','~nu_tauLbar',3*' ',  
1610      &'~chi_1-',2*' ','~chi_2-',' ','~d_Rbar','~u_Rbar','~s_Rbar',      
1611      &'~c_Rbar','~b_2bar','~t_2bar','~e_R+','~nu_eRbar','~mu_R+'/       
1612       DATA (CHAF(I,2),I= 326, 500)/'~nu_muRbar','~tau_2+',              
1613      &'~nu_tauRbar',' ','pi_tc-',3*' ','rho_tc-',8*' ','d*bar','u*bar', 
1614      &'e*bar+','nu*_ebar0',5*' ','W_R-','H_L--','H_R--',' ',            
1615      &'pi_diffr-',3*' ','n_diffrbar0','p_diffrbar-',7*' ','a_tc-',     
1616      &81*' ',
1617 C...UED
1618      &'d*_Sbar','u*_Sbar','s*_Sbar','c*_Sbar','b*_Sbar','t*_Sbar',
1619      &'d*_Dbar','u*_Dbar','s*_Dbar','c*_Dbar','b*_Dbar','t*_Dbar',
1620      &'e*_Sbar+','mu*_Sbar+','tau*_Sbar+',
1621      &'nu*_eDbar','e*_Dbar+',
1622      &'nu*_muDbar','mu*_Dbar+',
1623      &'nu*_tauDbar','tau*_Dbar+',
1624      &'g*','gamma*','Z*0','W*-',25*' '/            
1625  
1626 C...PYDATR, with initial values for the random number generator.
1627       DATA MRPY/19780503,0,0,97,33,0/
1628  
1629 C...Default values for allowed processes and kinematics constraints.
1630       DATA MSEL/1/
1631       DATA MSUB/500*0/
1632       DATA ((KFIN(I,J),J=-40,40),I=1,2)/16*0,4*1,4*0,6*1,5*0,5*1,0,
1633      &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,
1634      &6*1,4*0,4*1,16*0/
1635       DATA CKIN/
1636      &  2.0D0, -1.0D0,  0.0D0, -1.0D0,  1.0D0,
1637      &  1.0D0,  -10D0,   10D0,  -40D0,   40D0,
1638      1  -40D0,   40D0,  -40D0,   40D0,  -40D0,
1639      1   40D0, -1.0D0,  1.0D0, -1.0D0,  1.0D0,
1640      2  0.0D0,  1.0D0,  0.0D0,  1.0D0, -1.0D0,
1641      2  1.0D0, -1.0D0,  1.0D0,    0D0,    0D0,
1642      3  2.0D0, -1.0D0,    0D0,    0D0,  0.0D0,
1643      3 -1.0D0,  0.0D0, -1.0D0,  4.0D0, -1.0D0,
1644      4 12.0D0, -1.0D0, 12.0D0, -1.0D0, 12.0D0,
1645      4 -1.0D0, 12.0D0, -1.0D0,    0D0,    0D0,
1646      5  0.0D0, -1.0D0,  0.0D0, -1.0D0,  0.0D0,
1647      5 -1.0D0,    0D0,    0D0,    0D0,    0D0,
1648      6 0.0001D0, 0.99D0, 0.0001D0, 0.99D0,    0D0,
1649      6   -1D0,    0D0,   -1D0,    0D0,   -1D0,
1650      7    0D0,   -1D0, 0.0001D0, 0.99D0, 0.0001D0,
1651      7 0.99D0,    2D0,   -1D0,    0D0,    0D0,
1652      8  120*0D0/
1653  
1654 C...Default values for main switches and parameters. Reset information.
1655       DATA (MSTP(I),I=1,100)/
1656      &  3,    1,    2,    0,    0,    0,    0,    0,    0,    0,
1657      1  1,    0,    1,   30,    0,    1,    4,    3,    4,    3,
1658      2  1,    0,    1,    0,    0,    0,    0,    0,    0,    1,
1659      3  1,    8,    0,    1,    0,    2,    1,    5,    2,    0,
1660      4  2,    1,    3,    7,    3,    1,    1,    0,    1,    0,
1661      5  7,    1,    3,    1,    5,    1,    1,    5,    1,    7,
1662      6  2,    3,    2,    2,    1,    5,    2,    3,    0,    0,
1663      7  1,    1,    0,    0,    0,    0,    0,    0,    0,    0,
1664      8  1,    4,  100,    1,    1,    2,    4,    1,    1,    0,
1665      9  1,    3,    1,    3,    1,    0,    0,    0,    0,    0/
1666       DATA (MSTP(I),I=101,200)/
1667      &  3,    1,    0,    0,    0,    0,    0,    0,    0,    0,
1668      1  1,    1,    1,    0,    0,    0,    0,    0,    0,    0,
1669      2  0,    1,    2,    1,    1,  100,    0,    0,   10,    0,
1670      3  0,    4,    0,    1,    0,    0,    0,    0,    0,    0,
1671      4  0,    0,    0,    0,    0,    1,    0,    0,    0,    0,
1672      5  0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1673      6  0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1674      7  0,    2,    0,    0,    0,    0,    0,    0,    0,    0,
1675      8  6,  424, 2010,   10,   21,    0,    0,    0,    0,    0,
1676      9  0,    0,    0,    0,    0,    0,    0,    0,    0,    0/
1677       DATA (PARP(I),I=1,100)/
1678      &  0.25D0,  10D0, 8*0D0,
1679      1  0D0, 0D0, 1.0D0, 0.01D0, 0.5D0, 1.0D0, 1.0D0, 0.4D0, 2*0D0,
1680      2  10*0D0,
1681      3  1.5D0,2.0D0,0.075D0,1.0D0,0.2D0,0D0,1.0D0,0.70D0,0.006D0,0D0,
1682      4  0.02D0,2.0D0,0.10D0,1000D0,2054D0,123D0,246D0,50D0,0D0,0.054D0,
1683      5  10*0D0,
1684      6  0.25D0, 1.0D0,0.25D0, 1.0D0, 2.0D0,1D-3, 4.0D0,1D-3,2*0D0,
1685      7  4.0D0, 0.25D0, 5*0D0, 0.025D0, 2.0D0, 0.1D0,
1686      8  1.90D0, 2.0D0, 0.5D0, 0.4D0, 0.90D0,
1687      8  0.95D0, 0.7D0, 0.5D0, 1800D0, 0.25D0,
1688      9  2.0D0,0.40D0,5.0D0,1.0D0,0.0D0,3.0D0,1.0D0,0.75D0,1.0D0,5.0D0/
1689       DATA (PARP(I),I=101,200)/
1690      &  0.5D0, 0.28D0,  1.0D0, 0.8D0, 0D0, 0D0, 0D0, 0D0, 0D0, 1D0,
1691      1  2.0D0, 3*0D0, 1.5D0, 0.5D0, 0.6D0, 2.5D0, 2.0D0, 1.0D0,
1692      2  1.0D0,  0.4D0, 8*0D0,
1693      3  0.01D0, 9*0D0,
1694      4  1.16D0, 0.0119D0, 0.01D0, 0.01D0, 0.05D0, 
1695      4  9.28D0, 0.15D0, 0.02D0, 0.48D0, 0.09D0,
1696      5  0D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0,
1697      6  2.20D0, 23.6D0, 18.4D0, 11.5D0, 0.5D0, 0D0, 0D0, 0D0, 2*0D0,
1698      7  0D0,   0D0,   0D0,  1.0D0, 6*0D0,
1699      8  0.1D0, 0.01D0, 0.01D0, 0.01D0, 0.1D0, 0.01D0, 0.01D0, 0.01D0,
1700      8  0.3D0, 0.64D0,
1701      9  0.64D0, 5.0D0, 1.0D4, 1.0D4, 6*0D0/
1702       DATA MSTI/200*0/
1703       DATA PARI/200*0D0/
1704       DATA MINT/400*0/
1705       DATA VINT/400*0D0/
1706  
1707 C...Constants for the generation of the various processes.
1708       DATA (ISET(I),I=1,100)/
1709      &  1,    1,    1,   -1,    3,   -1,   -1,    3,   -2,    2,
1710      1  2,    2,    2,    2,    2,    2,   -1,    2,    2,    2,
1711      2 -1,    2,    2,    2,    2,    2,   -1,    2,    2,    2,
1712      3  2,    2,    2,    2,    2,    2,   -1,   -1,   -1,   -1,
1713      4 -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,
1714      5 -1,   -1,    2,    2,   -1,   -1,   -1,    2,   -1,   -1,
1715      6 -1,   -1,   -1,   -1,   -1,   -1,   -1,    2,    2,    2,
1716      7  4,    4,    4,   -1,   -1,    4,    4,   -1,   -1,    2,
1717      8  2,    2,    2,    2,    2,    2,    2,    2,    2,   -2,
1718      9  0,    0,    0,    0,    0,    9,   -2,   -2,    8,   -2/
1719       DATA (ISET(I),I=101,200)/
1720      & -1,    1,    1,    1,    1,    2,    2,    2,   -2,    2,
1721      1  2,    2,    2,    2,    2,   -1,   -1,   -1,   -2,   -2,
1722      2  5,    5,    5,    5,   -2,   -2,   -2,   -2,   -2,   -2,
1723      3  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1724      4  1,    1,    1,    1,    1,    1,    1,    1,    1,   -2,
1725      5  1,    1,    1,   -2,   -2,    1,    1,    1,   -2,   -2,
1726      6  2,    2,    2,    2,    2,    2,    2,    2,    2,   -2,
1727      7  2,    2,    5,    5,   -2,    2,    2,    5,    5,   -2,
1728      8  5,    5,    2,    2,    2,    5,    5,    2,    2,    2,
1729      9  1,    1,    1,    2,    2,   -2,   -2,   -2,   -2,   -2/
1730       DATA (ISET(I),I=201,300)/
1731      &  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1732      1  2,    2,    2,    2,   -2,    2,    2,    2,    2,    2,
1733      2  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1734      3  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1735      4  2,    2,    2,    2,   -1,    2,    2,    2,    2,    2,
1736      5  2,    2,    2,    2,   -1,    2,   -1,    2,    2,   -2,
1737      6  2,    2,    2,    2,    2,   -1,   -1,   -1,   -1,   -1,
1738      7  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1739      8  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1740      9  2,    2,    2,    2,    2,    2,    2,    2,    2,    2/
1741       DATA (ISET(I),I=301,500)/
1742      &  2, 9*-2, 9*2, 21*-2,
1743      4  1,    1,    2,    2,    2,    2,    2,    2,    2,    2,
1744      5  5,    5,    1,    1,   -1,   -1,   -1,   -1,   -1,   -1,
1745      6  2,    2,    2,    2,    2,    2,    2,    2,   -1,    2,
1746      7  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1747      8  2,    2,    2,    2,    2,    2,    2,    2,   -2,   -2,
1748      9  1,    1,    2,    2,    2, 5*-2,
1749      &  5,    5, 18*-2,
1750      2  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1751      3  2,    2,    2,    2,    2,    2,    2,    2,    2, 21*-2,
1752      6  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1753      7  2,    2,    2,    2,    2,    2,    2,    2,    2,   -2,
1754      8  2,    2,  18*-2/
1755       DATA ((KFPR(I,J),J=1,2),I=1,50)/
1756      &  23,    0,   24,    0,   25,    0,   24,    0,   25,    0,
1757      &  24,    0,   23,    0,   25,    0,    0,    0,    0,    0,
1758      1   0,    0,    0,    0,   21,   21,   21,   22,   21,   23,
1759      1  21,   24,   21,   25,   22,   22,   22,   23,   22,   24,
1760      2  22,   25,   23,   23,   23,   24,   23,   25,   24,   24,
1761      2  24,   25,   25,   25,    0,   21,    0,   22,    0,   23,
1762      3   0,   24,    0,   25,    0,   21,    0,   22,    0,   23,
1763      3   0,   24,    0,   25,    0,   21,    0,   22,    0,   23,
1764      4   0,   24,    0,   25,    0,   21,    0,   22,    0,   23,
1765      4   0,   24,    0,   25,    0,   21,    0,   22,    0,   23/
1766       DATA ((KFPR(I,J),J=1,2),I=51,100)/
1767      5   0,   24,    0,   25,    0,    0,    0,    0,    0,    0,
1768      5   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1769      6   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1770      6   0,    0,    0,    0,   21,   21,   24,   24,   23,   24,
1771      7  23,   23,   24,   24,   23,   24,   23,   25,   22,   22,
1772      7  23,   23,   24,   24,   24,   25,   25,   25,    0,  211,
1773      8   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1774      8 443,   21,10441,   21,20443,   21,  445,   21,    0,    0,
1775      9   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1776      9   0,    0,    0,    0,    0,    0,    0,    0,    0,    0/
1777       DATA ((KFPR(I,J),J=1,2),I=101,150)/
1778      &  23,    0,   25,    0,   25,    0,10441,    0,  445,    0,
1779      & 443,   22,  443,   21,  443,   22,    0,    0,   22,   25,
1780      1  21,   25,    0,   25,   21,   25,   22,   22,   21,   22,
1781      1  22,   23,   23,   23,   24,   24,    0,    0,    0,    0,
1782      2  25,    6,   25,    6,   25,    0,   25,    0,    0,    0,
1783      2   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1784      3   0,   21,    0,   21,    0,   22,    0,   22,    0,    0,
1785      3   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1786      4  32,    0,   34,    0,   37,    0,   41,    0,   42,    0,
1787      4 4000011, 0, 4000001, 0, 4000002, 0, 3000331, 0,   0,    0/
1788       DATA ((KFPR(I,J),J=1,2),I=151,200)/
1789      5  35,    0,   35,    0,   35,    0,    0,    0,    0,    0,
1790      5  36,    0,   36,    0,   36,    0,    0,    0,    0,    0,
1791      6   6,   37,   42,    0,   42,   42,   42,   42,   11,    0,
1792      6  11,    0, 0, 4000001, 0, 4000002, 0, 4000011,    0,    0,
1793      7  23,   35,   24,   35,   35,    0,   35,    0,    0,    0,
1794      7  23,   36,   24,   36,   36,    0,   36,    0,    0,    0,
1795      8  35,    6,   35,    6,   21,   35,    0,   35,   21,   35,
1796      8  36,    6,   36,    6,   21,   36,    0,   36,   21,   36,
1797      9  3000113, 0, 3000213, 0, 3000223, 0, 11,    0,   11,    0,
1798      9   0,    0,    0,    0,    0,    0,    0,    0,    0,    0/
1799       DATA ((KFPR(I,J),J=1,2),I=201,240)/
1800      &  1000011,   1000011,   2000011,   2000011,   1000011,
1801      &  2000011,   1000013,   1000013,   2000013,   2000013,
1802      &  1000013,   2000013,   1000015,   1000015,   2000015,
1803      &  2000015,   1000015,   2000015,   1000011,   1000012,
1804      1  1000015,   1000016,   2000015,   1000016,   1000012,
1805      1  1000012,   1000016,   1000016,         0,         0,
1806      1  1000022,   1000022,   1000023,   1000023,   1000025,
1807      1  1000025,   1000035,   1000035,   1000022,   1000023,
1808      2  1000022,   1000025,   1000022,   1000035,   1000023,
1809      2  1000025,   1000023,   1000035,   1000025,   1000035,
1810      2  1000024,   1000024,   1000037,   1000037,   1000024,
1811      2  1000037,   1000022,   1000024,   1000023,   1000024,
1812      3  1000025,   1000024,   1000035,   1000024,   1000022,
1813      3  1000037,   1000023,   1000037,   1000025,   1000037,
1814      3  1000035,   1000037,   1000021,   1000022,   1000021,
1815      3  1000023,   1000021,   1000025,   1000021,   1000035/
1816       DATA ((KFPR(I,J),J=1,2),I=241,280)/
1817      4  1000021,   1000024,   1000021,   1000037,   1000021,
1818      4  1000021,   1000021,   1000021,         0,         0,
1819      4  1000002,   1000022,   2000002,   1000022,   1000002,
1820      4  1000023,   2000002,   1000023,   1000002,   1000025,
1821      5  2000002,   1000025,   1000002,   1000035,   2000002,
1822      5  1000035,   1000001,   1000024,   2000005,   1000024,
1823      5  1000001,   1000037,   2000005,   1000037,   1000002,
1824      5  1000021,   2000002,   1000021,         0,         0,
1825      6  1000006,   1000006,   2000006,   2000006,   1000006,
1826      6  2000006,   1000006,   1000006,   2000006,   2000006,
1827      6        0,         0,         0,         0,         0,
1828      6        0,         0,         0,         0,         0,
1829      7  1000002,   1000002,   2000002,   2000002,   1000002,
1830      7  2000002,   1000002,   1000002,   2000002,   2000002,
1831      7  1000002,   2000002,   1000002,   1000002,   2000002,
1832      7  2000002,   1000002,   1000002,   2000002,   2000002/
1833       DATA ((KFPR(I,J),J=1,2),I=281,350)/
1834      8  1000005,   1000002,   2000005,   2000002,   1000005,
1835      8  2000002,   1000005,   1000002,   2000005,   2000002,
1836      8  1000005,   2000002,   1000005,   1000005,   2000005,
1837      8  2000005,   1000005,   1000005,   2000005,   2000005,
1838      9  1000005,   1000005,   2000005,   2000005,   1000005,
1839      9  2000005,   1000005,   1000021,   2000005,   1000021,
1840      9  1000005,   2000005,        37,        25,        37,
1841      9       35,        36,        25,        36,        35,
1842      &       37,        37,      18*0,
1843 C...UED: 311-319
1844      &  5100021,   5100021, 
1845      &  5100002,   5100021, 
1846      &  5100002,   5100001,
1847      &  5100002,  -5100002, 
1848      &  5100002,  -5100002,
1849      &  5100002,  -6100001,
1850      &  5100002,  -5100001,
1851      &  5100002,   6100001,
1852      &  5100001,  -5100001,
1853      &  42*0,
1854      4  9900041,         0,   9900042,         0,   9900041,
1855      4       11,   9900042,        11,   9900041,        13,
1856      4  9900042,        13,   9900041,        15,   9900042,
1857      4       15,   9900041,   9900041,   9900042,   9900042/
1858       DATA ((KFPR(I,J),J=1,2),I=351,400)/
1859      5  9900041,         0,   9900042,         0,   9900023,
1860      5        0,   9900024,         0,         0,         0,
1861      5        0,         0,         0,         0,         0,
1862      5        0,         0,         0,         0,         0,
1863      6       24,        24,        24,   3000211,   3000211,
1864      6  3000211,        22,   3000111,        22,   3000221,
1865      6       23,   3000111,        23,   3000221,        24,
1866      6  3000211,         0,         0,        24,        23,
1867      7       24,   3000111,   3000211,        23,   3000211,
1868      7  3000111,        22,   3000211,        23,   3000211,
1869      7       24,   3000111,        24,   3000221,        22,
1870      7       24,        22,        23,        23,        23,
1871      8   0,    0,    0,    0,   21,   21,    0,   21,    0,    0,
1872      8  21,   21,    0,    0,    0,    0,    0,    0,    0,    0,
1873      9  5000039,         0,   5000039,         0,        21,
1874      9  5000039,         0,   5000039,        21,   5000039,
1875      9     10*0/
1876       DATA ((KFPR(I,J),J=1,2),I=401,500)/
1877      &  37,    6,   37,    6,    36*0,
1878      2      443,        21,   9900443,        21,   9900441,
1879      2       21,   9910441,        21,         0,   9900443,
1880      2        0,   9900441,         0,   9910441,        21,
1881      2  9900443,        21,   9900441,        21,   9910441,
1882      3 10441, 21, 20443,  21,  445,   21,    0, 10441,   0, 20443,
1883      3   0,  445,   21, 10441,  21, 20443,  21,  445,  42*0,
1884      6      553,        21,   9900553,        21,   9900551,
1885      6       21,   9910551,        21,         0,   9900553,
1886      6        0,   9900551,         0,   9910551,        21,
1887      6  9900553,        21,   9900551,        21,   9910551,
1888      7 10551, 21, 20553,  21,  555,   21,    0, 10551,   0, 20553,
1889      7   0,  555,   21, 10551,  21, 20553,  21,  555, 42*0/
1890       DATA COEF/10000*0D0/
1891       DATA (((ICOL(I,J,K),K=1,2),J=1,4),I=1,40)/
1892      &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,
1893      &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,
1894      &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,
1895      &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,
1896      &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,
1897      &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,
1898      &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,
1899      &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,
1900      &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,
1901      &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/
1902  
1903 C...Treatment of resonances.
1904       DATA (MWID(I)  ,I=   1, 500)/5*0,3*1,8*0,1,5*0,3*1,6*0,1,0,4*1,   
1905      &3*0,2*1,254*0,19*2,0,7*2,0,2,0,2,0,26*1,7*0,6*2,2*1,
1906      &81*0,21*1,4*1,25*0/
1907  
1908 C...Character constants: name of processes.
1909       DATA PROC(0)/                    'All included subprocesses   '/
1910       DATA (PROC(I),I=1,20)/
1911      &'f + fbar -> gamma*/Z0       ',  'f + fbar'' -> W+/-           ',
1912      &'f + fbar -> h0              ',  'gamma + W+/- -> W+/-        ',
1913      &'Z0 + Z0 -> h0               ',  'Z0 + W+/- -> W+/-           ',
1914      &'                            ',  'W+ + W- -> h0               ',
1915      &'                            ',  'f + f'' -> f + f'' (QFD)      ',
1916      1'f + f'' -> f + f'' (QCD)      ','f + fbar -> f'' + fbar''      ',
1917      1'f + fbar -> g + g           ',  'f + fbar -> g + gamma       ',
1918      1'f + fbar -> g + Z0          ',  'f + fbar'' -> g + W+/-       ',
1919      1'f + fbar -> g + h0          ',  'f + fbar -> gamma + gamma   ',
1920      1'f + fbar -> gamma + Z0      ',  'f + fbar'' -> gamma + W+/-   '/
1921       DATA (PROC(I),I=21,40)/
1922      2'f + fbar -> gamma + h0      ',  'f + fbar -> Z0 + Z0         ',
1923      2'f + fbar'' -> Z0 + W+/-      ', 'f + fbar -> Z0 + h0         ',
1924      2'f + fbar -> W+ + W-         ',  'f + fbar'' -> W+/- + h0      ',
1925      2'f + fbar -> h0 + h0         ',  'f + g -> f + g              ',
1926      2'f + g -> f + gamma          ',  'f + g -> f + Z0             ',
1927      3'f + g -> f'' + W+/-          ', 'f + g -> f + h0             ',
1928      3'f + gamma -> f + g          ',  'f + gamma -> f + gamma      ',
1929      3'f + gamma -> f + Z0         ',  'f + gamma -> f'' + W+/-      ',
1930      3'f + gamma -> f + h0         ',  'f + Z0 -> f + g             ',
1931      3'f + Z0 -> f + gamma         ',  'f + Z0 -> f + Z0            '/
1932       DATA (PROC(I),I=41,60)/
1933      4'f + Z0 -> f'' + W+/-         ', 'f + Z0 -> f + h0            ',
1934      4'f + W+/- -> f'' + g          ', 'f + W+/- -> f'' + gamma      ',
1935      4'f + W+/- -> f'' + Z0         ', 'f + W+/- -> f'' + W+/-       ',
1936      4'f + W+/- -> f'' + h0         ', 'f + h0 -> f + g             ',
1937      4'f + h0 -> f + gamma         ',  'f + h0 -> f + Z0            ',
1938      5'f + h0 -> f'' + W+/-         ', 'f + h0 -> f + h0            ',
1939      5'g + g -> f + fbar           ',  'g + gamma -> f + fbar       ',
1940      5'g + Z0 -> f + fbar          ',  'g + W+/- -> f + fbar''       ',
1941      5'g + h0 -> f + fbar          ',  'gamma + gamma -> f + fbar   ',
1942      5'gamma + Z0 -> f + fbar      ',  'gamma + W+/- -> f + fbar''   '/
1943       DATA (PROC(I),I=61,80)/
1944      6'gamma + h0 -> f + fbar      ',  'Z0 + Z0 -> f + fbar         ',
1945      6'Z0 + W+/- -> f + fbar''      ', 'Z0 + h0 -> f + fbar         ',
1946      6'W+ + W- -> f + fbar         ',  'W+/- + h0 -> f + fbar''      ',
1947      6'h0 + h0 -> f + fbar         ',  'g + g -> g + g              ',
1948      6'gamma + gamma -> W+ + W-    ',  'gamma + W+/- -> Z0 + W+/-   ',
1949      7'Z0 + Z0 -> Z0 + Z0          ',  'Z0 + Z0 -> W+ + W-          ',
1950      7'Z0 + W+/- -> Z0 + W+/-      ',  'Z0 + Z0 -> Z0 + h0          ',
1951      7'W+ + W- -> gamma + gamma    ',  'W+ + W- -> Z0 + Z0          ',
1952      7'W+/- + W+/- -> W+/- + W+/-  ',  'W+/- + h0 -> W+/- + h0      ',
1953      7'h0 + h0 -> h0 + h0          ',  'q + gamma -> q'' + pi+/-     '/
1954       DATA (PROC(I),I=81,100)/
1955      8'q + qbar -> Q + Qbar, mass  ',  'g + g -> Q + Qbar, massive  ',
1956      8'f + q -> f'' + Q, massive    ', 'g + gamma -> Q + Qbar, mass ',
1957      8'gamma + gamma -> F + Fbar, m',  'g + g -> J/Psi + g          ',
1958      8'g + g -> chi_0c + g         ',  'g + g -> chi_1c + g         ',
1959      8'g + g -> chi_2c + g         ',  '                            ',
1960      9'Elastic scattering          ',  'Single diffractive (XB)     ',
1961      9'Single diffractive (AX)     ',  'Double  diffractive         ',
1962      9'Low-pT scattering           ',  'Semihard QCD 2 -> 2         ',
1963      9'                            ',  '                            ',
1964      9'q + gamma* -> q             ',  '                            '/
1965       DATA (PROC(I),I=101,120)/
1966      &'g + g -> gamma*/Z0          ',  'g + g -> h0                 ',
1967      &'gamma + gamma -> h0         ',  'g + g -> chi_0c             ',
1968      &'g + g -> chi_2c             ',  'g + g -> J/Psi + gamma      ',
1969      &'gamma + g -> J/Psi + g      ',  'gamma+gamma -> J/Psi + gamma',
1970      &'                            ',  'f + fbar -> gamma + h0      ',
1971      1'q + qbar -> g + h0          ',  'q + g -> q + h0             ',
1972      1'g + g -> g + h0             ',  'g + g -> gamma + gamma      ',
1973      1'g + g -> g + gamma          ',  'g + g -> gamma + Z0         ',
1974      1'g + g -> Z0 + Z0            ',  'g + g -> W+ + W-            ',
1975      1'                            ',  '                            '/
1976       DATA (PROC(I),I=121,140)/
1977      2'g + g -> Q + Qbar + h0      ',  'q + qbar -> Q + Qbar + h0   ',
1978      2'f + f'' -> f + f'' + h0       ',
1979      2'f + f'' -> f" + f"'' + h0     ',
1980      2'                            ',  '                            ',
1981      2'                            ',  '                            ',
1982      2'                            ',  '                            ',
1983      3'f + gamma*_T -> f + g       ',  'f + gamma*_L -> f + g       ',
1984      3'f + gamma*_T -> f + gamma   ',  'f + gamma*_L -> f + gamma   ',
1985      3'g + gamma*_T -> f + fbar    ',  'g + gamma*_L -> f + fbar    ',
1986      3'gamma*_T+gamma*_T -> f+fbar ',  'gamma*_T+gamma*_L -> f+fbar ',
1987      3'gamma*_L+gamma*_T -> f+fbar ',  'gamma*_L+gamma*_L -> f+fbar '/
1988       DATA (PROC(I),I=141,160)/
1989      4'f + fbar -> gamma*/Z0/Z''0   ', 'f + fbar'' -> W''+/-          ',
1990      4'f + fbar'' -> H+/-           ', 'f + fbar'' -> R              ',
1991      4'q + l -> LQ                 ',  'e + gamma -> e*             ',
1992      4'd + g -> d*                 ',  'u + g -> u*                 ',
1993      4'g + g -> eta_tc             ',  '                            ',
1994      5'f + fbar -> H0              ',  'g + g -> H0                 ',
1995      5'gamma + gamma -> H0         ',  '                            ',
1996      5'                            ',  'f + fbar -> A0              ',
1997      5'g + g -> A0                 ',  'gamma + gamma -> A0         ',
1998      5'                            ',  '                            '/
1999       DATA (PROC(I),I=161,180)/
2000      6'f + g -> f'' + H+/-          ', 'q + g -> LQ + lbar          ',
2001      6'g + g -> LQ + LQbar         ',  'q + qbar -> LQ + LQbar      ',
2002      6'f + fbar -> f'' + fbar'' (g/Z)',
2003      6'f +fbar'' -> f" + fbar"'' (W) ',
2004      6'q + q'' -> q" + d*           ',  'q + q'' -> q" + u*           ',
2005      6'q + qbar -> e + e*          ',  '                            ',
2006      7'f + fbar -> Z0 + H0         ', 'f + fbar'' -> W+/- + H0      ',
2007      7'f + f'' -> f + f'' + H0       ',
2008      7'f + f'' -> f" + f"'' + H0     ',
2009      7'                            ',  'f + fbar -> Z0 + A0         ',
2010      7'f + fbar'' -> W+/- + A0      ',
2011      7'f + f'' -> f + f'' + A0       ',
2012      7'f + f'' -> f" + f"'' + A0     ',
2013      7'                            '/
2014       DATA (PROC(I),I=181,200)/
2015      8'g + g -> Q + Qbar + H0      ',  'q + qbar -> Q + Qbar + H0   ',
2016      8'q + qbar -> g + H0          ',  'q + g -> q + H0             ',
2017      8'g + g -> g + H0             ',  'g + g -> Q + Qbar + A0      ',
2018      8'q + qbar -> Q + Qbar + A0   ',  'q + qbar -> g + A0          ',
2019      8'q + g -> q + A0             ',  'g + g -> g + A0             ',
2020      9'f + fbar -> rho_tc0         ',  'f + f'' -> rho_tc+/-         ',
2021      9'f + fbar -> omega_tc0      ',  'f+fbar -> f''+fbar'' (ETC)  ',
2022      9'f+fbar'' -> f"+fbar"'' (ETC)','                          ',
2023      9'                            ',  '                            ',
2024      9'                            ',  '                            '/
2025       DATA (PROC(I),I=201,220)/
2026      &'f + fbar -> ~e_L + ~e_Lbar  ',  'f + fbar -> ~e_R + ~e_Rbar  ',
2027      &'f + fbar -> ~e_R + ~e_Lbar  ',  'f + fbar -> ~mu_L + ~mu_Lbar',
2028      &'f + fbar -> ~mu_R + ~mu_Rbar',  'f + fbar -> ~mu_L + ~mu_Rbar',
2029      &'f+fbar -> ~tau_1 + ~tau_1bar',  'f+fbar -> ~tau_2 + ~tau_2bar',
2030      &'f+fbar -> ~tau_1 + ~tau_2bar',  'q + qbar'' -> ~l_L + ~nulbar ',
2031      1'q+qbar''-> ~tau_1 + ~nutaubar', 'q+qbar''-> ~tau_2 + ~nutaubar',
2032      1'f + fbar -> ~nul + ~nulbar  ',  'f+fbar -> ~nutau + ~nutaubar',
2033      1'                            ',  'f + fbar -> ~chi1 + ~chi1   ',
2034      1'f + fbar -> ~chi2 + ~chi2   ',  'f + fbar -> ~chi3 + ~chi3   ',
2035      1'f + fbar -> ~chi4 + ~chi4   ',  'f + fbar -> ~chi1 + ~chi2   '/
2036       DATA (PROC(I),I=221,240)/
2037      2'f + fbar -> ~chi1 + ~chi3   ',  'f + fbar -> ~chi1 + ~chi4   ',
2038      2'f + fbar -> ~chi2 + ~chi3   ',  'f + fbar -> ~chi2 + ~chi4   ',
2039      2'f + fbar -> ~chi3 + ~chi4   ',  'f+fbar -> ~chi+-1 + ~chi-+1 ',
2040      2'f+fbar -> ~chi+-2 + ~chi-+2 ',  'f+fbar -> ~chi+-1 + ~chi-+2 ',
2041      2'q + qbar'' -> ~chi1 + ~chi+-1', 'q + qbar'' -> ~chi2 + ~chi+-1',
2042      3'q + qbar'' -> ~chi3 + ~chi+-1', 'q + qbar'' -> ~chi4 + ~chi+-1',
2043      3'q + qbar'' -> ~chi1 + ~chi+-2', 'q + qbar'' -> ~chi2 + ~chi+-2',
2044      3'q + qbar'' -> ~chi3 + ~chi+-2', 'q + qbar'' -> ~chi4 + ~chi+-2',
2045      3'q + qbar -> ~chi1 + ~g      ',  'q + qbar -> ~chi2 + ~g      ',
2046      3'q + qbar -> ~chi3 + ~g      ',  'q + qbar -> ~chi4 + ~g      '/
2047       DATA (PROC(I),I=241,260)/
2048      4'q + qbar'' -> ~chi+-1 + ~g   ', 'q + qbar'' -> ~chi+-2 + ~g  ',
2049      4'q + qbar -> ~g + ~g         ',  'g + g -> ~g + ~g            ',
2050      4'                            ',  'qj + g -> ~qj_L + ~chi1     ',
2051      4'qj + g -> ~qj_R + ~chi1     ',  'qj + g -> ~qj_L + ~chi2     ',
2052      4'qj + g -> ~qj_R + ~chi2     ',  'qj + g -> ~qj_L + ~chi3     ',
2053      5'qj + g -> ~qj_R + ~chi3     ',  'qj + g -> ~qj_L + ~chi4     ',
2054      5'qj + g -> ~qj_R + ~chi4     ',  'qj + g -> ~qk_L + ~chi+-1   ',
2055      5'qj + g -> ~qk_R + ~chi+-1   ',  'qj + g -> ~qk_L + ~chi+-2   ',
2056      5'qj + g -> ~qk_R + ~chi+-2   ',  'qj + g -> ~qj_L + ~g        ',
2057      5'qj + g -> ~qj_R + ~g        ',  '                            '/
2058       DATA (PROC(I),I=261,300)/
2059      6'f + fbar -> ~t_1 + ~t_1bar  ',  'f + fbar -> ~t_2 + ~t_2bar  ',
2060      6'f + fbar -> ~t_1 + ~t_2bar  ',  'g + g -> ~t_1 + ~t_1bar     ',
2061      6'g + g -> ~t_2 + ~t_2bar     ',  '                            ',
2062      6'                            ',  '                            ',
2063      6'                            ',  '                            ',
2064      7'qi + qj -> ~qi_L + ~qj_L    ',  'qi + qj -> ~qi_R + ~qj_R    ',
2065      7'qi + qj -> ~qi_L + ~qj_R    ',  'qi+qjbar -> ~qi_L + ~qj_Lbar',
2066      7'qi+qjbar -> ~qi_R + ~qj_Rbar',  'qi+qjbar -> ~qi_L + ~qj_Rbar',
2067      7'f + fbar -> ~qi_L + ~qi_Lbar',  'f + fbar -> ~qi_R + ~qi_Rbar',
2068      7'g + g -> ~qi_L + ~qi_Lbar   ',  'g + g -> ~qi_R + ~qi_Rbar   ',
2069      8'b + qj -> ~b_1 + ~qj_L      ',  'b + qj -> ~b_2 + ~qj_R      ',
2070      8'b + qj -> ~b_1 + ~qj_R      ',  'b + qjbar -> ~b_1 + ~qj_Lbar',
2071      8'b + qjbar -> ~b_2 + ~qj_Rbar',  'b + qjbar -> ~b_1 + ~qj_Rbar',
2072      8'f + fbar -> ~b_1 + ~b_1bar  ',  'f + fbar -> ~b_2 + ~b_2bar  ',
2073      8'g + g -> ~b_1 + ~b_1bar     ',  'g + g -> ~b_2 + ~b_2bar     ',
2074      9'b + b -> ~b_1 + ~b_1        ',  'b + b -> ~b_2 + ~b_2        ',
2075      9'b + b -> ~b_1 + ~b_2        ',  'b + g -> ~b_1 + ~g          ',
2076      9'b + g -> ~b_2 + ~g          ',  'b + bbar -> ~b_1 + ~b_2bar  ',
2077      9'f + fbar'' -> H+/- + h0     ',  'f + fbar -> H+/- + H0       ',
2078      9'f + fbar -> A0 + h0         ',  'f + fbar -> A0 + H0         '/
2079       DATA (PROC(I),I=301,340)/
2080      &'f + fbar -> H+ + H-         ',
2081      &9*'                          ',  'g + g -> g* + g*            ',
2082      &'q + g -> q*_D + g*          ',  'qi + qj -> q*_Di + q*_Dj    ',
2083      &'g + g -> q*_D + q*_Dbar     ',  'q  + qbar -> q*_D + q*_Dbar ',
2084      &'qi + qbarj -> q*Di + q*Sbarj',  'qi + qjbar -> q*Di + q*Dbarj',
2085      &'qi + qj -> q*_Di + q*_Sj    ',  'qi + qibar -> q*Dj + q*Dbarj',
2086      &21*'                          '/
2087       DATA (PROC(I),I=341,380)/
2088      4'l + l -> H_L++/--           ',  'l + l -> H_R++/--           ',
2089      4'l + gamma -> H_L++/-- e-/+  ',  'l + gamma -> H_R++/-- e-/+  ',
2090      4'l + gamma -> H_L++/-- mu-/+ ',  'l + gamma -> H_R++/-- mu-/+ ',
2091      4'l + gamma -> H_L++/-- tau-/+',  'l + gamma -> H_R++/-- tau-/+',
2092      4'f + fbar -> H_L++ + H_L--   ',  'f + fbar -> H_R++ + H_R--   ',
2093      5'f + f -> f'' + f'' + H_L++/-- ',
2094      5'f + f -> f'' + f'' + H_R++/-- ','f + fbar -> Z_R0            ',
2095      5'f + fbar'' -> W_R+/-         ',5*'                            ',
2096      6'                            ',  'f + fbar -> W_L+ W_L-       ',
2097      6'f + fbar -> W_L+/- pi_T-/+  ',  'f + fbar -> pi_T+ pi_T-     ',
2098      6'f + fbar -> gamma pi_T0     ',  'f + fbar -> gamma pi_T0''    ',
2099      6'f + fbar -> Z0 pi_T0        ',  'f + fbar -> Z0 pi_T0''       ',
2100      6'f + fbar -> W+/- pi_T-/+    ',  '                            ',
2101      7'f + fbar'' -> W_L+/- Z_L0    ', 'f + fbar'' -> W_L+/- pi_T0   ',
2102      7'f + fbar'' -> pi_T+/- Z_L0   ', 'f + fbar'' -> pi_T+/- pi_T0  ',
2103      7'f + fbar'' -> gamma pi_T+/-  ', 'f + fbar'' -> Z0 pi_T+/-     ',
2104      7'f + fbar'' -> W+/- pi_T0     ',
2105      7'f + fbar'' -> W+/- pi_T0''    ',
2106      7'f + fbar'' -> gamma W+/-(ETC)','f + fbar -> gamma Z0 (ETC)',
2107      7'f + fbar -> Z0 Z0 (ETC)     '/
2108       DATA (PROC(I),I=381,420)/
2109      8'f + f'' -> f + f'' (ETC)      ','f + fbar -> f'' + fbar'' (ETC)',
2110      8'f + fbar -> g + g (ETC)     ',  'f + g -> f + g (ETC)        ',
2111      8'g + g -> f + fbar (ETC)     ',  'g + g -> g + g (ETC)        ',
2112      8'q + qbar -> Q + Qbar (ETC)  ',  'g + g -> Q + Qbar (ETC)     ',
2113      8'                            ',  '                            ',
2114      9'f + fbar -> G*              ',  'g + g -> G*                 ',
2115      9'q + qbar -> g + G*          ',  'q + g -> q + G*             ',
2116      9'g + g -> g + G*             ',  '                            ',
2117      9 4*'                         ',
2118      &'g + g -> t + b + H+/-       ',  'q + qbar -> t + b + H+/-    ',
2119      & 18*'                            '/
2120       DATA (PROC(I),I=421,460)/
2121      2'g + g  -> cc~[3S1(1)] + g   ',  'g + g  -> cc~[3S1(8)] + g   ',
2122      2'g + g  -> cc~[1S0(8)] + g   ',  'g + g  -> cc~[3PJ(8)] + g   ',
2123      2'g + q  -> q + cc~[3S1(8)]   ',  'g + q  -> q + cc~[1S0(8)]   ',
2124      2'g + q  -> q + cc~[3PJ(8)]   ',  'q + q~ -> g + cc~[3S1(8)]   ',
2125      2'q + q~ -> g + cc~[1S0(8)]   ',  'q + q~ -> g + cc~[3PJ(8)]   ',
2126      3'g + g  -> cc~[3P0(1)] + g   ',  'g + g  -> cc~[3P1(1)] + g   ',
2127      3'g + g  -> cc~[3P2(1)] + g   ',  'q + g  -> q + cc~[3P0(1)]   ',
2128      3'q + g  -> q + cc~[3P1(1)]   ',  'q + g  -> q + cc~[3P2(1)]   ',
2129      3'q + q~ -> g + cc~[3P0(1)]   ',  'q + q~ -> g + cc~[3P1(1)]   ',
2130      3'q + q~ -> g + cc~[3P2(1)]   ',
2131      3     21 *'                            '/
2132       DATA (PROC(I),I=461,500)/
2133      6'g + g  -> bb~[3S1(1)] + g   ',  'g + g  -> bb~[3S1(8)] + g   ',
2134      6'g + g  -> bb~[1S0(8)] + g   ',  'g + g  -> bb~[3PJ(8)] + g   ',
2135      6'g + q  -> q + bb~[3S1(8)]   ',  'g + q  -> q + bb~[1S0(8)]   ',
2136      6'g + q  -> q + bb~[3PJ(8)]   ',  'q + q~ -> g + bb~[3S1(8)]   ',
2137      6'q + q~ -> g + bb~[1S0(8)]   ',  'q + q~ -> g + bb~[3PJ(8)]   ',
2138      7'g + g  -> bb~[3P0(1)] + g   ',  'g + g  -> bb~[3P1(1)] + g   ',
2139      7'g + g  -> bb~[3P2(1)] + g   ',  'q + g  -> q + bb~[3P0(1)]   ',
2140      7'q + g  -> q + bb~[3P1(1)]   ',  'q + g  -> q + bb~[3P2(1)]   ',
2141      7'q + q~ -> g + bb~[3P0(1)]   ',  'q + q~ -> g + bb~[3P1(1)]   ',
2142      7'q + q~ -> g + bb~[3P2(1)]   ',
2143      7     21 *'                            '/
2144  
2145 C...Cross sections and slope offsets.
2146       DATA SIGT/294*0D0/
2147  
2148 C...Supersymmetry switches and parameters.
2149       DATA IMSS/0,
2150      &  0,  0,  0,  1,  0,  0,  0,  0,  0,  0,
2151      1  89*0/
2152       DATA RMSS/0D0,
2153      &  80D0,160D0,500D0,800D0,2D0,250D0,200D0,800D0,700D0,800D0,
2154      1  700D0,500D0,250D0,200D0,800D0,400D0,0D0,0.1D0,850D0,0.041D0,
2155      2   1D0,800D0,1D4,1D4,1D4,0D0,0D0,0D0,24D17,0D0,
2156      3  10*0D0,  
2157      4  0D0,1D0,8*0D0,  
2158      5  49*0D0/
2159 C...Initial values for R-violating SUSY couplings.
2160 C...Should not be changed here. See PYMSIN.
2161       DATA RVLAM/27*0D0/
2162       DATA RVLAMP/27*0D0/
2163       DATA RVLAMB/27*0D0/
2164  
2165 C...Technicolor switches and parameters
2166       DATA ITCM/0,
2167      &  4,  0,  0,  0,  0,  0,  0,  0,  0,  0,
2168      1  89*0/
2169       DATA RTCM/0D0,
2170      &  82D0,1.333D0,.333D0,0.408D0,1D0,1D0,.0182D0,1D0,0D0,1.333D0,
2171      1  .05D0,200D0,200D0,0D0,0D0,0D0,0D0,0D0,0D0,0D0,
2172      2  .283D0,.707D0,0D0,0D0,0D0,1.667D0,250D0,250D0,.707D0,0D0,
2173      3  .707D0,0D0,1D0,0D0,0D0,0D0,0D0,0D0,0D0,0D0,
2174      4  1000D0, 1D0, 1D0, 1D0, 1D0, 0D0, 1D0, 3*200D0,
2175      4  200D0, 48*0D0/
2176  
2177 C...UED switches and parameters.
2178 C... IUED(0) empty IUED vector element
2179 C... IUED(1) UED ON(=1)/OFF(=0) switch
2180 C... IUED(2) ON(=1)/OFF(=0) switch for gravity mediated decays
2181 C... IUED(3) NFLAVOURS Number of KK excitation quark flavours
2182 C... IUED(4) N the number of large extra dimensions
2183 C... IUED(5) Selects whether the code takes Lambda (=0)
2184 C...         or Lambda*R (=1) as input.
2185 C... IUED(6) With radiative corrections to the masses (=1)
2186 C...         or without (=0)
2187 C...
2188 C... RUED(0) empty RUED vector element
2189 C... RUED(1) RINV (1/R) the curvature of the extra dimension
2190 C... RUED(2) XMD the (4+N)-dimensional Planck scale
2191 C... RUED(3) LAMUED (Lambda cutoff scale)
2192 C... RUED(4) LAMUED/RINV (feasible values are order of 10-20)
2193 C...
2194       DATA IUED/0,0,0,5,6,0,1,93*0/
2195       DATA RUED/0.D0,1000D0,5000D0,20000.,20.,95*0D0/
2196 
2197 C...Data for histogramming routines.
2198       DATA IHIST/1000,20000,55,1/
2199       DATA INDX/1000*0/
2200 
2201 C...Data for SUSY Les Houches Accord.
2202       DATA CPRO/'PYTHIA      ','PYTHIA      '/
2203       DATA CVER/'6.4         ','6.4         '/
2204       DATA MODSEL/200*0/
2205       DATA PARMIN/100*0D0/
2206       DATA RMSOFT/101*0D0/
2207       DATA AU/9*0D0/
2208       DATA AD/9*0D0/
2209       DATA AE/9*0D0/
2210  
2211       END