EPOS/epos167/epos-dky-165.f

   1 c-----------------------------------------------------------------------
   2       subroutine decayall(n)
   3 c-----------------------------------------------------------------------
   4 c  decay of objects n to nptl, including their children
   5 c-----------------------------------------------------------------------
   6       include 'epos.inc'
   7       common/cttaus/tpro,zpro,ttar,ztar,ttaus,detap,detat
   8       double precision tpro,zpro,ttar,ztar,ttaus,detap,detat
   9       ttaus=1
  10       np1=n
  11  1    np2=nptl
  12       do ip=np1,np2
  13          if(istptl(ip).eq.0)then  !consider last generation particles
  14             call hdecas(ip,iret)
  15             if(iret.eq.1)stop'error in hdecas detected in decay'
  16          endif
  17       enddo
  18       np1=np2+1
  19       if(np1.le.nptl)goto1
  20       end
  21
  22
  23 c-----------------------------------------------------------------------
  24       subroutine hdecas(i,iret)
  25 c-----------------------------------------------------------------------
  26 c  decay of object i  (main decay routine)
  27 c-----------------------------------------------------------------------
  28
  29       include 'epos.inc'
  30       double precision tpro,zpro,ttar,ztar,ttaus,detap,detat,zor,tor
  31       common/cttaus/tpro,zpro,ttar,ztar,ttaus,detap,detat
  32       double precision ttaux,ttauz
  33       integer jcdu(nflav,2)
  34       common/cdfptl/idfptl(mxptl)
  35
  36       iret=0
  37       nptlb=nptl
  38
  39 c no last generation -> no decay
  40
  41       if(istptl(i).ne.0)return
  42
  43       if(nptl.gt.mxptl-10)then
  44         call alist('end&',1,nptl)
  45         call utstop('hdecas: mxptl too small&')
  46       endif
  47 c entry
  48
  49       call utpri('hdecas',ish,ishini,5)
  50       ttauz=ttaus
  51
  52
  53 c small droplet decay
  54
  55       if(iabs(idptl(i)).gt.1e8)then
  56         stop'hdecas: no longer supported (2).       '
  57       endif
  58
  59 c  ordinary decay
  60
  61       call idmass(111,amrho0)
  62       call idmass(221,amomeg)
  63       ioi=iorptl(i)
  64       if(ioi.gt.0)then
  65       if(.not.(iabs(idptl(ioi)).lt.10000
  66      *.and.jorptl(i).eq.0))then
  67       if(idptl(i).eq.111)idptl(i)=221
  68       if(idptl(i).eq.221.and.rangen().gt.0.5)idptl(i)=111
  69       endif
  70       endif
  71
  72       ida=iabs(idptl(i))
  73       if(.not.(iappl.eq.7.and.i.eq.1))then
  74       if(mod(ndecay        ,10).eq.1
  75      *.and.ida.ne.0.and.ida.lt.10000)goto1000
  76       if(mod(ndecay/10     ,10).eq.1.and.ida.eq.  20)goto1000
  77       if(mod(ndecay/100    ,10).eq.1.and.ida.eq.2130)goto1000
  78       if(mod(ndecay/1000   ,10).eq.1.and.ida.eq.1130)goto1000
  79       if(mod(ndecay/1000   ,10).eq.1.and.ida.eq.2230)goto1000
  80       if(mod(ndecay/10000  ,10).eq.1.and.ida.eq.2330)goto1000
  81       if(mod(ndecay/10000  ,10).eq.1.and.ida.eq.1330)goto1000
  82       if(mod(ndecay/100000 ,10).eq.1.and.ida.eq.3331)goto1000
  83       if(mod(ndecay/1000000,10).eq.1.and.ida.eq. 110)goto1000
  84
  85       if(nrnody.gt.0)then
  86       do nod=1,nrnody
  87       if(idptl(i).eq.nody(nod))goto1000
  88       enddo
  89       endif
  90
  91       endif
  92
  93       call hdecay(i,iret)
  94       if(iret.eq.1)goto1000
  95       if(nptl.le.nptlb)then
  96         iret=-1
  97         goto 1000
  98       endif
  99
 100 2000  continue ! ---successful decay---
 101
 102       istptl(i)=1
 103       ifrptl(1,i)=nptlb+1
 104       ifrptl(2,i)=nptl
 105
 106       t=tivptl(2,i)
 107       x=xorptl(1,i)+(t-xorptl(4,i))*pptl(1,i)/pptl(4,i)
 108       y=xorptl(2,i)+(t-xorptl(4,i))*pptl(2,i)/pptl(4,i)
 109       z=xorptl(3,i)+(t-xorptl(4,i))*pptl(3,i)/pptl(4,i)
 110       call jtaux(t,z,ttaux)
 111       ttaus=ttaux
 112       if( ttaus.gt.0d0 ) then
 113         call jtauin
 114         call jtaus(z,ttest,sz)
 115         if (abs(t-ttest).gt.1e-5.and.ish.ge.1) then
 116           call utmsg('hdecas')
 117           write(ifch,*)'*****  t /= ttest'
 118           write(ifch,*)t,ttest,i,z,t,xorptl(3,i),xorptl(4,i)
 119      $         ,pptl(3,i),pptl(4,i)
 120           call utmsgf
 121         endif
 122       endif
 123
 124 c loop over decay products
 125
 126       do 20 n=nptlb+1,nptl
 127       iorptl(n)=i
 128       jorptl(n)=0
 129       istptl(n)=0
 130       ifrptl(1,n)=0
 131       ifrptl(2,n)=0
 132       if(idfptl(i).eq.0)then
 133        idfptl(n)=0
 134       else
 135        idfptl(n)=1
 136       endif
 137       rad=0
 138       phi=0
 139       ti=t
 140       zi=z
 141       xorptl(1,n)=x + rad*cos(phi)
 142       xorptl(2,n)=y + rad*sin(phi)
 143       xorptl(3,n)=zi
 144       xorptl(4,n)=ti
 145       io=n
 146 1     io=iorptl(io)
 147       if(ish.ge.4)write(ifch,*)'io = ',io,'  origin: ',iorptl(io)
 148       if(iorptl(io).gt.0)goto1
 149       if(ish.ge.4)write(ifch,*)'origin: ',io,idptl(io)
 150       zor=xorptl(3,io)
 151       tor=xorptl(4,io)
 152       call idquac(io,nq,ndummy1,ndummy2,jcdu)
 153       r=rangen()
 154       tauran=-taurea*alog(r)
 155       call jtaix(n,tauran,zor,tor,zis,tis)
 156       tivptl(1,n)=amax1(ti,tis)
 157       call idtau(idptl(n),pptl(4,n),pptl(5,n),taugm)
 158       r=rangen()
 159       tivptl(2,n)=t+taugm*(-alog(r))
 160       ityptl(n)=ityptl(i)
 161       radptl(n)=0.
 162       dezptl(n)=0.
 163       itsptl(n)=itsptl(i)
 164       rinptl(n)=rinptl(i)
 165 20    continue
 166
 167       if(iabs(idptl(nptlb+1)).le.6) then
 168         call gakli2(0,0)
 169         write (*,*) 'nptlb+1,nptl:',nptlb+1,nptl
 170         istptl(nptlb+1)=1
 171         do n=nptlb+2,nptl
 172           istptl(n)=20
 173         enddo
 174         call gakfra(iret)
 175         if(iret.eq.1)goto1000
 176         call gakli2(0,0)
 177       endif
 178
 179 1000  continue
 180       ttaus=ttauz
 181       call jtauin
 182       call utprix('hdecas',ish,ishini,5)
 183       return
 184       end
 185
 186 c-----------------------------------------------------------------------
 187       subroutine hdecay(ip,iret)
 188 c-----------------------------------------------------------------------
 189 c  decays particle ip from /cptl/
 190 c  for ip being no resonance: call hnbaaa
 191 c  for ip being resonance: standard resonance decay  procedure
 192 c-----------------------------------------------------------------------
 193       include 'epos.inc'
 194       double precision tpro,zpro,ttar,ztar,ttaus,detap,detat
 195       common/cttaus/tpro,zpro,ttar,ztar,ttaus,detap,detat
 196       common/wco/wmass2,wgam2
 197       parameter (mxlook=10000,mxdky=2000)
 198       common/dkytab/look(mxlook),cbr(mxdky),mode(5,mxdky)
 199       dimension pgen(5,10),rnd(10),u(3),beta(3)
 200      1     ,reduce(10)
 201       dimension prest(4,10),kno(10)
 202       data reduce/1.,1.,2.,5.,15.,60.,250.,1500.,1.2E4,1.2E5/
 203       data twome/1.022006e-3/
 204
 205 c          fctn definitions
 206       dot(i1,i2)=prest(4,i1)*prest(4,i2)-prest(1,i1)*prest(1,i2)
 207      *-prest(2,i1)*prest(2,i2)-prest(3,i1)*prest(3,i2)
 208 c          charged w propagator.
 209       wprop(z)=(z-wmass2**2)**2+(wmass2*wgam2)**2
 210
 211       call utpri('hdecay',ish,ishini,5)
 212
 213       ipp=ip
 214       iret=0
 215       nptlb=nptl
 216
 217       if(ish.ge.5)write(ifch,*)'ip,id,mass: ',ip,idptl(ip),pptl(5,ip)
 218 c     no k_long decay
 219 c     ---------------
 220 c     if(idptl(ip).eq.-20)goto1000
 221
 222 c     select decay mode
 223 c     -----------------
 224       ntry=0
 225 2     ntry=ntry+1
 226            if(ntry.gt.100)then
 227       if(ish.ge.1)then
 228       call utmsg('hdecay')
 229       write(ifch,*)'*****  decay not possible. iret = 1.'
 230       call utmsgf
 231       endif
 232       iret=1
 233       goto1000
 234            endif
 235       idlv1=idptl(ip)
 236       amss=pptl(5,ip)
 237
 238 c Decay of deuteron
 239
 240       if(idlv1.eq.17)then
 241         amss=1.01*amss
 242         naddptl=2
 243         call idmass(1120,amnew)
 244         pptl(5,nptl+1)=amnew
 245         idptl(nptl+1)=1120
 246         sum=amnew
 247         call idmass(1220,amnew)
 248         pptl(5,nptl+2)=amnew
 249         idptl(nptl+2)=1220
 250         sum=sum+amnew
 251         goto 111
 252       endif
 253
 254 c Decay of triton
 255
 256       if(idlv1.eq.18)then
 257         amss=1.01*amss
 258         naddptl=3
 259         call idmass(1120,amnew)
 260         pptl(5,nptl+1)=amnew
 261         idptl(nptl+1)=1120
 262         sum=amnew
 263         call idmass(1220,amnew)
 264         pptl(5,nptl+2)=amnew
 265         idptl(nptl+2)=1220
 266         sum=sum+amnew
 267         call idmass(1220,amnew)
 268         pptl(5,nptl+3)=amnew
 269         idptl(nptl+3)=1220
 270         sum=sum+amnew
 271          goto 111
 272       endif
 273
 274 c Decay of alpha
 275
 276       if(idlv1.eq.19)then
 277         amss=1.01*amss
 278         naddptl=4
 279         call idmass(1120,amnew)
 280         pptl(5,nptl+1)=amnew
 281         idptl(nptl+1)=1120
 282         sum=amnew
 283         call idmass(1220,amnew)
 284         pptl(5,nptl+2)=amnew
 285         idptl(nptl+2)=1220
 286         sum=sum+amnew
 287         call idmass(1120,amnew)
 288         pptl(5,nptl+3)=amnew
 289         idptl(nptl+3)=1120
 290         sum=sum+amnew
 291         call idmass(1220,amnew)
 292         pptl(5,nptl+4)=amnew
 293         idptl(nptl+4)=1220
 294         sum=sum+amnew
 295         goto 111
 296       endif
 297
 298 1     continue             !select one of the decay channel
 299       ipoint=look(iabs(idlv1))-1
 300       if(idlv1.eq.-20)ipoint=look(320)-1
 301       if(ipoint.lt.0) goto1000
 302       try=rangen()
 303 100   ipoint=ipoint+1
 304       if(ish.ge.4)write(ifch,*)'ipoint,cbr,try',ipoint,cbr(ipoint),try
 305       if(try.gt.cbr(ipoint)) goto100
 306       naddptl=0
 307       sum=0.
 308 c      nstart=nptl+1                  !?????????????????unused
 309       do 110 i=1,5         !store id and mass of products
 310         if(mode(i,ipoint).eq.0) goto 110
 311         if(nptl+naddptl+1.gt.mxptl) goto 9999
 312         if(iabs( mode(1,ipoint)) .le. 6.and.i.eq.2)then   !decay into quark ???
 313           call vedi(mode(1,ipoint),mode(2,ipoint),k3,idlv1)
 314           idptl(new)=idlv1
 315           call idmass(idlv1,amnew)
 316           pptl(5,new)=amnew
 317           sum=pptl(5,new)
 318         else                                 !decay into particles
 319           naddptl=naddptl+1
 320           new=nptl+naddptl
 321           idptl(new)=mode(i,ipoint)
 322           idlv1=idptl(new)
 323           call idmass(idlv1,pptl(5,new))
 324           sum=sum+pptl(5,new)
 325         endif
 326  110  continue
 327  111  continue
 328       if(naddptl.ne.1.and.sum.ge.amss)goto 2
 329  112  naddptl1=naddptl-1
 330       do 120 j=1,5
 331       pgen(j,1)=pptl(j,ip)
 332 120   continue
 333       pgen(5,1)=amss !needed because of deuteron, triton and alpha decay and OK
 334
 335       pgen(5,naddptl)=pptl(5,nptl+naddptl)
 336       if(naddptl.eq.1) goto 700            !one body decay
 337       if(naddptl.eq.2) goto 400            !two body decay
 338
 339       if(ish.ge.4)write(ifch,*)'>= 3 body decay'
 340
 341 c     use kroll-wada distribution for pi0 and eta dalitz decays.
 342 c     ----------------------------------------------
 343       if(.not.((idptl(ip).eq.110.or.idptl(ip).eq.220).and.
 344      1iabs(idptl(nptl+2)).eq.12)) goto 130
 345       ntry=0             !decay of pi0 or eta into electron
 346 125   ntry=ntry+1
 347            if(ntry.gt.10)then
 348       if(ish.ge. 0)then
 349       call utmsg('hdecay')
 350       write(ifch,*)'*****  ntry > 10. iret = 1.'
 351       write(ifch,*)'***** amee,ree,wtee',amee,ree,wtee
 352       call utmsgf
 353       endif
 354       iret=1
 355       goto1000
 356            endif
 357       amee=twome*(pptl(5,ip)/twome)**rangen()
 358       ree=(twome/amee)**2
 359       wtee=(1.-(amee/pptl(5,ip))**2)**3*sqrt(1.-ree)*(1.+.5*ree)
 360       if(wtee.lt.rangen()) goto125
 361       pgen(5,2)=amee
 362       goto400
 363 130   continue
 364
 365 c     calculate maximum phase-space weight
 366 c     ------------------------------------
 367       wtmax=1./reduce(naddptl)
 368       sum1=pgen(5,1)
 369       sum2=sum-pptl(5,nptl+1)
 370       do 200 i=1,naddptl1
 371       wtmax=wtmax*utpcm(sum1,sum2,pptl(5,nptl+i))
 372       sum1=sum1-pptl(5,nptl+i)
 373       sum2=sum2-pptl(5,nptl+i+1)
 374 200   continue
 375
 376 c     generate uniform naddptl-body phase space
 377 c     --------------------------------------
 378       ntry=0
 379 300   ntry=ntry+1
 380            if(ntry.gt.10000)then
 381       if(ish.ge. 0)then
 382       call utmsg('hdecay')
 383       write(ifch,*)'*****  infinite loop (2). iret = 1.'
 384       write(ifch,*)'***** ip,idptl(ip),pptl(5,ip):'
 385      *,ip,idptl(ip),pptl(5,ip)
 386       write(ifch,*)'***** wt,wtmax:',wt,wtmax
 387       write(ifch,*)'***** i,pgen(5,i),pptl(5,nptl+i),idptl(nptl+i):'
 388       do i=1,naddptl
 389       write(ifch,*)i,pgen(5,i),pptl(5,nptl+i),idptl(nptl+i)
 390       enddo
 391       call utmsgf
 392       endif
 393       iret=1
 394       goto1000
 395            endif
 396       rnd(1)=1.
 397       do 310 i=2,naddptl1
 398       rnew=rangen()
 399       i1=i-1
 400       do 320 jj1=1,i1
 401       j=i-jj1
 402       jsave=j+1
 403       if(rnew.le.rnd(j)) goto310
 404       rnd(jsave)=rnd(j)
 405 320   continue
 406 310   rnd(jsave)=rnew
 407       rnd(naddptl)=0.
 408       wt=1.
 409       sum1=sum
 410       do 330 i=2,naddptl
 411       sum1=sum1-pptl(5,nptl+i-1)
 412       pgen(5,i)=sum1+rnd(i)*(pgen(5,1)-sum)
 413       a=pgen(5,i-1)
 414       b=pgen(5,i)
 415       c=pptl(5,nptl+i-1)
 416       wt=wt*utpcm(a,b,c)
 417 330   continue
 418       if(wt.lt.rangen()*wtmax) goto300
 419
 420 c     carry out two-body decays in pgen frames
 421 c     ----------------------------------------
 422 400   continue
 423       if(ish.ge.4)write(ifch,*)'2 body decay'
 424       do 410 i=1,naddptl1
 425       qcm=utpcm(pgen(5,i),pgen(5,i+1),pptl(5,nptl+i))
 426       u(3)=2.*rangen()-1.
 427       phi=2.*pi*rangen()
 428       u(1)=sqrt(1.-u(3)**2)*cos(phi)
 429       u(2)=sqrt(1.-u(3)**2)*sin(phi)
 430       do 420 j=1,3
 431       pptl(j,nptl+i)=qcm*u(j)
 432       pgen(j,i+1)=-pptl(j,nptl+i)
 433 420   continue
 434       pptl(4,nptl+i)=sqrt(qcm**2+pptl(5,nptl+i)**2)
 435       pgen(4,i+1)=sqrt(qcm**2+pgen(5,i+1)**2)
 436 410   continue
 437       do 430 j=1,4
 438       pptl(j,nptl+naddptl)=pgen(j,naddptl)
 439 430   continue
 440
 441 c     boost pgen frames to lab frame
 442 c          also save momenta in rest frame (last frame)
 443 c     -------------------------------------------------
 444       do 500 ii=1,naddptl1
 445       i=naddptl-ii
 446       do 510 j=1,3
 447       beta(j)=pgen(j,i)/pgen(4,i)
 448 510   continue
 449       gamma=pgen(4,i)/pgen(5,i)
 450       do 520 k=i,naddptl
 451       k1=nptl+k
 452       bp=beta(1)*pptl(1,k1)+beta(2)*pptl(2,k1)+beta(3)*pptl(3,k1)
 453       do 530 j=1,3
 454       prest(j,k)=pptl(j,k1)
 455       pptl(j,k1)=pptl(j,k1)+gamma*beta(j)*(pptl(4,k1)
 456      1+bp*gamma/(gamma+1.))
 457 530   continue
 458       prest(4,k)=pptl(4,k1)
 459       pptl(4,k1)=gamma*(pptl(4,k1)+bp)
 460       if(pptl(4,k1).lt.1.d-5)then
 461         pptl(4,k1)=sqrt(pptl(1,k1)*pptl(1,k1)+pptl(2,k1)*pptl(2,k1)
 462      &                 +pptl(3,k1)*pptl(3,k1))
 463       endif
 464 520   continue
 465 500   continue
 466
 467 c     matrix elements
 468 c     ---------------
 469         if(iabs(idptl(ip)).eq.14)then                  !muon decay
 470           goto 650
 471         elseif(naddptl.eq.3)then
 472           if(idptl(ip).eq.221.or.idptl(ip).eq.331)then  !omeg and phi decay
 473             goto 610
 474           elseif(iabs(idptl(ip)).eq.130.or.       !Kl and K decay
 475      1       idptl(ip).eq.-20)then
 476             if(iabs(idptl(nptl+2)).lt.20)then   !semi-leptonic
 477               goto 630
 478             else                                !hadronic
 479               goto 640
 480             endif
 481           elseif(iabs(idptl(nptl+1)).lt.20.and. !other semi-leptonic decay
 482      1       idptl(nptl+1).ne.10)then
 483             goto 620
 484           elseif(iabs(idptl(nptl+2)).le.6)then
 485             goto 605            !decay into quark
 486           else
 487             goto 800
 488           endif
 489         else
 490          goto 800
 491         endif
 492
 493  605    wt=pptl(5,ip)*pptl(5,nptl+1)*dot(2,3)
 494         IF(wt.LT.rangen()*pptl(5,ip)**4/16.) goto 300
 495         ams=sqrt(dot(2,2)+dot(3,3)+2.*dot(2,3))
 496         kno(1)=idptl(nptl+2)
 497         kno(2)=idptl(nptl+3)
 498         if(ammin(kno(1),kno(2)).gt.ams)then
 499           call vedi(kno(1),kno(2),iddum,idlv2)
 500           idptl(nptl+2)=idlv2
 501           call idmass(idlv2,amnew2)
 502           pptl(5,nptl+2)=amnew2
 503           naddptl=2
 504           goto 112
 505         endif
 506 c......multiplicity
 507         PS =sqrt(dot(2,2))
 508         psq=sqrt(dot(3,3))
 509 c        PSP=PS                  !!???????????????unused
 510         np=0                    !!!!?????
 511         nq=2
 512         CNDE=4.5*LOG(MAX((ams-PS-PSQ)/0.7,1.1))
 513 c        IF(MMAT.EQ.12) CNDE=CNDE+PARJ(63)
 514  769    NTRY=NTRY+1
 515         IF(NTRY.GT.1000) THEN
 516           write(*,*)'hdecay caught in infinite loop'
 517           goto 1000
 518         ENDIF
 519         GAUSS=SQRT(-2.*CNDE*LOG(MAX(1E-10,rangen())))*
 520      &       SIN(2.*pi*rangen())
 521         ND=0.5+0.5*NP+0.25*NQ+CNDE+GAUSS
 522         IF(ND.LT.NP+NQ/2.OR.ND.LT.2.OR.ND.GT.10) GOTO 769
 523
 524
 525 c......choose hadrons
 526
 527
 528  132    kno(3)=kno(1)
 529         kno(4)=kno(2)
 530
 531         CONTINUE
 532         IF(ND.EQ.NP+NQ/2) GOTO 773
 533         DO I=nptl+2,nptl+2+nd-nq/2-1
 534           JT=2+1+INT((NQ-1) * rangen() )
 535           CALL vedi(kno(JT),0,KFL2,idlv3)
 536           idptl(i)=idlv3
 537 c          IF(K(I,2).EQ.0) GOTO 769
 538           kno(JT)=-KFL2
 539         enddo
 540  773    CONTINUE
 541         CALL vedi(kno(3),kno(4),KFLDMP,idlv4)
 542         idptl(nptl+2+nd-nq/2)=idlv4
 543         sum=0.
 544         do i=nptl+2,nptl+2+nd-nq/2
 545           call idmass(idptl(i),am)
 546           pptl(5,i)=am
 547           sum=sum+am
 548         enddo
 549         if(sum.gt.ams) goto 769
 550 c......goto phase space dis....
 551         ip=nptl+2+nd-nq/2+1
 552         do j=1,4
 553           pptl(j,ip)=pptl(j,ipp)-pptl(j,nptl+1)
 554         enddo
 555         pptl(5,ip)=ams
 556         idptl(ip)=sign(80,idptl(ipp))
 557         nptl=nptl+1
 558         naddptl=nd
 559         goto 112
 560
 561
 562 c     omeg and phi decay
 563 c          use vectors in rest frame
 564 c     ------------------------------
 565 610   wt=(pptl(5,nptl+1)*pptl(5,nptl+2)*pptl(5,nptl+3))**2
 566      1-(pptl(5,nptl+1)*dot(2,3))**2
 567      2-(pptl(5,nptl+2)*dot(1,3))**2
 568      3-(pptl(5,nptl+3)*dot(1,2))**2
 569      4+2.*dot(1,2)*dot(2,3)*dot(1,3)
 570       if(wt.lt.rangen()*pptl(5,ip)**6/108.) goto300
 571       goto800
 572
 573 c     semileptonic and quark decays
 574 c          use vectors in rest frame, where ip has (m,0,0,0)
 575 c          include w propagator
 576 c     ------------------------------------------------------
 577 620   wt=(pptl(5,ip)*prest(4,2))*dot(1,3)
 578       s12=pptl(5,nptl+1)**2+pptl(5,nptl+2)**2+2.*dot(1,2)
 579       s12max=pptl(5,ip)**2
 580       wt=wt*wprop(s12max)/wprop(s12)
 581       if(wt.lt.rangen()*pptl(5,ip)**4/16.) goto 300
 582       goto 800
 583
 584 c     semileptonic kaon decays
 585 c          use vectors in rest frame, where ip has (m,0,0,0)
 586 c          include form factor FML
 587 c     ------------------------------------------------------
 588 630   if(iabs(idptl(ip)).eq.130)then
 589         if(iabs(idptl(nptl+2)).eq.12)then
 590           ncha=1          !K   -> Pi0 + e + Nu
 591         else
 592           ncha=2          !K   -> Pi0 + Mu + Nu
 593         endif
 594       else
 595         if(iabs(idptl(nptl+2)).eq.12)then
 596           ncha=3          !K0  -> Pi + e + Nu
 597         else
 598           ncha=4          !K0  -> Pi + Mu + Nu
 599         endif
 600       endif
 601
 602       wt=FML(ncha,pptl(5,ip),pptl(5,nptl+1),pptl(5,nptl+2)
 603      &       ,prest(4,1),prest(4,2),prest(4,3))
 604       if(wt.lt.rangen()) goto 300
 605       goto 800
 606
 607 c     hadronic kaon decays
 608 c          use vectors in rest frame, where ip has (m,0,0,0)
 609 c          include form factor FM
 610 c     ------------------------------------------------------
 611 640   if(iabs(idptl(ip)).eq.130)then
 612         if(iabs(idptl(nptl+3)).eq.120)then
 613           ncha=1          !K   -> 3 Pi
 614         else
 615           ncha=2          !K   ->  Pi + 2 Pi0
 616         endif
 617       else
 618         if(iabs(idptl(nptl+1)).eq.110)then
 619           ncha=3          !K0  -> 3 Pi0
 620         else
 621           ncha=4          !K0  -> 2 Pi + Pi0
 622         endif
 623       endif
 624       S0=(pptl(5,ip)**2+pptl(5,nptl+1)**2+pptl(5,nptl+2)**2
 625      &   +pptl(5,nptl+3)**2)/3.d0
 626       S1=pptl(5,ip)**2+pptl(5,nptl+1)**2-2.*prest(4,1)*pptl(5,ip)
 627       S2=pptl(5,ip)**2+pptl(5,nptl+2)**2-2.*prest(4,2)*pptl(5,ip)
 628       S3=pptl(5,ip)**2+pptl(5,nptl+3)**2-2.*prest(4,3)*pptl(5,ip)
 629       wt=FM(ncha,S0,S1,S2,S3)
 630       if(wt.lt.rangen()) goto 300
 631       goto 800
 632
 633 c     muon decays
 634 c          use vectors in rest frame, where ip has (m,0,0,0)
 635 c          include form factor FMU
 636 c     ------------------------------------------------------
 637 650   xxx=2.*prest(4,1)/pptl(5,ip)            !reduced energy of electron
 638       if(xxx.gt.1.) goto 300
 639       wt=FMU(xxx)
 640       rrr=rangen()
 641       if(wt.lt.rrr) goto 300
 642       goto 800
 643
 644 c     one-particle decays
 645 c     -------------------
 646 700   continue
 647       do 710 j=1,5
 648       pptl(j,nptl+1)=pptl(j,ip)
 649 710   continue
 650
 651 c     swap particles and antiparticles if idptl(ip)<0
 652 c     -----------------------------------------------
 653  800    continue
 654         if(iabs(idptl(ip)).eq.80)then
 655           nptl=nptl-1
 656           naddptl=naddptl+1
 657         endif
 658         if(idptl(ipp).ge.0.or.iabs(idptl(ipp)).eq.20) goto 900
 659         do 810 i=1,naddptl
 660           idabs=iabs(idptl(nptl+i))
 661           ifl1=idabs/1000
 662           ifl2=mod(idabs/100,10)
 663           ifl3=mod(idabs/10,10)
 664           if(ifl1.eq.0.and.ifl2.ne.0.and.ifl2.eq.ifl3) goto 810
 665           if(idabs.eq.9.or.idabs.eq.10.or.idabs.eq.20) goto 810
 666           if(idabs.eq.29.or.idabs.eq.30.or.idabs.eq.40) goto 810
 667           idptl(nptl+i)=-idptl(nptl+i)
 668  810    continue
 669
 670  900    continue
 671         nptl=nptl+naddptl
 672         if(nptl.gt.mxptl)call utstop('hdecay: nptl>mxptl&')
 673 c        nqk=0           !???????????????????unused
 674         if(iabs(idptl(nptl)).lt.10.or.mod(idptl(nptl),100).eq.0)then
 675 c          call utstop('hdecay: decay ptcl is parton&')
 676         endif
 677
 678 c     print
 679 c     -----
 680
 681       if(ish.ge.3)then
 682       write(ifch,140)sngl(ttaus)
 683   140 format(/' ----------------------------'/
 684      *'    decay  at tau =',f6.2/
 685      *' ----------------------------')
 686       write(ifch,*)'decaying object:'
 687       call alist('&',ip,ip)
 688       write(ifch,*)'decay products:'
 689       call alist('&',nptlb+1,nptl)
 690       endif
 691       if(ish.ge.5)then
 692       write(ifch,*)'momentum sum:'
 693       do kk=1,5
 694       pptl(kk,nptl+1)=0
 695       do ii=nptlb+1,nptl
 696       pptl(kk,nptl+1)=pptl(kk,nptl+1)+pptl(kk,ii)
 697       enddo
 698       enddo
 699       call alist('&',nptl+1,nptl+1)
 700       endif
 701
 702 c     exit
 703 c     ----
 704
 705  1000 continue
 706       ip=ipp
 707       if(iret.ne.0.and.ish.ge.1)write(ifmt,'(a)')'hdecay: redo event'
 708       call utprix('hdecay',ish,ishini,5)
 709       return
 710
 711  9999   call utstop('hdecay: mxptl too small&')
 712         end
 713
 714 c---------------------------------------------------------------------
 715       subroutine vedi(k1,k2,k3,id)
 716 c---------------------------------------------------------------------
 717       include 'epos.inc'
 718       if(k2.eq.0)then
 719         if(rangen().lt.pdiqua.and.iabs(k1).lt.6)then
 720           ifl1=int(rangen()/pud)+1
 721           ifl2=int(rangen()/pud)+1
 722           k3=-min(ifl1,ifl2)*1000-max(ifl1,ifl2)*100
 723         else
 724           k3=int(rangen()/pud)+1
 725         endif
 726         if(k1.gt.0.and.k1.le.6)k3=-k3
 727         if(k1.lt.-1000)k3=-k3
 728       else
 729         k3=k2
 730       endif
 731       id=idsp(k1,k3)
 732       if(iabs(id).le.999) then
 733         ids=max(mod(iabs(id)/100,10),mod(iabs(id)/10,10))
 734         if(ids.le.2)then
 735           idr=sign(iabs(id)+int(rangen()+0.5),id)
 736         elseif(ids.eq.3)then
 737           idr=sign(iabs(id)+int(rangen()+0.6),id)
 738         else
 739           idr=sign(iabs(id)+int(rangen()+0.75),id)
 740         endif
 741       else
 742         idr=sign(iabs(id)+int(0.5+rangen()),id)
 743       endif
 744       id=idr
 745       if(ish.ge.5)write(ifch,*) 'Flavor:',k1,k2,k3,id
 746       end
 747
 748 c-----------------------------------------------------------------------
 749       subroutine hdecin(lprint)
 750 c-----------------------------------------------------------------------
 751 c     sets up /dkytab/
 752 c-----------------------------------------------------------------------
 753       include 'epos.inc'
 754       common/wco/wmass2,wgam2
 755       dimension imode(6)
 756       character*8 idlabl,lmode(6),lres
 757       character*8 iblank
 758       logical lprint
 759       parameter (mxlook=10000,mxdky=2000)
 760       common/dkytab/look(mxlook),cbr(mxdky),mode(5,mxdky)
 761       common/nodcay/nodcay,noeta,nopi0,nonunu,noevol,nohadr
 762       logical nodcay,noeta,nopi0,nonunu,noevol,nohadr
 763       parameter (ndectb=1192)
 764       real dectab(7,ndectb)
 765
 766       data ((dectab(i,j),i=1,7),j=  1, 18)/
 767      *  110., .98850,  10.,  10.,   0.,   0.,   0.
 768      *, 110.,1.00000,  10.,  12., -12.,   0.,   0.
 769      *, 220., .38000,  10.,  10.,   0.,   0.,   0.
 770      *, 220., .71000, 110., 110., 110.,   0.,   0.
 771      *, 220., .94600, 120.,-120., 110.,   0.,   0.
 772      *, 220., .99500, 120.,-120.,  10.,   0.,   0.
 773      *, 220.,1.00000,  10.,  12., -12.,   0.,   0.
 774      *, 330., .44100, 220., 120.,-120.,   0.,   0.
 775      *, 330., .66100, 220., 110., 110.,   0.,   0.
 776      *, 330., .95900, 111.,  10.,   0.,   0.,   0.
 777      *, 330., .98000, 221.,  10.,   0.,   0.,   0.
 778      *, 330.,1.00000,  10.,  10.,   0.,   0.,   0.
 779      *, 121.,1.00000, 120., 110.,   0.,   0.,   0.
 780      *, 111., .99989, 120.,-120.,   0.,   0.,   0.
 781      *, 111., .99993,  12., -12.,   0.,   0.,   0.
 782      *, 111.,1.00000,  14., -14.,   0.,   0.,   0.
 783      *, 221., .89900, 120.,-120., 110.,   0.,   0.
 784      *, 221., .91200, 120.,-120.,   0.,   0.,   0./
 785       data ((dectab(i,j),i=1,7),j= 19, 36)/
 786      *  221., .99992, 110.,  10.,   0.,   0.,   0.
 787      *, 221.,1.00000,  12., -12.,   0.,   0.,   0.
 788      *, 331., .48600, 130.,-130.,   0.,   0.,   0.
 789      *, 331., .83700, 230.,-230.,   0.,   0.,   0.
 790      *, 331., .98400, 120.,-120., 110.,   0.,   0.
 791      *, 331., .99944, 220.,  10.,   0.,   0.,   0.
 792      *, 331., .99975,  12., -12.,   0.,   0.,   0.
 793      *, 331.,1.00000,  14., -14.,   0.,   0.,   0.
 794      *, 230., .50000,  20.,   0.,   0.,   0.,   0.
 795      *, 230.,1.00000, -20.,   0.,   0.,   0.,   0.
 796      *, 131., .66670, 230., 120.,   0.,   0.,   0.
 797      *, 131.,1.00000, 130., 110.,   0.,   0.,   0.
 798      *, 231., .66670, 130.,-120.,   0.,   0.,   0.
 799      *, 231.,1.00000, 230., 110.,   0.,   0.,   0.
 800      *, 240., .11000,  12., -11., 230.,   0.,   0.
 801      *, 240., .17000,  12., -11., 231.,   0.,   0.
 802      *, 240., .28000,  14., -13., 230.,   0.,   0.
 803      *, 240., .34000,  14., -13., 231.,   0.,   0./
 804       data ((dectab(i,j),i=1,7),j= 37, 54)/
 805      *  240., .37800, 230.,-120.,   0.,   0.,   0.
 806      *, 240., .56300, 230.,-121.,   0.,   0.,   0.
 807      *, 240., .60800, 231.,-120.,   0.,   0.,   0.
 808      *, 240., .62100, 230.,-120., 110.,   0.,   0.
 809      *, 240., .71000, 130.,-120.,-120.,   0.,   0.
 810      *, 240., .80100, 230.,-120.,-120., 120.,   0.
 811      *, 240., .87900, 130.,-120.,-120., 110.,   0.
 812      *, 240., .95400, 230.,-120., 110., 110.,   0.
 813      *, 240., .96600, 230.,-130.,   0.,   0.,   0.
 814      *, 240., .97600, 331.,-120.,   0.,   0.,   0.
 815      *, 240., .98800,-130., 231.,   0.,   0.,   0.
 816      *, 240.,1.00000,-131., 230.,   0.,   0.,   0.
 817      *, 140., .04500, -12.,  11., 130.,   0.,   0.
 818      *, 140., .07500, -12.,  11., 131.,   0.,   0.
 819      *, 140., .12000, -14.,  13., 130.,   0.,   0.
 820      *, 140., .15000, -14.,  13., 131.,   0.,   0.
 821      *, 140., .20300, 130.,-120.,   0.,   0.,   0.
 822      *, 140., .22700, 230., 110.,   0.,   0.,   0./
 823       data ((dectab(i,j),i=1,7),j= 55, 72)/
 824      *  140., .24700, 230., 220.,   0.,   0.,   0.
 825      *, 140., .28900, 230., 221.,   0.,   0.,   0.
 826      *, 140., .45100, 130.,-121.,   0.,   0.,   0.
 827      *, 140., .53600, 131.,-120.,   0.,   0.,   0.
 828      *, 140., .56200, 231., 110.,   0.,   0.,   0.
 829      *, 140., .57600, 230., 111.,   0.,   0.,   0.
 830      *, 140., .58700, 130.,-120., 110.,   0.,   0.
 831      *, 140., .60300, 230.,-120., 120.,   0.,   0.
 832      *, 140., .72700, 130.,-120.,-120., 120.,   0.
 833      *, 140., .87600, 230.,-120., 120., 110.,   0.
 834      *, 140., .96900, 130.,-120., 110., 110.,   0.
 835      *, 140.,1.00000, 230., 110., 110., 110.,   0.
 836      *, 340., .03250,  12., -11., 220.,   0.,   0.
 837      *, 340., .06500,  12., -11., 331.,   0.,   0.
 838      *, 340., .09750,  14., -13., 220.,   0.,   0.
 839      *, 340., .13000,  14., -13., 331.,   0.,   0.
 840      *, 340., .17900,-130., 230.,   0.,   0.,   0.
 841      *, 340., .22800,-120., 220.,   0.,   0.,   0./
 842       data ((dectab(i,j),i=1,7),j= 73, 90)/
 843      *  340., .33800,-131., 230.,   0.,   0.,   0.
 844      *, 340., .44800,-130., 231.,   0.,   0.,   0.
 845      *, 340., .55800,-120., 331.,   0.,   0.,   0.
 846      *, 340., .57500,-130., 230., 110.,   0.,   0.
 847      *, 340., .59200,-230., 230.,-120.,   0.,   0.
 848      *, 340., .69400,-130., 230.,-120., 120.,   0.
 849      *, 340., .79600,-130., 230., 110., 110.,   0.
 850      *, 340., .89800,-130., 130.,-120., 110.,   0.
 851      *, 340.,1.00000,-230., 230.,-120., 110.,   0.
 852      *, 241., .64000, 140.,-120.,   0.,   0.,   0.
 853      *, 241., .92000, 240., 110.,   0.,   0.,   0.
 854      *, 241.,1.00000, 240.,  10.,   0.,   0.,   0.
 855      *, 141., .55000, 140., 110.,   0.,   0.,   0.
 856      *, 141.,1.00000, 140.,  10.,   0.,   0.,   0.
 857      *, 341.,1.00000, 340.,  10.,   0.,   0.,   0.
 858      *, 441., .07400,  12., -12.,   0.,   0.,   0.
 859      *, 441., .14800,  14., -14.,   0.,   0.,   0.
 860      *, 441., .15210,-121., 120.,   0.,   0.,   0./
 861       data ((dectab(i,j),i=1,7),j= 91,108)/
 862      *  441., .15620, 111., 110.,   0.,   0.,   0.
 863      *, 441., .16020, 121.,-120.,   0.,   0.,   0.
 864      *, 441., .16300,-121., 111., 120.,   0.,   0.
 865      *, 441., .16580, 121.,-121., 110.,   0.,   0.
 866      *, 441., .16860, 121., 111.,-120.,   0.,   0.
 867      *, 441., .28740, 120.,-120., 130.,-130.,   0.
 868      *, 441., .40620, 110., 110., 130.,-130.,   0.
 869      *, 441., .52500, 120.,-120., 120.,-120.,   0.
 870      *, 441., .64380, 120.,-120., 110., 110.,   0.
 871      *, 441., .76260, 110., 110., 110., 110.,   0.
 872      *, 441., .88130, 120.,-120., 230.,-230.,   0.
 873      *, 441.,1.00000, 110., 110., 230., 230.,   0.
 874      *, 150., .06000, -12.,  11., 140.,   0.,   0.
 875      *, 150., .12000, -12.,  11., 141.,   0.,   0.
 876      *, 150., .18000, -14.,  13., 140.,   0.,   0.
 877      *, 150., .24000, -14.,  13., 141.,   0.,   0.
 878      *, 150., .25500, -16.,  15., 140.,   0.,   0.
 879      *, 150., .27000, -16.,  15., 141.,   0.,   0./
 880       data ((dectab(i,j),i=1,7),j=109,122)/
 881      *  150., .28050, 140., 120.,   0.,   0.,   0.
 882      *, 150., .29100, 140., 121.,   0.,   0.,   0.
 883      *, 150., .30150, 141., 120.,   0.,   0.,   0.
 884      *, 150., .31200, 141., 121.,   0.,   0.,   0.
 885      *, 150., .32650, 140.,-340.,   0.,   0.,   0.
 886      *, 150., .34100, 140.,-341.,   0.,   0.,   0.
 887      *, 150., .35550, 141.,-340.,   0.,   0.,   0.
 888      *, 150., .37000, 141.,-341.,   0.,   0.,   0.
 889      *, 150., 0.820  ,   1.,  -4.,   1.,  -2.,   0.
 890      *, 150., 0.920  ,   1.,  -2.,   1.,  -4.,   0.
 891      *, 150., 0.975  ,   1.,  -4.,   4.,  -3.,   0.
 892      *, 150., 0.985  ,   1.,  -3.,   4.,  -4.,   0.
 893      *, 150., 0.995  ,   1.,  -1.,   1.,  -2.,   0.
 894      *, 150., 1.     ,   1.,  -1.,   4.,  -3.,   0./
 895       data ((dectab(i,j),i=1,7),j=123,142)
 896      */ 250., .06000, -12.,  11., 240.,   0.,   0.
 897      *, 250., .12000, -12.,  11., 241.,   0.,   0.
 898      *, 250., .18000, -14.,  13., 240.,   0.,   0.
 899      *, 250., .24000, -14.,  13., 241.,   0.,   0.
 900      *, 250., .25500, -16.,  15., 240.,   0.,   0.
 901      *, 250., .27000, -16.,  15., 241.,   0.,   0.
 902      *, 250., .28050, 240., 120.,   0.,   0.,   0.
 903      *, 250., .29100, 240., 121.,   0.,   0.,   0.
 904      *, 250., .30150, 241., 120.,   0.,   0.,   0.
 905      *, 250., .31200, 241., 121.,   0.,   0.,   0.
 906      *, 250., .32650, 240.,-340.,   0.,   0.,   0.
 907      *, 250., .34100, 240.,-341.,   0.,   0.,   0.
 908      *, 250., .35550, 241.,-340.,   0.,   0.,   0.
 909      *, 250., .37000, 241.,-341.,   0.,   0.,   0.
 910      *, 250., 0.820  ,   2.,  -4.,   1.,  -2.,   0.
 911      *, 250., 0.920  ,   2.,  -2.,   1.,  -4.,   0.
 912      *, 250., 0.975  ,   2.,  -4.,   4.,  -3.,   0.
 913      *, 250., 0.985  ,   2.,  -3.,   4.,  -4.,   0.
 914      *, 250., 0.995  ,   2.,  -1.,   1.,  -2.,   0.
 915      *, 250., 1.     ,   2.,  -1.,   4.,  -3.,   0./
 916       data ((dectab(i,j),i=1,7),j=143,176)/
 917      *     238*1. /
 918       data ((dectab(i,j),i=1,7),j=177,190)
 919      * /350., .06000,  12., -11., 340.,   0.,   0.
 920      *, 350., .12000,  12., -11., 341.,   0.,   0.
 921      *, 350., .18000,  14., -13., 340.,   0.,   0.
 922      *, 350., .24000,  14., -13., 341.,   0.,   0.
 923      *, 350., .25500,  16., -15., 340.,   0.,   0.
 924      *, 350., .27000,  16., -15., 341.,   0.,   0.
 925      *, 350., .28050, 340., 120.,   0.,   0.,   0.
 926      *, 350., .29100, 340., 121.,   0.,   0.,   0.
 927      *, 350., .30150, 341., 120.,   0.,   0.,   0.
 928      *, 350., .31200, 341., 121.,   0.,   0.,   0.
 929      *, 350., .32650, 340.,-340.,   0.,   0.,   0.
 930      *, 350., .34100, 340.,-341.,   0.,   0.,   0.
 931      *, 350., .35550, 341.,-340.,   0.,   0.,   0.
 932      *, 350., .37000, 341.,-341.,   0.,   0.,   0./
 933       data ((dectab(i,j),i=1,7),j=191,196)/
 934      *  350., 0.820  ,   3.,  -4.,   1.,  -2.,   0.
 935      *, 350., 0.920  ,   3.,  -2.,   1.,  -4.,   0.
 936      *, 350., 0.975  ,   3.,  -4.,   4.,  -3.,   0.
 937      *, 350., 0.985  ,   3.,  -3.,   4.,  -4.,   0.
 938      *, 350., 0.995  ,   3.,  -1.,   1.,  -2.,   0.
 939      *, 350., 1.     ,   3.,  -1.,   4.,  -3.,   0./
 940       data ((dectab(i,j),i=1,7),j=197,244)/
 941      *     336*1. /
 942       data ((dectab(i,j),i=1,7),j=245,262)/
 943      *  160., .33330,  -1.,   2.,  -5.,   0.,   0.
 944      *, 160., .66660,  -4.,   3.,  -5.,   0.,   0.
 945      *, 160., .77770,  11., -12.,  -5.,   0.,   0.
 946      *, 160., .88880,  13., -14.,  -5.,   0.,   0.
 947      *, 160.,1.00000, -15.,  16.,  -5.,   0.,   0.
 948      *, 260., .33330,  -1.,   2.,  -5.,   0.,   0.
 949      *, 260., .66660,  -4.,   3.,  -5.,   0.,   0.
 950      *, 260., .77770, -11.,  12.,  -5.,   0.,   0.
 951      *, 260., .88880, -13.,  14.,  -5.,   0.,   0.
 952      *, 260.,1.00000, -15.,  16.,  -5.,   0.,   0.
 953      *, 360., .33330,  -1.,   2.,  -5.,   0.,   0.
 954      *, 360., .66660,  -4.,   3.,  -5.,   0.,   0.
 955      *, 360., .77770, -11.,  12.,  -5.,   0.,   0.
 956      *, 360., .88880, -13.,  14.,  -5.,   0.,   0.
 957      *, 360.,1.00000, -15.,  16.,  -5.,   0.,   0.
 958      *, 151.,1.00000, 150.,  10.,   0.,   0.,   0.
 959      *, 251.,1.00000, 250.,  10.,   0.,   0.,   0.
 960      *, 351.,1.00000, 350.,  10.,   0.,   0.,   0./
 961       data ((dectab(i,j),i=1,7),j=263,280)/
 962      *  161.,1.00000, 160.,  10.,   0.,   0.,   0.
 963      *, 261.,1.00000, 260.,  10.,   0.,   0.,   0.
 964      *, 361.,1.00000, 360.,  10.,   0.,   0.,   0.
 965      *,1230.,1.00000,2130.,  10.,   0.,   0.,   0.
 966      *,1111.,1.00000,1120., 120.,   0.,   0.,   0.
 967      *,1121., .66670,1120., 110.,   0.,   0.,   0.
 968      *,1121.,1.00000,1220., 120.,   0.,   0.,   0.
 969      *,1221., .66670,1220., 110.,   0.,   0.,   0.
 970      *,1221.,1.00000,1120.,-120.,   0.,   0.,   0.
 971      *,2221.,1.00000,1220.,-120.,   0.,   0.,   0.
 972      *,1131., .88000,2130., 120.,   0.,   0.,   0.
 973      *,1131., .94000,1130., 110.,   0.,   0.,   0.
 974      *,1131.,1.00000,1230., 120.,   0.,   0.,   0.
 975      *,1231., .88000,2130., 110.,   0.,   0.,   0.
 976      *,1231., .94000,1130.,-120.,   0.,   0.,   0.
 977      *,1231.,1.00000,2230., 120.,   0.,   0.,   0.
 978      *,2231., .88000,2130.,-120.,   0.,   0.,   0.
 979      *,2231., .94000,1230.,-120.,   0.,   0.,   0./
 980       data ((dectab(i,j),i=1,7),j=281,298)/
 981      * 2231.,1.00000,2230., 110.,   0.,   0.,   0.
 982      *,1331., .66670,2330., 120.,   0.,   0.,   0.
 983      *,1331.,1.00000,1330., 110.,   0.,   0.,   0.
 984      *,2331., .66670,1330.,-120.,   0.,   0.,   0.
 985      *,2331.,1.00000,2330., 110.,   0.,   0.,   0.
 986      *,  16., .18000,  12., -11.,  15.,   0.,   0.
 987      *,  16., .36000,  14., -13.,  15.,   0.,   0.
 988      *,  16., .45100,-120.,  15.,   0.,   0.,   0.
 989      *,  16., .66000,-121.,  15.,   0.,   0.,   0.
 990      *,  16., .78000, 110., 110.,-120.,  15.,   0.
 991      *,  16., .83600, 120.,-120.,-120.,  15.,   0.
 992      *,  16.,1.00000, 120., 110.,-120.,-120.,  15.
 993      *,2140., .03750, -12.,  11.,2130.,   0.,   0.
 994      *,2140., .07500, -12.,  11.,1231.,   0.,   0.
 995      *,2140., .11250, -14.,  13.,2130.,   0.,   0.
 996      *,2140., .15000, -14.,  13.,1231.,   0.,   0.
 997      *,2140., .18200,2130., 120.,   0.,   0.,   0.
 998      *,2140., .21300,1230., 120.,   0.,   0.,   0./
 999       data ((dectab(i,j),i=1,7),j=299,316)/
1000      * 2140., .24400,1120.,-230.,   0.,   0.,   0.
1001      *,2140., .29500,1131., 110.,   0.,   0.,   0.
1002      *,2140., .34600,1231., 120.,   0.,   0.,   0.
1003      *,2140., .39700,1121.,-230.,   0.,   0.,   0.
1004      *,2140., .44800,1111.,-130.,   0.,   0.,   0.
1005      *,2140., .49900,1130., 111.,   0.,   0.,   0.
1006      *,2140., .55000,1230., 121.,   0.,   0.,   0.
1007      *,2140., .60100,1120.,-231.,   0.,   0.,   0.
1008      *,2140., .65800,1120.,-230., 120.,-120.,   0.
1009      *,2140., .71500,1120.,-230., 110., 110.,   0.
1010      *,2140., .77200,1120.,-130., 120., 110.,   0.
1011      *,2140., .82900,1220.,-230., 120., 110.,   0.
1012      *,2140., .88600,1220.,-130., 120., 120.,   0.
1013      *,2140., .94300,2130., 120., 120.,-120.,   0.
1014      *,2140.,1.00000,2130., 120., 110., 110.,   0.
1015      *,1140.,1.00000,2140., 120.,   0.,   0.,   0.
1016      *,1240.,1.00000,2140., 110.,   0.,   0.,   0.
1017      *,2240.,1.00000,2140.,-120.,   0.,   0.,   0./
1018       data ((dectab(i,j),i=1,7),j=317,334)/
1019      * 1340., .03750, -12.,  11.,1330.,   0.,   0.
1020      *,1340., .07500, -12.,  11.,1331.,   0.,   0.
1021      *,1340., .11250, -14.,  13.,1330.,   0.,   0.
1022      *,1340., .15000, -14.,  13.,1331.,   0.,   0.
1023      *,1340., .19900,1330., 120.,   0.,   0.,   0.
1024      *,1340., .24800,1231., 130.,   0.,   0.,   0.
1025      *,1340., .28800,1330., 120.,   0.,   0.,   0.
1026      *,1340., .32800,1131.,-230.,   0.,   0.,   0.
1027      *,1340., .36800,1330., 121.,   0.,   0.,   0.
1028      *,1340., .40800,1130.,-230.,   0.,   0.,   0.
1029      *,1340., .44800,1330., 120., 110.,   0.,   0.
1030      *,1340., .48800,2330., 120., 120.,   0.,   0.
1031      *,1340., .52800,1130.,-130., 120.,   0.,   0.
1032      *,1340., .56800,1130.,-230., 110.,   0.,   0.
1033      *,1340., .60800,1230.,-230., 120.,   0.,   0.
1034      *,1340., .66400,2130.,-230., 120., 110.,   0.
1035      *,1340., .72000,2130.,-130., 120., 120.,   0.
1036      *,1340., .77600,1130.,-230., 120.,-120.,   0./
1037       data ((dectab(i,j),i=1,7),j=335,352)/
1038      * 1340., .83200,1130.,-230., 110., 110.,   0.
1039      *,1340., .88800,1330., 120., 120.,-120.,   0.
1040      *,1340., .94400,1330., 120., 110., 110.,   0.
1041      *,1340.,1.00000,2330., 120., 120., 110.,   0.
1042      *,3140., .03750, -12.,  11.,1330.,   0.,   0.
1043      *,3140., .07500, -12.,  11.,1331.,   0.,   0.
1044      *,3140., .11250, -14.,  13.,1330.,   0.,   0.
1045      *,3140., .15000, -14.,  13.,1331.,   0.,   0.
1046      *,3140., .19900,1330., 120.,   0.,   0.,   0.
1047      *,3140., .24800,1231., 130.,   0.,   0.,   0.
1048      *,3140., .28800,1330., 120.,   0.,   0.,   0.
1049      *,3140., .32800,1131.,-230.,   0.,   0.,   0.
1050      *,3140., .36800,1330., 121.,   0.,   0.,   0.
1051      *,3140., .40800,1130.,-230.,   0.,   0.,   0.
1052      *,3140., .44800,1330., 120., 110.,   0.,   0.
1053      *,3140., .48800,2330., 120., 120.,   0.,   0.
1054      *,3140., .52800,1130.,-130., 120.,   0.,   0.
1055      *,3140., .56800,1130.,-230., 110.,   0.,   0./
1056       data ((dectab(i,j),i=1,7),j=353,370)/
1057      * 3140., .60800,1230.,-230., 120.,   0.,   0.
1058      *,3140., .66400,2130.,-230., 120., 110.,   0.
1059      *,3140., .72000,2130.,-130., 120., 120.,   0.
1060      *,3140., .77600,1130.,-230., 120.,-120.,   0.
1061      *,3140., .83200,1130.,-230., 110., 110.,   0.
1062      *,3140., .88800,1330., 120., 120.,-120.,   0.
1063      *,3140., .94400,1330., 120., 110., 110.,   0.
1064      *,3140.,1.00000,2330., 120., 120., 110.,   0.
1065      *,2340., .03750, -12.,  11.,2330.,   0.,   0.
1066      *,2340., .07500, -12.,  11.,2331.,   0.,   0.
1067      *,2340., .11250, -14.,  13.,2330.,   0.,   0.
1068      *,2340., .15000, -14.,  13.,2331.,   0.,   0.
1069      *,2340., .17500,2330., 120.,   0.,   0.,   0.
1070      *,2340., .20000,1330., 110.,   0.,   0.,   0.
1071      *,2340., .22500,1130.,-130.,   0.,   0.,   0.
1072      *,2340., .25000,1230.,-230.,   0.,   0.,   0.
1073      *,2340., .29500,2331., 120.,   0.,   0.,   0.
1074      *,2340., .34000,1331., 110.,   0.,   0.,   0./
1075       data ((dectab(i,j),i=1,7),j=371,388)/
1076      * 2340., .38500,1131.,-130.,   0.,   0.,   0.
1077      *,2340., .43000,1231.,-230.,   0.,   0.,   0.
1078      *,2340., .47500,2330., 121.,   0.,   0.,   0.
1079      *,2340., .52000,1330., 111.,   0.,   0.,   0.
1080      *,2340., .56500,1130.,-131.,   0.,   0.,   0.
1081      *,2340., .61000,1230.,-231.,   0.,   0.,   0.
1082      *,2340., .64900,2130.,-230., 120.,-120.,   0.
1083      *,2340., .68800,2130.,-230., 110., 110.,   0.
1084      *,2340., .72700,2130.,-130., 120., 110.,   0.
1085      *,2340., .76600,1130.,-230.,-120., 110.,   0.
1086      *,2340., .80500,1130.,-130., 120.,-120.,   0.
1087      *,2340., .84400,1130.,-130., 110., 110.,   0.
1088      *,2340., .88300,1330., 120.,-120., 110.,   0.
1089      *,2340., .92200,1330., 110., 110., 110.,   0.
1090      *,2340., .96100,2330., 120., 120.,-120.,   0.
1091      *,2340.,1.00000,2330., 120., 110., 110.,   0.
1092      *,3240., .03750, -12.,  11.,2330.,   0.,   0.
1093      *,3240., .07500, -12.,  11.,2331.,   0.,   0./
1094       data ((dectab(i,j),i=1,7),j=389,406)/
1095      * 3240., .11250, -14.,  13.,2330.,   0.,   0.
1096      *,3240., .15000, -14.,  13.,2331.,   0.,   0.
1097      *,3240., .17500,2330., 120.,   0.,   0.,   0.
1098      *,3240., .20000,1330., 110.,   0.,   0.,   0.
1099      *,3240., .22500,1130.,-130.,   0.,   0.,   0.
1100      *,3240., .25000,1230.,-230.,   0.,   0.,   0.
1101      *,3240., .29500,2331., 120.,   0.,   0.,   0.
1102      *,3240., .34000,1331., 110.,   0.,   0.,   0.
1103      *,3240., .38500,1131.,-130.,   0.,   0.,   0.
1104      *,3240., .43000,1231.,-230.,   0.,   0.,   0.
1105      *,3240., .47500,2330., 121.,   0.,   0.,   0.
1106      *,3240., .52000,1330., 111.,   0.,   0.,   0.
1107      *,3240., .56500,1130.,-131.,   0.,   0.,   0.
1108      *,3240., .61000,1230.,-231.,   0.,   0.,   0.
1109      *,3240., .64900,2130.,-230., 120.,-120.,   0.
1110      *,3240., .68800,2130.,-230., 110., 110.,   0.
1111      *,3240., .72700,2130.,-130., 120., 110.,   0.
1112      *,3240., .76600,1130.,-230.,-120., 110.,   0./
1113       data ((dectab(i,j),i=1,7),j=407,424)/
1114      * 3240., .80500,1130.,-130., 120.,-120.,   0.
1115      *,3240., .84400,1130.,-130., 110., 110.,   0.
1116      *,3240., .88300,1330., 120.,-120., 110.,   0.
1117      *,3240., .92200,1330., 110., 110., 110.,   0.
1118      *,3240., .96100,2330., 120., 120.,-120.,   0.
1119      *,3240.,1.00000,2330., 120., 110., 110.,   0.
1120      *,3340., .07500, -12.,  11.,3331.,   0.,   0.
1121      *,3340., .15000, -12.,  11.,3331.,   0.,   0.
1122      *,3340., .25000,1330.,-230.,   0.,   0.,   0.
1123      *,3340., .31000,3331., 120.,   0.,   0.,   0.
1124      *,3340., .37000,1331.,-230.,   0.,   0.,   0.
1125      *,3340., .43000,1330.,-231.,   0.,   0.,   0.
1126      *,3340., .49000,2330.,-230., 120.,   0.,   0.
1127      *,3340., .55000,1330.,-230., 110.,   0.,   0.
1128      *,3340., .61000,1330.,-130., 120.,   0.,   0.
1129      *,3340., .67500,3331., 120., 120.,-120.,   0.
1130      *,3340., .74000,3331., 120., 110., 110.,   0.
1131      *,3340., .80500,1330.,-230., 120.,-120.,   0./
1132       data ((dectab(i,j),i=1,7),j=425,442)/
1133      * 3340., .87000,1330.,-230., 110., 110.,   0.
1134      *,3340., .93500,2330.,-230., 120., 110.,   0.
1135      *,3340.,1.00000,2330.,-130., 120., 120.,   0.
1136      *,1141.,1.00000,2140., 120.,   0.,   0.,   0.
1137      *,1241.,1.00000,2140., 110.,   0.,   0.,   0.
1138      *,2241.,1.00000,2140.,-120.,   0.,   0.,   0.
1139      *,1341., .66670,2340., 120.,   0.,   0.,   0.
1140      *,1341.,1.00000,1340., 110.,   0.,   0.,   0.
1141      *,2341., .66670,1340.,-120.,   0.,   0.,   0.
1142      *,2341.,1.00000,2340., 110.,   0.,   0.,   0.
1143      *,3341.,1.00000,3340., 110.,   0.,   0.,   0.
1144      *,1150., .06000,  12., -11.,1140.,   0.,   0.
1145      *,1150., .12000,  12., -11.,1141.,   0.,   0.
1146      *,1150., .18000,  14., -13.,1140.,   0.,   0.
1147      *,1150., .24000,  14., -13.,1141.,   0.,   0.
1148      *,1150., .25500,  16., -15.,1140.,   0.,   0.
1149      *,1150., .27000,  16., -15.,1141.,   0.,   0.
1150      *,1150., .28925,1140.,-120.,   0.,   0.,   0./
1151       data ((dectab(i,j),i=1,7),j=443,460)/
1152      * 1150., .30850,1140.,-121.,   0.,   0.,   0.
1153      *,1150., .32775,1141.,-120.,   0.,   0.,   0.
1154      *,1150., .34700,1141.,-121.,   0.,   0.,   0.
1155      *,1150., .35775,1140., 340.,   0.,   0.,   0.
1156      *,1150., .36850,1140., 341.,   0.,   0.,   0.
1157      *,1150., .37925,1141., 340.,   0.,   0.,   0.
1158      *,1150., .39000,1141., 341.,   0.,   0.,   0.
1159      *,1150., .42050,1140.,-120., 110.,   0.,   0.
1160      *,1150., .45100,1140.,-120., 220.,   0.,   0.
1161      *,1150., .48150,1140.,-120., 111.,   0.,   0.
1162      *,1150., .51200,1140.,-120., 221.,   0.,   0.
1163      *,1150., .54250,1140.,-121., 110.,   0.,   0.
1164      *,1150., .57300,1140.,-121., 220.,   0.,   0.
1165      *,1150., .60350,1140.,-121., 111.,   0.,   0.
1166      *,1150., .63400,1140.,-121., 221.,   0.,   0.
1167      *,1150., .66450,1141.,-120., 110.,   0.,   0.
1168      *,1150., .69500,1141.,-120., 220.,   0.,   0.
1169      *,1150., .72550,1141.,-120., 111.,   0.,   0./
1170       data ((dectab(i,j),i=1,7),j=461,478)/
1171      * 1150., .75600,1141.,-120., 221.,   0.,   0.
1172      *,1150., .78650,1141.,-121., 110.,   0.,   0.
1173      *,1150., .81700,1141.,-121., 220.,   0.,   0.
1174      *,1150., .84750,1141.,-121., 111.,   0.,   0.
1175      *,1150., .87800,1141.,-121., 221.,   0.,   0.
1176      *,1150., .89325,1140.,-130., 230.,   0.,   0.
1177      *,1150., .90850,1140.,-130., 231.,   0.,   0.
1178      *,1150., .92375,1140.,-131., 230.,   0.,   0.
1179      *,1150., .93900,1140.,-131., 231.,   0.,   0.
1180      *,1150., .95425,1141.,-130., 230.,   0.,   0.
1181      *,1150., .96950,1141.,-130., 231.,   0.,   0.
1182      *,1150., .98475,1141.,-131., 230.,   0.,   0.
1183      *,1150.,1.00000,1141.,-131., 231.,   0.,   0.
1184      *,1250., .06000,  12., -11.,1240.,   0.,   0.
1185      *,1250., .12000,  12., -11.,1241.,   0.,   0.
1186      *,1250., .18000,  14., -13.,1240.,   0.,   0.
1187      *,1250., .24000,  14., -13.,1241.,   0.,   0.
1188      *,1250., .25500,  16., -15.,1240.,   0.,   0./
1189       data ((dectab(i,j),i=1,7),j=479,496)/
1190      * 1250., .27000,  16., -15.,1241.,   0.,   0.
1191      *,1250., .28925,1240.,-120.,   0.,   0.,   0.
1192      *,1250., .30850,1240.,-121.,   0.,   0.,   0.
1193      *,1250., .32775,1241.,-120.,   0.,   0.,   0.
1194      *,1250., .34700,1241.,-121.,   0.,   0.,   0.
1195      *,1250., .35775,1240., 340.,   0.,   0.,   0.
1196      *,1250., .36850,1240., 341.,   0.,   0.,   0.
1197      *,1250., .37925,1241., 340.,   0.,   0.,   0.
1198      *,1250., .39000,1241., 341.,   0.,   0.,   0.
1199      *,1250., .42050,1240.,-120., 110.,   0.,   0.
1200      *,1250., .45100,1240.,-120., 220.,   0.,   0.
1201      *,1250., .48150,1240.,-120., 111.,   0.,   0.
1202      *,1250., .51200,1240.,-120., 221.,   0.,   0.
1203      *,1250., .54250,1240.,-121., 110.,   0.,   0.
1204      *,1250., .57300,1240.,-121., 220.,   0.,   0.
1205      *,1250., .60350,1240.,-121., 111.,   0.,   0.
1206      *,1250., .63400,1240.,-121., 221.,   0.,   0.
1207      *,1250., .66450,1241.,-120., 110.,   0.,   0./
1208       data ((dectab(i,j),i=1,7),j=497,514)/
1209      * 1250., .69500,1241.,-120., 220.,   0.,   0.
1210      *,1250., .72550,1241.,-120., 111.,   0.,   0.
1211      *,1250., .75600,1241.,-120., 221.,   0.,   0.
1212      *,1250., .78650,1241.,-121., 110.,   0.,   0.
1213      *,1250., .81700,1241.,-121., 220.,   0.,   0.
1214      *,1250., .84750,1241.,-121., 111.,   0.,   0.
1215      *,1250., .87800,1241.,-121., 221.,   0.,   0.
1216      *,1250., .89325,1240.,-130., 230.,   0.,   0.
1217      *,1250., .90850,1240.,-130., 231.,   0.,   0.
1218      *,1250., .92375,1240.,-131., 230.,   0.,   0.
1219      *,1250., .93900,1240.,-131., 231.,   0.,   0.
1220      *,1250., .95425,1241.,-130., 230.,   0.,   0.
1221      *,1250., .96950,1241.,-130., 231.,   0.,   0.
1222      *,1250., .98475,1241.,-131., 230.,   0.,   0.
1223      *,1250.,1.00000,1241.,-131., 231.,   0.,   0.
1224      *,1350., .06000,  12., -11.,1340.,   0.,   0.
1225      *,1350., .12000,  12., -11.,1341.,   0.,   0.
1226      *,1350., .18000,  14., -13.,1340.,   0.,   0./
1227       data ((dectab(i,j),i=1,7),j=515,532)/
1228      * 1350., .24000,  14., -13.,1341.,   0.,   0.
1229      *,1350., .25500,  16., -15.,1340.,   0.,   0.
1230      *,1350., .27000,  16., -15.,1341.,   0.,   0.
1231      *,1350., .28925,1340.,-120.,   0.,   0.,   0.
1232      *,1350., .30850,1340.,-121.,   0.,   0.,   0.
1233      *,1350., .32775,1341.,-120.,   0.,   0.,   0.
1234      *,1350., .34700,1341.,-121.,   0.,   0.,   0.
1235      *,1350., .35775,1340., 340.,   0.,   0.,   0.
1236      *,1350., .36850,1340., 341.,   0.,   0.,   0.
1237      *,1350., .37925,1341., 340.,   0.,   0.,   0.
1238      *,1350., .39000,1341., 341.,   0.,   0.,   0.
1239      *,1350., .42050,1340.,-120., 110.,   0.,   0.
1240      *,1350., .45100,1340.,-120., 220.,   0.,   0.
1241      *,1350., .48150,1340.,-120., 111.,   0.,   0.
1242      *,1350., .51200,1340.,-120., 221.,   0.,   0.
1243      *,1350., .54250,1340.,-121., 110.,   0.,   0.
1244      *,1350., .57300,1340.,-121., 220.,   0.,   0.
1245      *,1350., .60350,1340.,-121., 111.,   0.,   0./
1246       data ((dectab(i,j),i=1,7),j=533,550)/
1247      * 1350., .63400,1340.,-121., 221.,   0.,   0.
1248      *,1350., .66450,1341.,-120., 110.,   0.,   0.
1249      *,1350., .69500,1341.,-120., 220.,   0.,   0.
1250      *,1350., .72550,1341.,-120., 111.,   0.,   0.
1251      *,1350., .75600,1341.,-120., 221.,   0.,   0.
1252      *,1350., .78650,1341.,-121., 110.,   0.,   0.
1253      *,1350., .81700,1341.,-121., 220.,   0.,   0.
1254      *,1350., .84750,1341.,-121., 111.,   0.,   0.
1255      *,1350., .87800,1341.,-121., 221.,   0.,   0.
1256      *,1350., .89325,1340.,-130., 230.,   0.,   0.
1257      *,1350., .90850,1340.,-130., 231.,   0.,   0.
1258      *,1350., .92375,1340.,-131., 230.,   0.,   0.
1259      *,1350., .93900,1340.,-131., 231.,   0.,   0.
1260      *,1350., .95425,1341.,-130., 230.,   0.,   0.
1261      *,1350., .96950,1341.,-130., 231.,   0.,   0.
1262      *,1350., .98475,1341.,-131., 230.,   0.,   0.
1263      *,1350.,1.00000,1341.,-131., 231.,   0.,   0.
1264      *,2150., .06000,  12., -11.,2140.,   0.,   0./
1265       data ((dectab(i,j),i=1,7),j=551,568)/
1266      * 2150., .12000,  12., -11.,1241.,   0.,   0.
1267      *,2150., .18000,  14., -13.,2140.,   0.,   0.
1268      *,2150., .24000,  14., -13.,1241.,   0.,   0.
1269      *,2150., .25500,  16., -15.,2140.,   0.,   0.
1270      *,2150., .27000,  16., -15.,1241.,   0.,   0.
1271      *,2150., .28925,2140.,-120.,   0.,   0.,   0.
1272      *,2150., .30850,2140.,-121.,   0.,   0.,   0.
1273      *,2150., .32775,1241.,-120.,   0.,   0.,   0.
1274      *,2150., .34700,1241.,-121.,   0.,   0.,   0.
1275      *,2150., .35775,2140., 340.,   0.,   0.,   0.
1276      *,2150., .36850,2140., 341.,   0.,   0.,   0.
1277      *,2150., .37925,1241., 340.,   0.,   0.,   0.
1278      *,2150., .39000,1241., 341.,   0.,   0.,   0.
1279      *,2150., .42050,2140.,-120., 110.,   0.,   0.
1280      *,2150., .45100,2140.,-120., 220.,   0.,   0.
1281      *,2150., .48150,2140.,-120., 111.,   0.,   0.
1282      *,2150., .51200,2140.,-120., 221.,   0.,   0.
1283      *,2150., .54250,2140.,-121., 110.,   0.,   0./
1284       data ((dectab(i,j),i=1,7),j=569,586)/
1285      * 2150., .57300,2140.,-121., 220.,   0.,   0.
1286      *,2150., .60350,2140.,-121., 111.,   0.,   0.
1287      *,2150., .63400,2140.,-121., 221.,   0.,   0.
1288      *,2150., .66450,1241.,-120., 110.,   0.,   0.
1289      *,2150., .69500,1241.,-120., 220.,   0.,   0.
1290      *,2150., .72550,1241.,-120., 111.,   0.,   0.
1291      *,2150., .75600,1241.,-120., 221.,   0.,   0.
1292      *,2150., .78650,1241.,-121., 110.,   0.,   0.
1293      *,2150., .81700,1241.,-121., 220.,   0.,   0.
1294      *,2150., .84750,1241.,-121., 111.,   0.,   0.
1295      *,2150., .87800,1241.,-121., 221.,   0.,   0.
1296      *,2150., .89325,2140.,-130., 230.,   0.,   0.
1297      *,2150., .90850,2140.,-130., 231.,   0.,   0.
1298      *,2150., .92375,2140.,-131., 230.,   0.,   0.
1299      *,2150., .93900,2140.,-131., 231.,   0.,   0.
1300      *,2150., .95425,1241.,-130., 230.,   0.,   0.
1301      *,2150., .96950,1241.,-130., 231.,   0.,   0.
1302      *,2150., .98475,1241.,-131., 230.,   0.,   0./
1303       data ((dectab(i,j),i=1,7),j=587,604)/
1304      * 2150.,1.00000,1241.,-131., 231.,   0.,   0.
1305      *,2250., .06000,  12., -11.,2240.,   0.,   0.
1306      *,2250., .12000,  12., -11.,2241.,   0.,   0.
1307      *,2250., .18000,  14., -13.,2240.,   0.,   0.
1308      *,2250., .24000,  14., -13.,2241.,   0.,   0.
1309      *,2250., .25500,  16., -15.,2240.,   0.,   0.
1310      *,2250., .27000,  16., -15.,2241.,   0.,   0.
1311      *,2250., .28925,2240.,-120.,   0.,   0.,   0.
1312      *,2250., .30850,2240.,-121.,   0.,   0.,   0.
1313      *,2250., .32775,2241.,-120.,   0.,   0.,   0.
1314      *,2250., .34700,2241.,-121.,   0.,   0.,   0.
1315      *,2250., .35775,2240., 340.,   0.,   0.,   0.
1316      *,2250., .36850,2240., 341.,   0.,   0.,   0.
1317      *,2250., .37925,2241., 340.,   0.,   0.,   0.
1318      *,2250., .39000,2241., 341.,   0.,   0.,   0.
1319      *,2250., .42050,2240.,-120., 110.,   0.,   0.
1320      *,2250., .45100,2240.,-120., 220.,   0.,   0.
1321      *,2250., .48150,2240.,-120., 111.,   0.,   0./
1322       data ((dectab(i,j),i=1,7),j=605,622)/
1323      * 2250., .51200,2240.,-120., 221.,   0.,   0.
1324      *,2250., .54250,2240.,-121., 110.,   0.,   0.
1325      *,2250., .57300,2240.,-121., 220.,   0.,   0.
1326      *,2250., .60350,2240.,-121., 111.,   0.,   0.
1327      *,2250., .63400,2240.,-121., 221.,   0.,   0.
1328      *,2250., .66450,2241.,-120., 110.,   0.,   0.
1329      *,2250., .69500,2241.,-120., 220.,   0.,   0.
1330      *,2250., .72550,2241.,-120., 111.,   0.,   0.
1331      *,2250., .75600,2241.,-120., 221.,   0.,   0.
1332      *,2250., .78650,2241.,-121., 110.,   0.,   0.
1333      *,2250., .81700,2241.,-121., 220.,   0.,   0.
1334      *,2250., .84750,2241.,-121., 111.,   0.,   0.
1335      *,2250., .87800,2241.,-121., 221.,   0.,   0.
1336      *,2250., .89325,2240.,-130., 230.,   0.,   0.
1337      *,2250., .90850,2240.,-130., 231.,   0.,   0.
1338      *,2250., .92375,2240.,-131., 230.,   0.,   0.
1339      *,2250., .93900,2240.,-131., 231.,   0.,   0.
1340      *,2250., .95425,2241.,-130., 230.,   0.,   0./
1341       data ((dectab(i,j),i=1,7),j=623,640)/
1342      * 2250., .96950,2241.,-130., 231.,   0.,   0.
1343      *,2250., .98475,2241.,-131., 230.,   0.,   0.
1344      *,2250.,1.00000,2241.,-131., 231.,   0.,   0.
1345      *,2350., .06000,  12., -11.,2340.,   0.,   0.
1346      *,2350., .12000,  12., -11.,2341.,   0.,   0.
1347      *,2350., .18000,  14., -13.,2340.,   0.,   0.
1348      *,2350., .24000,  14., -13.,2341.,   0.,   0.
1349      *,2350., .25500,  16., -15.,2340.,   0.,   0.
1350      *,2350., .27000,  16., -15.,2341.,   0.,   0.
1351      *,2350., .28925,2340.,-120.,   0.,   0.,   0.
1352      *,2350., .30850,2340.,-121.,   0.,   0.,   0.
1353      *,2350., .32775,2341.,-120.,   0.,   0.,   0.
1354      *,2350., .34700,2341.,-121.,   0.,   0.,   0.
1355      *,2350., .35775,2340., 340.,   0.,   0.,   0.
1356      *,2350., .36850,2340., 341.,   0.,   0.,   0.
1357      *,2350., .37925,2341., 340.,   0.,   0.,   0.
1358      *,2350., .39000,2341., 341.,   0.,   0.,   0.
1359      *,2350., .42050,2340.,-120., 110.,   0.,   0./
1360       data ((dectab(i,j),i=1,7),j=641,658)/
1361      * 2350., .45100,2340.,-120., 220.,   0.,   0.
1362      *,2350., .48150,2340.,-120., 111.,   0.,   0.
1363      *,2350., .51200,2340.,-120., 221.,   0.,   0.
1364      *,2350., .54250,2340.,-121., 110.,   0.,   0.
1365      *,2350., .57300,2340.,-121., 220.,   0.,   0.
1366      *,2350., .60350,2340.,-121., 111.,   0.,   0.
1367      *,2350., .63400,2340.,-121., 221.,   0.,   0.
1368      *,2350., .66450,2341.,-120., 110.,   0.,   0.
1369      *,2350., .69500,2341.,-120., 220.,   0.,   0.
1370      *,2350., .72550,2341.,-120., 111.,   0.,   0.
1371      *,2350., .75600,2341.,-120., 221.,   0.,   0.
1372      *,2350., .78650,2341.,-121., 110.,   0.,   0.
1373      *,2350., .81700,2341.,-121., 220.,   0.,   0.
1374      *,2350., .84750,2341.,-121., 111.,   0.,   0.
1375      *,2350., .87800,2341.,-121., 221.,   0.,   0.
1376      *,2350., .89325,2340.,-130., 230.,   0.,   0.
1377      *,2350., .90850,2340.,-130., 231.,   0.,   0.
1378      *,2350., .92375,2340.,-131., 230.,   0.,   0./
1379       data ((dectab(i,j),i=1,7),j=659,720)/
1380      * 434*1./
1381       data ((dectab(i,j),i=1,7),j=721,738)/
1382      * 2350., .93900,2340.,-131., 231.,   0.,   0.
1383      *,2350., .95425,2341.,-130., 230.,   0.,   0.
1384      *,2350., .96950,2341.,-130., 231.,   0.,   0.
1385      *,2350., .98475,2341.,-131., 230.,   0.,   0.
1386      *,2350.,1.00000,2341.,-131., 231.,   0.,   0.
1387      *,3150., .06000,  12., -11.,3140.,   0.,   0.
1388      *,3150., .12000,  12., -11.,1341.,   0.,   0.
1389      *,3150., .18000,  14., -13.,3140.,   0.,   0.
1390      *,3150., .24000,  14., -13.,1341.,   0.,   0.
1391      *,3150., .25500,  16., -15.,3140.,   0.,   0.
1392      *,3150., .27000,  16., -15.,1341.,   0.,   0.
1393      *,3150., .28925,3140.,-120.,   0.,   0.,   0.
1394      *,3150., .30850,3140.,-121.,   0.,   0.,   0.
1395      *,3150., .32775,1341.,-120.,   0.,   0.,   0.
1396      *,3150., .34700,1341.,-121.,   0.,   0.,   0.
1397      *,3150., .35775,3140., 340.,   0.,   0.,   0.
1398      *,3150., .36850,3140., 341.,   0.,   0.,   0.
1399      *,3150., .37925,1341., 340.,   0.,   0.,   0./
1400       data ((dectab(i,j),i=1,7),j=739,756)/
1401      * 3150., .39000,1341., 341.,   0.,   0.,   0.
1402      *,3150., .42050,3140.,-120., 110.,   0.,   0.
1403      *,3150., .45100,3140.,-120., 220.,   0.,   0.
1404      *,3150., .48150,3140.,-120., 111.,   0.,   0.
1405      *,3150., .51200,3140.,-120., 221.,   0.,   0.
1406      *,3150., .54250,3140.,-121., 110.,   0.,   0.
1407      *,3150., .57300,3140.,-121., 220.,   0.,   0.
1408      *,3150., .60350,3140.,-121., 111.,   0.,   0.
1409      *,3150., .63400,3140.,-121., 221.,   0.,   0.
1410      *,3150., .66450,1341.,-120., 110.,   0.,   0.
1411      *,3150., .69500,1341.,-120., 220.,   0.,   0.
1412      *,3150., .72550,1341.,-120., 111.,   0.,   0.
1413      *,3150., .75600,1341.,-120., 221.,   0.,   0.
1414      *,3150., .78650,1341.,-121., 110.,   0.,   0.
1415      *,3150., .81700,1341.,-121., 220.,   0.,   0.
1416      *,3150., .84750,1341.,-121., 111.,   0.,   0.
1417      *,3150., .87800,1341.,-121., 221.,   0.,   0.
1418      *,3150., .89325,3140.,-130., 230.,   0.,   0./
1419       data ((dectab(i,j),i=1,7),j=757,774)/
1420      * 3150., .90850,3140.,-130., 231.,   0.,   0.
1421      *,3150., .92375,3140.,-131., 230.,   0.,   0.
1422      *,3150., .93900,3140.,-131., 231.,   0.,   0.
1423      *,3150., .95425,1341.,-130., 230.,   0.,   0.
1424      *,3150., .96950,1341.,-130., 231.,   0.,   0.
1425      *,3150., .98475,1341.,-131., 230.,   0.,   0.
1426      *,3150.,1.00000,1341.,-131., 231.,   0.,   0.
1427      *,3250., .06000,  12., -11.,3240.,   0.,   0.
1428      *,3250., .12000,  12., -11.,2341.,   0.,   0.
1429      *,3250., .18000,  14., -13.,3240.,   0.,   0.
1430      *,3250., .24000,  14., -13.,2341.,   0.,   0.
1431      *,3250., .25500,  16., -15.,3240.,   0.,   0.
1432      *,3250., .27000,  16., -15.,2341.,   0.,   0.
1433      *,3250., .28925,3240.,-120.,   0.,   0.,   0.
1434      *,3250., .30850,3240.,-121.,   0.,   0.,   0.
1435      *,3250., .32775,2341.,-120.,   0.,   0.,   0.
1436      *,3250., .34700,2341.,-121.,   0.,   0.,   0.
1437      *,3250., .35775,3240., 340.,   0.,   0.,   0./
1438       data ((dectab(i,j),i=1,7),j=775,792)/
1439      * 3250., .36850,3240., 341.,   0.,   0.,   0.
1440      *,3250., .37925,2341., 340.,   0.,   0.,   0.
1441      *,3250., .39000,2341., 341.,   0.,   0.,   0.
1442      *,3250., .42050,3240.,-120., 110.,   0.,   0.
1443      *,3250., .45100,3240.,-120., 220.,   0.,   0.
1444      *,3250., .48150,3240.,-120., 111.,   0.,   0.
1445      *,3250., .51200,3240.,-120., 221.,   0.,   0.
1446      *,3250., .54250,3240.,-121., 110.,   0.,   0.
1447      *,3250., .57300,3240.,-121., 220.,   0.,   0.
1448      *,3250., .60350,3240.,-121., 111.,   0.,   0.
1449      *,3250., .63400,3240.,-121., 221.,   0.,   0.
1450      *,3250., .66450,2341.,-120., 110.,   0.,   0.
1451      *,3250., .69500,2341.,-120., 220.,   0.,   0.
1452      *,3250., .72550,2341.,-120., 111.,   0.,   0.
1453      *,3250., .75600,2341.,-120., 221.,   0.,   0.
1454      *,3250., .78650,2341.,-121., 110.,   0.,   0.
1455      *,3250., .81700,2341.,-121., 220.,   0.,   0.
1456      *,3250., .84750,2341.,-121., 111.,   0.,   0./
1457       data ((dectab(i,j),i=1,7),j=793,810)/
1458      * 3250., .87800,2341.,-121., 221.,   0.,   0.
1459      *,3250., .89325,3240.,-130., 230.,   0.,   0.
1460      *,3250., .90850,3240.,-130., 231.,   0.,   0.
1461      *,3250., .92375,3240.,-131., 230.,   0.,   0.
1462      *,3250., .93900,3240.,-131., 231.,   0.,   0.
1463      *,3250., .95425,2341.,-130., 230.,   0.,   0.
1464      *,3250., .96950,2341.,-130., 231.,   0.,   0.
1465      *,3250., .98475,2341.,-131., 230.,   0.,   0.
1466      *,3250.,1.00000,2341.,-131., 231.,   0.,   0.
1467      *,3350., .06000,  12., -11.,3340.,   0.,   0.
1468      *,3350., .12000,  12., -11.,3341.,   0.,   0.
1469      *,3350., .18000,  14., -13.,3340.,   0.,   0.
1470      *,3350., .24000,  14., -13.,3341.,   0.,   0.
1471      *,3350., .25500,  16., -15.,3340.,   0.,   0.
1472      *,3350., .27000,  16., -15.,3341.,   0.,   0.
1473      *,3350., .28925,3340.,-120.,   0.,   0.,   0.
1474      *,3350., .30850,3340.,-121.,   0.,   0.,   0.
1475      *,3350., .32775,3341.,-120.,   0.,   0.,   0./
1476       data ((dectab(i,j),i=1,7),j=811,828)/
1477      * 3350., .34700,3341.,-121.,   0.,   0.,   0.
1478      *,3350., .35775,3340., 340.,   0.,   0.,   0.
1479      *,3350., .36850,3340., 341.,   0.,   0.,   0.
1480      *,3350., .37925,3341., 340.,   0.,   0.,   0.
1481      *,3350., .39000,3341., 341.,   0.,   0.,   0.
1482      *,3350., .42050,3340.,-120., 110.,   0.,   0.
1483      *,3350., .45100,3340.,-120., 220.,   0.,   0.
1484      *,3350., .48150,3340.,-120., 111.,   0.,   0.
1485      *,3350., .51200,3340.,-120., 221.,   0.,   0.
1486      *,3350., .54250,3340.,-121., 110.,   0.,   0.
1487      *,3350., .57300,3340.,-121., 220.,   0.,   0.
1488      *,3350., .60350,3340.,-121., 111.,   0.,   0.
1489      *,3350., .63400,3340.,-121., 221.,   0.,   0.
1490      *,3350., .66450,3341.,-120., 110.,   0.,   0.
1491      *,3350., .69500,3341.,-120., 220.,   0.,   0.
1492      *,3350., .72550,3341.,-120., 111.,   0.,   0.
1493      *,3350., .75600,3341.,-120., 221.,   0.,   0.
1494      *,3350., .78650,3341.,-121., 110.,   0.,   0./
1495       data ((dectab(i,j),i=1,7),j=829,846)/
1496      * 3350., .81700,3341.,-121., 220.,   0.,   0.
1497      *,3350., .84750,3341.,-121., 111.,   0.,   0.
1498      *,3350., .87800,3341.,-121., 221.,   0.,   0.
1499      *,3350., .89325,3340.,-130., 230.,   0.,   0.
1500      *,3350., .90850,3340.,-130., 231.,   0.,   0.
1501      *,3350., .92375,3340.,-131., 230.,   0.,   0.
1502      *,3350., .93900,3340.,-131., 231.,   0.,   0.
1503      *,3350., .95425,3341.,-130., 230.,   0.,   0.
1504      *,3350., .96950,3341.,-130., 231.,   0.,   0.
1505      *,3350., .98475,3341.,-131., 230.,   0.,   0.
1506      *,3350.,1.00000,3341.,-131., 231.,   0.,   0.
1507      *,1160., .33300,   1.,  -2.,1500.,   0.,   0.
1508      *,1160., .66700,   4.,  -3.,1500.,   0.,   0.
1509      *,1160., .77800, -12.,  11.,1500.,   0.,   0.
1510      *,1160., .88900, -14.,  13.,1500.,   0.,   0.
1511      *,1160.,1.00000, -16.,  15.,1500.,   0.,   0.
1512      *,1260., .33300,   1.,  -2.,2500.,   0.,   0.
1513      *,1260., .66700,   4.,  -3.,2500.,   0.,   0./
1514       data ((dectab(i,j),i=1,7),j=847,864)/
1515      * 1260., .77800, -12.,  11.,2500.,   0.,   0.
1516      *,1260., .88900, -14.,  13.,2500.,   0.,   0.
1517      *,1260.,1.00000, -16.,  15.,2500.,   0.,   0.
1518      *,2260., .33300,   1.,  -2.,2500.,   0.,   0.
1519      *,2260., .66700,   4.,  -3.,2500.,   0.,   0.
1520      *,2260., .77800, -12.,  11.,2500.,   0.,   0.
1521      *,2260., .88900, -14.,  13.,2500.,   0.,   0.
1522      *,2260.,1.00000, -16.,  15.,2500.,   0.,   0.
1523      *,2160., .33300,   1.,  -2.,1500.,   0.,   0.
1524      *,2160., .66700,   4.,  -3.,1500.,   0.,   0.
1525      *,2160., .77800, -12.,  11.,1500.,   0.,   0.
1526      *,2160., .88900, -14.,  13.,1500.,   0.,   0.
1527      *,2160.,1.00000, -16.,  15.,1500.,   0.,   0.
1528      *,1360., .33300,   1.,  -2.,3500.,   0.,   0.
1529      *,1360., .66700,   4.,  -3.,3500.,   0.,   0.
1530      *,1360., .77800, -12.,  11.,3500.,   0.,   0.
1531      *,1360., .88900, -14.,  13.,3500.,   0.,   0.
1532      *,1360.,1.00000, -16.,  15.,3500.,   0.,   0./
1533       data ((dectab(i,j),i=1,7),j=865,882)/
1534      * 2360., .33300,   1.,  -2.,3500.,   0.,   0.
1535      *,2360., .66700,   4.,  -3.,3500.,   0.,   0.
1536      *,2360., .77800, -12.,  11.,3500.,   0.,   0.
1537      *,2360., .88900, -14.,  13.,3500.,   0.,   0.
1538      *,2360.,1.00000, -16.,  15.,3500.,   0.,   0.
1539      *,3360., .33300,   1.,  -2.,3500.,   0.,   0.
1540      *,3360., .66700,   4.,  -3.,3500.,   0.,   0.
1541      *,3360., .77800, -12.,  11.,3500.,   0.,   0.
1542      *,3360., .88900, -14.,  13.,3500.,   0.,   0.
1543      *,3360.,1.00000, -16.,  15.,3500.,   0.,   0.
1544      *,1151.,1.00000,1150.,  10.,   0.,   0.,   0.
1545      *,1251.,1.00000,1250.,  10.,   0.,   0.,   0.
1546      *,2251.,1.00000,2250.,  10.,   0.,   0.,   0.
1547      *,1351.,1.00000,1350.,  10.,   0.,   0.,   0.
1548      *,2351.,1.00000,2350.,  10.,   0.,   0.,   0.
1549      *,3351.,1.00000,3350.,  10.,   0.,   0.,   0.
1550      *,1161.,1.00000,1160.,  10.,   0.,   0.,   0.
1551      *,1261.,1.00000,1260.,  10.,   0.,   0.,   0./
1552       data ((dectab(i,j),i=1,7),j=883,886)/
1553      * 2261.,1.00000,2260.,  10.,   0.,   0.,   0.
1554      *,1361.,1.00000,1360.,  10.,   0.,   0.,   0.
1555      *,2361.,1.00000,2360.,  10.,   0.,   0.,   0.
1556      *,3361.,1.00000,3360.,  10.,   0.,   0.,   0./
1557 c    *---------------------------------------------
1558 c    *    delta++ resonances
1559 c    *---------------------------------------------
1560       data ((dectab(i,j),i=1,7),j=887,900)/
1561 c    *--dl++(1620)---------------------------------
1562      * 1112., .30000,1120., 120.,   0.,   0.,   0.
1563      *,1112., .66000,1111., 110.,   0.,   0.,   0.
1564      *,1112., .90000,1121., 120.,   0.,   0.,   0.
1565      *,1112.,1.00000,1120., 120., 110.,   0.,   0.
1566 c    *--dl++(1700)---------------------------------
1567      *,1113., .15000,1120., 120.,   0.,   0.,   0.
1568      *,1113., .51000,1111., 110.,   0.,   0.,   0.
1569      *,1113., .75000,1121., 120.,   0.,   0.,   0.
1570      *,1113.,1.00000,1120., 120., 110.,   0.,   0.
1571 c    *--dl++(1925)---------------------------------
1572      *,1114., .28000,1120., 120.,   0.,   0.,   0.
1573      *,1114., .40600,1111., 110.,   0.,   0.,   0.
1574      *,1114., .49000,1121., 120.,   0.,   0.,   0.
1575      *,1114., .69000,1120., 121.,   0.,   0.,   0.
1576      *,1114., .70000,1130., 130.,   0.,   0.,   0.
1577      *,1114.,1.00000,1122., 120.,   0.,   0.,   0./
1578 c    *---------------------------------------------
1579 c    *    delta- resonances
1580 c    *---------------------------------------------
1581       data ((dectab(i,j),i=1,7),j=901,914)/
1582 c    *--dl-(1620)----------------------------------
1583      * 2222., .30000,1220.,-120.,   0.,   0.,   0.
1584      *,2222., .66000,2221., 110.,   0.,   0.,   0.
1585      *,2222., .90000,1221.,-120.,   0.,   0.,   0.
1586      *,2222.,1.00000,1220., 110.,-120.,   0.,   0.
1587 c    *--dl-(1700)----------------------------------
1588      *,2223., .15000,1220.,-120.,   0.,   0.,   0.
1589      *,2223., .51000,2221., 110.,   0.,   0.,   0.
1590      *,2223., .75000,1221.,-120.,   0.,   0.,   0.
1591      *,2223.,1.00000,1220., 110.,-120.,   0.,   0.
1592 c    *--dl-(1925)----------------------------------
1593      *,2224., .28000,1220.,-120.,   0.,   0.,   0.
1594      *,2224., .40600,2221., 110.,   0.,   0.,   0.
1595      *,2224., .49000,1221.,-120.,   0.,   0.,   0.
1596      *,2224., .69000,1220.,-121.,   0.,   0.,   0.
1597      *,2224., .70000,2230., 230.,   0.,   0.,   0.
1598      *,2224.,1.00000,1222.,-120.,   0.,   0.,   0./
1599 c    *---------------------------------------------
1600 c    *    n*+ resonances + delta+ resonances
1601 c    *---------------------------------------------
1602       data ((dectab(i,j),i=1,7),j=915,931)/
1603 c    *--n*+(1440)----------------------------------
1604      * 1122., .20000,1120., 110.,   0.,   0.,   0.
1605      *,1122., .60000,1220., 120.,   0.,   0.,   0.
1606      *,1122., .68000,1111.,-120.,   0.,   0.,   0.
1607      *,1122., .73000,1121., 110.,   0.,   0.,   0.
1608      *,1122., .76000,1221., 120.,   0.,   0.,   0.
1609      *,1122., .84000,1120., 120.,-120.,   0.,   0.
1610      *,1122., .87000,1120., 110., 110.,   0.,   0.
1611      *,1122.,1.00000,1220., 120., 110.,   0.,   0.
1612 c    *--n*+(1530)----------------------------------
1613      *,1123., .17000,1120., 110.,   0.,   0.,   0.
1614      *,1123., .51000,1220., 120.,   0.,   0.,   0.
1615      *,1123., .57000,1111.,-120.,   0.,   0.,   0.
1616      *,1123., .61000,1121., 110.,   0.,   0.,   0.
1617      *,1123., .63000,1221., 120.,   0.,   0.,   0.
1618      *,1123., .67000,1120., 120.,-120.,   0.,   0.
1619      *,1123., .68000,1120., 110., 110.,   0.,   0.
1620      *,1123., .75000,1220., 120., 110.,   0.,   0.
1621      *,1123.,1.00000,1120., 220.,   0.,   0.,   0./
1622       data ((dectab(i,j),i=1,7),j=932,948)/
1623 c    *--dl+(1620)----------------------------------
1624      * 1124., .20000,1120., 110.,   0.,   0.,   0.
1625      *,1124., .30000,1220., 120.,   0.,   0.,   0.
1626      *,1124., .54000,1111.,-120.,   0.,   0.,   0.
1627      *,1124., .58000,1121., 110.,   0.,   0.,   0.
1628      *,1124., .90000,1221., 120.,   0.,   0.,   0.
1629      *,1124., .96000,1120., 120.,-120.,   0.,   0.
1630      *,1124.,1.00000,1220., 120., 110.,   0.,   0.
1631 c    *--n*+(1665)----------------------------------
1632      *,1125., .16700,1120., 110.,   0.,   0.,   0.
1633      *,1125., .49970,1220., 120.,   0.,   0.,   0.
1634      *,1125., .62470,1111.,-120.,   0.,   0.,   0.
1635      *,1125., .70800,1121., 110.,   0.,   0.,   0.
1636      *,1125., .74970,1221., 120.,   0.,   0.,   0.
1637      *,1125., .82080,1120., 120.,-120.,   0.,   0.
1638      *,1125., .85190,1120., 110., 110.,   0.,   0.
1639      *,1125., .96300,1220., 120., 110.,   0.,   0.
1640      *,1125., .97300,1120., 220.,   0.,   0.,   0.
1641      *,1125.,1.00000,2130., 130.,   0.,   0.,   0./
1642       data ((dectab(i,j),i=1,7),j=949,955)/
1643 c    *--dl+(1700)----------------------------------
1644      * 1126., .10000,1120., 110.,   0.,   0.,   0.
1645      *,1126., .15000,1220., 120.,   0.,   0.,   0.
1646      *,1126., .39000,1111.,-120.,   0.,   0.,   0.
1647      *,1126., .43000,1121., 110.,   0.,   0.,   0.
1648      *,1126., .75000,1221., 120.,   0.,   0.,   0.
1649      *,1126., .91500,1120., 120.,-120.,   0.,   0.
1650      *,1126.,1.00000,1220., 120., 110.,   0.,   0./
1651       data ((dectab(i,j),i=1,7),j=956,969)/
1652 c    *--n*+(1710)----------------------------------
1653      * 1127., .04430,1120., 110.,   0.,   0.,   0.
1654      *,1127., .13290,1220., 120.,   0.,   0.,   0.
1655      *,1127., .23790,1111.,-120.,   0.,   0.,   0.
1656      *,1127., .30790,1121., 110.,   0.,   0.,   0.
1657      *,1127., .34290,1221., 120.,   0.,   0.,   0.
1658      *,1127., .41190,1120., 120.,-120.,   0.,   0.
1659      *,1127., .48090,1120., 110., 110.,   0.,   0.
1660      *,1127., .54990,1220., 120., 110.,   0.,   0.
1661      *,1127., .66070,1120., 220.,   0.,   0.,   0.
1662      *,1127., .72800,2130., 130.,   0.,   0.,   0.
1663      *,1127., .74930,1230., 130.,   0.,   0.,   0.
1664      *,1127., .76000,1130., 230.,   0.,   0.,   0.
1665      *,1127., .84000,1120., 111.,   0.,   0.,   0.
1666      *,1127.,1.00000,1220., 121.,   0.,   0.,   0./
1667       data ((dectab(i,j),i=1,7),j=970,980)/
1668 c    *--dl+(1925)----------------------------------
1669      * 1128., .18700,1120., 110.,   0.,   0.,   0.
1670      *,1128., .28000,1220., 120.,   0.,   0.,   0.
1671      *,1128., .36400,1111.,-120.,   0.,   0.,   0.
1672      *,1128., .37800,1121., 110.,   0.,   0.,   0.
1673      *,1128., .49000,1221., 120.,   0.,   0.,   0.
1674      *,1128., .62300,1120., 111.,   0.,   0.,   0.
1675      *,1128., .69000,1220., 121.,   0.,   0.,   0.
1676      *,1128., .69350,1130., 230.,   0.,   0.,   0.
1677      *,1128., .69900,1230., 130.,   0.,   0.,   0.
1678      *,1128., .89900,1122., 110.,   0.,   0.,   0.
1679      *,1128.,1.00000,1222., 120.,   0.,   0.,   0./
1680 c    *---------------------------------------------
1681 c    *    n*0  resonances + delta0 resonances
1682 c    *---------------------------------------------
1683       data ((dectab(i,j),i=1,7),j=981,997)/
1684 c    *----------n*0(1440)--------------------------
1685      * 1222., .20000,1220., 110.,   0.,   0.,   0.
1686      *,1222., .60000,1120.,-120.,   0.,   0.,   0.
1687      *,1222., .68000,2221., 120.,   0.,   0.,   0.
1688      *,1222., .73000,1221., 110.,   0.,   0.,   0.
1689      *,1222., .76000,1121.,-120.,   0.,   0.,   0.
1690      *,1222., .84000,1220., 120.,-120.,   0.,   0.
1691      *,1222., .87000,1220., 110., 110.,   0.,   0.
1692      *,1222.,1.00000,1120.,-120., 110.,   0.,   0.
1693 c    *----------n*0(1530)--------------------------
1694      *,1223., .17000,1220., 110.,   0.,   0.,   0.
1695      *,1223., .51000,1120.,-120.,   0.,   0.,   0.
1696      *,1223., .57000,2221., 120.,   0.,   0.,   0.
1697      *,1223., .61000,1221., 110.,   0.,   0.,   0.
1698      *,1223., .63000,1121.,-120.,   0.,   0.,   0.
1699      *,1223., .67000,1220., 120.,-120.,   0.,   0.
1700      *,1223., .68000,1220., 110., 110.,   0.,   0.
1701      *,1223., .75000,1120.,-120., 110.,   0.,   0.
1702      *,1223.,1.00000,1220., 220.,   0.,   0.,   0./
1703       data ((dectab(i,j),i=1,7),j=998,1014)/
1704 c    *----------dl0(1620)--------------------------
1705      * 1224., .20000,1220., 110.,   0.,   0.,   0.
1706      *,1224., .30000,1120.,-120.,   0.,   0.,   0.
1707      *,1224., .54000,2221., 120.,   0.,   0.,   0.
1708      *,1224., .58000,1221., 110.,   0.,   0.,   0.
1709      *,1224., .90000,1121.,-120.,   0.,   0.,   0.
1710      *,1224., .96500,1220., 120.,-120.,   0.,   0.
1711      *,1224.,1.00000,1120.,-120., 110.,   0.,   0.
1712 c    *----------n*0(1665)--------------------------
1713      *,1225., .16700,1220., 110.,   0.,   0.,   0.
1714      *,1225., .49970,1120.,-120.,   0.,   0.,   0.
1715      *,1225., .62470,2221., 120.,   0.,   0.,   0.
1716      *,1225., .70800,1221., 110.,   0.,   0.,   0.
1717      *,1225., .74970,1121.,-120.,   0.,   0.,   0.
1718      *,1225., .82080,1220., 120.,-120.,   0.,   0.
1719      *,1225., .85190,1220., 110., 110.,   0.,   0.
1720      *,1225., .96300,1120.,-120., 110.,   0.,   0.
1721      *,1225., .97300,1220., 220.,   0.,   0.,   0.
1722      *,1225.,1.00000,2130., 230.,   0.,   0.,   0./
1723       data ((dectab(i,j),i=1,7),j=1015,1021)/
1724 c    *----------dl0(1700)--------------------------
1725      * 1226., .10000,1220., 110.,   0.,   0.,   0.
1726      *,1226., .15000,1120.,-120.,   0.,   0.,   0.
1727      *,1226., .39000,2221., 120.,   0.,   0.,   0.
1728      *,1226., .43000,1221., 110.,   0.,   0.,   0.
1729      *,1226., .75000,1121.,-120.,   0.,   0.,   0.
1730      *,1226., .91500,1220., 120.,-120.,   0.,   0.
1731      *,1226.,1.00000,1120.,-120., 110.,   0.,   0./
1732       data ((dectab(i,j),i=1,7),j=1022,1035)/
1733 c    *----------n*0(1710)--------------------------
1734      * 1227., .04430,1220., 110.,   0.,   0.,   0.
1735      *,1227., .13290,1120.,-120.,   0.,   0.,   0.
1736      *,1227., .23790,2221., 120.,   0.,   0.,   0.
1737      *,1227., .30790,1221., 110.,   0.,   0.,   0.
1738      *,1227., .34290,1121.,-120.,   0.,   0.,   0.
1739      *,1227., .41190,1220., 120.,-120.,   0.,   0.
1740      *,1227., .48090,1220., 110., 110.,   0.,   0.
1741      *,1227., .54990,1120.,-120., 110.,   0.,   0.
1742      *,1227., .66070,1220., 220.,   0.,   0.,   0.
1743      *,1227., .72800,2130., 230.,   0.,   0.,   0.
1744      *,1227., .73870,1230., 230.,   0.,   0.,   0.
1745      *,1227., .76000,2230., 130.,   0.,   0.,   0.
1746      *,1227., .92000,1120.,-121.,   0.,   0.,   0.
1747      *,1227.,1.00000,1220., 111.,   0.,   0.,   0./
1748       data ((dectab(i,j),i=1,7),j=1036,1046)/
1749 c    *----------dl0(1925)--------------------------
1750      * 1228., .18700,1220., 110.,   0.,   0.,   0.
1751      *,1228., .28000,1120.,-120.,   0.,   0.,   0.
1752      *,1228., .36400,2221., 120.,   0.,   0.,   0.
1753      *,1228., .37800,1221., 110.,   0.,   0.,   0.
1754      *,1228., .49000,1121.,-120.,   0.,   0.,   0.
1755      *,1228., .55700,1220., 111.,   0.,   0.,   0.
1756      *,1228., .69000,1120.,-121.,   0.,   0.,   0.
1757      *,1228., .69350,2230., 130.,   0.,   0.,   0.
1758      *,1228., .70000,1230., 230.,   0.,   0.,   0.
1759      *,1228., .80000,1122.,-120.,   0.,   0.,   0.
1760      *,1228.,1.00000,1222., 110.,   0.,   0.,   0./
1761 c    *---------------------------------------------
1762 c    *   lambda resonances + sigma0 resonances
1763 c    *---------------------------------------------
1764       data ((dectab(i,j),i=1,7),j=1047,1059)/
1765 c    *----------lambda(1405)-----------------------
1766      * 1233., .33000,1230., 110.,   0.,   0.,   0.
1767      *,1233., .66000,2230., 120.,   0.,   0.,   0.
1768      *,1233.,1.00000,1130.,-120.,   0.,   0.,   0.
1769 c    *----------lambda(1520)-----------------------
1770      *,1234., .22500,1120.,-130.,   0.,   0.,   0.
1771      *,1234., .48000,1220.,-230.,   0.,   0.,   0.
1772      *,1234., .62000,1230., 110.,   0.,   0.,   0.
1773      *,1234., .76000,2230., 120.,   0.,   0.,   0.
1774      *,1234., .90000,1130.,-120.,   0.,   0.,   0.
1775      *,1234., .96000,2130., 120.,-120.,   0.,   0.
1776      *,1234., .99000,2130., 110., 110.,   0.,   0.
1777      *,1234., .99330,1130.,-120., 110.,   0.,   0.
1778      *,1234., .99660,2230., 120., 110.,   0.,   0.
1779      *,1234.,1.00000,1230., 120.,-120.,   0.,   0./
1780       data ((dectab(i,j),i=1,7),j=1060,1075)/
1781 c    *----------lambda(1645)-----------------------
1782      * 1235., .10000,1120.,-130.,   0.,   0.,   0.
1783      *,1235., .20000,1220.,-230.,   0.,   0.,   0.
1784      *,1235., .35000,1230., 110.,   0.,   0.,   0.
1785      *,1235., .50000,2230., 120.,   0.,   0.,   0.
1786      *,1235., .65000,1130.,-120.,   0.,   0.,   0.
1787      *,1235., .75000,2130., 120.,-120.,   0.,   0.
1788      *,1235., .80000,2130., 110., 110.,   0.,   0.
1789      *,1235., .84500,1130.,-120., 110.,   0.,   0.
1790      *,1235., .89000,2230., 120., 110.,   0.,   0.
1791      *,1235., .93500,1230., 120.,-120.,   0.,   0.
1792      *,1235.,1.00000,2130., 220.,   0.,   0.,   0.
1793 c    *----------sigma0(1665)-----------------------
1794      *,1236., .10000,1120.,-130.,   0.,   0.,   0.
1795      *,1236., .20000,1220.,-230.,   0.,   0.,   0.
1796      *,1236., .40000,2230., 120.,   0.,   0.,   0.
1797      *,1236., .60000,1130.,-120.,   0.,   0.,   0.
1798      *,1236.,1.00000,2130., 110.,   0.,   0.,   0./
1799       data ((dectab(i,j),i=1,7),j=1076,1084)/
1800 c    *----------sigma0(1776)-----------------------
1801      * 1237., .17500,1120.,-130.,   0.,   0.,   0.
1802      *,1237., .35000,1220.,-230.,   0.,   0.,   0.
1803      *,1237., .38750,2230., 120.,   0.,   0.,   0.
1804      *,1237., .42500,1130.,-120.,   0.,   0.,   0.
1805      *,1237., .57500,2130., 110.,   0.,   0.,   0.
1806      *,1237., .60000,2231., 120.,   0.,   0.,   0.
1807      *,1237., .62500,1131.,-120.,   0.,   0.,   0.
1808      *,1237., .75000,1234., 110.,   0.,   0.,   0.
1809      *,1237.,1.00000,1230., 220.,   0.,   0.,   0./
1810       data ((dectab(i,j),i=1,7),j=1085,1094)/
1811 c    *----------lambda(1845)-----------------------
1812      * 1238., .17000,1120.,-130.,   0.,   0.,   0.
1813      *,1238., .34000,1220.,-230.,   0.,   0.,   0.
1814      *,1238., .44000,1230., 110.,   0.,   0.,   0.
1815      *,1238., .54000,2230., 120.,   0.,   0.,   0.
1816      *,1238., .64000,1130.,-120.,   0.,   0.,   0.
1817      *,1238., .70000,1231., 110.,   0.,   0.,   0.
1818      *,1238., .76000,2231., 120.,   0.,   0.,   0.
1819      *,1238., .82000,1131.,-120.,   0.,   0.,   0.
1820      *,1238., .91000,1120.,-131.,   0.,   0.,   0.
1821      *,1238.,1.00000,1220.,-231.,   0.,   0.,   0./
1822       data ((dectab(i,j),i=1,7),j=1095,1106)/
1823 c    *----------sigma0(1930)-----------------------
1824      * 1239., .07500,1120.,-130.,   0.,   0.,   0.
1825      *,1239., .15000,1220.,-230.,   0.,   0.,   0.
1826      *,1239., .20000,1121.,-130.,   0.,   0.,   0.
1827      *,1239., .25000,1221.,-230.,   0.,   0.,   0.
1828      *,1239., .32500,1120.,-131.,   0.,   0.,   0.
1829      *,1239., .40000,1220.,-231.,   0.,   0.,   0.
1830      *,1239., .47500,2230., 120.,   0.,   0.,   0.
1831      *,1239., .55000,1130.,-120.,   0.,   0.,   0.
1832      *,1239., .70000,2130., 110.,   0.,   0.,   0.
1833      *,1239., .77500,2231., 120.,   0.,   0.,   0.
1834      *,1239., .85000,1131.,-120.,   0.,   0.,   0.
1835      *,1239.,1.00000,1234., 110.,   0.,   0.,   0./
1836 c    *---------------------------------------------
1837 c    *            sigma+ resonances
1838 c    *---------------------------------------------
1839       data ((dectab(i,j),i=1,7),j=1107,1118)/
1840 c    *----------sigma+(1665)-----------------------
1841      * 1132., .20000,1120.,-230.,   0.,   0.,   0.
1842      *,1132., .40000,1130., 110.,   0.,   0.,   0.
1843      *,1132., .60000,1230., 120.,   0.,   0.,   0.
1844      *,1132.,1.00000,2130., 120.,   0.,   0.,   0.
1845 c    *----------sigma+(1776)-----------------------
1846      *,1133., .35000,1120.,-230.,   0.,   0.,   0.
1847      *,1133., .38750,1130., 110.,   0.,   0.,   0.
1848      *,1133., .42500,1230., 120.,   0.,   0.,   0.
1849      *,1133., .57500,2130., 120.,   0.,   0.,   0.
1850      *,1133., .60000,1131., 110.,   0.,   0.,   0.
1851      *,1133., .62500,1231., 120.,   0.,   0.,   0.
1852      *,1133., .75000,1234., 120.,   0.,   0.,   0.
1853      *,1133.,1.00000,1130., 220.,   0.,   0.,   0./
1854       data ((dectab(i,j),i=1,7),j=1119,1128)/
1855 c    *----------sigma+(1930)-----------------------
1856      * 1134., .15000,1120.,-230.,   0.,   0.,   0.
1857      *,1134., .22500,1111.,-130.,   0.,   0.,   0.
1858      *,1134., .25000,1121.,-230.,   0.,   0.,   0.
1859      *,1134., .40000,1120.,-231.,   0.,   0.,   0.
1860      *,1134., .47500,1130., 110.,   0.,   0.,   0.
1861      *,1134., .55000,1230., 120.,   0.,   0.,   0.
1862      *,1134., .70000,2130., 120.,   0.,   0.,   0.
1863      *,1134., .77500,1131., 110.,   0.,   0.,   0.
1864      *,1134., .85000,1231., 120.,   0.,   0.,   0.
1865      *,1134.,1.00000,1234., 120.,   0.,   0.,   0./
1866 c    *---------------------------------------------
1867 c    *            sigma- resonances
1868 c    *---------------------------------------------
1869       data ((dectab(i,j),i=1,7),j=1129,1140)/
1870 c    *----------sigma-(1665)-----------------------
1871      * 2232., .20000,1220.,-130.,   0.,   0.,   0.
1872      *,2232., .40000,2230., 110.,   0.,   0.,   0.
1873      *,2232., .60000,1230.,-120.,   0.,   0.,   0.
1874      *,2232.,1.00000,2130.,-120.,   0.,   0.,   0.
1875 c    *----------sigma-(1776)-----------------------
1876      *,2233., .35000,1220.,-130.,   0.,   0.,   0.
1877      *,2233., .38750,2230., 110.,   0.,   0.,   0.
1878      *,2233., .42500,1230.,-120.,   0.,   0.,   0.
1879      *,2233., .57500,2130.,-120.,   0.,   0.,   0.
1880      *,2233., .60000,2231., 110.,   0.,   0.,   0.
1881      *,2233., .62500,1231.,-120.,   0.,   0.,   0.
1882      *,2233., .75000,1234.,-120.,   0.,   0.,   0.
1883      *,2233.,1.00000,2230., 220.,   0.,   0.,   0./
1884       data ((dectab(i,j),i=1,7),j=1141,1150)/
1885 c    *----------sigma-(1930)-----------------------
1886      * 2234., .15000,1220.,-130.,   0.,   0.,   0.
1887      *,2234., .17500,1221.,-130.,   0.,   0.,   0.
1888      *,2234., .25000,2221.,-230.,   0.,   0.,   0.
1889      *,2234., .40000,1220.,-131.,   0.,   0.,   0.
1890      *,2234., .47500,2230., 110.,   0.,   0.,   0.
1891      *,2234., .55000,1230.,-120.,   0.,   0.,   0.
1892      *,2234., .70000,2130.,-120.,   0.,   0.,   0.
1893      *,2234., .77500,2231., 110.,   0.,   0.,   0.
1894      *,2234., .85000,1231.,-120.,   0.,   0.,   0.
1895      *,2234.,1.00000,1234.,-120.,   0.,   0.,   0./
1896 c    *---------------------------------------------
1897 c    *      additional mesonresonances
1898 c    *---------------------------------------------
1899       data ((dectab(i,j),i=1,7),j=1151,1159)/
1900 c    *-----------f0(975)---------------------------
1901      *  332., .50000, 120.,-120.,   0.,   0.,   0.
1902      *, 332., .75000, 110., 110.,   0.,   0.,   0.
1903      *, 332., .87500, 130.,-130.,   0.,   0.,   0.
1904      *, 332.,1.00000, 230.,-230.,   0.,   0.,   0.
1905 c    *-----------a0(980)---------------------------
1906      *, 112., .56000, 110., 220.,   0.,   0.,   0.
1907      *, 112., .78000, 130.,-130.,   0.,   0.,   0.
1908      *, 112.,1.00000, 230.,-230.,   0.,   0.,   0.
1909 c    *-----------a+(980)---------------------------
1910      *, 122., .60000, 120., 220.,   0.,   0.,   0.
1911      *, 122.,1.00000, 130.,-230.,   0.,   0.,   0./
1912 c    *---------------------------------------------
1913 c    *      weak baryon decays
1914 c    *---------------------------------------------
1915       data ((dectab(i,j),i=1,7),j=1160,1169)/
1916 c    *-----------lambda(1116)----------------------
1917      * 2130.,0.64200,1120.,-120.,   0.,   0.,   0.
1918      *,2130.,1.00000,1220., 110.,   0.,   0.,   0.
1919 c    *-----------sigma+(1180)----------------------
1920      *,1130.,0.51580,1120., 110.,   0.,   0.,   0.
1921      *,1130.,1.00000,1220., 120.,   0.,   0.,   0.
1922 c    *-----------sigma-(1180)----------------------
1923      *,2230.,1.00000,1220.,-120.,   0.,   0.,   0.
1924 c    *---------kaskade-(1360)----------------------
1925      *,2330.,1.00000,2130.,-120.,   0.,   0.,   0.
1926 c    *---------kaskade0(1360)----------------------
1927      *,1330.,1.00000,2130., 110.,   0.,   0.,   0.
1928 c    *---------omega-(1680)------------------------
1929      *,3331.,0.68000,2130.,-130.,   0.,   0.,   0.
1930      *,3331.,0.82000,1330.,-120.,   0.,   0.,   0.
1931      *,3331.,1.00000,2330., 110.,   0.,   0.,   0./
1932 c    *---------------------------------------------
1933 c    *      weak meson decays
1934 c    *---------------------------------------------
1935       data ((dectab(i,j),i=1,7),j=1170,1171)/
1936 c    *-----------k0s()--------------------------
1937      *   20., .68610, 120.,-120.,   0.,   0.,   0.
1938      *,  20.,1.00000, 110., 110.,   0.,   0.,   0./
1939       data ((dectab(i,j),i=1,7),j=1172,1192)/
1940 c    *-----------k0l-------------------------------
1941      *  320., .2113, 110., 110., 110.,   0.,   0.
1942      *, 320., .2113, 110., 110., 110.,   0.,   0.
1943      *, 320., .2120, 110., 110., 110.,   0.,   0.
1944      *, 320., .3380, 120.,-120., 110.,   0.,   0.
1945      *, 320., .4744, 120.,  14., -13.,   0.,   0.
1946      *, 320., .6108,-120., -14.,  13.,   0.,   0.
1947      *, 320., .8054, 120.,  12., -11.,   0.,   0.
1948      *, 320.,1.0000,-120., -12.,  11.,   0.,   0.
1949 c    *-----------k+-------------------------------
1950      *, 130., .6352 , -14.,  13.,   0.,   0.,   0.
1951      *, 130., .8468 , 120., 110.,   0.,   0.,   0.
1952      *, 130., .9027 , 120., 120.,-120.,   0.,   0.
1953      *, 130., .92   , 120., 110., 110.,   0.,   0.
1954      *, 130., .9518 , 110., -14.,  13.,   0.,   0.
1955      *, 130.,1.     , 110., -12.,  11.,   0.,   0.
1956 c    *-----------pi+------------------------------
1957      *, 120., 1.    , -14.,  13.,   0.,   0.,   0.
1958 c    *-----------mu-------------------------------
1959      *,  14., 1.    ,  12., -11.,  13.,   0.,   0.
1960 c    *-----------etac-------------------------------
1961      *, 440.,  .32  , 230.,-230., 110.,   0.,   0.
1962      *, 440.,  .64  , 220., 110., 110.,   0.,   0.
1963      *, 440.,  .76  , 120.,-120., 130.,-130.,   0.
1964      *, 440.,  .88  , 120.,-120., 120.,-120.,   0.
1965      *, 440., 1.    , 130.,-130., 130.,-130.,   0./
1966
1967       call idresi
1968
1969 c     determine wmass2,wgam2
1970 c     ----------------------
1971       alfa=1./137.036
1972       gf=1.16570e-5
1973       sin2w=.215
1974       sinw=sqrt(sin2w)
1975 c      cosw=sqrt(1.-sin2w)           !?????????????????unused
1976       amw=sqrt(pi*alfa/(.9304*sqrt(2.)*gf))/sinw
1977       wmass2=amw
1978       call idmass(5,amlep5)
1979       call idmass(6,amlep6)
1980       ngam=12
1981       if(amlep5+amlep6.gt.amw) ngam=9
1982       wgam2=gf*amw**3/(6.*pi*sqrt(2.))*ngam
1983
1984       data iblank/' '/
1985       ird=0
1986       do 1 i=1,mxlook
1987 1     look(i)=0
1988       do 2 i=1,mxdky
1989       do 3 j=1,5
1990 3     mode(j,i)=0
1991 2     cbr(i)=0.
1992       nodcay=.false.
1993       noeta=.false.
1994       nopi0=.false.
1995       nonunu=.false.
1996       noevol=.false.
1997       nohadr=.false.
1998       if(lprint) write(ifch,10)
1999 10    format('1',30('*')/' *',28x,'*'/
2000      1' *',5x,'isajet decay table',5x,'*'/
2001      2' *',28x,'*'/' ',30('*')//
2002      36x,'part',18x,'decay mode',19x,'cum br',15x,'ident',17x,
2003      4'decay ident')
2004       loop=0
2005       iold=0
2006       if(nodcay) return
2007
2008 200   loop=loop+1
2009       if(loop.gt.mxdky) goto9999
2010 220   do 210 i=1,5
2011       imode(i)=0
2012       lmode(i)=iblank
2013 210   continue
2014       ird=ird+1
2015       if(ird.gt.ndectb)return
2016 c      if(ird.gt.1171)return   ! ??????????????????????????
2017       ires=nint(dectab(1,ird))
2018       br=dectab(2,ird)
2019       do 215 i=1,5
2020 215   imode(i)=nint(dectab(2+i,ird))
2021       if(nopi0.and.ires.eq.110) goto220
2022       if(noeta.and.ires.eq.220) goto220
2023       if(ires.eq.iold) goto230
2024       if(ires.lt.0.or.ires.gt.mxlook)
2025      *call utstop('hdecin: ires out of range&')
2026       look(ires)=loop
2027 230   iold=ires
2028       cbr(loop)=br
2029       do 240 i=1,5
2030       mode(i,loop)=imode(i)
2031       if(imode(i).ne.0) lmode(i)=idlabl(imode(i))
2032 240   continue
2033       lres=idlabl(ires)
2034       if(lprint) write(ifch,20) lres,(lmode(k),k=1,5),
2035      1br,ires,(imode(k),k=1,5)
2036 20    format(6x,a5,6x,5(a5,2x),3x,f8.5,15x,i5,4x,5(i5,2x))
2037       goto200
2038
2039 9999  write(ifch,*)'loop=', loop
2040       call utstop('hdecin: loop > mxdky&')
2041
2042       end
2043
2044 C -----------------------------------------------
2045       FUNCTION FM(NQ,S0,S1,S2,S3)
2046 C -----------------------------------------------
2047 C Normalized TRANSITION MATRIX FOR THE DALIZT PLOT DISTRI.
2048 C OF K -> 3 PIONS. PARAMETRIZATION OF WEINBERG
2049 C AS DESCRIBE IN PARTICLE DATA BOOK.
2050 C G IS THE LINEAR COEFFICIENT (SLOPE g)
2051 C H IS THE QUADRATIC COEFFICIENT h
2052 C D IS THE QUADRATIC COEFFICIENT k
2053 C Amax is the maximum of this amplitude (taken from Corsika by D. Heck)
2054 C NQ is the decay channel :
2055 C   1 - K -> 3 Pi
2056 C   2 - K -> Pi + 2 Pi0
2057 C   3 - K0 -> 3 Pi0
2058 C   4 - K0 -> 2 Pi + Pi0
2059 C -----------------------------------------------
2060       DIMENSION G(4),H(4),D(4),Amax(4)
2061       PARAMETER (PIM=139.57E-3)
2062       DATA G/-0.2154,0.594,0.,0.67/
2063       DATA H/0.01,0.035,0.,0.079/
2064       DATA D/-0.01,0.,0.,0.0098/
2065       DATA Amax/1.27,1.84,1.,2.22/
2066
2067       FM=1.+G(NQ)*(S3-S0)/(PIM*PIM)+H(NQ)*((S3-S0)/(PIM*PIM))**2
2068      *+D(NQ)*((S2-S1)/(PIM*PIM))**2
2069       FM=FM/Amax(NQ)
2070
2071       RETURN
2072       END
2073 C -----------------------------------------------
2074       FUNCTION FML(N,AM,RM1,RM2,E1S,E2S,E3S)
2075 C -----------------------------------------------
2076 C Normalized DALITZ PLOT DENSITY (RHO)
2077 C OF K -> 1 PION + 2 LEPTONS
2078 C AS DESCRIBE IN PARTICLE DATA BOOK.
2079 C CLP IS THE LAMBDA + FORM FACTOR COEFFICIENT
2080 C CLN IS THE LAMBDA 0 FORM FACTOR COEFFICIENT
2081 C EEP IS E'pion
2082 C GP IS THE F+(t) FORM FACTOR (t=AM*AM+SM1-2.D0*AM*E1S)
2083 C H IS EPS(t)=F-(t)/F+(t) WHERE F- IS CALCULATED FROM F0
2084 C Amax is the maximum of this density (taken from Corsika by D. Heck)
2085 C N is the decay channel :
2086 C   1 - K -> Pi0 + e + Nu
2087 C   2 - K -> Pi0 + Mu + Nu
2088 C   3 - K0 -> Pi + e + Nu
2089 C   4 - K0 -> Pi + Mu + Nu
2090 C -----------------------------------------------
2091       DIMENSION CLP(4),CLN(4),Amax(4)
2092       DATA CLP/0.0276,0.031,0.0288,0.034/
2093       DATA CLN/0.0,0.006,0.,0.025/
2094       DATA Amax/1.28e-2,1.194e-2,1.31e-2,1.241e-2/
2095
2096       SM1=RM1*RM1
2097       SM2=RM2*RM2
2098       EEP=0.5D0*(AM*AM+SM1-SM2)/AM-E1S
2099       GP=1.+CLP(N)*(AM*AM+SM1-2.*AM*E1S)/SM1
2100       H=(AM*AM-SM1)/SM1*(CLN(N)-CLP(N))/GP
2101       FML=GP*GP*(AM*(2.*E2S*E3S-AM*EEP)+
2102      *SM2*(0.25*EEP-E3S)+H*SM2*(E3S-0.5*EEP)+
2103      *0.25*H*H*SM2*EEP)
2104       FML=FML/Amax(N)
2105       RETURN
2106       END
2107 C -----------------------------------------------
2108       FUNCTION FMU(X)
2109 C -----------------------------------------------
2110 C PROBABILITY DISTRI. FOR ELECTRON ENERGY FROM MUON DECAY :
2111 C MU -> 2NU + E. DESCRIBE IN PARTICLE DATA BOOK.
2112 C (SIMPLIFY DIFFERENTIAL DECAY RATE INTEGRATED)
2113 C X REDUCED ENERGY OF PARTICLE
2114 C -----------------------------------------------
2115
2116       FMU=2.*(3.-2.*X)*X*X
2117
2118       RETURN
2119       END