HERWIG/src/hwudpr.f

   1
   2 CDECK  ID>, HWUDPR.
   3
   4 *CMZ :-        -27/07/99  13.33.03  by  Mike Seymour
   5
   6 *-- Author :    Ian Knowles, Bryan Webber & Kosuke Odagiri
   7
   8 C-----------------------------------------------------------------------
   9
  10       SUBROUTINE HWUDPR
  11
  12 C-----------------------------------------------------------------------
  13
  14 C     Prints out particle properies/decay tables in a number of formats:
  15
  16 C     If (PRNDEF) ASCII to stout
  17
  18 C     If (PRNTEX) LaTeX to the file HW_decays.tex
  19
  20 C                 Paper size and offsets as set in HWUEPR
  21
  22 C                 Uses the package longtable.sty
  23
  24 C                 Designed to be printed as landscape
  25
  26 C     If (PRNWEB) HTML  to the file HW_decays/index.html
  27
  28 C                                            /PART0000001.html etc.
  29
  30 C-----------------------------------------------------------------------
  31
  32       INCLUDE 'HERWIG61.INC'
  33
  34       INTEGER MMWIDE,MMLONG,MMHOFF,MMVOFF,IUNITT,IUNTW1,IUNTW2,I,NM,J,K,
  35
  36      & L,M
  37
  38       CHARACTER*1  Z
  39
  40       CHARACTER*2  ZZ,ACHRG
  41
  42       CHARACTER*3  ASPIN(0:10)
  43
  44       CHARACTER*6  BGCOLS(5),TBCOLS(3)
  45
  46       CHARACTER*7  HWUNST,TMPNME
  47
  48       CHARACTER*17 FNAMEP
  49
  50       CHARACTER*33 FNAMEW
  51
  52       COMMON/PAPER/MMWIDE,MMLONG,MMHOFF,MMVOFF
  53
  54       EXTERNAL HWUNST
  55
  56       DATA BGCOLS/'ffffff','0000aa','aa0000','00aa00','aa00ff'/
  57
  58       DATA TBCOLS/'ccccff','9966ff','ffff00'/
  59
  60       DATA ASPIN/' 0 ','1/2',' 1 ','3/2',' 2 ','5/2',' 3 ','7/2',
  61
  62      &           ' 4 ','9/2',' 5 '/
  63
  64 C
  65
  66       Z=CHAR(92)
  67
  68       ZZ=Z//Z
  69
  70 C
  71
  72       IUNITT=50
  73
  74       IUNTW1=51
  75
  76       IUNTW2=52
  77
  78 C Open and write out file header information for index file
  79
  80       IF (PRNDEF) THEN
  81
  82         IF (NPRFMT.LE.1) THEN
  83
  84           WRITE (6,10) NRES
  85
  86         ELSE
  87
  88           WRITE (6,20) NRES
  89
  90         END IF
  91
  92       END IF
  93
  94       IF (PRNTEX) THEN
  95
  96         OPEN(IUNITT,STATUS='UNKNOWN',FILE='HW_decays.tex')
  97
  98         IF (NPRFMT.LE.1) THEN
  99
 100           WRITE(IUNITT,30) Z,Z,Z,MMWIDE,Z,MMLONG,Z,MMHOFF,Z,MMVOFF,
 101
 102      &     Z,Z,Z,Z,Z,ZZ,Z,Z,Z,Z,Z,Z,NRES,ZZ,Z,Z,ZZ,Z,Z
 103
 104         ELSE
 105
 106           WRITE(IUNITT,40) Z,Z,Z,MMLONG,Z,MMWIDE,Z,MMHOFF,Z,MMVOFF,
 107
 108      &     Z,Z,Z,Z,Z,Z,ZZ,Z,Z,Z,Z,Z,Z,NRES,ZZ,Z,Z,Z,ZZ,Z,Z
 109
 110         END IF
 111
 112       ENDIF
 113
 114       IF (PRNWEB) THEN
 115
 116         OPEN(IUNTW1,STATUS='UNKNOWN',FILE='HW_decays/index.html')
 117
 118         WRITE(IUNTW1,50) BGCOLS,TBCOLS,NRES,((TBCOLS(I),I=2,3),J=1,7)
 119
 120       ENDIF
 121
 122    10 FORMAT(1H1//15X,'TABLE OF PROPERTIES OF',I4,' PARTICLES USED'/)
 123
 124    20 FORMAT(1H1//30X,'TABLE OF PROPERTIES OF',I4,' PARTICLES USED'//
 125
 126      & 5X,'Name        IDPDG   Mass   Chg Spn Lifetime Modes ',
 127
 128      & ' Branching fractions ME codes and decay products')
 129
 130    30 FORMAT(A1,'documentclass{article}'/A1,'usepackage{longtable}'/
 131
 132      & A1,'textwidth ',I4,'mm ',A1,'textheight ',I4,'mm'/
 133
 134      & A1,'hoffset   ',I4,'mm ',A1,'voffset    ',I4,'mm'/
 135
 136      & A1,'pagestyle{empty}'/A1,'begin{document}'/
 137
 138      & A1,'begin{center}'/A1,'begin{longtable}{|r|c|r|r|r|r|r|r|}'/
 139
 140      & A1,'hline'/'Id HW & Name & Id PDG & Mass & Charge & Spin ',
 141
 142      & '& Lifetime & Modes ',A2/A1,'hline'/
 143
 144      & A1,'endhead'/A1,'hline'/A1,'endfoot'/A1,'hline'/
 145
 146      & A1,'multicolumn{8}{|c|}{HERWIG 6.0: Table of properties',
 147
 148      & ' of the ',I3,' particles used} ',A2/A1,'hline',A1,'hline'/
 149
 150      & 'Id HW & Name & Id PDG & Mass & Charge & Spin & ',
 151
 152      & 'Lifetime & Modes ',A2/A1,'hline'/A1,'endfirsthead')
 153
 154    40 FORMAT(A1,'documentclass{article}'/A1,'usepackage{longtable}'/
 155
 156      & A1,'textwidth ',I4,'mm ',A1,'textheight ',I4,'mm'/
 157
 158      & A1,'hoffset   ',I4,'mm ',A1,'voffset    ',I4,'mm'/
 159
 160      & A1,'pagestyle{empty}'/A1,'begin{document}'/A1,'begin{center}'/
 161
 162      & A1,'begin{longtable}{|r|c|r|r|r|r|r|r|c|r|ccccc|}'/
 163
 164      & A1,'hline'/'Id HW & Name & Id PDG & Mass & Charge & Spin ',
 165
 166      & '& Lifetime & Modes & B.R. & M.E. & ' /
 167
 168      & A1,'multicolumn{5}{|c|}{Decay Products} ',A2/A1,'hline'/
 169
 170      & A1,'endhead'/A1,'hline'/A1,'endfoot'/A1,'hline'/
 171
 172      & A1,'multicolumn{15}{|c|}{HERWIG 6.0: Table of properties',
 173
 174      & ' of the ',I3,' particles used} ',A2/A1,'hline',A1,'hline'/
 175
 176      & 'Id HW & Name & Id PDG & Mass & Charge & Spin & ',
 177
 178      & 'Lifetime & Modes & B.R. & M.E. & '/
 179
 180      & A1,'multicolumn{5}{|c|}{Decay Products} ',A2/A1,'hline'/
 181
 182      & A1,'endfirsthead')
 183
 184    50 FORMAT('<!-- No Less productions -->'/'<HTML>'/'<HEAD>'/
 185
 186      & '<TITLE>HERWIG 6.0 Particle Properties</TITLE>'/'</HEAD>'/
 187
 188      & '<BODY BGCOLOR=#',A6,' TEXT=#',A6,' LINK=#',A6,
 189
 190      & ' ALINK=#',A6,' VLINK=#',A6,'>'/'<CENTER>'/
 191
 192      & '<TABLE ALIGN="CENTER" BGCOLOR=#',A6,'>',
 193
 194      & '<TR>'/'<TH COLSPAN=8 BGCOLOR=#',A6,' ALIGN="CENTER">',
 195
 196      & '<A HREF=="http://hepwww.rl.ac.uk/theory/seymour/herwig/">',
 197
 198      & 'HERWIG 6.0:</A><FONT COLOR=#',A6,'> Table of properties of',
 199
 200      & ' the ',I3,' particles used</FONT></TH>'/'<TR>'/'<TH></TH>'/
 201
 202      & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Name</FONT></TH>'/
 203
 204      & '<TH BGCOLOR=#',A6,' ALIGN="CENTER"><FONT COLOR=#',A6,'>',
 205
 206      & 'Id PDG</FONT></TH>'/
 207
 208      & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Mass</FONT></TH>'/
 209
 210      & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Charge</FONT></TH>'/
 211
 212      & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Spin</FONT></TH>'/
 213
 214      & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Lifetime</FONT></TH>'/
 215
 216      & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Modes</FONT></TH>'/
 217
 218      & '</TR>')
 219
 220 C Loop through resonances
 221
 222       DO 260 I=1,NRES
 223
 224 C Skip particles that can't be produced or blank lines
 225
 226       IF ((VTOCDK(I).AND.VTORDK(I)).OR.
 227
 228      &    (RNAME(I).EQ.'        ')) GOTO 260
 229
 230 C Open and write out header information for particle file
 231
 232       IF (PRNWEB) THEN
 233
 234         TMPNME = HWUNST(I)
 235
 236         WRITE(FNAMEP,'(A5,A7,A5)') 'PART_',TMPNME,'.html'
 237
 238         WRITE(FNAMEW,'(A,A17)') 'HW_decays/',FNAMEP
 239
 240         OPEN(IUNTW2,STATUS='UNKNOWN',FILE=FNAMEW)
 241
 242         WRITE(IUNTW2,60) RNAME(I),BGCOLS
 243
 244         WRITE(IUNTW2,70) TBCOLS,((TBCOLS(L),L=2,3),M=1,6)
 245
 246       ENDIF
 247
 248    60 FORMAT('<!-- No Less productions -->'/'<HTML>'/'<HEAD>'/
 249
 250      & '<TITLE>HERWIG 6.0: ',A8,' properties</TITLE>'/'</HEAD>'/
 251
 252      & '<BODY BGCOLOR=#',A6,' TEXT=#',A6,' LINK=#',A6,
 253
 254      & ' ALINK=#',A6,' VLINK=#',A6,'>'/'<CENTER>')
 255
 256    70 FORMAT('<TABLE ALIGN="CENTER" BGCOLOR=#',A6,'>'/
 257
 258      & '<TR>'/'<TH></TH>'/
 259
 260      & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Name</FONT></TH>'/
 261
 262      & '<TH BGCOLOR=#',A6,' ALIGN="CENTER"><FONT COLOR=#',A6,
 263
 264      & '>Id PDG</FONT></TH>'/
 265
 266      & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Mass</FONT></TH>'/
 267
 268      & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Charge</FONT></TH>'/
 269
 270      & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Spin</FONT></TH>'/
 271
 272      & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Lifetime</FONT></TH>'/
 273
 274      & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Modes</FONT></TH>'/
 275
 276      & '</TR>')
 277
 278 C Trick to output charge in fractions for di/s - quarks
 279
 280       IF ((I.GE.  1.AND.I.LE. 12).OR.(I.GE.109.AND.I.LE.120).OR.
 281
 282      &    (I.GE.209.AND.I.LE.218).OR.(I.GE.401.AND.I.LE.424)) THEN
 283
 284         ACHRG='/3'
 285
 286       ELSE
 287
 288         ACHRG='  '
 289
 290       ENDIF
 291
 292 C Write out special particles with no decay modes
 293
 294       IF (NMODES(I).EQ.0) THEN
 295
 296         IF (PRNDEF) THEN
 297
 298           IF (NPRFMT.LE.1) THEN
 299
 300             WRITE(6,80) I,RNAME(I),IDPDG(I),RMASS(I),ICHRG(I),
 301
 302      &                  ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),0
 303
 304           ELSE
 305
 306             WRITE(6,90) I,RNAME(I),IDPDG(I),RMASS(I),ICHRG(I),
 307
 308      &                  ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),0
 309
 310           ENDIF
 311
 312         ENDIF
 313
 314 C Add particle to LaTeX file
 315
 316         IF (PRNTEX) THEN
 317
 318           IF (NPRFMT.LE.1) THEN
 319
 320             WRITE(IUNITT,100) Z,I,TXNAME(1,I),IDPDG(I),RMASS(I),
 321
 322      &       ICHRG(I),ACHRG,ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),0,ZZ
 323
 324           ELSE
 325
 326             WRITE(IUNITT,110) Z,I,TXNAME(1,I),IDPDG(I),RMASS(I),
 327
 328      &       ICHRG(I),ACHRG,ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),0,Z,ZZ
 329
 330           ENDIF
 331
 332         ENDIF
 333
 334         IF (PRNWEB) THEN
 335
 336 C Add properties to Web index
 337
 338           WRITE(IUNTW1,120) TBCOLS(2),TBCOLS(3),I,FNAMEP,TXNAME(2,I),
 339
 340      &                      IDPDG(I),RMASS(I),ICHRG(I),ACHRG,
 341
 342      &                      ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),0
 343
 344 C Add properties to Web particle file
 345
 346           WRITE(IUNTW2,130) TBCOLS(2),TBCOLS(3),I,TXNAME(2,I),
 347
 348      &                      IDPDG(I),RMASS(I),ICHRG(I),ACHRG,
 349
 350      &                      ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),0
 351
 352         ENDIF
 353
 354    80   FORMAT(/1X,I3,1X,A8,' IDPDG=',I8,', M=',F8.3,', Q=',I2,',  J=',
 355
 356      &   A3,', T=',1P,E9.3,',',I3,' Modes')
 357
 358    90   FORMAT(/1X,I3,1X,A8,1X,I8,1X,F8.3,1X,I2,1X,A3,1X,1P,E9.3,1X,I3)
 359
 360   100   FORMAT(A1,'hline',I4,' & ',A37,' & $',I8,'$ & ',F8.3,' & $',I2,
 361
 362      &   A2,'$ & ',A3,' & $',1P,E9.3,'$ & ',I3,' ',A2)
 363
 364   110   FORMAT(A1,'cline{1-8}'/
 365
 366      &   I4,' & ',A37,' & $',I8,'$ & ',F8.3,' & $',I2,A2,'$ & ',A3,
 367
 368      &   ' & $',1P,E9.3,'$ & ',I3,' & ',A1,'multicolumn{7}{|c|}{} ',A2)
 369
 370   120   FORMAT('<TR>'/
 371
 372      &   '<TD ALIGN="RIGHT" BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>',I3,
 373
 374      &   '</FONT></TD>'/
 375
 376      &   '<TD ALIGN="CENTER"><A HREF="',A17,'">',A37,'</A></TD>'/
 377
 378      &   '<TD ALIGN="RIGHT">',I8,'</TD>'/
 379
 380      &   '<TD ALIGN="RIGHT">',F8.3,'</TD>'/
 381
 382      &   '<TD ALIGN="RIGHT">',I2,A2,'</TD>'/
 383
 384      &   '<TD ALIGN="RIGHT">',A3,'</TD>'/
 385
 386      &   '<TD ALIGN="RIGHT">',1P,E9.3,'</TD>'/
 387
 388      &   '<TD ALIGN="RIGHT">',I3,'</TD>'/'</TR>')
 389
 390   130   FORMAT('<TR>'/
 391
 392      &   '<TD ALIGN="RIGHT" BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>',I3,
 393
 394      &   '</FONT></TD>'/
 395
 396      &   '<TD ALIGN="CENTER">',A37,'</TD>'/
 397
 398      &   '<TD ALIGN="RIGHT">',I8,'</TD>'/
 399
 400      &   '<TD ALIGN="RIGHT">',F8.3,'</TD>'/
 401
 402      &   '<TD ALIGN="RIGHT">',I2,A2,'</TD>'/
 403
 404      &   '<TD ALIGN="RIGHT">',A3,'</TD>'/
 405
 406      &   '<TD ALIGN="RIGHT">',1P,E9.3,'</TD>'/
 407
 408      &   '<TD ALIGN="RIGHT">',I3,'</TD>'/'</TR>'/'</TABLE>'/'<P>')
 409
 410       ELSE
 411
 412 C Particle with decay modes
 413
 414         IF (RSTAB(I)) THEN
 415
 416           NM=0
 417
 418         ELSEIF (VTOCDK(I)) THEN
 419
 420           NM=-NMODES(I)
 421
 422         ELSE
 423
 424           NM=NMODES(I)
 425
 426         ENDIF
 427
 428         K=LSTRT(I)
 429
 430 C Write out properties and first decay mode
 431
 432         IF (PRNDEF) THEN
 433
 434           IF (NPRFMT.LE.1) THEN
 435
 436             WRITE(6, 80) I,RNAME(I),IDPDG(I),RMASS(I),ICHRG(I),
 437
 438      &       ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),NM
 439
 440             WRITE(6,140) (RNAME(IDKPRD(L,K)),L=1,5),BRFRAC(K),NME(K)
 441
 442           ELSE
 443
 444             WRITE(6,150) I,RNAME(I),IDPDG(I),RMASS(I),ICHRG(I),
 445
 446      &       ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),NM,BRFRAC(K),NME(K),
 447
 448      &       (RNAME(IDKPRD(L,K)),L=1,5)
 449
 450           ENDIF
 451
 452         ENDIF
 453
 454         IF (PRNTEX) THEN
 455
 456           IF (NPRFMT.LE.1) THEN
 457
 458             WRITE(IUNITT,160) Z,I,TXNAME(1,I),IDPDG(I),RMASS(I),
 459
 460      &       ICHRG(I),ACHRG,ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),NM,ZZ,Z
 461
 462             WRITE(IUNITT,170) Z,Z,(TXNAME(1,IDKPRD(L,K)),L=1,5),Z,
 463
 464      &       BRFRAC(K),Z,NME(K),ZZ
 465
 466           ELSE
 467
 468             WRITE(IUNITT,180) Z,I,TXNAME(1,I),IDPDG(I),RMASS(I),
 469
 470      &       ICHRG(I),ACHRG,ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),NM,
 471
 472      &       BRFRAC(K),NME(K),(TXNAME(1,IDKPRD(L,K)),L=1,5),ZZ,Z
 473
 474           ENDIF
 475
 476         END IF
 477
 478         IF (PRNWEB) THEN
 479
 480 C Add properties to index
 481
 482           WRITE(IUNTW1,120) TBCOLS(2),TBCOLS(3),I,FNAMEP,TXNAME(2,I),
 483
 484      &     IDPDG(I),RMASS(I),ICHRG(I),ACHRG,ASPIN(INT(TWO*RSPIN(I))),
 485
 486      &     RLTIM(I),NM
 487
 488 C Add properties to Web particle file
 489
 490           WRITE(IUNTW2,130) TBCOLS(2),TBCOLS(3),I,TXNAME(2,I),IDPDG(I),
 491
 492      &     RMASS(I),ICHRG(I),ACHRG,ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),NM
 493
 494           WRITE(IUNTW2,190) TBCOLS,TXNAME(2,I),
 495
 496      &     ((TBCOLS(L),L=2,3),M=1,3)
 497
 498           WRITE(IUNTW2,200) TBCOLS(2),TBCOLS(3),1,BRFRAC(K),NME(K),
 499
 500      &     (TXNAME(2,IDKPRD(L,K)),L=1,5)
 501
 502         ENDIF
 503
 504   140   FORMAT(5X,'BR[ -->',5(1X,A6),']=',F5.3,', ME code',I4)
 505
 506   150   FORMAT(/1X,I3,1X,A8,1X,I8,1X,F8.3,1X,I2,1X,A3,1X,1P,E9.3,1X,I3,
 507
 508      &   2X,F5.3,1X,I3,5(1X,A8))
 509
 510   160   FORMAT(A1,'hline',
 511
 512      &   I4,' & ',A37,' & $',I8,'$ & ',F8.3,' & $',I2,A2,'$ & ',
 513
 514      &   A3,' & $',1P,E9.3,'$ & ',I3,' ',A2/A1,'cline{2-8}')
 515
 516   170   FORMAT(' & & ',A1,'multicolumn{2}{l}{$',A1,'longrightarrow$'/
 517
 518      &   5(A37,' '),'}'/' & ',A1,'multicolumn{2}{l}{BR = ',F5.3,'} & ',
 519
 520      &   A1,'multicolumn{2}{l|}{ME code = ',I3,'} ',A2)
 521
 522   180   FORMAT(A1,'hline'/
 523
 524      &   I4,' & ',A37,' & $',I8,'$ & ',F8.3,' & $',I2,A2,'$ & ',
 525
 526      &   A3,' & $',1P,E9.3,'$ & ',I3,' & ',F5.3,' & ',I3,
 527
 528      &   5(' & ',A37), ' ',A2/A1,'cline{2-8}')
 529
 530   190   FORMAT('<TABLE ALIGN="CENTER" BGCOLOR=#',A6,'>'/'<TR>'/
 531
 532      &   '<TH COLSPAN=8 BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>',A37,
 533
 534      &   ' Decay Modes</FONT></TH>'/'</TR>'/'<TR>'/'<TH></TH>',
 535
 536      &   '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>B.R.</FONT></TH>'/
 537
 538      &   '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>M.E.</FONT></TH>'/
 539
 540      &   '<TH BGCOLOR=#',A6,' ALIGN="CENTER" COLSPAN=5>',
 541
 542      &   '<FONT COLOR=#',A6,'>Decay products</FONT></TH>'/'</TR>')
 543
 544   200   FORMAT('<TR>'/
 545
 546      &   '<TD ALIGN="RIGHT" BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>',
 547
 548      &   I3,'</FONT></TD>'/
 549
 550      &   '<TD ALIGN="RIGHT">',F5.3,'</TD>'/
 551
 552      &   '<TD ALIGN="RIGHT">',I3,'</TD>'/
 553
 554      &   5('<TD ALIGN="CENTER">',A37,'</TD>'/),'</TR>')
 555
 556 C Write out additional decay modes
 557
 558         IF (NMODES(I).GE.2) THEN
 559
 560           DO 210 J=2,NMODES(I)
 561
 562           K=LNEXT(K)
 563
 564           IF (PRNDEF) THEN
 565
 566             IF (NPRFMT.LE.1) THEN
 567
 568               WRITE(6,140) (RNAME(IDKPRD(L,K)),L=1,5),BRFRAC(K),NME(K)
 569
 570             ELSE
 571
 572               WRITE(6,220) BRFRAC(K),NME(K),(RNAME(IDKPRD(L,K)),L=1,5)
 573
 574             END IF
 575
 576           END IF
 577
 578           IF (PRNTEX) THEN
 579
 580             IF (NPRFMT.LE.1) THEN
 581
 582               WRITE(IUNITT,170) Z,Z,(TXNAME(1,IDKPRD(L,K)),L=1,5),Z,
 583
 584      &         BRFRAC(K),Z,NME(K),ZZ
 585
 586             ELSE
 587
 588               WRITE(IUNITT,230) Z,BRFRAC(K),NME(K),
 589
 590      &         (TXNAME(1,IDKPRD(L,K)),L=1,5),ZZ
 591
 592             ENDIF
 593
 594           ENDIF
 595
 596           IF (PRNWEB) WRITE(IUNTW2,200) TBCOLS(2),TBCOLS(3),J,
 597
 598      &     BRFRAC(K),NME(K),(TXNAME(2,IDKPRD(L,K)),L=1,5)
 599
 600   210     CONTINUE
 601
 602           IF (PRNTEX.AND.NPRFMT.EQ.2.AND.NMODES(I+1).EQ.0)
 603
 604      &     WRITE(IUNITT,240) Z
 605
 606   220     FORMAT(54X,F5.3,1X,I3,5(1X,A8))
 607
 608   230     FORMAT(' & ',A1,'multicolumn{7}{|c|}{} & ',F5.3,' & ',I3,
 609
 610      &     5(' & ',A37),' ',A2)
 611
 612   240     FORMAT(A1,'hline')
 613
 614         ENDIF
 615
 616       ENDIF
 617
 618 C Close Web particle file
 619
 620       IF (PRNWEB) THEN
 621
 622         WRITE(IUNTW2,250)
 623
 624         CLOSE(IUNTW2)
 625
 626       ENDIF
 627
 628   250 FORMAT('</TABLE>'/'</CENTER>'/'<P>'/
 629
 630      & 'Main particle <A HREF="index.html">index</A>'/
 631
 632      & '</BODY>'/'</HTML>')
 633
 634   260 CONTINUE
 635
 636 C Close the LaTeX file
 637
 638       IF (PRNTEX) THEN
 639
 640         WRITE(IUNITT,270) Z,Z,Z
 641
 642         CLOSE(IUNITT)
 643
 644       ENDIF
 645
 646 C Close the index file
 647
 648       IF (PRNWEB) THEN
 649
 650         WRITE(IUNTW1,280)
 651
 652         CLOSE(IUNTW1)
 653
 654       ENDIF
 655
 656   270 FORMAT(A1,'end{longtable}'/A1,'end{center}'/A1,'end{document}')
 657
 658   280 FORMAT('</TABLE>'/'</CENTER>'/'</BODY>'/'</HTML>')
 659
 660       RETURN
 661
 662       END