Call to LHAPDF initialisation via pdfset_herwig

[u/mrichter/AliRoot.git] / HERWIG / herwig6510.f

diff --git a/HERWIG/herwig6510.f b/HERWIG/herwig6510.f
index fd8b4743a83084a078bddc61f90f180e32893977..11113180f8f6977f028ec0798932b20de203ef5d 100644 (file)
--- a/HERWIG/herwig6510.f
+++ b/HERWIG/herwig6510.f
@@ -152,7 +152,7 @@ C used by TAUOLA
 C
 C written by Martin W. Gruenewald (91/01/28)
 C ----------------------------------------------------------------------
 C
 C written by Martin W. Gruenewald (91/01/28)
 C ----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       LOGICAL QEDRAD
       COMMON /PHORAD/ QEDRAD(NMXHEP)
       INTEGER N,IHEP,IST,ID,JMO1,JMO2,JDA1,JDA2,I,IP
       LOGICAL QEDRAD
       COMMON /PHORAD/ QEDRAD(NMXHEP)
       INTEGER N,IHEP,IST,ID,JMO1,JMO2,JDA1,JDA2,I,IP


@@ -256,7 +256,7 @@ C-----------------------------------------------------------------------
 C     Azimuthal correlation functions for Collins' algorithm,
 C     see I.G.Knowles, Comp. Phys. Comm. 58 (90) 271 for notation.
 C-----------------------------------------------------------------------
 C     Azimuthal correlation functions for Collins' algorithm,
 C     see I.G.Knowles, Comp. Phys. Comm. 58 (90) 271 for notation.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION Z1,Z2,DOT12,DOT23,DOT31,TR,FN(7),VEC1(2),VEC2(2),
      & VEC3(2),VEC(2)
       INTEGER IPAR,JPAR
       DOUBLE PRECISION Z1,Z2,DOT12,DOT23,DOT31,TR,FN(7),VEC1(2),VEC2(2),
      & VEC3(2),VEC(2)
       INTEGER IPAR,JPAR


@@ -330,7 +330,7 @@ C     MODIFIED 11/01/01 BY PR  FOR SPIN CORRELATIONS(PROBLEM WITH ORDER
 C                                                    OF DECAYS)
 C     NEW VARAIBLE BACK TO ALLOW CODE TO SEARCH DOWN CHAIN
 C-----------------------------------------------------------------------
 C                                                    OF DECAYS)
 C     NEW VARAIBLE BACK TO ALLOW CODE TO SEARCH DOWN CHAIN
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IHEP,IST,ID,JC,KC,JD,JHEP,LHEP,ID2,NTRY,KHEP
       LOGICAL BACK
       IF (IERROR.NE.0) RETURN
       INTEGER IHEP,IST,ID,JC,KC,JD,JHEP,LHEP,ID2,NTRY,KHEP
       LOGICAL BACK
       IF (IERROR.NE.0) RETURN


@@ -508,7 +508,7 @@ C
 C********MODIFIED 13/11/00 BY BRW TO ALLOW MULTIPLE APPLICATION IN
 C*******SAME EVENT (FOR WW AND ZZ) N.B. NO CLEANUP CALLS FOR THESE!
 C-----------------------------------------------------------------------
 C********MODIFIED 13/11/00 BY BRW TO ALLOW MULTIPLE APPLICATION IN
 C*******SAME EVENT (FOR WW AND ZZ) N.B. NO CLEANUP CALLS FOR THESE!
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWBVMC,HWRGEN,HWUALF,HWUSQR,X(3),W,WMAX,WSUM,
      & X1MIN,X1MAX,X2MIN,X2MAX,QSCALE,GAMFAC,GLUFAC,R(3,3),CS,SN,M(3),
      & E(3),LAMBDA,A,B,C,PTSQ,EM,P1(5),P2(5),PVRT(4),EPS,MASDEP
       DOUBLE PRECISION HWBVMC,HWRGEN,HWUALF,HWUSQR,X(3),W,WMAX,WSUM,
      & X1MIN,X1MAX,X2MIN,X2MAX,QSCALE,GAMFAC,GLUFAC,R(3,3),CS,SN,M(3),
      & E(3),LAMBDA,A,B,C,PTSQ,EM,P1(5),P2(5),PVRT(4),EPS,MASDEP


@@ -779,7 +779,7 @@ C     FILL MISSING AREA OF DIS PHASE-SPACE WITH 2+1-JET EVENTS
 C     IF (IOPT.EQ.1) SET UP EVENT RECORD
 C     IF (IOPT.EQ.2) CLEAN UP EVENT RECORD AFTER SHOWERING
 C-----------------------------------------------------------------------
 C     IF (IOPT.EQ.1) SET UP EVENT RECORD
 C     IF (IOPT.EQ.2) CLEAN UP EVENT RECORD AFTER SHOWERING
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWBVMC,HWUALF,HWULDO,P1(5),P2(5),P3(5),
      & PCMF(5),L(5),R(3,3),Q,XBJ,RN,XPMIN,XPMAX,XP,ZPMIN,ZPMAX,ZP,FAC,
      & X1,X2,XTSQ,XT,PTSQ,SIN1,SIN2,W1,W2,CFAC,PDFOLD(13),PDFNEW(13),
       DOUBLE PRECISION HWRGEN,HWBVMC,HWUALF,HWULDO,P1(5),P2(5),P3(5),
      & PCMF(5),L(5),R(3,3),Q,XBJ,RN,XPMIN,XPMAX,XP,ZPMIN,ZPMAX,ZP,FAC,
      & X1,X2,XTSQ,XT,PTSQ,SIN1,SIN2,W1,W2,CFAC,PDFOLD(13),PDFNEW(13),


@@ -1252,7 +1252,7 @@ C-----------------------------------------------------------------------
       SUBROUTINE HWBDYP(IOPT)
 C     MATRIX ELEMENT CORRECTIONS TO DRELL-YAN PROCESSES
 C-----------------------------------------------------------------------
       SUBROUTINE HWBDYP(IOPT)
 C     MATRIX ELEMENT CORRECTIONS TO DRELL-YAN PROCESSES
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWBVMC,HWRGEN,HWUALF,HWUSQR,PMODK,AZ,CZ,
      & T,U,S,EM,TMIN,TMAX,PMOD2,GLUFAC,SMIN,SMAX,SZ,TEST,
      & JAC,M(3),W1,W,PMOD3,SCAPR,CPHI,SPHI,SCALE,XI1,XI2,
       DOUBLE PRECISION HWBVMC,HWRGEN,HWUALF,HWUSQR,PMODK,AZ,CZ,
      & T,U,S,EM,TMIN,TMAX,PMOD2,GLUFAC,SMIN,SMAX,SZ,TEST,
      & JAC,M(3),W1,W,PMOD3,SCAPR,CPHI,SPHI,SCALE,XI1,XI2,


@@ -1862,7 +1862,7 @@ C-----------------------------------------------------------------------
 C     DELETES INTERNAL LINES FROM SHOWER, MAKES COLOUR CONNECTION INDEX
 C     AND COPIES INTO /HEPEVT/ IN COLOUR ORDER.
 C-----------------------------------------------------------------------
 C     DELETES INTERNAL LINES FROM SHOWER, MAKES COLOUR CONNECTION INDEX
 C     AND COPIES INTO /HEPEVT/ IN COLOUR ORDER.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IHEP,ID,IJET,KHEP,IPAR,JPAR,NXPAR,IP,JP
       IF (IERROR.NE.0) RETURN
 C---SAVE VIRTUAL PARTON DATA
       INTEGER IHEP,ID,IJET,KHEP,IPAR,JPAR,NXPAR,IP,JP
       IF (IERROR.NE.0) RETURN
 C---SAVE VIRTUAL PARTON DATA


@@ -1979,7 +1979,7 @@ C     BRANCHING GENERATOR WITH INTERFERING GLUONS
 C     HWBGEN EVOLVES QCD JETS ACCORDING TO THE METHOD OF
 C     G.MARCHESINI & B.R.WEBBER, NUCL. PHYS. B238(1984)1
 C-----------------------------------------------------------------------
 C     HWBGEN EVOLVES QCD JETS ACCORDING TO THE METHOD OF
 C     G.MARCHESINI & B.R.WEBBER, NUCL. PHYS. B238(1984)1
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWULDO,HWRGAU,EINHEP,ERTXI,RTXI,XF
       INTEGER NTRY,LASHEP,IHEP,NRHEP,ID,IST,JHEP,KPAR,I,J,IRHEP(NMXJET),
      & IRST(NMXJET),JPR
       DOUBLE PRECISION HWULDO,HWRGAU,EINHEP,ERTXI,RTXI,XF
       INTEGER NTRY,LASHEP,IHEP,NRHEP,ID,IST,JHEP,KPAR,I,J,IRHEP(NMXJET),
      & IRST(NMXJET),JPR


@@ -2194,7 +2194,7 @@ C     Makes the colour connections and performs the parton shower
 C     for events read in from the GUPI (Generic User Process Interface)
 C     event common block
 C----------------------------------------------------------------------
 C     for events read in from the GUPI (Generic User Process Interface)
 C     event common block
 C----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER MAXNUP
       PARAMETER (MAXNUP=500)
       INTEGER NUP,IDPRUP,IDUP,ISTUP,MOTHUP,ICOLUP
       INTEGER MAXNUP
       PARAMETER (MAXNUP=500)
       INTEGER NUP,IDPRUP,IDUP,ISTUP,MOTHUP,ICOLUP


@@ -2386,7 +2386,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     COMBINES JETS WITH REQUIRED KINEMATICS
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     COMBINES JETS WITH REQUIRED KINEMATICS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWULDO,EPS,PTX,PTY,PF,PTINF,PTCON,CN,CP,SP,PP0,
      & PM0,ET0,DET,ECM,EMJ,EMP,EMS,DMS,ES,DPF,ALF,AL(2),ET(2),PP(2),
      & PT(3),PA(5),PB(5),PC(5),PQ(5),PR(5),PS(5),RR(3,3),RS(3,3),ETC,
       DOUBLE PRECISION HWULDO,EPS,PTX,PTY,PF,PTINF,PTCON,CN,CP,SP,PP0,
      & PM0,ET0,DET,ECM,EMJ,EMP,EMS,DMS,ES,DPF,ALF,AL(2),ET(2),PP(2),
      & PT(3),PA(5),PB(5),PC(5),PQ(5),PR(5),PS(5),RR(3,3),RS(3,3),ETC,


@@ -2956,7 +2956,7 @@ C
 C     On output PPAR(1-5,*) (if TMPAR(*)), containts:
 C     (P-trans,Xi or Xilast,P-long,E,M)
 C-----------------------------------------------------------------------
 C     On output PPAR(1-5,*) (if TMPAR(*)), containts:
 C     (P-trans,Xi or Xilast,P-long,E,M)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWUSQR,EXI,PISQ,PJPK,EJEK,PTSQ,Z,ZMIN,ZMAX,
      $     EMI,EMJ,EMK,C,NQ,HWBVMC,RHO,POLD,PNEW,EOLD,ENEW,A,B
       INTEGER IPAR,JPAR,KPAR,MPAR,I,J,K
       DOUBLE PRECISION HWUSQR,EXI,PISQ,PJPK,EJEK,PTSQ,Z,ZMIN,ZMAX,
      $     EMI,EMJ,EMK,C,NQ,HWBVMC,RHO,POLD,PNEW,EOLD,ENEW,A,B
       INTEGER IPAR,JPAR,KPAR,MPAR,I,J,K


@@ -3070,7 +3070,7 @@ C-----------------------------------------------------------------------
 C     BRANCHES TIMELIKE PARTON KPAR INTO TWO, PUTS PRODUCTS
 C     INTO NPAR+1 AND NPAR+2, AND INCREASES NPAR BY TWO
 C-----------------------------------------------------------------------
 C     BRANCHES TIMELIKE PARTON KPAR INTO TWO, PUTS PRODUCTS
 C     INTO NPAR+1 AND NPAR+2, AND INCREASES NPAR BY TWO
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWBVMC,HWRGEN,HWUALF,HWUTAB,HWRUNI,HWULDO,PMOM,
      & QNOW,QLST,QKTHR,RN,QQBAR,DQQ,QGTHR,SNOW,QSUD,ZMIN,ZMAX,ZRAT,WMIN,
      & QLAM,Z1,Z2,ETEST,ZTEST,ENOW,XI,XIPREV,EPREV,QMAX,QGAM,SLST,SFNL,
       DOUBLE PRECISION HWBVMC,HWRGEN,HWUALF,HWUTAB,HWRUNI,HWULDO,PMOM,
      & QNOW,QLST,QKTHR,RN,QQBAR,DQQ,QGTHR,SNOW,QSUD,ZMIN,ZMAX,ZRAT,WMIN,
      & QLAM,Z1,Z2,ETEST,ZTEST,ENOW,XI,XIPREV,EPREV,QMAX,QGAM,SLST,SFNL,


@@ -3537,7 +3537,7 @@ C-----------------------------------------------------------------------
 C     SUBROUTINE TO REPLACE HWBCON IN RPARITY VIOLATING SUSY
 C     BASED ON HWBCON BY BRW
 C-----------------------------------------------------------------------
 C     SUBROUTINE TO REPLACE HWBCON IN RPARITY VIOLATING SUSY
 C     BASED ON HWBCON BY BRW
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IHEP,IST,ID,JC,JD,JHEP,IDP,IDM,IDM2,
      &        RHEP,IST2,ANTC,XHEP,IP,COLP
       LOGICAL BVVUSE,BVVHRD,BVDEC1,BVDEC2,COLRD,ACOLRD,BVHRD,BVHRD2,
       INTEGER IHEP,IST,ID,JC,JD,JHEP,IDP,IDM,IDM2,
      &        RHEP,IST2,ANTC,XHEP,IP,COLP
       LOGICAL BVVUSE,BVVHRD,BVDEC1,BVDEC2,COLRD,ACOLRD,BVHRD,BVHRD2,


@@ -4233,7 +4233,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C--Function to find the right daugther of a decaying gluino
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C--Function to find the right daugther of a decaying gluino
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER ID,JHEP,KC,JC
       LOGICAL COL,IFGO
 C---N.B. WILL NEED MODS WHEN SUSY PARTICLES CAN SHOWER
       INTEGER ID,JHEP,KC,JC
       LOGICAL COL,IFGO
 C---N.B. WILL NEED MODS WHEN SUSY PARTICLES CAN SHOWER


@@ -4275,7 +4275,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C--Function to search in the jet for the particle
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C--Function to search in the jet for the particle
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER JC,JD,QHEP,LHEP,IHEP,JHEP,IDM,NCOUNT,ID,IP,IDM2,COLP
       LOGICAL CON,BVVUSE,FLA,AFLA,BVVHRD
       FLA(IP)  = (IP.LE.6.OR.(IP.GE.115.AND.IP.LE.120).
       INTEGER JC,JD,QHEP,LHEP,IHEP,JHEP,IDM,NCOUNT,ID,IP,IDM2,COLP
       LOGICAL CON,BVVUSE,FLA,AFLA,BVVHRD
       FLA(IP)  = (IP.LE.6.OR.(IP.GE.115.AND.IP.LE.120).


@@ -4494,7 +4494,7 @@ C     jet started by parton no.2 interference partner 1 and spin density
 C     DECPAR(2). RHOPAR(2) gives the jet spin density matrix.
 C     See I.G. Knowles, Comp. Phys. Comm. 58 (90) 271.
 C-----------------------------------------------------------------------
 C     DECPAR(2). RHOPAR(2) gives the jet spin density matrix.
 C     See I.G. Knowles, Comp. Phys. Comm. 58 (90) 271.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,DMIN,PT,EIKON,EISCR,EINUM,EIDEN1,EIDEN2,
      & WT,SPIN,Z1,Z2,TR,PRMAX,CX,SX,CAZ,ROHEP(3),RMAT(3,3),ZERO2(2)
       INTEGER IPAR,JPAR,KPAR,LPAR,MPAR,JSTR,LSTR,MSTR
       DOUBLE PRECISION HWRGEN,DMIN,PT,EIKON,EISCR,EINUM,EIDEN1,EIDEN2,
      & WT,SPIN,Z1,Z2,TR,PRMAX,CX,SX,CAZ,ROHEP(3),RMAT(3,3),ZERO2(2)
       INTEGER IPAR,JPAR,KPAR,LPAR,MPAR,JSTR,LSTR,MSTR


@@ -4629,7 +4629,7 @@ C     Constructs appropriate spin density/decay matrix for parton
 C     in hard subprocess, otherwise zero. Assignments based upon
 C     Comp. Phys. Comm. 58 (1990) 271.
 C-----------------------------------------------------------------------
 C     in hard subprocess, otherwise zero. Assignments based upon
 C     Comp. Phys. Comm. 58 (1990) 271.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION C,V12,V23,V13,TR,C1,C2,C3,R1(2),R2(2)
       INTEGER IST
       SAVE R1,R2,V12
       DOUBLE PRECISION C,V12,V23,V13,TR,C1,C2,C3,R1(2),R2(2)
       INTEGER IST
       SAVE R1,R2,V12


@@ -4736,7 +4736,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     COMPUTES (OR READS) TABLES OF SUDAKOV FORM FACTORS
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     COMPUTES (OR READS) TABLES OF SUDAKOV FORM FACTORS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWUGAU,HWBVMC,HWBSUG,HWBSU1,HWBSU2,G1,G2,QRAT,
      & QLAM,POWER,AFAC,QMIN,QFAC,QNOW,ZMIN,ZMAX,Q1,QCOLD,VGOLD,VQOLD,
      & RMOLD(6),ACOLD,ZLO,ZHI
       DOUBLE PRECISION HWUGAU,HWBVMC,HWBSUG,HWBSU1,HWBSU2,G1,G2,QRAT,
      & QLAM,POWER,AFAC,QMIN,QFAC,QNOW,ZMIN,ZMAX,Q1,QCOLD,VGOLD,VQOLD,
      & RMOLD(6),ACOLD,ZLO,ZHI


@@ -4900,7 +4900,7 @@ C     LOGARITHMIC PART OF INTEGRAND IN EXPONENT OF SUDAKOV FORM FACTOR.
 C     THE SECOND ORDER ALPHAS CASE COMES FROM CONVERTING INTEGRAL OVER
 C     Q^2 INTO ONE OVER ALPHAS, WITH FLAVOUR THRESHOLDS.
 C-----------------------------------------------------------------------
 C     THE SECOND ORDER ALPHAS CASE COMES FROM CONVERTING INTEGRAL OVER
 C     Q^2 INTO ONE OVER ALPHAS, WITH FLAVOUR THRESHOLDS.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWBSUL,HWUALF,Z,QRAT,QLAM,U,AL,BL,QNOW,QMIN,
      & BET(6),BEP(6),MUMI(6),MUMA(6),ALMI(6),ALMA(6),FINT(6),ALFINT,
      & MUMIN,MUMAX,ALMIN,ALMAX
       DOUBLE PRECISION HWBSUL,HWUALF,Z,QRAT,QLAM,U,AL,BL,QNOW,QMIN,
      & BET(6),BEP(6),MUMI(6),MUMA(6),ALMI(6),ALMA(6),FINT(6),ALFINT,
      & MUMIN,MUMAX,ALMIN,ALMAX


@@ -4982,7 +4982,7 @@ C     Includes azimuthal angular correlations between branching planes
 C     due to spin (if AZSPIN) using the algorithm of Knowles & Collins.
 C     Ses Nucl. Phys. B304 (1988) 794 & Comp. Phys. Comm. 58 (1990) 271.
 C-----------------------------------------------------------------------
 C     due to spin (if AZSPIN) using the algorithm of Knowles & Collins.
 C     Ses Nucl. Phys. B304 (1988) 794 & Comp. Phys. Comm. 58 (1990) 271.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,DMIN,PT,EIKON,EINUM,EIDEN1,EIDEN2,EISCR,
      & WT,SPIN,Z1,Z2,PRMAX,CAZ,CX,SX,ROHEP(3),RMAT(3,3),ZERO2(2)
       INTEGER INITBR,INTERF,IPAR,JPAR,KPAR,LPAR,MPAR,NTRY,JOLD
       DOUBLE PRECISION HWRGEN,DMIN,PT,EIKON,EINUM,EIDEN1,EIDEN2,EISCR,
      & WT,SPIN,Z1,Z2,PRMAX,CAZ,CX,SX,ROHEP(3),RMAT(3,3),ZERO2(2)
       INTEGER INITBR,INTERF,IPAR,JPAR,KPAR,LPAR,MPAR,NTRY,JOLD


@@ -5136,7 +5136,7 @@ CDECK  ID>, HWBTOP.
 C-----------------------------------------------------------------------
       SUBROUTINE HWBTOP
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
       SUBROUTINE HWBTOP
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWBVMC,HWRGEN,HWUALF,HWUSQR,X(3),W,
      & X3MIN,X3MAX,X1MIN,X1MAX,QSCALE,GLUFAC,R(3,3),M(3),
      & E(3),AW,PTSQ,EM,EPS,MASDEP,A,B,C,GAMDEP,LAMBDA,
       DOUBLE PRECISION HWBVMC,HWRGEN,HWUALF,HWUSQR,X(3),W,
      & X3MIN,X3MAX,X1MIN,X1MAX,QSCALE,GLUFAC,R(3,3),M(3),
      & E(3),AW,PTSQ,EM,EPS,MASDEP,A,B,C,GAMDEP,LAMBDA,


@@ -5297,7 +5297,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     VIRTUAL MASS CUTOFF FOR PARTON TYPE ID
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     VIRTUAL MASS CUTOFF FOR PARTON TYPE ID
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWBVMC
       INTEGER ID
       IF (ID.EQ.13) THEN
       DOUBLE PRECISION HWBVMC
       INTEGER ID
       IF (ID.EQ.13) THEN


@@ -5319,7 +5319,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to split a baryonic cluster containing two heavy quarks
 C  Based on HWCCUT
 C-----------------------------------------------------------------------
 C  Subroutine to split a baryonic cluster containing two heavy quarks
 C  Based on HWCCUT
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWUPCM,HWRGEN,HWVDOT,EMC,QM1,QM2,QM3,QM4,
      &                 PXY,PCX,PCY,RCM,PCL(5),AX(5),PA(5),PB(5),PC(5),
      &                 VCLUS(4),DQM,EMX,EMY,SKAPPA,RKAPPA,VTMP(4),
       DOUBLE PRECISION HWUPCM,HWRGEN,HWVDOT,EMC,QM1,QM2,QM3,QM4,
      &                 PXY,PCX,PCY,RCM,PCL(5),AX(5),PA(5),PB(5),PC(5),
      &                 VCLUS(4),DQM,EMX,EMY,SKAPPA,RKAPPA,VTMP(4),


@@ -5471,7 +5471,7 @@ C-----------------------------------------------------------------------
 C FINDS UNPAIRED PARTONS AFTER BARYON-NUMBER VIOLATION
 C  MODIFIED FOR RPARITY VIOLATING SUSY
 C-----------------------------------------------------------------------
 C FINDS UNPAIRED PARTONS AFTER BARYON-NUMBER VIOLATION
 C  MODIFIED FOR RPARITY VIOLATING SUSY
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       COMMON/HWBVIC/NBV,IBV(18)
       DOUBLE PRECISION HWRGEN,PDQ(5)
       INTEGER NBV,IBV,JBV,KBV,LBV,IHEP,IP1,IP2,IP3,JP1,JP2,JP3,
       COMMON/HWBVIC/NBV,IBV(18)
       DOUBLE PRECISION HWRGEN,PDQ(5)
       INTEGER NBV,IBV,JBV,KBV,LBV,IHEP,IP1,IP2,IP3,JP1,JP2,JP3,


@@ -5704,7 +5704,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C  Function to find the baryon number violating vertex a parton came from
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C  Function to find the baryon number violating vertex a parton came from
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER HWCBVT,IP,JP(2),KP,I,J,ID,TYPE,IDM,IDM2,IDM3,IDM4
       JP(1) = IP
       ID = IDHW(IP)
       INTEGER HWCBVT,IP,JP(2),KP,I,J,ID,TYPE,IDM,IDM2,IDM3,IDM4
       JP(1) = IP
       ID = IDHW(IP)


@@ -5769,7 +5769,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C  Subroutine to correct colour connections after the gluon splitting
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C  Subroutine to correct colour connections after the gluon splitting
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IHEP,STFSPT,LHEP,MHEP,RHEP
       IF(IERROR.NE.0) RETURN
 C--Find the first particle in the event record with status 150
       INTEGER IHEP,STFSPT,LHEP,MHEP,RHEP
       IF(IERROR.NE.0) RETURN
 C--Find the first particle in the event record with status 150


@@ -5826,7 +5826,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Cuts into 2 the cluster, momentum PCL, made of partons JHEP & KHEP
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Cuts into 2 the cluster, momentum PCL, made of partons JHEP & KHEP
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWREXQ,HWUPCM,HWRGEN,HWVDOT,EMC,QM1,QM2,EMX,EMY,
      & QM3,PXY,PCX,PCY,RCM,PCL(5),AX(5),PA(5),PB(5),PC(5),SKAPPA,DELTM,
      & VSCA,VTMP(4),RKAPPA,VCLUS
       DOUBLE PRECISION HWREXQ,HWUPCM,HWRGEN,HWVDOT,EMC,QM1,QM2,EMX,EMY,
      & QM3,PXY,PCX,PCY,RCM,PCL(5),AX(5),PA(5),PB(5),PC(5),SKAPPA,DELTM,
      & VSCA,VTMP(4),RKAPPA,VCLUS


@@ -5978,7 +5978,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     DECAYS CLUSTERS INTO PRIMARY HADRONS
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     DECAYS CLUSTERS INTO PRIMARY HADRONS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER JCL,KCL,IP,JP,KP,IST,ID1,ID2,ID3
       IF (IERROR.NE.0) RETURN
       IF (IPRO/10.EQ.9.OR.IPRO/10.EQ.5) THEN
       INTEGER JCL,KCL,IP,JP,KP,IST,ID1,ID2,ID3
       IF (IERROR.NE.0) RETURN
       IF (IPRO/10.EQ.9.OR.IPRO/10.EQ.5) THEN


@@ -6053,7 +6053,7 @@ C     Modified by IGK to include BRW's colour rearrangement and
 C     MHS's cluster vertices
 C     MODIFIED 16/10/97 BY BRW FOR SUSY PROCESSES
 C-----------------------------------------------------------------------
 C     MHS's cluster vertices
 C     MODIFIED 16/10/97 BY BRW FOR SUSY PROCESSES
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWULDO,HWVDOT,HWRGEN,HWUPCM,DCL0,DCL(4),DCL1,
      & DFAC,DISP1(4),DISP2(4),DMAX,PCL(5),DOT1,DOT2,FAC,VCLUS,SCA1,SCA2,
      & EM0,EM1,EM2,PC0,PC1
       DOUBLE PRECISION HWULDO,HWVDOT,HWRGEN,HWUPCM,DCL0,DCL(4),DCL1,
      & DFAC,DISP1(4),DISP2(4),DMAX,PCL(5),DOT1,DOT2,FAC,VCLUS,SCA1,SCA2,
      & EM0,EM1,EM2,PC0,PC1


@@ -6394,7 +6394,7 @@ C-----------------------------------------------------------------------
 C     SPLITS ANY TIMELIKE GLUONS REMAINING AFTER PERTURBATIVE
 C     BRANCHING INTO LIGHT (I.E. U OR D) Q-QBAR PAIRS
 C-----------------------------------------------------------------------
 C     SPLITS ANY TIMELIKE GLUONS REMAINING AFTER PERTURBATIVE
 C     BRANCHING INTO LIGHT (I.E. U OR D) Q-QBAR PAIRS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,PF
       INTEGER HWRINT,IHEP,JHEP,KHEP,LHEP,MHEP,ID,J,IST
       EXTERNAL HWRGEN,HWRINT
       DOUBLE PRECISION HWRGEN,PF
       INTEGER HWRINT,IHEP,JHEP,KHEP,LHEP,MHEP,ID,J,IST
       EXTERNAL HWRGEN,HWRINT


@@ -6555,7 +6555,7 @@ C     (IN SPECIAL CLUSTER CODE - SEE HWCFLA)
 C
 C MODIFIED 15/11/99 TO SMEAR POSITIONS OF HADRONS BY 1/(CLUSTER MASS)
 C-----------------------------------------------------------------------
 C
 C MODIFIED 15/11/99 TO SMEAR POSITIONS OF HADRONS BY 1/(CLUSTER MASS)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRGAU,HWVDOT,EM0,EM1,EM2,EMADU,EMSQ,
      & PCMAX,PCM,PTEST,PCQK,PP(5),EMLOW,RMAT(3,3),CT,ST,CX,SX,HPSMR
       INTEGER HWRINT,JCL,ID1,ID2,ID3,ID,IR1,IR2,NTRY,IDMIN,IMAX,I,MHEP,
       DOUBLE PRECISION HWRGEN,HWRGAU,HWVDOT,EM0,EM1,EM2,EMADU,EMSQ,
      & PCMAX,PCM,PTEST,PCQK,PP(5),EMLOW,RMAT(3,3),CT,ST,CX,SX,HPSMR
       INTEGER HWRINT,JCL,ID1,ID2,ID3,ID,IR1,IR2,NTRY,IDMIN,IMAX,I,MHEP,


@@ -6761,7 +6761,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Computes the width and maximum weight for a two body mode
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Computes the width and maximum weight for a two body mode
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IMODE,I
       DOUBLE PRECISION A(2),M(3),PCM,E1,E2,HWUPCM,PHS,WGT,MWGT,PCM2,
      &     M2(3)
       INTEGER IMODE,I
       DOUBLE PRECISION A(2),M(3),PCM,E1,E2,HWUPCM,PHS,WGT,MWGT,PCM2,
      &     M2(3)


@@ -6887,7 +6887,7 @@ C-----------------------------------------------------------------------
 C     Subroutine to perform the three body decays for spin correlations
 C     and SUSY three body modes
 C-----------------------------------------------------------------------
 C     Subroutine to perform the three body decays for spin correlations
 C     and SUSY three body modes
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IMODE,I,J,ID,IDP(4+NDIAGR),ITYPE,NDIA,ID1,ID2,
      &     DRTYPE(NDIAGR),NTRY,IDSPIN,NCTHRE,DRCF(NDIAGR)
       DOUBLE PRECISION A,B,MS,MWD,MA,MB,MA2,MB2,WGT,M342,HWRUNI,
       INTEGER IMODE,I,J,ID,IDP(4+NDIAGR),ITYPE,NDIA,ID1,ID2,
      &     DRTYPE(NDIAGR),NTRY,IDSPIN,NCTHRE,DRCF(NDIAGR)
       DOUBLE PRECISION A,B,MS,MWD,MA,MB,MA2,MB2,WGT,M342,HWRUNI,


@@ -7115,7 +7115,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Subroutine to calculate the matrix element for a given mode
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Subroutine to calculate the matrix element for a given mode
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER I,J,P0,P1,P2,P3,P0P,IB,ID,IDP(4+NDIAGR),IDSPIN,NDIA,
      &     DRTYPE(NDIAGR),NCTHRE,DRCF(NDIAGR)
       DOUBLE PRECISION A,B,MS,MWD,MA,MB,MA2,MB2,WGT,FJAC,M342,HWRUNI,
       INTEGER I,J,P0,P1,P2,P3,P0P,IB,ID,IDP(4+NDIAGR),IDSPIN,NDIA,
      &     DRTYPE(NDIAGR),NCTHRE,DRCF(NDIAGR)
       DOUBLE PRECISION A,B,MS,MWD,MA,MB,MA2,MB2,WGT,FJAC,M342,HWRUNI,


@@ -7395,7 +7395,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for the three body
 C  gauge boson exchange diagram
 C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for the three body
 C  gauge boson exchange diagram
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,C(2,2),E(2,2),ZI,APP(2,2),APM(2,2),
      &     AMP(2,2),AMM(2,2),F1M(2,2,8),F3(2,2,8)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,C(2,2),E(2,2),ZI,APP(2,2),APM(2,2),
      &     AMP(2,2),AMM(2,2),F1M(2,2,8),F3(2,2,8)


@@ -7473,7 +7473,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for the three body
 C  Higgs boson exchange diagram
 C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for the three body
 C  Higgs boson exchange diagram
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
      &     F3(2,2,8)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
      &     F3(2,2,8)


@@ -7523,7 +7523,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for the three body
 C  antisfermion exchange diagram
 C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for the three body
 C  antisfermion exchange diagram
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,V1(2,2),V2(2,2),ZI,F1M(2,2,8),
      &     F3(2,2,8)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,V1(2,2),V2(2,2),ZI,F1M(2,2,8),
      &     F3(2,2,8)


@@ -7572,7 +7572,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for the three body
 C  sfermion exchange diagram
 C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for the three body
 C  sfermion exchange diagram
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
      &     F3(2,2,8)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
      &     F3(2,2,8)


@@ -7621,7 +7621,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for the three body
 C  gauge boson exchange diagram (antiparticle decay)
 C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for the three body
 C  gauge boson exchange diagram (antiparticle decay)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,C(2,2),E(2,2),ZI,APP(2,2),APM(2,2),
      &     AMP(2,2),AMM(2,2),F1M(2,2,8),F3(2,2,8),F01(2,2,8,8)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,C(2,2),E(2,2),ZI,APP(2,2),APM(2,2),
      &     AMP(2,2),AMM(2,2),F1M(2,2,8),F3(2,2,8),F01(2,2,8,8)


@@ -7699,7 +7699,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for the three body
 C  gauge boson exchange diagram
 C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for the three body
 C  gauge boson exchange diagram
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,C(2,2),ZI,APP(2,2),APM(2,2),
      &     AMP(2,2),AMM(2,2),F1M(2,2,8),F3(2,2,8),F01(2,2,8,8)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,C(2,2),ZI,APP(2,2),APM(2,2),
      &     AMP(2,2),AMM(2,2),F1M(2,2,8),F3(2,2,8),F01(2,2,8,8)


@@ -7771,7 +7771,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for the three body
 C  decay fermion --> gravitino fermion antifermion (via gauge boson)
 C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for the three body
 C  decay fermion --> gravitino fermion antifermion (via gauge boson)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,F1M(2,2,8),F3(2,2,8)
       DOUBLE PRECISION A,B,MS,MWD,MA,MB,MA2,MB2,M342,M232,M242,MR,
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,F1M(2,2,8),F3(2,2,8)
       DOUBLE PRECISION A,B,MS,MWD,MA,MB,MA2,MB2,M342,M232,M242,MR,


@@ -7814,7 +7814,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for 1st 3 body RPV
 C  diagram f--> fbar fbar f
 C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for 1st 3 body RPV
 C  diagram f--> fbar fbar f
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
      &     F3(2,2,8)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
      &     F3(2,2,8)


@@ -7863,7 +7863,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for 2nd 3 body RPV
 C  diagram f --> fbar fbar f
 C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for 2nd 3 body RPV
 C  diagram f --> fbar fbar f
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,V1(2,2),V2(2,2),ZI,F1M(2,2,8),
      &     F3(2,2,8)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,V1(2,2),V2(2,2),ZI,F1M(2,2,8),
      &     F3(2,2,8)


@@ -7912,7 +7912,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for 3rd 3 body RPV
 C  diagram f --> fbar fbar f
 C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for 3rd 3 body RPV
 C  diagram f --> fbar fbar f
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
      &     F3(2,2,8)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
      &     F3(2,2,8)


@@ -7961,7 +7961,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for 4th 3 body RPV
 C  diagram f --> f f f
 C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for 4th 3 body RPV
 C  diagram f --> f f f
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
      &     F3(2,2,8)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
      &     F3(2,2,8)


@@ -8010,7 +8010,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for 5th 3 body RPV
 C  diagram f --> f f f
 C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for 5th 3 body RPV
 C  diagram f --> f f f
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,V1(2,2),V2(2,2),ZI,F1M(2,2,8),
      &     F3(2,2,8)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,V1(2,2),V2(2,2),ZI,F1M(2,2,8),
      &     F3(2,2,8)


@@ -8059,7 +8059,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for 6th 3 body RPV
 C  diagram f --> f f f
 C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for 6th 3 body RPV
 C  diagram f --> f f f
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
      &     F3(2,2,8)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
      &     F3(2,2,8)


@@ -8108,7 +8108,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for 7th 3 body RPV
 C  diagram f --> fbar fbar fbar
 C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for 7th 3 body RPV
 C  diagram f --> fbar fbar fbar
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
      &     F3(2,2,8)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
      &     F3(2,2,8)


@@ -8157,7 +8157,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for 8th 3 body RPV
 C  diagram f --> fbar fbar fbar
 C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for 8th 3 body RPV
 C  diagram f --> fbar fbar fbar
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,V1(2,2),V2(2,2),ZI,F1M(2,2,8),
      &     F3(2,2,8)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,V1(2,2),V2(2,2),ZI,F1M(2,2,8),
      &     F3(2,2,8)


@@ -8206,7 +8206,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for 9th 3 body RPV
 C  diagram f --> fbar fbar fbar
 C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for 9th 3 body RPV
 C  diagram f --> fbar fbar fbar
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
      &     F3(2,2,8)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
      &     F3(2,2,8)


@@ -8255,7 +8255,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for the three body
 C  Higgs boson exchange diagram antifermion decay
 C-----------------------------------------------------------------------
 C  Subroutine to calculate the helicity amplitudes for the three body
 C  Higgs boson exchange diagram antifermion decay
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
      &     F3(2,2,8)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
      &     F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
      &     F3(2,2,8)


@@ -8304,7 +8304,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Subroutine to perform the four body Higgs decays
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Subroutine to perform the four body Higgs decays
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IMODE,I,J,ID,IDP(4+NDIAGR),ITYPE(2),NTRY,ITYPE1,ITYPE2
       DOUBLE PRECISION A,B,MS,MWD,M,M2,WGT,HWRUNI,BRW(6),BRZ(12),
      &     HWUPCM,WMAX,WSUM,WSSUM,MR,PRE,TEMP,HWRGEN,WTMAX,P(5,5)
       INTEGER IMODE,I,J,ID,IDP(4+NDIAGR),ITYPE(2),NTRY,ITYPE1,ITYPE2
       DOUBLE PRECISION A,B,MS,MWD,M,M2,WGT,HWRUNI,BRW(6),BRZ(12),
      &     HWUPCM,WMAX,WSUM,WSSUM,MR,PRE,TEMP,HWRGEN,WTMAX,P(5,5)


@@ -8414,7 +8414,7 @@ C-----------------------------------------------------------------------
 C     Subroutine to calculate the matrix element for a given four body
 C     decay mode
 C-----------------------------------------------------------------------
 C     Subroutine to calculate the matrix element for a given four body
 C     decay mode
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER I,J,P0,P1,P2,P3,ID,O(2),IDP(4+NDIAGR),II,P4
       DOUBLE PRECISION A,B,MS,MWD,M,M2,WGT,HWRUNI,
      &     M23,PCMA,PCMB(2),HWUPCM,PHS,N(3),HWVDOT,PP,HWULDO,EPS,
       INTEGER I,J,P0,P1,P2,P3,ID,O(2),IDP(4+NDIAGR),II,P4
       DOUBLE PRECISION A,B,MS,MWD,M,M2,WGT,HWRUNI,
      &     M23,PCMA,PCMB(2),HWUPCM,PHS,N(3),HWVDOT,PP,HWULDO,EPS,


@@ -8644,7 +8644,7 @@ C     IF BOSON CAME FROM HIGGS DECAY, GIVE BOTH THE SAME HELICITY (EPR)
 C     IF BOSON CAME FROM W+1JET, GIVE IT THE CORRECT DECAY CORRELATIONS
 C--BRW FIX 20/07/04: ADD FULL DECAY CORRELATIONS FOR W/Z+HIGGS
 C-----------------------------------------------------------------------
 C     IF BOSON CAME FROM W+1JET, GIVE IT THE CORRECT DECAY CORRELATIONS
 C--BRW FIX 20/07/04: ADD FULL DECAY CORRELATIONS FOR W/Z+HIGGS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,HWUPCM,HWULDO,R(3,3),CV,CA,BR,PCM,
      & PBOS(5),PMAX,PROB,RRLL,RLLR
       INTEGER HWRINT,IBOS,IBOSON,IPAIR,ICMF,IOPT,IHEL,IMOTH,
       DOUBLE PRECISION HWRGEN,HWRUNI,HWUPCM,HWULDO,R(3,3),CV,CA,BR,PCM,
      & PBOS(5),PMAX,PROB,RRLL,RLLR
       INTEGER HWRINT,IBOS,IBOSON,IPAIR,ICMF,IOPT,IHEL,IMOTH,


@@ -8849,7 +8849,7 @@ C-----------------------------------------------------------------------
 C     CHOOSE DECAY MODE OF BOSON
 C     IOPT=2 TO RESET COUNTERS, 1 FOR BOSON PAIR, 0 FOR ANY OTHERS
 C-----------------------------------------------------------------------
 C     CHOOSE DECAY MODE OF BOSON
 C     IOPT=2 TO RESET COUNTERS, 1 FOR BOSON PAIR, 0 FOR ANY OTHERS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,BRMODE(12,3),CV,CA,BR,BRLST,BRCOM,FACZ,
      & FACW
       INTEGER HWRINT,IDBOS,IDEC,IDMODE(2,12,3),IFER,IANT,IOPT,I1,I2,
       DOUBLE PRECISION HWRGEN,BRMODE(12,3),CV,CA,BR,BRLST,BRCOM,FACZ,
      & FACW
       INTEGER HWRINT,IDBOS,IDEC,IDMODE(2,12,3),IFER,IANT,IOPT,I1,I2,


@@ -9012,7 +9012,7 @@ C     IOPT=2 TO RESET COUNTERS, 1 FOR BOSON PAIR, 0 FOR ANY OTHERS
 C     IDENTICAL TO HWDBOZ BUT REQUIRES DECAY MODE ACCESSIBLE FOR GIVEN
 C     MASS
 C-----------------------------------------------------------------------
 C     IDENTICAL TO HWDBOZ BUT REQUIRES DECAY MODE ACCESSIBLE FOR GIVEN
 C     MASS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,BRMODE(12,3),CV,CA,BR,BRLST,BRCOM,FACZ,
      & FACW,MSMODE(12,3),MASS
       INTEGER HWRINT,IDBOS,IDEC,IDMODE(2,12,3),IFER,IANT,IOPT,I1,I2,
       DOUBLE PRECISION HWRGEN,BRMODE(12,3),CV,CA,BR,BRLST,BRCOM,FACZ,
      & FACW,MSMODE(12,3),MASS
       INTEGER HWRINT,IDBOS,IDEC,IDMODE(2,12,3),IFER,IANT,IOPT,I1,I2,


@@ -9189,7 +9189,7 @@ C-----------------------------------------------------------------------
 C     Checks line L of decay table is compatible with decay of particle
 C     IDKY, tidies up the line and sets NPRODS.
 C-----------------------------------------------------------------------
 C     Checks line L of decay table is compatible with decay of particle
 C     IDKY, tidies up the line and sets NPRODS.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION EPS,QS,Q,DM
       INTEGER IDKY,L,IFAULT,I,ID,J
       LOGICAL IFGO
       DOUBLE PRECISION EPS,QS,Q,DM
       INTEGER IDKY,L,IFAULT,I,ID,J
       LOGICAL IFGO


@@ -9267,7 +9267,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C INTERFACE TO QQ-CLEO MONTE CARLO (LS 11/12/91)
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C INTERFACE TO QQ-CLEO MONTE CARLO (LS 11/12/91)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IHEP,IIHEP,NHEPHF,QQLMAT
       LOGICAL QQLERR
       CHARACTER*8 NAME
       INTEGER IHEP,IIHEP,NHEPHF,QQLMAT
       LOGICAL QQLERR
       CHARACTER*8 NAME


@@ -9435,7 +9435,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C INTERFACE TO EURODEC PACKAGE (LS 10/29/91)
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C INTERFACE TO EURODEC PACKAGE (LS 10/29/91)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IHEP,IIHEP,NHEPHF,IEUPDG,IPDGEU
       CHARACTER*8 NAME
 C---EURODEC COMMON'S : INITIAL INPUT
       INTEGER IHEP,IIHEP,NHEPHF,IEUPDG,IPDGEU
       CHARACTER*8 NAME
 C---EURODEC COMMON'S : INITIAL INPUT


@@ -9538,7 +9538,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Generates 4-body decay 0->1+2+3+4 using pure phase space
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Generates 4-body decay 0->1+2+3+4 using pure phase space
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,P0(5),P1(5),P2(5),P3(5),P4(5),B,C,AA,BB,
      & CC,DD,EE,TT,S1,RS1,FF,S2,PP,QQ,RR,P1CM,P234(5),P2CM,P34(5),P3CM
       INTEGER NTRY
       DOUBLE PRECISION HWRGEN,P0(5),P1(5),P2(5),P3(5),P4(5),B,C,AA,BB,
      & CC,DD,EE,TT,S1,RS1,FF,S2,PP,QQ,RR,P1CM,P234(5),P2CM,P34(5),P3CM
       INTEGER NTRY


@@ -9589,7 +9589,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Generates 5-body decay 0->1+2+3+4+5 using pure phase space
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Generates 5-body decay 0->1+2+3+4+5 using pure phase space
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,P0(5),P1(5),P2(5),P3(5),P4(5),P5(5),B,C,
      & AA,BB,CC,DD,EE,FF,TT,S1,RS1,GG,S2,RS2,HH,S3,PP,QQ,RR,SS,P1CM,
      & P2345(5),P2CM,P345(5),P3CM,P45(5),P4CM
       DOUBLE PRECISION HWRGEN,P0(5),P1(5),P2(5),P3(5),P4(5),P5(5),B,C,
      & AA,BB,CC,DD,EE,FF,TT,S1,RS1,GG,S2,RS2,HH,S3,PP,QQ,RR,SS,P1CM,
      & P2345(5),P2CM,P345(5),P3CM,P45(5),P4CM


@@ -9651,7 +9651,7 @@ C-----------------------------------------------------------------------
 C     GENERATES DECAYS OF UNSTABLE HADRONS AND LEPTONS
 C     Modified for TAUOLA interface 16/10/01 PR
 C-----------------------------------------------------------------------
 C     GENERATES DECAYS OF UNSTABLE HADRONS AND LEPTONS
 C     Modified for TAUOLA interface 16/10/01 PR
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       COMMON/FFS/TB,BT
       COMMON/SFF/IT1,IB1,IT2,IB2
       DOUBLE PRECISION TB,BT
       COMMON/FFS/TB,BT
       COMMON/SFF/IT1,IB1,IT2,IB2
       DOUBLE PRECISION TB,BT


@@ -9928,7 +9928,7 @@ C-----------------------------------------------------------------------
 C  CALCULATE THE COMPLEX FUNCTION F OF HHG eq 2.18
 C  FOR USE IN H-->GAMMGAMM DECAYS
 C-----------------------------------------------------------------------
 C  CALCULATE THE COMPLEX FUNCTION F OF HHG eq 2.18
 C  FOR USE IN H-->GAMMGAMM DECAYS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION TAU,FNREAL,FNIMAG,FNLOG,FNSQR
       IF (TAU.GT.ONE) THEN
         FNREAL=(ASIN(1/SQRT(TAU)))**2
       DOUBLE PRECISION TAU,FNREAL,FNIMAG,FNLOG,FNSQR
       IF (TAU.GT.ONE) THEN
         FNREAL=(ASIN(1/SQRT(TAU)))**2


@@ -9952,7 +9952,7 @@ C-----------------------------------------------------------------------
 C  CALCULATE THE DOUBLE BREIT-WIGNER INTEGRAL
 C  X=(EMV/EMH)**2 , Y=EMV*GAMV/EMH**2
 C-----------------------------------------------------------------------
 C  CALCULATE THE DOUBLE BREIT-WIGNER INTEGRAL
 C  X=(EMV/EMH)**2 , Y=EMV*GAMV/EMH**2
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWDHGF,X,Y,CHANGE,X1,X2,FAC1,FAC2,TH1,TH2,TH1HI,
      & TH1LO,TH2HI,TH2LO,X2MAX,SQFAC
       INTEGER NBIN,IBIN1,IBIN2
       DOUBLE PRECISION HWDHGF,X,Y,CHANGE,X1,X2,FAC1,FAC2,TH1,TH2,TH1HI,
      & TH1LO,TH2HI,TH2LO,X2MAX,SQFAC
       INTEGER NBIN,IBIN1,IBIN2


@@ -10017,7 +10017,7 @@ C     A) FOR GAMinp=0 FIND AND DECAY HIGGS
 C     B) FOR GAMinp>0 CALCULATE TOTAL HIGGS WIDTH
 C                     FOR EMH=GAMINP. STORE RESULT IN GAMINP.
 C-----------------------------------------------------------------------
 C     B) FOR GAMinp>0 CALCULATE TOTAL HIGGS WIDTH
 C                     FOR EMH=GAMINP. STORE RESULT IN GAMINP.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWDHGF,HWRGEN,HWRUNI,HWUSQR,HWUPCM,GAMINP,EMH,
      & EMF,COLFAC,ENF,K1,K0,BET0,BET1,GAM0,GAM1,SCLOG,CFAC,XF,EM,GAMLIM,
      & GAM,XW,EMW,XZ,EMZ,YW,YZ,EMI,TAUT,TAUW,WIDHIG,VECDEC,EMB,GAMB,
       DOUBLE PRECISION HWDHGF,HWRGEN,HWRUNI,HWUSQR,HWUPCM,GAMINP,EMH,
      & EMF,COLFAC,ENF,K1,K0,BET0,BET1,GAM0,GAM1,SCLOG,CFAC,XF,EM,GAMLIM,
      & GAM,XW,EMW,XZ,EMZ,YW,YZ,EMI,TAUT,TAUW,WIDHIG,VECDEC,EMB,GAMB,


@@ -10260,7 +10260,7 @@ C   MODIFIED TO INCLUDE R-PARITY VIOLATING SUSY PR 9/4/99
 C   MODIFIED TO CALL A NUMBER OF ROUTINES TO DO THE VARIOUS BITS OF
 C   THE PROCESS
 C-----------------------------------------------------------------------
 C   MODIFIED TO CALL A NUMBER OF ROUTINES TO DO THE VARIOUS BITS OF
 C   THE PROCESS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION PW(5)
       INTEGER IHEP,IS,ID,IM,KHEP,LHEP,MHEP,NPR,CLSAVE(2),NHEPST
       LOGICAL FOUND
       DOUBLE PRECISION PW(5)
       INTEGER IHEP,IS,ID,IM,KHEP,LHEP,MHEP,NPR,CLSAVE(2),NHEPST
       LOGICAL FOUND


@@ -10315,7 +10315,7 @@ C   Subroutine to perform the first part of the heavy object decays
 C   IE to select the decay mode
 C   was part of HWDHOB
 C-----------------------------------------------------------------------
 C   IE to select the decay mode
 C   was part of HWDHOB
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWUMBW,HWRGEN,SDKM,RN,BF
       INTEGER IST(3),IHEP,ID,IM,I,JHEP,LHEP,MHEP,NPR,MTRY,NTRY,IS
       EXTERNAL HWRGEN
       DOUBLE PRECISION HWUMBW,HWRGEN,SDKM,RN,BF
       INTEGER IST(3),IHEP,ID,IM,I,JHEP,LHEP,MHEP,NPR,MTRY,NTRY,IS
       EXTERNAL HWRGEN


@@ -10447,7 +10447,7 @@ C   Subroutine to perform the second part of the heavy object decays
 C   IE generate the kinematics for the decay
 C   was part of HWDHOB
 C-----------------------------------------------------------------------
 C   IE generate the kinematics for the decay
 C   was part of HWDHOB
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       COMMON/FFS/TB,BT
       COMMON/SFF/IT1,IB1,IT2,IB2
       DOUBLE PRECISION TB,BT
       COMMON/FFS/TB,BT
       COMMON/SFF/IT1,IB1,IT2,IB2
       DOUBLE PRECISION TB,BT


@@ -10598,7 +10598,7 @@ C   Subroutine to perform the third part of the heavy object decays
 C   IE setup the colour connections
 C   was part of HWDHOB
 C-----------------------------------------------------------------------
 C   IE setup the colour connections
 C   was part of HWDHOB
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER ID,IM,KHEP,LHEP,MHEP,NPR,CLSAVE(2)
       IF (IERROR.NE.0) RETURN
 C Colour connections
       INTEGER ID,IM,KHEP,LHEP,MHEP,NPR,CLSAVE(2)
       IF (IERROR.NE.0) RETURN
 C Colour connections


@@ -10793,7 +10793,7 @@ C   Subroutine to perform the fourth part of the heavy object decays
 C   IE parton-showers with special treatment for top
 C   was part of HWDHOB
 C-----------------------------------------------------------------------
 C   IE parton-showers with special treatment for top
 C   was part of HWDHOB
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION PW(5),PDW(5,3)
       INTEGER IHEP,ID,IM,I,KHEP,LHEP,MHEP,NPR,NTRY,WHEP,SHEP
       DOUBLE COMPLEX RHOIN(2,2)
       DOUBLE PRECISION PW(5),PDW(5,3)
       INTEGER IHEP,ID,IM,I,KHEP,LHEP,MHEP,NPR,NTRY,WHEP,SHEP
       DOUBLE COMPLEX RHOIN(2,2)


@@ -10903,7 +10903,7 @@ C   Subroutine to perform the fifth part of the heavy object decays
 C   IE sort out RPV colour connections
 C   was part of HWDHOB
 C-----------------------------------------------------------------------
 C   IE sort out RPV colour connections
 C   was part of HWDHOB
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER ID,LHEP,MHEP,IDM,IDM2,THEP,CLSAVE(2)
       IF (IERROR.NE.0) RETURN
 C--New to correct colour connections in Rslash
       INTEGER ID,LHEP,MHEP,IDM,IDM2,THEP,CLSAVE(2)
       IF (IERROR.NE.0) RETURN
 C--New to correct colour connections in Rslash


@@ -10968,7 +10968,7 @@ C-----------------------------------------------------------------------
 C   Subroutine to perform the final part of the heavy object decays
 C   IE sort out any colour connection problems
 C-----------------------------------------------------------------------
 C   Subroutine to perform the final part of the heavy object decays
 C   IE sort out any colour connection problems
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IHEP,IM,JHEP,ISM,JCM
       IF (IERROR.NE.0) RETURN
 C Fix any SUSY colour disconnections
       INTEGER IHEP,IM,JHEP,ISM,JCM
       IF (IERROR.NE.0) RETURN
 C Fix any SUSY colour disconnections


@@ -11008,7 +11008,7 @@ C     Performs partonic decays of hadrons containing heavy quark(s):
 C     either, meson/baryon spectator model weak decays;
 C     or, quarkonia -> 2-gluons, q-qbar, 3-gluons, or 2-gluons + photon.
 C-----------------------------------------------------------------------
 C     either, meson/baryon spectator model weak decays;
 C     or, quarkonia -> 2-gluons, q-qbar, 3-gluons, or 2-gluons + photon.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       COMMON/FFS/TB,BT
       COMMON/SFF/IT1,IB1,IT2,IB2
       DOUBLE PRECISION TB,BT
       COMMON/FFS/TB,BT
       COMMON/SFF/IT1,IB1,IT2,IB2
       DOUBLE PRECISION TB,BT


@@ -11250,7 +11250,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Sets the colour connections in Baryon number violating decays
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Sets the colour connections in Baryon number violating decays
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IHEP,MHEP,ID,ID2,IDM2,IDM3,COLCON(2,2,3),FLACON(2,3),JHEP,
      &        DECAY,COLANT,KHEP,IDM,IDMB,IDMB2,IDMB3,IDMB4,QHEP,IDM4,
      &        CLSAVE(2),XHEP,I,HWRINT,THEP
       INTEGER IHEP,MHEP,ID,ID2,IDM2,IDM3,COLCON(2,2,3),FLACON(2,3),JHEP,
      &        DECAY,COLANT,KHEP,IDM,IDMB,IDMB2,IDMB3,IDMB4,QHEP,IDM4,
      &        CLSAVE(2),XHEP,I,HWRINT,THEP


@@ -11409,7 +11409,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SUBROUTINE TO IMPLEMENT ALL RPARITY DECAY MATRIX ELEMENTS
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SUBROUTINE TO IMPLEMENT ALL RPARITY DECAY MATRIX ELEMENTS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION SM(6),SW(6),HWULDO,INFCOL,AM, M12SQ,M23SQ,MSGN,
      &                 M13SQ,A(6),B(6),SWEAK,MW,DECMOM(5),TEST(3),EPS,
      &                 M12SQT(6),M23SQT(6),M13SQT(6),LIMIT,M(4),RAND,
       DOUBLE PRECISION SM(6),SW(6),HWULDO,INFCOL,AM, M12SQ,M23SQ,MSGN,
      &                 M13SQ,A(6),B(6),SWEAK,MW,DECMOM(5),TEST(3),EPS,
      &                 M12SQT(6),M23SQT(6),M13SQT(6),LIMIT,M(4),RAND,


@@ -12141,7 +12141,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C  Subroutine to perform decays including spin correlations
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C  Subroutine to perform decays including spin correlations
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION PW(5)
       INTEGER IDEC,IP,IS,IHEP,ID,IM,LHEP,MHEP,NPR,KHEP,CLSAVE(2),NTRY,
      &     ID1
       DOUBLE PRECISION PW(5)
       INTEGER IDEC,IP,IS,IHEP,ID,IM,LHEP,MHEP,NPR,KHEP,CLSAVE(2),NTRY,
      &     ID1


@@ -12234,7 +12234,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C  Subroutine to check a vertex and decide which branch to treat
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C  Subroutine to check a vertex and decide which branch to treat
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IDEC,I,IPICK(5),IP,HWRINT,P1,P2,P3,P4,P3P,P4P,NPR,P0,P0P,
      &     P1P,P2P,IF1,IF2,P5,P5P
       DOUBLE PRECISION NORM
       INTEGER IDEC,I,IPICK(5),IP,HWRINT,P1,P2,P3,P4,P3P,P4P,NPR,P0,P0P,
      &     P1P,P2P,IF1,IF2,P5,P5P
       DOUBLE PRECISION NORM


@@ -12684,7 +12684,7 @@ C   IE generate the kinematics for the decay
 C   including spin correlations
 C   was part of HWDHOB
 C-----------------------------------------------------------------------
 C   including spin correlations
 C   was part of HWDHOB
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,PW(5),HWDPWT,HWDWWT,PCM,HWUPCM
       INTEGER IHEP,IM,KHEP,MHEP,NPR,ISN,RHEP
       EXTERNAL HWRGEN,HWDPWT,HWDWWT,HWUPCM
       DOUBLE PRECISION HWRGEN,PW(5),HWDPWT,HWDWWT,PCM,HWUPCM
       INTEGER IHEP,IM,KHEP,MHEP,NPR,ISN,RHEP
       EXTERNAL HWRGEN,HWDPWT,HWDWWT,HWUPCM


@@ -12804,7 +12804,7 @@ C     averages spin if not using TAUOLA
 C     if using TAUOLA selects the spin and uses TAUOLA to perform the
 C     decay
 C-----------------------------------------------------------------------
 C     if using TAUOLA selects the spin and uses TAUOLA to perform the
 C     decay
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IP,IHEP,ID1,ID,NTRY
       DOUBLE PRECISION PPOL,HWRGEN,POL
       EXTERNAL HWRGEN
       INTEGER IP,IHEP,ID1,ID,NTRY
       DOUBLE PRECISION PPOL,HWRGEN,POL
       EXTERNAL HWRGEN


@@ -12870,7 +12870,7 @@ C-----------------------------------------------------------------------
 C  Subroutine to calculate the two body matrix element for spin
 C  correlations
 C-----------------------------------------------------------------------
 C  Subroutine to calculate the two body matrix element for spin
 C  correlations
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IOUT1,IOUT2,IMODE,IDSPIN,ID,I,J,IDP(3),P0,P1,P2,O(2),P0P,
      &     NTRY
       DOUBLE PRECISION XMASS,PLAB,PRW,PCM,PREF(5),P(5,3),PM(5,3),PCMA,
       INTEGER IOUT1,IOUT2,IMODE,IDSPIN,ID,I,J,IDP(3),P0,P1,P2,O(2),P0P,
      &     NTRY
       DOUBLE PRECISION XMASS,PLAB,PRW,PCM,PREF(5),P(5,3),PM(5,3),PCMA,


@@ -13195,7 +13195,7 @@ C-----------------------------------------------------------------------
 C     Master subroutine for three body SUSY and spin ME's
 C     Uses HWD3ME to generate the momenta etc
 C-----------------------------------------------------------------------
 C     Master subroutine for three body SUSY and spin ME's
 C     Uses HWD3ME to generate the momenta etc
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX F0(2,2,8),F1(2,2,8),F1M(2,2,8),F3(2,2,8),
      &     F0M(2,2,8),F2(2,2,8),RHOIN(2,2),F01(2,2,8,8)
       DOUBLE PRECISION A,B,MS,MWD,MA,MB,MA2,MB2,M342,M232,M242,MR,
       DOUBLE COMPLEX F0(2,2,8),F1(2,2,8),F1M(2,2,8),F3(2,2,8),
      &     F0M(2,2,8),F2(2,2,8),RHOIN(2,2),F01(2,2,8,8)
       DOUBLE PRECISION A,B,MS,MWD,MA,MB,MA2,MB2,M342,M232,M242,MR,


@@ -13475,7 +13475,7 @@ C     Subroutine to perform the four body decays
 C     IOPT = 1 select decay mode and generate momenta
 C     IOPT = 2 enter first decays and perform parton shower
 C-----------------------------------------------------------------------
 C     IOPT = 1 select decay mode and generate momenta
 C     IOPT = 2 enter first decays and perform parton shower
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IOPT,ID,IOUT1,IOUT2,IB(2),I,IDF(4),ITYPE(2),IMODE,
      &     IDP(4+NDIAGR),ID1,ID2,J
       DOUBLE PRECISION CV,CA,A,B,MS,MWD,MR,M,M2,P(5,5),PW(5,2),BR
       INTEGER IOPT,ID,IOUT1,IOUT2,IB(2),I,IDF(4),ITYPE(2),IMODE,
      &     IDP(4+NDIAGR),ID1,ID2,J
       DOUBLE PRECISION CV,CA,A,B,MS,MWD,MR,M,M2,P(5,5),PW(5,2),BR


@@ -13586,7 +13586,7 @@ C     IOPT =-1 initialises
 C     IOPT = 1 performs decay
 C     IOPT = 2 write outs final TAUOLA information
 C-----------------------------------------------------------------------
 C     IOPT = 1 performs decay
 C     IOPT = 2 write outs final TAUOLA information
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IOPT,IHEP,ID,ITAU,I,IMO,NHEPPO
       DOUBLE PRECISION POL
       REAL POL1(4)
       INTEGER IOPT,IHEP,ID,ITAU,I,IMO,NHEPPO
       DOUBLE PRECISION POL
       REAL POL1(4)


@@ -13729,7 +13729,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     DECIDES WHETHER TO DO TOP QUARK DECAY BEFORE HADRONIZATION
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     DECIDES WHETHER TO DO TOP QUARK DECAY BEFORE HADRONIZATION
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       LOGICAL DECAY
       DECAY=RMASS(6).GT.130D0
       END
       LOGICAL DECAY
       DECAY=RMASS(6).GT.130D0
       END


@@ -13792,7 +13792,7 @@ C-----------------------------------------------------------------------
 C     MATRIX ELEMENT SQUARED FOR
 C     ((V-A)*TB1+(V+A)*CT1)*((V-A)*TB2+(V+A)*CT2)) WEAK DECAY
 C-----------------------------------------------------------------------
 C     MATRIX ELEMENT SQUARED FOR
 C     ((V-A)*TB1+(V+A)*CT1)*((V-A)*TB2+(V+A)*CT2)) WEAK DECAY
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       COMMON/FFS/TB,BT
       COMMON/SFF/IT1,IB1,IT2,IB2
       DOUBLE PRECISION TB,BT
       COMMON/FFS/TB,BT
       COMMON/SFF/IT1,IB1,IT2,IB2
       DOUBLE PRECISION TB,BT


@@ -13820,7 +13820,7 @@ C-----------------------------------------------------------------------
 C     Sets STAB=.TRUE. if DKVRTX lies outside the specified region.
 C  Revised 05/09/00 by BRW to put parameters in common
 C-----------------------------------------------------------------------
 C     Sets STAB=.TRUE. if DKVRTX lies outside the specified region.
 C  Revised 05/09/00 by BRW to put parameters in common
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION DKVRTX(4),RR
       LOGICAL STAB
       STAB=.FALSE.
       DOUBLE PRECISION DKVRTX(4),RR
       LOGICAL STAB
       STAB=.FALSE.


@@ -13863,7 +13863,7 @@ C     TERMINAL CALCULATIONS ON ELEMENTARY PROCESS
 C     Modified 28/03/01 by BRW to handle negative weights
 C     Modified 15/07/02 by PR for Les Houches Accord
 C-----------------------------------------------------------------------
 C     Modified 28/03/01 by BRW to handle negative weights
 C     Modified 15/07/02 by PR for Les Houches Accord
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER I
       DOUBLE PRECISION RNWGT,SPWGT,ERWGT
 C--Les Houches Common Block
       INTEGER I
       DOUBLE PRECISION RNWGT,SPWGT,ERWGT
 C--Les Houches Common Block


@@ -13986,7 +13986,7 @@ C-----------------------------------------------------------------------
 C     GENERATES A PHOTON IN WEIZSACKER-WILLIAMS (WWA=.TRUE.) OR
 C     ELSE EQUIVALENT PHOTON APPROX FROM INCOMING E+, E-, MU+ OR MU-
 C-----------------------------------------------------------------------
 C     GENERATES A PHOTON IN WEIZSACKER-WILLIAMS (WWA=.TRUE.) OR
 C     ELSE EQUIVALENT PHOTON APPROX FROM INCOMING E+, E-, MU+ OR MU-
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,EGMIN,ZMIN,ZMAX,ZGAM,SS,ZMI,ZMA,
      & PPL,PMI,QT2,Q2,QQMIN,QQMAX,S0,A
       INTEGER IHEP,IHADIS
       DOUBLE PRECISION HWRGEN,HWRUNI,EGMIN,ZMIN,ZMAX,ZGAM,SS,ZMI,ZMA,
      & PPL,PMI,QT2,Q2,QQMIN,QQMAX,S0,A
       INTEGER IHEP,IHADIS


@@ -14109,7 +14109,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     FIND MINIMUM INVARIANT MASS SQUARED NEEDED FOR HARD PROCESS, S0
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     FIND MINIMUM INVARIANT MASS SQUARED NEEDED FOR HARD PROCESS, S0
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION S0,RPM(2)
       INTEGER HQ,I
       IF (IPRO.EQ.13.OR.IPRO.EQ.14) THEN
       DOUBLE PRECISION S0,RPM(2)
       INTEGER HQ,I
       IF (IPRO.EQ.13.OR.IPRO.EQ.14) THEN


@@ -14263,7 +14263,7 @@ C-----------------------------------------------------------------------
 C     INITIALISES ELEMENTARY PROCESS
 C     Modified 28/03/01 by BRW to handle negative weights
 C-----------------------------------------------------------------------
 C     INITIALISES ELEMENTARY PROCESS
 C     Modified 28/03/01 by BRW to handle negative weights
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRSET,DUMMY,SAFETY
       EXTERNAL HWRSET
       PARAMETER (SAFETY=1.001)
       DOUBLE PRECISION HWRSET,DUMMY,SAFETY
       EXTERNAL HWRSET
       PARAMETER (SAFETY=1.001)


@@ -14328,7 +14328,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     GENERATES AN ISR PHOTON FROM INCOMING E+, E-, MU+ OR MU-
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     GENERATES AN ISR PHOTON FROM INCOMING E+, E-, MU+ OR MU-
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION CIRCKP(2)
       COMMON /HWCIR2/CIRCKP
       DOUBLE PRECISION HWRGEN,QSQMAX,QSQMIN,A,B,B1,B2,B3,B4,B5,B6,B7,B8,
       DOUBLE PRECISION CIRCKP(2)
       COMMON /HWCIR2/CIRCKP
       DOUBLE PRECISION HWRGEN,QSQMAX,QSQMIN,A,B,B1,B2,B3,B4,B5,B6,B7,B8,


@@ -14562,7 +14562,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SETS UP 2->1 (COLOUR SINGLET) HARD SUBPROCESS
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SETS UP 2->1 (COLOUR SINGLET) HARD SUBPROCESS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION PA
       INTEGER ICMF,I,IBM,IHEP
 C---INCOMING LINES
       DOUBLE PRECISION PA
       INTEGER ICMF,I,IBM,IHEP
 C---INCOMING LINES


@@ -14617,7 +14617,7 @@ C     WHEN NEVHEP=0, CHOOSES X VALUES AND FINDS WEIGHT FOR PROCESS IPROC
 C     OTHERWISE, CHOOSES AND LOADS ALL VARIABLES FOR HARD PROCESS
 C     modifications for Les Houches accord by PR (7/15/02)
 C-----------------------------------------------------------------------
 C     OTHERWISE, CHOOSES AND LOADS ALL VARIABLES FOR HARD PROCESS
 C     modifications for Les Houches accord by PR (7/15/02)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION CIRCKP(2)
       COMMON /HWCIR2/CIRCKP
       DOUBLE PRECISION Z1,Z2,C1,C2,B1,B2,CIRCEE,CIRCGG,RS,MISS,ETA,
       DOUBLE PRECISION CIRCKP(2)
       COMMON /HWCIR2/CIRCKP
       DOUBLE PRECISION Z1,Z2,C1,C2,B1,B2,CIRCEE,CIRCGG,RS,MISS,ETA,


@@ -15122,7 +15122,7 @@ C-----------------------------------------------------------------------
 C     SETS UP 2->2 HARD SUBPROCESS
 c BRW change 18/8/04: BW smearing of mass i only if SMRi is true
 C-----------------------------------------------------------------------
 C     SETS UP 2->2 HARD SUBPROCESS
 c BRW change 18/8/04: BW smearing of mass i only if SMRi is true
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWUMBW,HWUPCM,PA,PCM
       INTEGER ICMF,IBM,I,J,K,IHEP,NTRY
       LOGICAL SMR3,SMR4
       DOUBLE PRECISION HWUMBW,HWUPCM,PA,PCM
       INTEGER ICMF,IBM,I,J,K,IHEP,NTRY
       LOGICAL SMR3,SMR4


@@ -15217,7 +15217,7 @@ C...It factorises (PIFAC*ALPHA/SWEIN/RMW/RMW/SQRT(2.))**2
 C
 C...First release:  1-APR-1998 by Stefano Moretti
 C-----------------------------------------------------------------------
 C
 C...First release:  1-APR-1998 by Stefano Moretti
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER I
       DOUBLE PRECISION P1(0:3),P2(0:3),P3(0:3),P4(0:3)
       DOUBLE PRECISION P(0:3)
       INTEGER I
       DOUBLE PRECISION P1(0:3),P2(0:3),P3(0:3),P4(0:3)
       DOUBLE PRECISION P(0:3)


@@ -15310,7 +15310,7 @@ C     SECTION ROUTINE
 C     I-L are the particles (all outgoing)
 C     Z1 and Z2 are the decay products of the Z
 C-----------------------------------------------------------------------
 C     I-L are the particles (all outgoing)
 C     Z1 and Z2 are the decay products of the Z
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER ND,I,J,K,L,Z1,Z2
       DOUBLE COMPLEX HWH2DD,ZI,S,D,F
       PARAMETER(ZI=(0.0D0,1.0D0))
       INTEGER ND,I,J,K,L,Z1,Z2
       DOUBLE COMPLEX HWH2DD,ZI,S,D,F
       PARAMETER(ZI=(0.0D0,1.0D0))


@@ -15353,7 +15353,7 @@ C
 C...First release:  01-APR-1998 by Kosuke Odagiri
 C...First modified: 12-APR-1998 by Stefano Moretti
 C-----------------------------------------------------------------------
 C...First release:  01-APR-1998 by Kosuke Odagiri
 C...First modified: 12-APR-1998 by Stefano Moretti
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER MU,IRES,IFL
       DOUBLE PRECISION P1(0:3),P2(0:3),P3(0:3),P4(0:3),P5(0:3)
       DOUBLE PRECISION EMB,EMT,EMW,EMH,EMH01,EMH02,EMH03
       INTEGER MU,IRES,IFL
       DOUBLE PRECISION P1(0:3),P2(0:3),P3(0:3),P4(0:3),P5(0:3)
       DOUBLE PRECISION EMB,EMT,EMW,EMH,EMH01,EMH02,EMH03


@@ -15453,7 +15453,7 @@ C-----------------------------------------------------------------------
 C     Subroutine to implement the F function of Eijk and Kliess
 C     fixed first momenta and all second momenta
 C-----------------------------------------------------------------------
 C     Subroutine to implement the F function of Eijk and Kliess
 C     fixed first momenta and all second momenta
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION P(5),MQ,PM(5),XMASS,PLAB,PRW,PCM,HWULDO,PDOT,EPS
       DOUBLE COMPLEX F(2,2,8),S,D,SIP(2),SJP(2)
       INTEGER I,J,NP
       DOUBLE PRECISION P(5),MQ,PM(5),XMASS,PLAB,PRW,PCM,HWULDO,PDOT,EPS
       DOUBLE COMPLEX F(2,2,8),S,D,SIP(2),SJP(2)
       INTEGER I,J,NP


@@ -15497,7 +15497,7 @@ C-----------------------------------------------------------------------
 C     Subroutine to implement the F function of Eijk and Kliess
 C     fixed second momenta and all first momenta
 C-----------------------------------------------------------------------
 C     Subroutine to implement the F function of Eijk and Kliess
 C     fixed second momenta and all first momenta
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION P(5),MQ,PM(5),XMASS,PLAB,PRW,PCM,HWULDO,PDOT,EPS
       DOUBLE COMPLEX F(2,2,8),S,D,SIP(2),SJP(2)
       INTEGER I,J,NP
       DOUBLE PRECISION P(5),MQ,PM(5),XMASS,PLAB,PRW,PCM,HWULDO,PDOT,EPS
       DOUBLE COMPLEX F(2,2,8),S,D,SIP(2),SJP(2)
       INTEGER I,J,NP


@@ -15541,7 +15541,7 @@ C-----------------------------------------------------------------------
 C     Subroutine to implement the F function of Eijk and Kliess
 C     All first and second momenta
 C-----------------------------------------------------------------------
 C     Subroutine to implement the F function of Eijk and Kliess
 C     All first and second momenta
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION P(5),MQ,PM(5),XMASS,PLAB,PRW,PCM,HWULDO,PDOT,EPS
       DOUBLE COMPLEX F(2,2,8,8),SIP(2),SJP(2),S,D
       INTEGER I,J,NP
       DOUBLE PRECISION P(5),MQ,PM(5),XMASS,PLAB,PRW,PCM,HWULDO,PDOT,EPS
       DOUBLE COMPLEX F(2,2,8,8),SIP(2),SJP(2),S,D
       INTEGER I,J,NP


@@ -15786,7 +15786,7 @@ C-----------------------------------------------------------------------
 C     Massless matrix elements for gg-->qqZ and qq-->qqZ
 C     using the matrix elements given in Nucl. Phys. B262 (1985) 235-242
 C-----------------------------------------------------------------------
 C     Massless matrix elements for gg-->qqZ and qq-->qqZ
 C     using the matrix elements given in Nucl. Phys. B262 (1985) 235-242
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IQ,I,J,OZ(2,2),IDZ,P1,P2,P3,P4,IQI,ID(2),K
       DOUBLE PRECISION MG(2),MQ(2,5),G(12,2),FLOW(3,3),CQFC,CQIFC,
      &     CGFC,CGIFC
       INTEGER IQ,I,J,OZ(2,2),IDZ,P1,P2,P3,P4,IQI,ID(2),K
       DOUBLE PRECISION MG(2),MQ(2,5),G(12,2),FLOW(3,3),CQFC,CQIFC,
      &     CGFC,CGIFC


@@ -15933,7 +15933,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Massive matrix elements for gg --> qqbarZ and qqbar --> qqbarZ
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Massive matrix elements for gg --> qqbarZ and qqbar --> qqbarZ
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IQ,I,IDZ,P1,P2,PL,PB,PBB,O(2),J,IQI
       DOUBLE PRECISION MG(2),MQ(2,5),G(12,2),CQFC,CQIFC,CGFC,CGIFC,
      &     PTMP(5,10),XMASS,PLAB,PRW,PCM,HWULDO,QBL,QBBL,Q2B,Q1B,Q2BB,
       INTEGER IQ,I,IDZ,P1,P2,PL,PB,PBB,O(2),J,IQI
       DOUBLE PRECISION MG(2),MQ(2,5),G(12,2),CQFC,CQIFC,CGFC,CGIFC,
      &     PTMP(5,10),XMASS,PLAB,PRW,PCM,HWULDO,QBL,QBBL,Q2B,Q1B,Q2BB,


@@ -16329,7 +16329,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Phase Space for vector boson plus 2 jets
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Phase Space for vector boson plus 2 jets
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION WEIGHT,XMASS,PLAB,PRW,PCM,Y(3),Y35,Y34,Y45,RAND,
      &     HWRGEN,HWRUNI,M35,M35S,G(IMAXCH),DEM,MT(3),PT(3),MJAC,ETOT,
      &     STOT,MQ(3),MQ2(3),PS35,HWUPCM,TWOPI2,MT35,PTJ(3),MT2(3),A,C,
       DOUBLE PRECISION WEIGHT,XMASS,PLAB,PRW,PCM,Y(3),Y35,Y34,Y45,RAND,
      &     HWRGEN,HWRUNI,M35,M35S,G(IMAXCH),DEM,MT(3),PT(3),MJAC,ETOT,
      &     STOT,MQ(3),MQ2(3),PS35,HWUPCM,TWOPI2,MT35,PTJ(3),MT2(3),A,C,


@@ -16801,7 +16801,7 @@ C     Subroutine to select virtual quark mass for HWH2PS
 C     IOPT=1 return the function at M2
 C     IOPT=2 calculate M2
 C-----------------------------------------------------------------------
 C     IOPT=1 return the function at M2
 C     IOPT=2 calculate M2
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IOPT
       DOUBLE PRECISION FJAC,MPOW,MMN,MQ2,M2,A1,A01,RPOW,QPOW,HWRGEN,MMX
       EXTERNAL HWRGEN
       INTEGER IOPT
       DOUBLE PRECISION FJAC,MPOW,MMN,MQ2,M2,A1,A01,RPOW,QPOW,HWRGEN,MMX
       EXTERNAL HWRGEN


@@ -16841,7 +16841,7 @@ C     Subroutine to select virtual quark mass for HWH2PS
 C     IOPT=1 return the function at M2
 C     IOPT=2 calculate M2
 C-----------------------------------------------------------------------
 C     IOPT=1 return the function at M2
 C     IOPT=2 calculate M2
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IOPT
       DOUBLE PRECISION FJAC,MPOW,A1,A01,RPOW,QPOW,HWRGEN,PT2,
      &     PPOW,PTMN2,PTMX2,Z
       INTEGER IOPT
       DOUBLE PRECISION FJAC,MPOW,A1,A01,RPOW,QPOW,HWRGEN,PT2,
      &     PPOW,PTMN2,PTMX2,Z


@@ -17143,7 +17143,7 @@ C-----------------------------------------------------------------------
 C     Subroutine to calculate the spinor products in the notation of
 C     Kleiss and Strirling S(1) is S and S(2) is T
 C-----------------------------------------------------------------------
 C     Subroutine to calculate the spinor products in the notation of
 C     Kleiss and Strirling S(1) is S and S(2) is T
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION WRN(2),K(5),KK(5),P(5,2),Q1,Q2,EPS,QTI,PTI,
      &     PT,QT,DPM,DMP,QP,QM,P1,P2,PP,PM
       DOUBLE COMPLEX S(2),ZI,Z1,ZT,ZQ,ZQS,ZPS,ZP,ZDPM,ZDMP
       DOUBLE PRECISION WRN(2),K(5),KK(5),P(5,2),Q1,Q2,EPS,QTI,PTI,
      &     PT,QT,DPM,DMP,QP,QM,P1,P2,PP,PM
       DOUBLE COMPLEX S(2),ZI,Z1,ZT,ZQ,ZQS,ZPS,ZP,ZDPM,ZDMP


@@ -17212,7 +17212,7 @@ C     I-L are the particles
 C     Z1 and Z2 are the decay products of the Z
 C     P1 is the polarization of the line I,J
 C-----------------------------------------------------------------------
 C     Z1 and Z2 are the decay products of the Z
 C     P1 is the polarization of the line I,J
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX HWH2T1,S,D
       INTEGER I,J,K,L,Z1,Z2,P1
       COMMON/HWHEWS/S(8,8,2),D(8,8)
       DOUBLE COMPLEX HWH2T1,S,D
       INTEGER I,J,K,L,Z1,Z2,P1
       COMMON/HWHEWS/S(8,8,2),D(8,8)


@@ -17236,7 +17236,7 @@ C     Z1 and Z2 are the decay products of the Z
 C     P1 is the polarization of the line I,J
 C     P2 is the polarization of the gluon K
 C-----------------------------------------------------------------------
 C     P1 is the polarization of the line I,J
 C     P2 is the polarization of the gluon K
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX HWH2T2,S,D
       INTEGER I,J,K,L,Z1,Z2,P1,P2
       DOUBLE PRECISION B(6)
       DOUBLE COMPLEX HWH2T2,S,D
       INTEGER I,J,K,L,Z1,Z2,P1,P2
       DOUBLE PRECISION B(6)


@@ -17269,7 +17269,7 @@ C     Z1 and Z2 are the decay products of the Z
 C     P1 is the polarization of the line I,J
 C     P2 is the polarization of the gluon K
 C-----------------------------------------------------------------------
 C     P1 is the polarization of the line I,J
 C     P2 is the polarization of the gluon K
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX HWH2T3,S,D
       INTEGER I,J,K,L,Z1,Z2,P1,P2
       DOUBLE PRECISION B(6)
       DOUBLE COMPLEX HWH2T3,S,D
       INTEGER I,J,K,L,Z1,Z2,P1,P2
       DOUBLE PRECISION B(6)


@@ -17300,7 +17300,7 @@ C     Z1 and Z2 are the decay products of the Z
 C     P1 is the polarization of the line I,J
 C     P2 is the polarization of the line K,L
 C-----------------------------------------------------------------------
 C     P1 is the polarization of the line I,J
 C     P2 is the polarization of the line K,L
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX HWH2T4,AP,AM,S,D
       INTEGER I,J,K,L,Z1,Z2,J1,J2,J3,J4,P1,P2
       DOUBLE PRECISION B(6)
       DOUBLE COMPLEX HWH2T4,AP,AM,S,D
       INTEGER I,J,K,L,Z1,Z2,J1,J2,J3,J4,P1,P2
       DOUBLE PRECISION B(6)


@@ -17335,7 +17335,7 @@ C     Z1 and Z2 are the decay products of the Z
 C     P1 is the polarization of the line I,J
 C     P2 is the polarization of the line K,L
 C-----------------------------------------------------------------------
 C     P1 is the polarization of the line I,J
 C     P2 is the polarization of the line K,L
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX HWH2T5,AP,AM,S,D
       INTEGER I,J,K,L,Z1,Z2,J1,J2,J3,J4,P1,P2
       DOUBLE PRECISION B(6)
       DOUBLE COMPLEX HWH2T5,AP,AM,S,D
       INTEGER I,J,K,L,Z1,Z2,J1,J2,J3,J4,P1,P2
       DOUBLE PRECISION B(6)


@@ -17371,7 +17371,7 @@ C     P1 is the polarization of the line I,J
 C     P2 is the polarization of the gluon K
 C     P3 is the polarization of the gluon L
 C-----------------------------------------------------------------------
 C     P2 is the polarization of the gluon K
 C     P3 is the polarization of the gluon L
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX HWH2T6,S,D
       INTEGER I,J,K,L,Z1,Z2,J1,J2,P1,P2,P3
       DOUBLE PRECISION B(6)
       DOUBLE COMPLEX HWH2T6,S,D
       INTEGER I,J,K,L,Z1,Z2,J1,J2,P1,P2,P3
       DOUBLE PRECISION B(6)


@@ -17419,7 +17419,7 @@ C     P1 is the polarization of the line I,J
 C     P2 is the polarization of the gluon K
 C     P3 is the polarization of the gluon L
 C-----------------------------------------------------------------------
 C     P2 is the polarization of the gluon K
 C     P3 is the polarization of the gluon L
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX HWH2T7,S,D
       INTEGER I,J,K,L,Z1,Z2,J1,J2,P1,P2,P3
       DOUBLE PRECISION B(6)
       DOUBLE COMPLEX HWH2T7,S,D
       INTEGER I,J,K,L,Z1,Z2,J1,J2,P1,P2,P3
       DOUBLE PRECISION B(6)


@@ -17468,7 +17468,7 @@ C     P1 is the polarization of the line I,J
 C     P2 is the polarization of the gluon K
 C     P3 is the polarization of the gluon L
 C-----------------------------------------------------------------------
 C     P2 is the polarization of the gluon K
 C     P3 is the polarization of the gluon L
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX HWH2T8,S,D
       INTEGER I,J,K,L,Z1,Z2,J1,J2,P1,P2,P3
       DOUBLE PRECISION B(6)
       DOUBLE COMPLEX HWH2T8,S,D
       INTEGER I,J,K,L,Z1,Z2,J1,J2,P1,P2,P3
       DOUBLE PRECISION B(6)


@@ -17517,7 +17517,7 @@ C     P1 is the polarization of the line I,J
 C     P2 is the polarization of the gluon K
 C     P3 is the polarization of the gluon L
 C-----------------------------------------------------------------------
 C     P2 is the polarization of the gluon K
 C     P3 is the polarization of the gluon L
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX HWH2T9,S,D
       INTEGER I,J,K,L,Z1,Z2,J1,J2,P1,P2,P3
       DOUBLE PRECISION B(6)
       DOUBLE COMPLEX HWH2T9,S,D
       INTEGER I,J,K,L,Z1,Z2,J1,J2,P1,P2,P3
       DOUBLE PRECISION B(6)


@@ -17556,7 +17556,7 @@ C     P1 is the polarization of the line I,J
 C     P2 is the polarization of the gluon K
 C     P3 is the polarization of the gluon L
 C-----------------------------------------------------------------------
 C     P2 is the polarization of the gluon K
 C     P3 is the polarization of the gluon L
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX HWH2T0,S,D
       INTEGER I,J,K,L,Z1,Z2,J1,J2,P1,P2,P3
       DOUBLE PRECISION B(6)
       DOUBLE COMPLEX HWH2T0,S,D
       INTEGER I,J,K,L,Z1,Z2,J1,J2,P1,P2,P3
       DOUBLE PRECISION B(6)


@@ -17662,7 +17662,7 @@ C     Matrix elements based on Ellis Ross & Terrano and Catani & Seymour
 C
 C     WARNING:  Phase space factor inaccurate for JADE y_cut > 0.14.
 C-----------------------------------------------------------------------
 C
 C     WARNING:  Phase space factor inaccurate for JADE y_cut > 0.14.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER LM,LP,IQK,I,J,IDMN,IDMX,ID1,ID2,IST(4)
       DOUBLE PRECISION HWRGEN,HWUALF,HWUAEM,HWULDO,HWH4J1,HWH4J2,
      & HWH4J4,HWH4J5,HWH4J6,HWH4J7,QNOW,Q2NOW,QLST,SCUT,PSFAC,FACT,
       INTEGER LM,LP,IQK,I,J,IDMN,IDMX,ID1,ID2,IST(4)
       DOUBLE PRECISION HWRGEN,HWUALF,HWUAEM,HWULDO,HWH4J1,HWH4J2,
      & HWH4J4,HWH4J5,HWH4J6,HWH4J7,QNOW,Q2NOW,QLST,SCUT,PSFAC,FACT,


@@ -18276,7 +18276,7 @@ C
 C      Mean EVWGT = cross section in nanoBarn
 C
 C-----------------------------------------------------------------------
 C      Mean EVWGT = cross section in nanoBarn
 C
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,Y,Q2,SHAT,Z,PHI,AJACOB,DSIGMA,ME,MP,
      & ML,MREMIF(18),MFIN1(18),MFIN2(18),RS,SMA,W2,RSHAT,FSIGMA(18),
      & SIGSUM,PROB,PRAN,PVRT(4),X
       DOUBLE PRECISION HWRGEN,Y,Q2,SHAT,Z,PHI,AJACOB,DSIGMA,ME,MP,
      & ML,MREMIF(18),MFIN1(18),MFIN2(18),RS,SMA,W2,RSHAT,FSIGMA(18),
      & SIGSUM,PROB,PRAN,PVRT(4),X


@@ -18602,7 +18602,7 @@ C     gives the fourmomenta in the laboratory system for the particles
 C     of the hard 2-->3 subprocess, to match with HERWIG routines of
 C     jet evolution.
 C----------------------------------------------------------------------
 C     of the hard 2-->3 subprocess, to match with HERWIG routines of
 C     jet evolution.
 C----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWUECM,HWUPCM,HWUSQR,Y,Q2,SHAT,Z,PHI,AJACOB,
      & DSIGMA,ME,MP,ML,MREMIF(18),MFIN1(18),MFIN2(18),RS,SMA,W2,RSHAT,
      & PGAMMA(5),SG,MF1,MF2,EP,PP,EL,PL,E1,E2,Q1,COSBET,SINBET,COSTHE,
       DOUBLE PRECISION HWUECM,HWUPCM,HWUSQR,Y,Q2,SHAT,Z,PHI,AJACOB,
      & DSIGMA,ME,MP,ML,MREMIF(18),MFIN1(18),MFIN2(18),RS,SMA,W2,RSHAT,
      & PGAMMA(5),SG,MF1,MF2,EP,PP,EL,PL,E1,E2,Q1,COSBET,SINBET,COSTHE,


@@ -18804,7 +18804,7 @@ C     corresponding Jacobian factor AJACOB
 C     Fill the logical vector INSIDE to tag contributing subprocesses
 C     to the cross-section
 C-----------------------------------------------------------------------
 C     Fill the logical vector INSIDE to tag contributing subprocesses
 C     to the cross-section
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       LOGICAL IFGO
       DOUBLE PRECISION HWRUNI,HWRGEN,HWUPCM,Y,Q2,SHAT,Z,PHI,AJACOB,
      & DSIGMA,ME,MP,ML,MREMIF(18),MFIN1(18),MFIN2(18),RS,SMA,W2,RSHAT,
       LOGICAL IFGO
       DOUBLE PRECISION HWRUNI,HWRGEN,HWUPCM,Y,Q2,SHAT,Z,PHI,AJACOB,
      & DSIGMA,ME,MP,ML,MREMIF(18),MFIN1(18),MFIN2(18),RS,SMA,W2,RSHAT,


@@ -19148,7 +19148,7 @@ C----------------------------------------------------------------------
 C     Returns differential cross section DSIGMA in (Y,Q2,ETA,Z,PHI)
 C     Scale for structure functions and alpha_s selected by BGSHAT
 C----------------------------------------------------------------------
 C     Returns differential cross section DSIGMA in (Y,Q2,ETA,Z,PHI)
 C     Scale for structure functions and alpha_s selected by BGSHAT
 C----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWUALF,HWUAEM,Y,Q2,SHAT,Z,PHI,AJACOB,DSIGMA,
      & ME,MP,ML,MREMIF(18),MFIN1(18),MFIN2(18),RS,SMA,W2,RSHAT,
      & SFUN(13),ALPHA,LDSIG,DLQ(7),SG,XG,MF1,MF2,MSUM,MDIF,MPRO,FFUN,
       DOUBLE PRECISION HWUALF,HWUAEM,Y,Q2,SHAT,Z,PHI,AJACOB,DSIGMA,
      & ME,MP,ML,MREMIF(18),MFIN1(18),MFIN2(18),RS,SMA,W2,RSHAT,
      & SFUN(13),ALPHA,LDSIG,DLQ(7),SG,XG,MF1,MF2,MSUM,MDIF,MPRO,FFUN,


@@ -19331,7 +19331,7 @@ C----------------------------------------------------------------------
 C----------------------------------------------------------------------
 C     DEEP INELASTIC LEPTON-HADRON SCATTERING: MEAN EVWGT = SIGMA IN NB
 C----------------------------------------------------------------------
 C----------------------------------------------------------------------
 C     DEEP INELASTIC LEPTON-HADRON SCATTERING: MEAN EVWGT = SIGMA IN NB
 C----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,HWUPCM,PRAN,PROB,SAMP,SIG,Q2,
      & XBJ,Y,W,S,MLEP,MHAD,MLSCAT,YMIN,YMAX,XXMAX,Q2JAC,XXJAC,
      & JACOBI,A1,A2,A3,B1,B2,PCM,PCMEP,PCMLW,PCMEQ,PCMLQ,COSPHI,PA,
       DOUBLE PRECISION HWRGEN,HWRUNI,HWUPCM,PRAN,PROB,SAMP,SIG,Q2,
      & XBJ,Y,W,S,MLEP,MHAD,MLSCAT,YMIN,YMAX,XXMAX,Q2JAC,XXJAC,
      & JACOBI,A1,A2,A3,B1,B2,PCM,PCMEP,PCMLW,PCMEQ,PCMLQ,COSPHI,PA,


@@ -19619,7 +19619,7 @@ C
 C     Modified 16/01/01 by BRW to implement Peter Richardson's
 C     fix for bug in lepton mass effects on branching ratio
 C-----------------------------------------------------------------------
 C     Modified 16/01/01 by BRW to implement Peter Richardson's
 C     fix for bug in lepton mass effects on branching ratio
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,HWUAEM,EPS,C1,C2,C3,EMSQZ,EMGMZ,
      & EMSQZP,EMGMZP,CQF(7,6,16),QPOW,RPOW,A01,A1,A02,A2,A03,A3,CRAN,
      & EMJ1,EMJ2,EMJ3,EMJAC,FACT,QSQ,HCS,FACTR,RCS,EXTRA,PMAX,PTHETA
       DOUBLE PRECISION HWRGEN,HWRUNI,HWUAEM,EPS,C1,C2,C3,EMSQZ,EMGMZ,
      & EMSQZP,EMGMZP,CQF(7,6,16),QPOW,RPOW,A01,A1,A02,A2,A03,A3,CRAN,
      & EMJ1,EMJ2,EMJ3,EMJAC,FACT,QSQ,HCS,FACTR,RCS,EXTRA,PMAX,PTHETA


@@ -19809,7 +19809,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Drell-Yan production with a q qbar pair
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Drell-Yan production with a q qbar pair
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER I,MAP(12),ORD,IFL,IDP(6),IFLOW,QCFL(2,2),GCFL(2),IDZ,IQ
       DOUBLE PRECISION HCS,RCS,MQ(2,5),HWRGEN,G(12,2),DIST(2),MG(2)
       LOGICAL FSTCLL,MASS
       INTEGER I,MAP(12),ORD,IFL,IDP(6),IFLOW,QCFL(2,2),GCFL(2),IDZ,IQ
       DOUBLE PRECISION HCS,RCS,MQ(2,5),HWRGEN,G(12,2),DIST(2),MG(2)
       LOGICAL FSTCLL,MASS


@@ -19884,7 +19884,7 @@ C     AFTER CUTS ON PT AND MASS OF CENTRE-OF-MASS SYSTEM
 C     AND COS(THETA) IN CENTRE-OF-MASS SYSTEM
 C     AND TIMES BRANCHING FRACTION IF WW
 C-----------------------------------------------------------------------
 C     AND COS(THETA) IN CENTRE-OF-MASS SYSTEM
 C     AND TIMES BRANCHING FRACTION IF WW
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWULDO,EMSQ,BETA,S,T,U,TMIN,TMAX,TRAT,
      & DSDT,PROB,X,Z(2),ZMIN,ZMAX,PCMIN,PCMAX,PCFAC,PLOGMI,PLOGMA,PTCMF,
      & Q,PC,BLOG,EMCMIN,EMCMAX,EMLMIN,EMLMAX,WGT(6),RWGT,CV,CA,BR,QT(2),
       DOUBLE PRECISION HWRGEN,HWULDO,EMSQ,BETA,S,T,U,TMIN,TMAX,TRAT,
      & DSDT,PROB,X,Z(2),ZMIN,ZMAX,PCMIN,PCMAX,PCFAC,PLOGMI,PLOGMA,PTCMF,
      & Q,PC,BLOG,EMCMIN,EMCMAX,EMLMIN,EMLMAX,WGT(6),RWGT,CV,CA,BR,QT(2),


@@ -20155,7 +20155,7 @@ C----------------------------------------------------------------------
 C     W + GAMMA --> FF'BAR :  MEAN EVWGT = CROSS SECTION IN NANOBARN
 C     BASED ON BOSON GLUON FUSION OF ABBIENDI AND STANCO
 C-----------------------------------------------------------------------
 C     W + GAMMA --> FF'BAR :  MEAN EVWGT = CROSS SECTION IN NANOBARN
 C     BASED ON BOSON GLUON FUSION OF ABBIENDI AND STANCO
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,GMASS,EV(3),RV,Y,Q2,SHAT,Z,PHI,AJACOB,
      & DSIGMA,ME,MP,ML,MREMIF(18),MFIN1(18),MFIN2(18),RS,SMA,W2,RSHAT
       INTEGER LEP
       DOUBLE PRECISION HWRGEN,GMASS,EV(3),RV,Y,Q2,SHAT,Z,PHI,AJACOB,
      & DSIGMA,ME,MP,ML,MREMIF(18),MFIN1(18),MFIN2(18),RS,SMA,W2,RSHAT
       INTEGER LEP


@@ -20359,7 +20359,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     COMPUTES DIFFERENTIAL CROSS SECTION DSIGMA IN (Y,Q2,ETA,Z,PHI)
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     COMPUTES DIFFERENTIAL CROSS SECTION DSIGMA IN (Y,Q2,ETA,Z,PHI)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION TMAX,TMIN,A1,A2,B1,B2,I0,I1,I2,I3,I4,I5,MUSQ,
      & MDSQ,ETA,Q1,COSTHE,S,G,T,U,C1,C2,D1,D2,F1,F2,COSBET,WPROP,D(4,4),
      & C(4,4),QU,QD,QE,QW,PHOTON,EMWSQ,EMSSQ,CFAC,Y,Q2,SHAT,Z,PHI,
       DOUBLE PRECISION TMAX,TMIN,A1,A2,B1,B2,I0,I1,I2,I3,I4,I5,MUSQ,
      & MDSQ,ETA,Q1,COSTHE,S,G,T,U,C1,C2,D1,D2,F1,F2,COSBET,WPROP,D(4,4),
      & C(4,4),QU,QD,QE,QW,PHOTON,EMWSQ,EMSSQ,CFAC,Y,Q2,SHAT,Z,PHI,


@@ -20480,7 +20480,7 @@ C     If IPROC=107: --> gg, distributed as sum of light quarks.
 C     If fermion flavour specified mass effects fully included.
 C     EVWGT=sig(e+e- --> ffbar) in nb
 C-----------------------------------------------------------------------
 C     If fermion flavour specified mass effects fully included.
 C     EVWGT=sig(e+e- --> ffbar) in nb
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,HWUPCM,HWUAEM,Q2NOW,Q2LST,FACTR,
      & VF2,VF,CLF(7),PRAN,PQWT,PMAX,PTHETA,SINTH2,CPHI,SPHI,C2PHI,S2PHI,
      & PPHI,SINTH,PCM,PP(5),EWGT
       DOUBLE PRECISION HWRGEN,HWRUNI,HWUPCM,HWUAEM,Q2NOW,Q2LST,FACTR,
      & VF2,VF,CLF(7),PRAN,PQWT,PMAX,PTHETA,SINTH2,CPHI,SPHI,C2PHI,S2PHI,
      & PPHI,SINTH,PCM,PP(5),EWGT


@@ -20625,7 +20625,7 @@ c     scheme, y_cut=1.-THMAX.
 C     If flavour specified mass effects fully included.
 C     EVWGT=sig(e^-e^+ --> qqbar g) in nb
 C-----------------------------------------------------------------------
 C     If flavour specified mass effects fully included.
 C     EVWGT=sig(e^-e^+ --> qqbar g) in nb
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWUALF,HWUAEM,HWULDO,HWDPWT,Q2NOW,Q2LST,
      & PHASP,QGMAX,QGMIN,FACTR,QM2,CLF(7),ORDER,PRAN,PQWT,QQG,QBG,SUM,
      & RUT,QQLM,QQLP,QBLM,QBLP,DYN1,DYN2,DYN3,DYN4,DYN5,DYN6,XQ2,X2SUM,
       DOUBLE PRECISION HWRGEN,HWUALF,HWUAEM,HWULDO,HWDPWT,Q2NOW,Q2LST,
      & PHASP,QGMAX,QGMIN,FACTR,QM2,CLF(7),ORDER,PRAN,PQWT,QQG,QBG,SUM,
      & RUT,QQLM,QQLP,QBLM,QBLP,DYN1,DYN2,DYN3,DYN4,DYN5,DYN6,XQ2,X2SUM,


@@ -20799,7 +20799,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SUSY E+E- -> 2 SLEPTON PROCESSES
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SUSY E+E- -> 2 SLEPTON PROCESSES
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWUAEM,EPS,HCS,RCS,S,PF,QPE,HWUPCM,PCM,
      & FACTR,SN2TH,MZ,ME2(2,2,6),EMSC2,HWUMBW,HWRUNI,T,SQPE
       INTEGER ID1,ID2,IL,IL1,IL2,I,J,IG,IG1,IHEP,NTRY,IDL,ILP,IDLR(2),
       DOUBLE PRECISION HWRGEN,HWUAEM,EPS,HCS,RCS,S,PF,QPE,HWUPCM,PCM,
      & FACTR,SN2TH,MZ,ME2(2,2,6),EMSC2,HWUMBW,HWRUNI,T,SQPE
       INTEGER ID1,ID2,IL,IL1,IL2,I,J,IG,IG1,IHEP,NTRY,IDL,ILP,IDLR(2),


@@ -21026,7 +21026,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SUSY E+E- -> 2 GAUGINO PROCESSES
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SUSY E+E- -> 2 GAUGINO PROCESSES
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWUAEM,HCS,RCS,MNU(4),MNU2(4),HWRUNI,
      &                 FACA,M1(4,4),S2W,XC(4),XD(4),MSNU,
      &                 MZ,HWHSS2,U,T,QPE,SQPE,MSL,MSL2,MSR,MSR2,
       DOUBLE PRECISION HWRGEN,HWUAEM,HCS,RCS,MNU(4),MNU2(4),HWRUNI,
      &                 FACA,M1(4,4),S2W,XC(4),XD(4),MSNU,
      &                 MZ,HWHSS2,U,T,QPE,SQPE,MSL,MSL2,MSR,MSR2,


@@ -21275,7 +21275,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SUSY E+E- -> 2 SPARTICLE PROCESSES
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SUSY E+E- -> 2 SPARTICLE PROCESSES
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION SAVWT(3),RANWT,HWRGEN
       EXTERNAL HWRGEN
       SAVE SAVWT
       DOUBLE PRECISION SAVWT(3),RANWT,HWRGEN
       EXTERNAL HWRGEN
       SAVE SAVWT


@@ -21317,7 +21317,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SUSY E+E- -> 2 SQUARK PROCESSES
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SUSY E+E- -> 2 SQUARK PROCESSES
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWUAEM,EPS,HCS,RCS,S,PF,QPE,HWUPCM,PCM,
      & FACTR,SN2TH,MZ,ME2(2,2,6),EMSC2,HWUMBW,HWRUNI,SQPE
       INTEGER ID1,ID2,IQ,IQ1,IQ2,I,J,IHEP,IDL,IDLR(2),IDSQU(2),NTRY
       DOUBLE PRECISION HWRGEN,HWUAEM,EPS,HCS,RCS,S,PF,QPE,HWUPCM,PCM,
      & FACTR,SN2TH,MZ,ME2(2,2,6),EMSC2,HWUMBW,HWRUNI,SQPE
       INTEGER ID1,ID2,IQ,IQ1,IQ2,I,J,IHEP,IDL,IDLR(2),IDSQU(2),NTRY


@@ -21478,7 +21478,7 @@ CDECK  ID>, HWHEW0.
 C-----------------------------------------------------------------------
       SUBROUTINE HWHEW0(IP,ETOT,XM,PR,WEIGHT,CR)
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
       SUBROUTINE HWHEW0(IP,ETOT,XM,PR,WEIGHT,CR)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,ETOT,XM(2),PR(5,2),WEIGHT,CR,XM1,XM2,S,
      & D1,PABS,D,CX,C,E,F,SC,G
       INTEGER IP,I
       DOUBLE PRECISION HWRGEN,ETOT,XM(2),PR(5,2),WEIGHT,CR,XM1,XM2,S,
      & D1,PABS,D,CX,C,E,F,SC,G
       INTEGER IP,I


@@ -21645,7 +21645,7 @@ C NOTE: EXTERNAL FACTOR OF COLOR AVERAGE AND SPIN AVERAGE AND
 C COUPLING (E**8/4/9) MUST BE INCLUDED AS WELL AS COMPENSATION
 C FOR ON POLE APPROXIMATION AS DESIRED.
 C-----------------------------------------------------------------------
 C COUPLING (E**8/4/9) MUST BE INCLUDED AS WELL AS COMPENSATION
 C FOR ON POLE APPROXIMATION AS DESIRED.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX HWHEW4,ZH,ZCH,ZD,ZAMP1,ZAMP3,DWW,CWW,BWW,AWW,
      & AWWM,AWWP,AMPTEM,ZTWO,ZHALF
       DOUBLE PRECISION XW,ZMASS,T3,EQ1,RR,RL,ZM2,AMP2,RKW,COLFAC(4),
       DOUBLE COMPLEX HWHEW4,ZH,ZCH,ZD,ZAMP1,ZAMP3,DWW,CWW,BWW,AWW,
      & AWWM,AWWP,AMPTEM,ZTWO,ZHALF
       DOUBLE PRECISION XW,ZMASS,T3,EQ1,RR,RL,ZM2,AMP2,RKW,COLFAC(4),


@@ -21762,7 +21762,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     E+E- -> W+W-/Z0Z0 (BASED ON ZOLTAN KUNSZT'S PROGRAM)
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     E+E- -> W+W-/Z0Z0 (BASED ON ZOLTAN KUNSZT'S PROGRAM)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ZH,ZCH,ZD
       DOUBLE PRECISION HWUAEM,HWRGEN,HWUPCM,ETOT,STOT,FLUXW,GAMM,GIMM,
      & WM2,WXMIN,WX1MAX,WX2MAX,FJAC1,FJAC2,WX1,WX2,WMM1,WMM2,XXM,W2BO,
       DOUBLE COMPLEX ZH,ZCH,ZD
       DOUBLE PRECISION HWUAEM,HWRGEN,HWUPCM,ETOT,STOT,FLUXW,GAMM,GIMM,
      & WM2,WXMIN,WX1MAX,WX2MAX,FJAC1,FJAC2,WX1,WX2,WMM1,WMM2,XXM,W2BO,


@@ -21955,7 +21955,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Hadron-Hadron to WW/WZ/ZZ (BASED ON ZOLTAN KUNSZT'S PROGRAM)
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Hadron-Hadron to WW/WZ/ZZ (BASED ON ZOLTAN KUNSZT'S PROGRAM)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX ZH,ZCH,ZD
       DOUBLE PRECISION HWUAEM,HWRGEN,HWUPCM,FLUXW,CSW,XMASS,
      &     PLAB,PRW,PCM,HWRUNI,P(5,10),AMPWW,DIST(4),MW2,CFAC1,AMP,
       DOUBLE COMPLEX ZH,ZCH,ZD
       DOUBLE PRECISION HWUAEM,HWRGEN,HWUPCM,FLUXW,CSW,XMASS,
      &     PLAB,PRW,PCM,HWRUNI,P(5,10),AMPWW,DIST(4),MW2,CFAC1,AMP,


@@ -22208,7 +22208,7 @@ C     ICH returns the channel used if OPTM=.FALSE.
 C     ICH specifies the channel to be used if OPTM=.TRUE.
 C     This is used in optimising the weights for the different channels
 C-----------------------------------------------------------------------
 C     ICH specifies the channel to be used if OPTM=.TRUE.
 C     This is used in optimising the weights for the different channels
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER ICH,IDBOS,ISM(2,IMAXCH),I,J,IB(2),IDRES,IDP,IOPT,IPRC,ID1
       DOUBLE PRECISION XMASS,PLAB,PRW,PCM,RAND,HWRGEN,BMS2(2),TJAC,PLM,
      &     MJAC(2),TWOPI2,SJAC,STOT,THAT,UHAT,TMIN,TMAX,UMIN,UMAX,PS(2),
       INTEGER ICH,IDBOS,ISM(2,IMAXCH),I,J,IB(2),IDRES,IDP,IOPT,IPRC,ID1
       DOUBLE PRECISION XMASS,PLAB,PRW,PCM,RAND,HWRGEN,BMS2(2),TJAC,PLM,
      &     MJAC(2),TWOPI2,SJAC,STOT,THAT,UHAT,TMIN,TMAX,UMIN,UMAX,PS(2),


@@ -22428,7 +22428,7 @@ C     ISM=2 select according to power law  for IDBOZ
 C     IOPT=1 return the function at MBOS2
 C     IOPT=2 calculate MBOS2
 C-----------------------------------------------------------------------
 C     IOPT=1 return the function at MBOS2
 C     IOPT=2 calculate MBOS2
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IDBOZ,ISM,IOPT
       DOUBLE PRECISION MBOZ,FJAC,GBOZ,GMBOZ,MPOW,MMIN,
      &     MBOS2,A1,A2,A01,A02,RPOW,QPOW,HWRGEN,MMAX,EMSQ
       INTEGER IDBOZ,ISM,IOPT
       DOUBLE PRECISION MBOZ,FJAC,GBOZ,GMBOZ,MPOW,MMIN,
      &     MBOS2,A1,A2,A01,A02,RPOW,QPOW,HWRGEN,MMAX,EMSQ


@@ -22491,7 +22491,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     WW cross section in hadron hadron
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     WW cross section in hadron hadron
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HCS,RCS,HWRGEN,DIST(2),CFAC,WAMP(2),S34,S56,KWW2,
      &     MW2,MZ2,GMW,GMZ,G(4,2),EE(4),CKM2(12),RF(2),LF(2),TAUI(2),
      &     CSW,CFAC1
       DOUBLE PRECISION HCS,RCS,HWRGEN,DIST(2),CFAC,WAMP(2),S34,S56,KWW2,
      &     MW2,MZ2,GMW,GMZ,G(4,2),EE(4),CKM2(12),RF(2),LF(2),TAUI(2),
      &     CSW,CFAC1


@@ -22571,7 +22571,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     ZZ cross section in hadron hadron
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     ZZ cross section in hadron hadron
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION AMP(2),RCS,HCS,HWRGEN,DIST(2),S34,S56,CFAC,
      &     MW2,MZ2,GMW,GMZ,G(4,2),EE(4),CKM2(12),RF(2),LF(2),TAUI(2),
      &     CSW,CFAC1
       DOUBLE PRECISION AMP(2),RCS,HCS,HWRGEN,DIST(2),S34,S56,CFAC,
      &     MW2,MZ2,GMW,GMZ,G(4,2),EE(4),CKM2(12),RF(2),LF(2),TAUI(2),
      &     CSW,CFAC1


@@ -22684,7 +22684,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     WZ cross section in hadron hadron
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     WZ cross section in hadron hadron
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION AMP(2),HCS,RCS,HWRGEN,W34,DIST(2),S34,S56,CFAC,
      &     TCS,S12,MW2,MZ2,GMW,GMZ,G(4,2),EE(4),CKM2(12),RF(2),LF(2),
      &     TAUI(2),CSW,CFAC1
       DOUBLE PRECISION AMP(2),HCS,RCS,HWRGEN,W34,DIST(2),S34,S56,CFAC,
      &     TCS,S12,MW2,MZ2,GMW,GMZ,G(4,2),EE(4),CKM2(12),RF(2),LF(2),
      &     TAUI(2),CSW,CFAC1


@@ -22802,7 +22802,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Subroutine to select t or u for HWHGBP
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Subroutine to select t or u for HWHGBP
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IOPT
       DOUBLE PRECISION FJAC,TPOW,TMIN,T,A1,A01,RPOW,QPOW,HWRGEN,TMAX,TN,
      &     TX,MT
       INTEGER IOPT
       DOUBLE PRECISION FJAC,TPOW,TMIN,T,A1,A01,RPOW,QPOW,HWRGEN,TMAX,TN,
      &     TX,MT


@@ -22845,7 +22845,7 @@ C     viz L = - G(mu,nu) T(mu,nu) / GRVLAM
 C     If GAMGRV is zero, it is revaluated during the first run
 C     MEAN EVWGT = SIGMA IN NB
 C-----------------------------------------------------------------------
 C     If GAMGRV is zero, it is revaluated during the first run
 C     MEAN EVWGT = SIGMA IN NB
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,EPS,EMSQG,
      & EMGMG,S,CC,SS,SS2,M1(16),M2(16),M3,M4,M5(3),M6(3),
      & RNGLU,FACT,HCS,FACTR,RCS,A2,A02,QPE,SQPE,RGRV
       DOUBLE PRECISION HWRGEN,HWRUNI,EPS,EMSQG,
      & EMGMG,S,CC,SS,SS2,M1(16),M2(16),M3,M4,M5(3),M6(3),
      & RNGLU,FACT,HCS,FACTR,RCS,A2,A02,QPE,SQPE,RGRV


@@ -23208,7 +23208,7 @@ C----------------------------------------------------------------------
 C     Use the GUPI (Generic User Process Interface) event common block
 C     as the hard process for HERWIG
 C----------------------------------------------------------------------
 C     Use the GUPI (Generic User Process Interface) event common block
 C     as the hard process for HERWIG
 C----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

 C--Les Houches Common Block
       INTEGER MAXPUP
       PARAMETER(MAXPUP=100)
 C--Les Houches Common Block
       INTEGER MAXPUP
       PARAMETER(MAXPUP=100)


@@ -23550,7 +23550,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     QCD HEAVY FLAVOUR PRODUCTION: MEAN EVWGT = SIGMA IN NB
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     QCD HEAVY FLAVOUR PRODUCTION: MEAN EVWGT = SIGMA IN NB
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,EPS,RCS,Z1,Z2,ET,EJ,
      & QM2,QPE,FACTR,S,T,U,ST,TU,US,TUS,UST,EN,RN,AF,ASTU,
      & AUST,CF,CN,CS,CSTU,CSUT,CTSU,CTUS,HCS,UT,SU,GT,DIST,KK,KK2,
       DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,EPS,RCS,Z1,Z2,ET,EJ,
      & QM2,QPE,FACTR,S,T,U,ST,TU,US,TUS,UST,EN,RN,AF,ASTU,
      & AUST,CF,CN,CS,CSTU,CSUT,CTSU,CTUS,HCS,UT,SU,GT,DIST,KK,KK2,


@@ -23830,7 +23830,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     HIGGS + HEAVY QUARK (BOTTOM & TOP) PRODUCTION (2HDM)
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     HIGGS + HEAVY QUARK (BOTTOM & TOP) PRODUCTION (2HDM)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN, HWUALF, HWUAEM, EPS, HCS, RCS,
      & DIST, SM, DM, QPE, PF, SQPE, EMSC2, FACTR, S, T3, U4,
      & SN2TH, ME2(0:4), MW, XWEIN, PT2MIN, PT2, GQH(0:4), G1, RMMIN,
       DOUBLE PRECISION HWRGEN, HWUALF, HWUAEM, EPS, HCS, RCS,
      & DIST, SM, DM, QPE, PF, SQPE, EMSC2, FACTR, S, T3, U4,
      & SN2TH, ME2(0:4), MW, XWEIN, PT2MIN, PT2, GQH(0:4), G1, RMMIN,


@@ -24095,7 +24095,7 @@ C
 C-----------------------------------------------------------------------
 C     ASSOCIATE PRODUCTION W+H- FROM QUARK FUSION (2HDM)
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     ASSOCIATE PRODUCTION W+H- FROM QUARK FUSION (2HDM)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER I,J,IHEL
       DOUBLE PRECISION EMH,EMHWT,RMW,EMW
       DOUBLE PRECISION RMH
       INTEGER I,J,IHEL
       DOUBLE PRECISION EMH,EMHWT,RMW,EMW
       DOUBLE PRECISION RMH


@@ -24368,7 +24368,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     PRODUCTION OF MSSM CHARGED HIGGSES FROM B-QUARK+LIGHT-QUARK FUSION
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     PRODUCTION OF MSSM CHARGED HIGGSES FROM B-QUARK+LIGHT-QUARK FUSION
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER I,J,K,L,M,N
       INTEGER II,JJ,ITMP
       INTEGER IFL,IRES
       INTEGER I,J,K,L,M,N
       INTEGER II,JJ,ITMP
       INTEGER IFL,IRES


@@ -24715,7 +24715,7 @@ C     IAPHIG (set in HWIGIN)=0: zero mass approximation =1: exact result
 C                           =2: infinite mass limit.
 C     Only top loop included. A factor (alpha_s**3*alpha_W) is extracted
 C-----------------------------------------------------------------------
 C                           =2: infinite mass limit.
 C     Only top loop included. A factor (alpha_s**3*alpha_W) is extracted
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX HWHIGB,HWHIGC,HWHIGD,HWHIG5,HWHIG1,HWHIG2,BI(4),
      & CI(7),DI(3),EPSI,TAMP(7)
       DOUBLE PRECISION S,T,U,EMH2,WTQQ,WTQG,WTGQ,WTGG,EMW2,RNGLU,RNQRK,
       DOUBLE COMPLEX HWHIGB,HWHIGC,HWHIGD,HWHIG5,HWHIG1,HWHIG2,BI(4),
      & CI(7),DI(3),EPSI,TAMP(7)
       DOUBLE PRECISION S,T,U,EMH2,WTQQ,WTQG,WTGQ,WTGG,EMW2,RNGLU,RNQRK,


@@ -24801,7 +24801,7 @@ C-----------------------------------------------------------------------
 C     One loop scalar integrals, used in HWHIGJ.
 C     If NOMASS=.TRUE. use a small mass approx. for particle in loop.
 C-----------------------------------------------------------------------
 C     One loop scalar integrals, used in HWHIGJ.
 C     If NOMASS=.TRUE. use a small mass approx. for particle in loop.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX HWHIGB,HWUCI2,HWULI2,EPSI,PII
       DOUBLE PRECISION S,T,EQ2,EH2,RAT
       LOGICAL NOMASS
       DOUBLE COMPLEX HWHIGB,HWUCI2,HWULI2,EPSI,PII
       DOUBLE PRECISION S,T,EQ2,EH2,RAT
       LOGICAL NOMASS


@@ -24837,7 +24837,7 @@ C-----------------------------------------------------------------------
 C     One loop scalar integrals, used in HWHIGJ.
 C     If NOMASS=.TRUE. use a small mass approx. for particle in loop.
 C-----------------------------------------------------------------------
 C     One loop scalar integrals, used in HWHIGJ.
 C     If NOMASS=.TRUE. use a small mass approx. for particle in loop.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX HWHIGC,HWUCI2,HWULI2,EPSI,PII
       DOUBLE PRECISION S,T,EQ2,EH2,RAT,COSH
       LOGICAL NOMASS
       DOUBLE COMPLEX HWHIGC,HWUCI2,HWULI2,EPSI,PII
       DOUBLE PRECISION S,T,EQ2,EH2,RAT,COSH
       LOGICAL NOMASS


@@ -24873,7 +24873,7 @@ C-----------------------------------------------------------------------
 C     One loop scalar integrals, used in HWHIGJ.
 C     If NOMASS=.TRUE. use a small mass approx. for particle in loop.
 C-----------------------------------------------------------------------
 C     One loop scalar integrals, used in HWHIGJ.
 C     If NOMASS=.TRUE. use a small mass approx. for particle in loop.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE COMPLEX HWHIGD,HWUCI2,HWULI2,EPSI,PII,Z1,Z2
       DOUBLE PRECISION S,T,EQ2,EH2,DLS,DLT,DLM,RZ12,DL1,DL2,
      & ST,ROOT,XP,XM
       DOUBLE COMPLEX HWHIGD,HWUCI2,HWULI2,EPSI,PII,Z1,Z2
       DOUBLE PRECISION S,T,EQ2,EH2,DLS,DLT,DLM,RZ12,DL1,DL2,
      & ST,ROOT,XP,XM


@@ -24934,7 +24934,7 @@ C
 C--------------------------------------------------------------------------
 C     LEPTOPRODUCTION OF MS(SM) HIGGSES IN ASSOCIATION WITH HEAVY QUARK PAIRS
 C--------------------------------------------------------------------------
 C--------------------------------------------------------------------------
 C     LEPTOPRODUCTION OF MS(SM) HIGGSES IN ASSOCIATION WITH HEAVY QUARK PAIRS
 C--------------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER JHIGGS
       INTEGER I,L,M,N,NN
       INTEGER IH,IQ,JQ,IIQ,JJQ
       INTEGER JHIGGS
       INTEGER I,L,M,N,NN
       INTEGER IH,IQ,JQ,IIQ,JJQ


@@ -25454,7 +25454,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     DRELL-YAN 2 PARTON -> 2 HIGGS PAIR (2HDM)
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     DRELL-YAN 2 PARTON -> 2 HIGGS PAIR (2HDM)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN, HWUAEM, EPS, HCS, RCS, DIST, S, PF, QPE,
      & FACTR, SN2TH, MZ, MW, MNN(2,2), MCC(2), MCN(3), EMSC2, GW2, GZ2,
      & GHH(4), XWEIN, S2W, ECM_MAX, X(3), XL(3),
       DOUBLE PRECISION HWRGEN, HWUAEM, EPS, HCS, RCS, DIST, S, PF, QPE,
      & FACTR, SN2TH, MZ, MW, MNN(2,2), MCC(2), MCN(3), EMSC2, GW2, GZ2,
      & GHH(4), XWEIN, S2W, ECM_MAX, X(3), XL(3),


@@ -25746,7 +25746,7 @@ C     QCD Higgs plus jet production; mean EVWGT = Sigma in nb*Higgs B.R.
 C     Adapted from the program of U. Baur and E.W.N. Glover
 C     See: Nucl. Phys. B339 (1990) 38
 C-----------------------------------------------------------------------
 C     Adapted from the program of U. Baur and E.W.N. Glover
 C     See: Nucl. Phys. B339 (1990) 38
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,HWUAEM,EPS,RCS,EMH,EMHWT,
      & EMHTMP,BR,CV,CA,EMH2,ET,EJ,PT,EMT,EMAX,YMAX,YHINF,YHSUP,EXYH,
      & YMIN,YJINF,YJSUP,EXYJ,S,T,U,FACT,AMPQQ,AMPQG,AMPGQ,AMPGG,HCS,
       DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,HWUAEM,EPS,RCS,EMH,EMHWT,
      & EMHTMP,BR,CV,CA,EMH2,ET,EJ,PT,EMT,EMAX,YMAX,YHINF,YHSUP,EXYH,
      & YMIN,YJINF,YJSUP,EXYJ,S,T,U,FACT,AMPQQ,AMPQG,AMPGQ,AMPGG,HCS,


@@ -25905,7 +25905,7 @@ C                   /  (EM**2-EMH**2)**2 + (GAM(EM)*EMH**2/EM)**2
 C       AS SUGGESTED IN M.H.SEYMOUR, PHYS.LETT.B354(1995)409.
 C     ENDIF
 C-----------------------------------------------------------------------
 C       AS SUGGESTED IN M.H.SEYMOUR, PHYS.LETT.B354(1995)409.
 C     ENDIF
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRUNI,EM,WEIGHT,EMH,DIF,FUN,THETA,T,EMHLST,W0,
      & W1,EMM,GAMEM,T0,TMIN,TMAX,THEMIN,THEMAX,ZMIN,ZMAX,Z,F,GAMOFS
       INTEGER I
       DOUBLE PRECISION HWRUNI,EM,WEIGHT,EMH,DIF,FUN,THETA,T,EMHLST,W0,
      & W1,EMM,GAMEM,T0,TMIN,TMAX,THEMIN,THEMAX,ZMIN,ZMAX,Z,F,GAMOFS
       INTEGER I


@@ -25996,7 +25996,7 @@ C
 C-----------------------------------------------------------------------
 C     PRODUCTION OF MSSM HIGGSES IN ASSOCIATION WITH HEAVY QUARK PAIRS
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     PRODUCTION OF MSSM HIGGSES IN ASSOCIATION WITH HEAVY QUARK PAIRS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER JHIGGS
       INTEGER I,J,K,L,M,N
       INTEGER IS,IH,IQ,JQ,IIQ,JJQ,IQMIN,IQMAX,IGG,IQQ
       INTEGER JHIGGS
       INTEGER I,J,K,L,M,N
       INTEGER IS,IH,IQ,JQ,IIQ,JJQ,IQMIN,IQMAX,IGG,IQQ


@@ -26521,7 +26521,7 @@ C-----------------------------------------------------------------------
 C     HIGGS PRODUCTION VIA GLUON OR QUARK FUSION
 C     MEAN EVWGT = HIGGS PRODN C-S * BRANCHING FRACTION IN NB
 C-----------------------------------------------------------------------
 C     HIGGS PRODUCTION VIA GLUON OR QUARK FUSION
 C     MEAN EVWGT = HIGGS PRODN C-S * BRANCHING FRACTION IN NB
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWUALF,HWHIGT,HWRGEN,HWUSQR,HWUAEM,BRHIGQ,EMH,
      & CSFAC(13),EVSUM(13),EMFAC,CV,CA,BR,RWGT,E1,E2,EMQ,GFACTR,RQM(6)
       INTEGER IDEC,I,J,ID1,ID2
       DOUBLE PRECISION HWUALF,HWHIGT,HWRGEN,HWUSQR,HWUAEM,BRHIGQ,EMH,
      & CSFAC(13),EVSUM(13),EMFAC,CV,CA,BR,RWGT,E1,E2,EMQ,GFACTR,RQM(6)
       INTEGER IDEC,I,J,ID1,ID2


@@ -26619,7 +26619,7 @@ C  CALCULATE MOD SQUARED I DEFINED AS IN BARGER & PHILLIPS p433
 C  WARNING: THIS IS A FACTOR OF 3 GREATER THAN EHLQ'S ETA FUNCTION
 C  PARITY=+1 FOR SCALAR AND -1 FOR PSEUDOSCALAR
 C-----------------------------------------------------------------------
 C  WARNING: THIS IS A FACTOR OF 3 GREATER THAN EHLQ'S ETA FUNCTION
 C  PARITY=+1 FOR SCALAR AND -1 FOR PSEUDOSCALAR
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWHIGT,RATIO,RAT2,EMH,FREAL,FIMAG,ETALOG,AIREAL,
      & AIIMAG
       INTEGER I,J,K,L
       DOUBLE PRECISION HWHIGT,RATIO,RAT2,EMH,FREAL,FIMAG,ETALOG,AIREAL,
      & AIIMAG
       INTEGER I,J,K,L


@@ -26694,7 +26694,7 @@ C-----------------------------------------------------------------------
 C     MSSM NEUTRAL HIGGS PRODUCTION IN ASSOCIATION WITH GAUGE BOSON
 C--BRW fix 27/8/04: corrected off-shell gauge boson mass dependence
 C-----------------------------------------------------------------------
 C     MSSM NEUTRAL HIGGS PRODUCTION IN ASSOCIATION WITH GAUGE BOSON
 C--BRW fix 27/8/04: corrected off-shell gauge boson mass dependence
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER I,J,K,L,M,N
       INTEGER IV,IDEC
       INTEGER ID1,ID2
       INTEGER I,J,K,L,M,N
       INTEGER IV,IDEC
       INTEGER ID1,ID2


@@ -26980,7 +26980,7 @@ C-----------------------------------------------------------------------
 C     HIGGS PRODUCTION VIA W/Z BOSON FUSION
 C     MEAN EVWGT = HIGGS PRODN C-S * BRANCHING FRACTION IN NB
 C-----------------------------------------------------------------------
 C     HIGGS PRODUCTION VIA W/Z BOSON FUSION
 C     MEAN EVWGT = HIGGS PRODN C-S * BRANCHING FRACTION IN NB
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWULDO,HWRUNI,HWRGEN,HWUAEM,K1MAX2,K1MIN2,K12,
      & K2MAX2,K2MIN2,K22,EMW2,EMW,ROOTS,EMH2,EMH,ROOTS2,P1,PHI1,PHI2,
      & COSPHI,COSTH1,SINTH1,COSTH2,SINTH2,P2,WEIGHT,TAU,TAULN,CSFAC,
       DOUBLE PRECISION HWULDO,HWRUNI,HWRGEN,HWUAEM,K1MAX2,K1MIN2,K12,
      & K2MAX2,K2MIN2,K22,EMW2,EMW,ROOTS,EMH2,EMH,ROOTS2,P1,PHI1,PHI2,
      & COSPHI,COSTH1,SINTH1,COSTH2,SINTH2,P2,WEIGHT,TAU,TAULN,CSFAC,


@@ -27378,7 +27378,7 @@ C     USES ALGORITHM OF BERENDS AND KLEISS: NUCL.PHYS. B260(1985)32
 C
 C     MEAN EVWGT = CROSS-SECTION (IN NB) * HIGGS BRANCHING FRACTION
 C-----------------------------------------------------------------------
 C
 C     MEAN EVWGT = CROSS-SECTION (IN NB) * HIGGS BRANCHING FRACTION
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWUAEM,HWHIGY,HWRUNI,HWRGEN,HWULDO,EMZ,CVE,CAE,
      & POL1,POL2,CE1,CE2,CE3,PMAX,EMZ2,S,B,FACTR,EMH,EMFAC,EMH2,A,XP,
      & CV,CA,BRHIGQ,BR,X1,X2,FAC1,FAC2,XPP,XPPSQ,COEF,X,XSQ,PROB,C1,C2,
       DOUBLE PRECISION HWUAEM,HWHIGY,HWRUNI,HWRGEN,HWULDO,EMZ,CVE,CAE,
      & POL1,POL2,CE1,CE2,CE3,PMAX,EMZ2,S,B,FACTR,EMH,EMFAC,EMH2,A,XP,
      & CV,CA,BRHIGQ,BR,X1,X2,FAC1,FAC2,XPP,XPPSQ,COEF,X,XSQ,PROB,C1,C2,


@@ -27589,7 +27589,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     PRODUCTION OF MSSM HIGGS PAIRS IN L+L- (L=E,MU)
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     PRODUCTION OF MSSM HIGGS PAIRS IN L+L- (L=E,MU)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN, HWUAEM, HCS, RCS, S, PF, QPE,
      & FACTR, SN2TH, MZ, MNN(2), MCC, EMSC2, GZ2,
      & GHH(4), XWEIN, S2W, X(1), XL(1),
       DOUBLE PRECISION HWRGEN, HWUAEM, HCS, RCS, S, PF, QPE,
      & FACTR, SN2TH, MZ, MNN(2), MCC, EMSC2, GZ2,
      & GHH(4), XWEIN, S2W, X(1), XL(1),


@@ -27747,7 +27747,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     PRODUCTION OF MSSM HIGGSES IN ASSOCIATION WITH B,T-SQUARK PAIRS
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     PRODUCTION OF MSSM HIGGSES IN ASSOCIATION WITH B,T-SQUARK PAIRS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       COMMON/SQSQH/JHIGGS,ILBL,JH,IF1MIN,IF1MAX,IF2MIN,IF2MAX
       INTEGER      JHIGGS,ILBL,JH,IF1MIN,IF1MAX,IF2MIN,IF2MAX
       INTEGER I,J,K,L,M,N
       COMMON/SQSQH/JHIGGS,ILBL,JH,IF1MIN,IF1MAX,IF2MIN,IF2MAX
       INTEGER      JHIGGS,ILBL,JH,IF1MIN,IF1MAX,IF2MIN,IF2MAX
       INTEGER I,J,K,L,M,N


@@ -28141,7 +28141,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     QQD direct photon pair production: mean EVWGT = sigma in nb
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     QQD direct photon pair production: mean EVWGT = sigma in nb
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,HWHPPB,EPS,RCS,ET,EJ,KK,KK2,
      & YJ1INF,YJ1SUP,Z1,YJ2INF,YJ2SUP,Z2,FACT,FACTR,RS,S,T,U,CSTU,TQSQ,
      & DSTU,HCS
       DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,HWHPPB,EPS,RCS,ET,EJ,KK,KK2,
      & YJ1INF,YJ1SUP,Z1,YJ2INF,YJ2SUP,Z2,FACT,FACTR,RS,S,T,U,CSTU,TQSQ,
      & DSTU,HCS


@@ -28235,7 +28235,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     QCD DIRECT PHOTON + JET PRODUCTION: MEAN EVWGT = SIGMA IN NB
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     QCD DIRECT PHOTON + JET PRODUCTION: MEAN EVWGT = SIGMA IN NB
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,HWHPPB,EPS,RCS,ET,EJ,KK,KK2,
      & YJ1INF,YJ1SUP,Z1,YJ2INF,YJ2SUP,Z2,FACT,FACTR,FACTF,RS,S,T,U,CF,
      & AF,CSTU,CTSU,CUST,DSTU,HCS,TQCH
       DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,HWHPPB,EPS,RCS,ET,EJ,KK,KK2,
      & YJ1INF,YJ1SUP,Z1,YJ2INF,YJ2SUP,Z2,FACT,FACTR,FACTF,RS,S,T,U,CF,
      & AF,CSTU,CTSU,CUST,DSTU,HCS,TQCH


@@ -28393,7 +28393,7 @@ C-----------------------------------------------------------------------
 C     point-like photon/QCD heavy flavour single excitation, using exact
 C     massive lightcone kinematics, mean EVWGT = sigma in nb.
 C-----------------------------------------------------------------------
 C     point-like photon/QCD heavy flavour single excitation, using exact
 C     massive lightcone kinematics, mean EVWGT = sigma in nb.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,EPS,PP1,PP2,QM2,FACTR,
      & PT,PJ,PT2,PTM,EXY,T,CC,EXY2,S,U,C,SIGE,HCS,RCS
       INTEGER IQ1,IQ2,ID1,ID2,IHAD1,IHAD2
       DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,EPS,PP1,PP2,QM2,FACTR,
      & PT,PJ,PT2,PTM,EXY,T,CC,EXY2,S,U,C,SIGE,HCS,RCS
       INTEGER IQ1,IQ2,ID1,ID2,IHAD1,IHAD2


@@ -28475,7 +28475,7 @@ C-----------------------------------------------------------------------
 C     Point-like photon/gluon heavy flavour pair production, with
 C     exact lightcone massive kinematics, mean EVWGT = sigma in nb.
 C-----------------------------------------------------------------------
 C     Point-like photon/gluon heavy flavour pair production, with
 C     exact lightcone massive kinematics, mean EVWGT = sigma in nb.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRUNI,HWUALF,EPS,PP1,PP2,QM2,FACTR,ET,EJ,ET2,
      & EXY,EXY2,S,T,U,C
       INTEGER IQ1,IHAD1,IHAD2
       DOUBLE PRECISION HWRUNI,HWUALF,EPS,PP1,PP2,QM2,FACTR,ET,EJ,ET2,
      & EXY,EXY2,S,T,U,C
       INTEGER IQ1,IHAD1,IHAD2


@@ -28545,7 +28545,7 @@ C     Point-like photon/QCD direct meson production
 C     See M. Benayoun, et al., Nucl. Phys. B282 (1987) 653 for details.
 C     mean EVWGT = sigma in nb
 C-----------------------------------------------------------------------
 C     See M. Benayoun, et al., Nucl. Phys. B282 (1987) 653 for details.
 C     mean EVWGT = sigma in nb
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,EPS,PP1,PP2,ET,EJ,EXY,EXY2,
      & FACT,FACTR,S,T,U,REDS,DELT(3,3),C1STU,C3STU,HCS,RCS,CMIX,SMIX,
      & C1WVFN,FPI,FETA8,FETA1,FRHO,FPHI8,FPHI1,FPI2,FETA2(3),FETAP2(3),
       DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,EPS,PP1,PP2,ET,EJ,EXY,EXY2,
      & FACT,FACTR,S,T,U,REDS,DELT(3,3),C1STU,C3STU,HCS,RCS,CMIX,SMIX,
      & C1WVFN,FPI,FETA8,FETA1,FRHO,FPHI8,FPHI1,FPI2,FETA2(3),FETAP2(3),


@@ -28795,7 +28795,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     point-like photon/QCD di-jet production: mean EVWGT = sigma in nb
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     point-like photon/QCD di-jet production: mean EVWGT = sigma in nb
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,EPS,RCS,PP1,PP2,ET,EJ,
      & EXY,EXY2,FACTR,RS,S,T,U,CSTU,CTSU,HCS
       INTEGER ID1,ID2,ID3,ID4,IHAD1,IHAD2
       DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,EPS,RCS,PP1,PP2,ET,EJ,
      & EXY,EXY2,FACTR,RS,S,T,U,CSTU,CTSU,HCS
       INTEGER ID1,ID2,ID3,ID4,IHAD1,IHAD2


@@ -28881,7 +28881,7 @@ C-----------------------------------------------------------------------
 C     Compton scattering of point-like photon and (anti)quark
 C     mean EVWGT = sigma in nb
 C-----------------------------------------------------------------------
 C     Compton scattering of point-like photon and (anti)quark
 C     mean EVWGT = sigma in nb
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,EPS,RCS,PP1,PP2,ET,EJ,EXY,EXY2,
      & FACTR,S,T,U,CTSU,HCS
       INTEGER ID1,ID2,IHAD1,IHAD2
       DOUBLE PRECISION HWRGEN,HWRUNI,EPS,RCS,PP1,PP2,ET,EJ,EXY,EXY2,
      & FACTR,S,T,U,CTSU,HCS
       INTEGER ID1,ID2,IHAD1,IHAD2


@@ -28951,7 +28951,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     QCD HARD 2->2 PROCESSES: MEAN EVWGT = SIGMA IN NB
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     QCD HARD 2->2 PROCESSES: MEAN EVWGT = SIGMA IN NB
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,RS,EPS,HF,RCS,Z1,Z2,ET,EJ,
      & FACTR,S,T,U,ST,TU,US,STU,TUS,UST,EN,RN,GFLA,AF,ASTU,ASUT,AUST,
      & BF,BSTU,BSUT,BUST,BUTS,CF,CSTU,CSUT,CTSU,CTUS,DF,DSTU,DTSU,DUTS,
       DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,RS,EPS,HF,RCS,Z1,Z2,ET,EJ,
      & FACTR,S,T,U,ST,TU,US,STU,TUS,UST,EN,RN,GFLA,AF,ASTU,ASUT,AUST,
      & BF,BSTU,BSUT,BUST,BUTS,CF,CSTU,CSUT,CTSU,CTUS,DF,DSTU,DTSU,DUTS,


@@ -29305,7 +29305,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     IDENTIFIES HARD SUBPROCESS
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     IDENTIFIES HARD SUBPROCESS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER HWRINT,ID3,ID4,IPERM,IHPR,ND3
       EXTERNAL HWRINT
       IHPRO=IHPR
       INTEGER HWRINT,ID3,ID4,IPERM,IHPR,ND3
       EXTERNAL HWRINT
       IHPRO=IHPR


@@ -29337,7 +29337,7 @@ C-----------------------------------------------------------------------
 C     HARD PROCESS: GAMGAM --> QQBAR/LLBAR/W+W-
 C     MEAN EVENT WEIGHT = CROSS-SECTION IN NB AFTER CUTS ON PT
 C-----------------------------------------------------------------------
 C     HARD PROCESS: GAMGAM --> QQBAR/LLBAR/W+W-
 C     MEAN EVENT WEIGHT = CROSS-SECTION IN NB AFTER CUTS ON PT
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION RCS,HCS,RS,S,EMSQ,BE,TMIN,TMAX,T,U,FACTR,Q,CFAC,
      $     HWRGEN
       INTEGER IHAD1,IHAD2,HQ,ID3,ID4,I1,I2
       DOUBLE PRECISION RCS,HCS,RS,S,EMSQ,BE,TMIN,TMAX,T,U,FACTR,Q,CFAC,
      $     HWRGEN
       INTEGER IHAD1,IHAD2,HQ,ID3,ID4,I1,I2


@@ -29424,7 +29424,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C  Subroutine for 2 parton -> 2 parton via UDD resonant squarks
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C  Subroutine for 2 parton -> 2 parton via UDD resonant squarks
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HCS,S,RCS,HWRGEN,MQ1,MQ2,TAU,LOWTLM,UPPTLM,RTAB,
      &                 SQSH,MATELM,SCF(12),CHANPB(2),HWRUNI,PCM,MIX(12),
      &                 ME(2,3,3,3,3),WD,MS(12),SWD(12),RAND,TAUA,
       DOUBLE PRECISION HCS,S,RCS,HWRGEN,MQ1,MQ2,TAU,LOWTLM,UPPTLM,RTAB,
      &                 SQSH,MATELM,SCF(12),CHANPB(2),HWRUNI,PCM,MIX(12),
      &                 ME(2,3,3,3,3),WD,MS(12),SWD(12),RAND,TAUA,


@@ -29657,7 +29657,7 @@ C-----------------------------------------------------------------------
 C  Subroutine for 2 parton -> parton SUSY particle via UDD resonant
 C  squarks.
 C-----------------------------------------------------------------------
 C  Subroutine for 2 parton -> parton SUSY particle via UDD resonant
 C  squarks.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HCS,S,RCS,HWRGEN,ME(4),CW,MER(6),MZ,TAU,TAUA,
      &                 TAUB,LOWTLM,UPPTLM,HWRUNI,SH,SQSH,SCF(12),MW2,
      &                 LAMC(3),CHANPB(2),PCM,ECM,RAND,MEN(7,6,3,3),
       DOUBLE PRECISION HCS,S,RCS,HWRGEN,ME(4),CW,MER(6),MZ,TAU,TAUA,
      &                 TAUB,LOWTLM,UPPTLM,HWRUNI,SH,SQSH,SCF(12),MW2,
      &                 LAMC(3),CHANPB(2),PCM,ECM,RAND,MEN(7,6,3,3),


@@ -30592,7 +30592,7 @@ C-----------------------------------------------------------------------
 C     SUSY E+E- --> SM PARTICLES VIA RPV
 C     MODIFIED TO INCLUDE BEAM POLARIZATION EFFECTS BY PETER RICHARDSON
 C-----------------------------------------------------------------------
 C     SUSY E+E- --> SM PARTICLES VIA RPV
 C     MODIFIED TO INCLUDE BEAM POLARIZATION EFFECTS BY PETER RICHARDSON
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,HWUPCM,HWUMBW,HWUAEM,HCS,RCS,FACA,
      &                 S,T,PCM,MQ1,MQ2,SP,TP,TPZ,TPN,TPN2,MSL2(3),MZ,
      &                 MZ2,MSU2(3,2),MWD(3),GL,GR,GLP,GRP,EC,EE,THTMIN,
       DOUBLE PRECISION HWRGEN,HWRUNI,HWUPCM,HWUMBW,HWUAEM,HCS,RCS,FACA,
      &                 S,T,PCM,MQ1,MQ2,SP,TP,TPZ,TPN,TPN2,MSL2(3),MZ,
      &                 MZ2,MSU2(3,2),MWD(3),GL,GR,GLP,GRP,EC,EE,THTMIN,


@@ -30861,7 +30861,7 @@ C-----------------------------------------------------------------------
 C     IDENTIFY THE REMNANTS OF THE HARD SCATTERING
 C     AND BREAK THEIR COLOUR CONNECTION IF NECESSARY
 C-----------------------------------------------------------------------
 C     IDENTIFY THE REMNANTS OF THE HARD SCATTERING
 C     AND BREAK THEIR COLOUR CONNECTION IF NECESSARY
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION PCL(5),
      $     P1P2,P1SQ,P2SQ,S,M1SQ,M2SQ,TMP1,TMP2,A,B,C,D,PTOT(4),HWULDO
       INTEGER IBEAM,ITARG,IHEP,NTEMP,I,ICOL,IANT
       DOUBLE PRECISION PCL(5),
      $     P1P2,P1SQ,P2SQ,S,M1SQ,M2SQ,TMP1,TMP2,A,B,C,D,PTOT(4),HWULDO
       INTEGER IBEAM,ITARG,IHEP,NTEMP,I,ICOL,IANT


@@ -30964,7 +30964,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SUSY E+E- RPV PRODUCTION
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SUSY E+E- RPV PRODUCTION
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       IF(IPROC.GE.800.AND.IPROC.LE.850) THEN
         CALL HWHRES
       ELSEIF(IPROC.GE.860.AND.IPROC.LT.890) THEN
       IF(IPROC.GE.800.AND.IPROC.LE.850) THEN
         CALL HWHRES
       ELSEIF(IPROC.GE.860.AND.IPROC.LT.890) THEN


@@ -30983,7 +30983,7 @@ C-----------------------------------------------------------------------
 C     SUSY E+E- --> RPV SINGLE SPARTICLE PRODUCTION
 C     POLARZATION EFFECTS ADDED 5/4/02 BY PETER RICHARDSON
 C-----------------------------------------------------------------------
 C     SUSY E+E- --> RPV SINGLE SPARTICLE PRODUCTION
 C     POLARZATION EFFECTS ADDED 5/4/02 BY PETER RICHARDSON
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWUAEM,HWRUNI,HWUPCM,HWUMBW,HCS,RCS,FACA,
      &                 FACB,FACC,FACD,FACE,M1(4,4),M2(2,4),M3(8,2),
      &                 MW,MZ,MSCL(2,2),MSCL2(2,2),MZ2,MSL2,MSR2,MSNU2,
       DOUBLE PRECISION HWRGEN,HWUAEM,HWRUNI,HWUPCM,HWUMBW,HCS,RCS,FACA,
      &                 FACB,FACC,FACD,FACE,M1(4,4),M2(2,4),M3(8,2),
      &                 MW,MZ,MSCL(2,2),MSCL2(2,2),MZ2,MSL2,MSR2,MSNU2,


@@ -31631,7 +31631,7 @@ C  Subroutine for resonant sleptons to standard model particles
 C  slepton mass and mass*width added to save statement to
 C  avoid problems with Linux by Peter Richardson
 C-----------------------------------------------------------------------
 C  slepton mass and mass*width added to save statement to
 C  avoid problems with Linux by Peter Richardson
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HCS,S,RCS,HWRGEN,FAC,ECM,TH,PCM,CFAC,CHANPB,SH,
      &                 TAU,TAUA,TAUB,LOWTLM,UPPTLM,HWRUNI,MSL(12),
      &                 SQSH,MET(2),SCF(12),MIX(12),ME(4,3,3,3,3,2),
       DOUBLE PRECISION HCS,S,RCS,HWRGEN,FAC,ECM,TH,PCM,CFAC,CHANPB,SH,
      &                 TAU,TAUA,TAUB,LOWTLM,UPPTLM,HWRUNI,MSL(12),
      &                 SQSH,MET(2),SCF(12),MIX(12),ME(4,3,3,3,3,2),


@@ -31872,7 +31872,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C  Subroutine for 2 parton -> sparticle + X via LQD
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C  Subroutine for 2 parton -> sparticle + X via LQD
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HCS,A(6,12),B(6,12),S,RCS,HWRGEN,CW,FAC2,EC,ME2,
      &               MW,G,TAU,TAUA,TAUB,LOWTLM,UPPTLM,HWRUNI,SW,SQSH,LC,
      &               SH,MSL(12),MSU(12),MST(6),C(2,6,12),D(2,6,12),UH,
       DOUBLE PRECISION HCS,A(6,12),B(6,12),S,RCS,HWRGEN,CW,FAC2,EC,ME2,
      &               MW,G,TAU,TAUA,TAUB,LOWTLM,UPPTLM,HWRUNI,SW,SQSH,LC,
      &               SH,MSL(12),MSU(12),MST(6),C(2,6,12),D(2,6,12),UH,


@@ -32677,7 +32677,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Subroutine for all hadron-hadron Rparity violating processes
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Subroutine for all hadron-hadron Rparity violating processes
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       IF(MOD(IPROC,10000).GE.4000.AND.MOD(IPROC,10000).LT.4060) THEN
 C--SINGLE SPARTICLE VIA LQD
         CALL HWHRLS
       IF(MOD(IPROC,10000).GE.4000.AND.MOD(IPROC,10000).LT.4060) THEN
 C--SINGLE SPARTICLE VIA LQD
         CALL HWHRLS


@@ -32703,7 +32703,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     IDENTIDY HARD R-PARITY VIOLATING PROCESS
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     IDENTIDY HARD R-PARITY VIOLATING PROCESS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER ID3, ID4, R4, IPERM,TYPE,ID1,ID2,NEUTD1(8),SLEPID(8),
      &        NEUTD2(13),SQUID(6),SGN,HWUANT,SQUID2(12),SLPID2(12),
      &        GAGID1(6),GAGID2(8)
       INTEGER ID3, ID4, R4, IPERM,TYPE,ID1,ID2,NEUTD1(8),SLEPID(8),
      &        NEUTD2(13),SQUID(6),SGN,HWUANT,SQUID2(12),SLPID2(12),
      &        GAGID1(6),GAGID2(8)


@@ -32803,7 +32803,7 @@ C     JMUEO IS THE UNDERLYING EVENT OPTION: 1=>VETO EVENTS WITH M
 C     SCATTERS ABOVE PTMIN WITH PROBABILITY 1/(M+1)
 C     PTJIM IS THE MINIMUM TRANSVERSE MOMENTUM FOR ADDITIONAL SCATTERS
 C-----------------------------------------------------------------------
 C     SCATTERS ABOVE PTMIN WITH PROBABILITY 1/(M+1)
 C     PTJIM IS THE MINIMUM TRANSVERSE MOMENTUM FOR ADDITIONAL SCATTERS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRGET,HWRSET,WGT,PBOOST(5),RBOOST(3,3),
      $     WJMAX,PT,PTJIM,DUMMY,HWUPCM
       INTEGER IHEP,IBM,ITG,IBMN,ITGN,IBMT,ITGT,I,REPORT,NHARD,
       DOUBLE PRECISION HWRGEN,HWRGET,HWRSET,WGT,PBOOST(5),RBOOST(3,3),
      $     WJMAX,PT,PTJIM,DUMMY,HWUPCM
       INTEGER IHEP,IBM,ITG,IBMN,ITGN,IBMT,ITGT,I,REPORT,NHARD,


@@ -33023,7 +33023,7 @@ C     SWAP THE HARD PROCESS GENERATION PARAMETERS,
 C     CALL HWHQCD, AND SWAP BACK
 C     WGT IS THE OUTPUT EVENT WEIGHT
 C-----------------------------------------------------------------------
 C     CALL HWHQCD, AND SWAP BACK
 C     WGT IS THE OUTPUT EVENT WEIGHT
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION WGT,PTJIM,XMIN,XMAX,XPOW,
      $     TMPXMN,TMPXMX,TMPXPW,TMPWGT
       LOGICAL FIRST
       DOUBLE PRECISION WGT,PTJIM,XMIN,XMAX,XPOW,
      $     TMPXMN,TMPXMX,TMPXPW,TMPWGT
       LOGICAL FIRST


@@ -33066,7 +33066,7 @@ C     MEAN EVWGT = SIGMA IN NB
 C     TREATS ALL PARTONS ON EQUAL FOOTING WITH HWHSNM(ID1,ID2,S,T)
 C     PROVIDING THE MATRIX ELEMENT SQUARED FOR PARTON TYPES ID1 AND ID2
 C-----------------------------------------------------------------------
 C     TREATS ALL PARTONS ON EQUAL FOOTING WITH HWHSNM(ID1,ID2,S,T)
 C     PROVIDING THE MATRIX ELEMENT SQUARED FOR PARTON TYPES ID1 AND ID2
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER ID1,ID2
       DOUBLE PRECISION HWRGEN,HWRUNI,HWHSNM,EPS,RCS,ET,EJ,KK,KK2,
      & YJ1INF,YJ1SUP,Z1,YJ2INF,YJ2SUP,Z2,FACT,S,T,U,HCS
       INTEGER ID1,ID2
       DOUBLE PRECISION HWRGEN,HWRUNI,HWHSNM,EPS,RCS,ET,EJ,KK,KK2,
      & YJ1INF,YJ1SUP,Z1,YJ2INF,YJ2SUP,Z2,FACT,S,T,U,HCS


@@ -33134,7 +33134,7 @@ C     FOR PHOTON EXCHANGE, INTERFERENCE WITH U-CHANNEL CONTRIBUTION IS
 C     INCLUDED FOR IDENTICAL QUARKS AND LIKEWISE S-CHANNEL CONTRIBUTION
 C     FOR IDENTICAL QUARK-ANTIQUARK PAIRS.
 C-----------------------------------------------------------------------
 C     INCLUDED FOR IDENTICAL QUARKS AND LIKEWISE S-CHANNEL CONTRIBUTION
 C     FOR IDENTICAL QUARK-ANTIQUARK PAIRS.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWHSNM,HWUAEM,HWUALF,S,T,ASQ,AINU,AINS,Y,SOLD,
      $ TOLD,QQ(13,13),ZETA3
       INTEGER ID1,ID2
       DOUBLE PRECISION HWHSNM,HWUAEM,HWUALF,S,T,ASQ,AINU,AINS,Y,SOLD,
      $ TOLD,QQ(13,13),ZETA3
       INTEGER ID1,ID2


@@ -33199,7 +33199,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Calculates the spin correlations for the hard process
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Calculates the spin correlations for the hard process
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ZI,S,D,ME(2,2,2,2,NCFMAX),MED(2,2,2,2),F3(2,2,8),
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ZI,S,D,ME(2,2,2,2,NCFMAX),MED(2,2,2,2),F3(2,2,8),


@@ -35882,7 +35882,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section f fbar --> gauge boson --> fermion fermion
 C  This diagram 1 from DAMTP-2001-83 with opposite sign of P4
 C-----------------------------------------------------------------------
 C  section f fbar --> gauge boson --> fermion fermion
 C  This diagram 1 from DAMTP-2001-83 with opposite sign of P4
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI


@@ -35922,7 +35922,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section  f fbar ---> fermion fermion via t-channel scalar exchange
 C  This diagram 2 from DAMTP-2001-83 with opposite sign of P4
 C-----------------------------------------------------------------------
 C  section  f fbar ---> fermion fermion via t-channel scalar exchange
 C  This diagram 2 from DAMTP-2001-83 with opposite sign of P4
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE


@@ -35955,7 +35955,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section  f fbar ---> fermion fermion via u-channel scalar exchange
 C  This diagram 3 from DAMTP-2001-83 with opposite sign of P4
 C-----------------------------------------------------------------------
 C  section  f fbar ---> fermion fermion via u-channel scalar exchange
 C  This diagram 3 from DAMTP-2001-83 with opposite sign of P4
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,


@@ -35988,7 +35988,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section f fbar --> gauge boson --> fermion antifermion
 C  This diagram 1 from DAMTP-2001-83
 C-----------------------------------------------------------------------
 C  section f fbar --> gauge boson --> fermion antifermion
 C  This diagram 1 from DAMTP-2001-83
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI


@@ -36029,7 +36029,7 @@ C  section gluon gluon --> fermion antifermion (1st colour flow)
 C  N.B. a gauge choice has been made to simplify the triple gluon vertex
 C  This diagram 4 from DAMTP-2001-83 with the gauge choice L1=2 L2=1
 C-----------------------------------------------------------------------
 C  N.B. a gauge choice has been made to simplify the triple gluon vertex
 C  This diagram 4 from DAMTP-2001-83 with the gauge choice L1=2 L2=1
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI


@@ -36067,7 +36067,7 @@ C  section gluon gluon --> fermion antifermion (2st colour flow)
 C  N.B. a gauge choice has been made to simplify the triple gluon vertex
 C  This diagram 5 from DAMTP-2001-83 with the gauge choice L1=2 L2=1
 C-----------------------------------------------------------------------
 C  N.B. a gauge choice has been made to simplify the triple gluon vertex
 C  This diagram 5 from DAMTP-2001-83 with the gauge choice L1=2 L2=1
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI


@@ -36105,7 +36105,7 @@ C  section gluon gluon --> fermion antifermion (triple gluon piece)
 C  N.B. a gauge choice has been made to simplify the triple gluon vertex
 C  This diagram 6 from DAMTP-2001-83 with the gauge choice L1=2 L2=1
 C-----------------------------------------------------------------------
 C  N.B. a gauge choice has been made to simplify the triple gluon vertex
 C  This diagram 6 from DAMTP-2001-83 with the gauge choice L1=2 L2=1
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI


@@ -36146,7 +36146,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section quark gluon --> fermion sfermion
 C  This diagram 7 from DAMTP-2001-83 with the gauge choice L2=1
 C-----------------------------------------------------------------------
 C  section quark gluon --> fermion sfermion
 C  This diagram 7 from DAMTP-2001-83 with the gauge choice L2=1
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI


@@ -36183,7 +36183,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section antiquark gluon --> fermion antisfermion
 C  This diagram 10 from DAMTP-2001-83 with the gauge choice L2=1
 C-----------------------------------------------------------------------
 C  section antiquark gluon --> fermion antisfermion
 C  This diagram 10 from DAMTP-2001-83 with the gauge choice L2=1
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI


@@ -36220,7 +36220,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section quark gluon --> fermion antisfermion (s-channel quark)
 C  This is diagram 8 from DAMTP-2001-83 with the gauge choice L2=1
 C-----------------------------------------------------------------------
 C  section quark gluon --> fermion antisfermion (s-channel quark)
 C  This is diagram 8 from DAMTP-2001-83 with the gauge choice L2=1
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI


@@ -36263,7 +36263,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section quark gluon --> fermion antisfermion (s-channel quark)
 C  This is diagram 11 from DAMTP-2001-83 with the gauge choice L2=1
 C-----------------------------------------------------------------------
 C  section quark gluon --> fermion antisfermion (s-channel quark)
 C  This is diagram 11 from DAMTP-2001-83 with the gauge choice L2=1
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI


@@ -36304,7 +36304,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section quark gluon --> fermion antisfermion (s-channel quark)
 C  This is diagram 9 from DAMTP-2001-83 with the gauge choice L2=1
 C-----------------------------------------------------------------------
 C  section quark gluon --> fermion antisfermion (s-channel quark)
 C  This is diagram 9 from DAMTP-2001-83 with the gauge choice L2=1
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI


@@ -36342,7 +36342,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section quark gluon --> fermion antisfermion (s-channel quark)
 C  This is diagram 12 from DAMTP-2001-83 with the gauge choice L2=1
 C-----------------------------------------------------------------------
 C  section quark gluon --> fermion antisfermion (s-channel quark)
 C  This is diagram 12 from DAMTP-2001-83 with the gauge choice L2=1
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI


@@ -36382,7 +36382,7 @@ C  N.B. a gauge choice has been made to simplify the triple gluon vertex
 C  This diagram 4 from DAMTP-2001-83 with opposite helicity for 4
 C  and gauge choice L1=2 L2=1
 C-----------------------------------------------------------------------
 C  This diagram 4 from DAMTP-2001-83 with opposite helicity for 4
 C  and gauge choice L1=2 L2=1
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI,
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI,


@@ -36422,7 +36422,7 @@ C  N.B. a gauge choice has been made to simplify the triple gluon vertex
 C  This diagram 5 from DAMTP-2001-83 with opposite helicity for 4
 C  and gauge choice L1=2 L2=1
 C-----------------------------------------------------------------------
 C  This diagram 5 from DAMTP-2001-83 with opposite helicity for 4
 C  and gauge choice L1=2 L2=1
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI,
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI,


@@ -36462,7 +36462,7 @@ C  N.B. a gauge choice has been made to simplify the triple gluon vertex
 C  This diagram 6 from DAMTP-2001-83 with opposite helicity for 4
 C  and gauge choice L1=2 L2=1
 C-----------------------------------------------------------------------
 C  This diagram 6 from DAMTP-2001-83 with opposite helicity for 4
 C  and gauge choice L1=2 L2=1
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI


@@ -36504,7 +36504,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section fermion fermion --> fermion fermion (t-channel boson)
 C  This diagram 13 from DAMTP-2001-83
 C-----------------------------------------------------------------------
 C  section fermion fermion --> fermion fermion (t-channel boson)
 C  This diagram 13 from DAMTP-2001-83
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI,
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI,


@@ -36547,7 +36547,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section fermion antifermion --> fermion antifermion (t-channel boson)
 C  This diagram 14 from DAMTP-2001-83
 C-----------------------------------------------------------------------
 C  section fermion antifermion --> fermion antifermion (t-channel boson)
 C  This diagram 14 from DAMTP-2001-83
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI,
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI,


@@ -36590,7 +36590,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section antifermion fermion --> antifermion fermion (t-channel boson)
 C  This diagram 15 from DAMTP-2001-83
 C-----------------------------------------------------------------------
 C  section antifermion fermion --> antifermion fermion (t-channel boson)
 C  This diagram 15 from DAMTP-2001-83
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI,
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI,


@@ -36633,7 +36633,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section antifermion fermion --> antifermion fermion (t-channel boson)
 C  This diagram 16 from DAMTP-2001-83
 C-----------------------------------------------------------------------
 C  section antifermion fermion --> antifermion fermion (t-channel boson)
 C  This diagram 16 from DAMTP-2001-83
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI,
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI,


@@ -36676,7 +36676,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section  f fbar ---> f fbar via s-channel scalar exchange
 C  This is diagram 1 from RPV notes
 C-----------------------------------------------------------------------
 C  section  f fbar ---> f fbar via s-channel scalar exchange
 C  This is diagram 1 from RPV notes
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE


@@ -36711,7 +36711,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section  f fbar ---> f fbar via t-channel scalar exchange
 C  This is diagram 2 from RPV notes
 C-----------------------------------------------------------------------
 C  section  f fbar ---> f fbar via t-channel scalar exchange
 C  This is diagram 2 from RPV notes
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE


@@ -36744,7 +36744,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section  f fbar ---> fermion fermion via t-channel scalar exchange
 C  This is diagram 3 from RPV notes
 C-----------------------------------------------------------------------
 C  section  f fbar ---> fermion fermion via t-channel scalar exchange
 C  This is diagram 3 from RPV notes
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE


@@ -36777,7 +36777,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section  f fbar ---> f f via s-channel scalar exchange
 C  This is diagram 4 from RPV notes
 C-----------------------------------------------------------------------
 C  section  f fbar ---> f f via s-channel scalar exchange
 C  This is diagram 4 from RPV notes
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE


@@ -36812,7 +36812,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section  f fbar ---> f f via u-channel scalar exchange
 C  This is diagram 5 from RPV notes
 C-----------------------------------------------------------------------
 C  section  f fbar ---> f f via u-channel scalar exchange
 C  This is diagram 5 from RPV notes
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE


@@ -36845,7 +36845,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section  f fbar ---> f f via t-channel scalar exchange
 C  This is diagram 6 from RPV notes
 C-----------------------------------------------------------------------
 C  section  f fbar ---> f f via t-channel scalar exchange
 C  This is diagram 6 from RPV notes
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE


@@ -36878,7 +36878,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section  f f ---> f fbar via s-channel scalar exchange
 C  This is diagram 7 from RPV notes
 C-----------------------------------------------------------------------
 C  section  f f ---> f fbar via s-channel scalar exchange
 C  This is diagram 7 from RPV notes
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE


@@ -36913,7 +36913,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section  f f ---> f fbar via t-channel scalar exchange
 C  This is diagram 8 from RPV notes
 C-----------------------------------------------------------------------
 C  section  f f ---> f fbar via t-channel scalar exchange
 C  This is diagram 8 from RPV notes
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE


@@ -36947,7 +36947,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section  f f ---> f fbar via u-channel scalar exchange
 C  This is diagram 9 from RPV notes
 C-----------------------------------------------------------------------
 C  section  f f ---> f fbar via u-channel scalar exchange
 C  This is diagram 9 from RPV notes
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE


@@ -36981,7 +36981,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section  fbar fbar ---> f f via s-channel scalar exchange
 C  This is diagram 10 from RPV notes
 C-----------------------------------------------------------------------
 C  section  fbar fbar ---> f f via s-channel scalar exchange
 C  This is diagram 10 from RPV notes
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE


@@ -37016,7 +37016,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section  fbar fbar ---> f f via t-channel scalar exchange
 C  This is diagram 11 from RPV notes
 C-----------------------------------------------------------------------
 C  section  fbar fbar ---> f f via t-channel scalar exchange
 C  This is diagram 11 from RPV notes
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE


@@ -37050,7 +37050,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section  fbar fbar ---> f f via u-channel scalar exchange
 C  This is diagram 12 from RPV notes
 C-----------------------------------------------------------------------
 C  section  fbar fbar ---> f f via u-channel scalar exchange
 C  This is diagram 12 from RPV notes
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE


@@ -37084,7 +37084,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section  f f ---> f f via s-channel scalar exchange
 C  This is diagram 13 from RPV
 C-----------------------------------------------------------------------
 C  section  f f ---> f f via s-channel scalar exchange
 C  This is diagram 13 from RPV
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE


@@ -37119,7 +37119,7 @@ C  Subroutine to calculate the helicity amplitudes for the 2-to-2 cross
 C  section  fbar fbar ---> fbar fbar via t-channel scalar exchange
 C  This is diagram 14 from RPV notes
 C-----------------------------------------------------------------------
 C  section  fbar fbar ---> fbar fbar via t-channel scalar exchange
 C  This is diagram 14 from RPV notes
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
       INTEGER NDIAHD
       PARAMETER(NDIAHD=10)
       DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE


@@ -37168,7 +37168,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     LL(BAR) -> GAUGINOS (including beam polarization)
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     LL(BAR) -> GAUGINOS (including beam polarization)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWHSS2, S, T, U, M3, M4, SGN
       DOUBLE COMPLEX CLL, CLR, CRL, CRR
       HWHSS2 =
       DOUBLE PRECISION HWHSS2, S, T, U, M3, M4, SGN
       DOUBLE COMPLEX CLL, CLR, CRL, CRR
       HWHSS2 =


@@ -37192,7 +37192,7 @@ C-----------------------------------------------------------------------
 C     SUSY 2 PARTON -> 2 GAUGINOS PROCESSES        (1 - 3)
 C                   -> GAUGINO + SPARTON PROCESSES (4 - 7)
 C-----------------------------------------------------------------------
 C     SUSY 2 PARTON -> 2 GAUGINOS PROCESSES        (1 - 3)
 C                   -> GAUGINO + SPARTON PROCESSES (4 - 7)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN, HWUALF, HWUAEM, EPS, HCS, RCS, DIST,
      & ML(6), ML2(6), MR(6), MR2(6), MCH(2), MCH2(2), MNU(4), MNU2(4),
      & MSQK, MG, MG2, SM, DM, DAB, QPE, SGN, PF, SQPE, EMSC2,
       DOUBLE PRECISION HWRGEN, HWUALF, HWUAEM, EPS, HCS, RCS, DIST,
      & ML(6), ML2(6), MR(6), MR2(6), MCH(2), MCH2(2), MNU(4), MNU2(4),
      & MSQK, MG, MG2, SM, DM, DAB, QPE, SGN, PF, SQPE, EMSC2,


@@ -37884,7 +37884,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SUSY 2 PARTON -> 2 SLEPTON PROCESSES
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SUSY 2 PARTON -> 2 SLEPTON PROCESSES
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN, HWUAEM, EPS, HCS, RCS, DIST, S, PF, QPE,
      & FACTR, SN2TH, MZ, MW, ME2(2,2,6,2), ME2W(2,3), EMSC2, GW2
       INTEGER IQ, IQ1, IQ2, ID1, ID2, IL, IL1, IL2, I, J
       DOUBLE PRECISION HWRGEN, HWUAEM, EPS, HCS, RCS, DIST, S, PF, QPE,
      & FACTR, SN2TH, MZ, MW, ME2(2,2,6,2), ME2W(2,3), EMSC2, GW2
       INTEGER IQ, IQ1, IQ2, ID1, ID2, IL, IL1, IL2, I, J


@@ -38102,7 +38102,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SUSY HARD 2 PARTON -> 2 SPARTON PROCESSES
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SUSY HARD 2 PARTON -> 2 SPARTON PROCESSES
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN, HWUALF, EPS, HCS, RCS, DIST, NC, NC2,
      & NC2C, ML2(6), ML4(6), MR2(6), MR4(6), MG2, SM, DM, QPE,
      & SQPE, FACTR, AFAC, AF, BONE, CFAC, CFC2, CFC3, CONE,
       DOUBLE PRECISION HWRGEN, HWUALF, EPS, HCS, RCS, DIST, NC, NC2,
      & NC2C, ML2(6), ML4(6), MR2(6), MR4(6), MG2, SM, DM, QPE,
      & SQPE, FACTR, AFAC, AF, BONE, CFAC, CFC2, CFC3, CONE,


@@ -39013,7 +39013,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SUSY HARD 2 PARTON -> 2 SPARTON/GAUGINO/SLEPTON PROCESSES
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SUSY HARD 2 PARTON -> 2 SPARTON/GAUGINO/SLEPTON PROCESSES
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION SAVWT(3),RANWT,HWRGEN,HWRUNI,Z1,Z2,ET,EJ,
      & QPE,S,T,U,KK,KK2,YJ1INF,YJ1SUP,YJ2INF,YJ2SUP,SVEMSC
       INTEGER ISP
       DOUBLE PRECISION SAVWT(3),RANWT,HWRGEN,HWRUNI,Z1,Z2,ET,EJ,
      & QPE,S,T,U,KK,KK2,YJ1INF,YJ1SUP,YJ2INF,YJ2SUP,SVEMSC
       INTEGER ISP


@@ -39090,7 +39090,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     IDENTIFIES HARD SUSY SUBPROCESS
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     IDENTIFIES HARD SUSY SUBPROCESS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER ID3, R3, ID4, R4, IPERM, IHPR, SSL
       PARAMETER (SSL = 400)
        IHPRO  = 3000 + IHPR
       INTEGER ID3, R3, ID4, R4, IPERM, IHPR, SSL
       PARAMETER (SSL = 400)
        IHPRO  = 3000 + IHPR


@@ -39111,7 +39111,7 @@ C   V + 1 JET PRODUCTION, WHERE V=W (IHPRO.LT.5) OR Z (IHPRO.GE.5).
 C   USES CROSS-SECTIONS OF EHLQ FOR ANNIHILATION AND COMPTON SCATTERING
 C   IHPRO=0 FOR BOTH, 1 FOR ANNIHILATION, AND 2 FOR COMPTON.
 C-----------------------------------------------------------------------
 C   USES CROSS-SECTIONS OF EHLQ FOR ANNIHILATION AND COMPTON SCATTERING
 C   IHPRO=0 FOR BOTH, 1 FOR ANNIHILATION, AND 2 FOR COMPTON.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,DISFAC(2,12,2),EMV2,DISMAX,S,T,U,
      & SHAT,THAT,UHAT,Z,HWUALF,PT,EMT,GFACTR,SIGANN,SIGCOM(2),CSFAC,ET,
      & EJ,YMIN,YMAX,VYMIN,VYMAX,EMAX,CV,CA,BR,EMV,GAMV,HWUAEM,TMIN,TMAX
       DOUBLE PRECISION HWRGEN,HWRUNI,DISFAC(2,12,2),EMV2,DISMAX,S,T,U,
      & SHAT,THAT,UHAT,Z,HWUALF,PT,EMT,GFACTR,SIGANN,SIGCOM(2),CSFAC,ET,
      & EJ,YMIN,YMAX,VYMIN,VYMAX,EMAX,CV,CA,BR,EMV,GAMV,HWUAEM,TMIN,TMAX


@@ -39308,7 +39308,7 @@ C     Vector Boson production with two hard jets
 C     Master subroutine for all vector boson + 2 jet processes
 C     Currently implemented qqbar Z only
 C-----------------------------------------------------------------------
 C     Master subroutine for all vector boson + 2 jet processes
 C     Currently implemented qqbar Z only
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER I,J,K,IDBS,IPRC,IDP(6),ORD,IB,ICMF,IHEP,IFLOW,IZ,IBRAD,
      &     ICOL(5),IDZ,IQ
       DOUBLE PRECISION HWRGEN,HWRUNI,XMASS,PLAB,PRW,PCM,HWUAEM,BR,FLUX,
       INTEGER I,J,K,IDBS,IPRC,IDP(6),ORD,IB,ICMF,IHEP,IFLOW,IZ,IBRAD,
      &     ICOL(5),IDZ,IQ
       DOUBLE PRECISION HWRGEN,HWRUNI,XMASS,PLAB,PRW,PCM,HWUAEM,BR,FLUX,


@@ -39596,7 +39596,7 @@ C     UNLESS USER SPECIFIES OTHERWISE BY MOD(IPROC,100)=1-8 RESPECTIVELY
 C---DSDCOS HOLDS THE CROSS-SECTIONS FOR THE PROCESSES LISTED ABOVE
 C   (1-8) ARE WITH B FROM BEAM 1, (9-16) ARE WITH B FROM BEAM 2.
 C-----------------------------------------------------------------------
 C---DSDCOS HOLDS THE CROSS-SECTIONS FOR THE PROCESSES LISTED ABOVE
 C   (1-8) ARE WITH B FROM BEAM 1, (9-16) ARE WITH B FROM BEAM 2.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,DSDCOS(16),EMT2,EMT,EMW2,EMW,
      & CMFMIN,TAUMIN,TAUMLN,S,T,U,ROOTS,DSMAX
       INTEGER HWRINT,IDHWEX(2,16),I
       DOUBLE PRECISION HWRGEN,HWRUNI,DSDCOS(16),EMT2,EMT,EMW2,EMW,
      & CMFMIN,TAUMIN,TAUMLN,S,T,U,ROOTS,DSMAX
       INTEGER HWRINT,IDHWEX(2,16),I


@@ -39677,7 +39677,7 @@ C-----------------------------------------------------------------------
 C     W+/- PRODUCTION AND DECAY VIA DRELL-YAN PROCESS
 C     MEAN EVWGT IS SIG(W+/-)*(BRANCHING FRACTION) IN NB
 C-----------------------------------------------------------------------
 C     W+/- PRODUCTION AND DECAY VIA DRELL-YAN PROCESS
 C     MEAN EVWGT IS SIG(W+/-)*(BRANCHING FRACTION) IN NB
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRUNI,HWUPCM,PRAN,PROB,COEF,CSFAC,EMW,
      & FTQK,PTOP,ETOP,EBOT,PMAX,FHAD,FTOT,BRAF,FLEP,TMIN,HWUAEM,TMAX
       INTEGER HWRINT,ICH,IC,IL,ID,IDEC,JDEC,IWP(2,16)
       DOUBLE PRECISION HWRGEN,HWRUNI,HWUPCM,PRAN,PROB,COEF,CSFAC,EMW,
      & FTQK,PTOP,ETOP,EBOT,PMAX,FHAD,FTOT,BRAF,FLEP,TMIN,HWUAEM,TMAX
       INTEGER HWRINT,ICH,IC,IL,ID,IDEC,JDEC,IWP(2,16)


@@ -39818,7 +39818,7 @@ CDECK  ID>, HWICHK.
 C-----------------------------------------------------------------------
       SUBROUTINE HWICHK
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
       SUBROUTINE HWICHK
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       IF(RMASS(1).LT.0.1.OR.RMASS(1).GT.1.0.OR.
      &   FMRS(1,1,20,1).LT.0.1.OR.FMRS(1,1,20,1).GT.1.0) THEN
         STOP 'Block data hwudat not loaded, stop execution'
       IF(RMASS(1).LT.0.1.OR.RMASS(1).GT.1.0.OR.
      &   FMRS(1,1,20,1).LT.0.1.OR.FMRS(1,1,20,1).GT.1.0) THEN
         STOP 'Block data hwudat not loaded, stop execution'


@@ -39837,7 +39837,7 @@ C     =1 PDG; =2 HERWIG numeric; =3 HERWIG character name.
 C     When reading in if IME =1 matrix element codes >= 100 are accepted
 C                             0                            are set zero.
 C-----------------------------------------------------------------------
 C     When reading in if IME =1 matrix element codes >= 100 are accepted
 C                             0                            are set zero.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IUNIT,IOPT,IME,JUNIT,I,J,K,L,IDKY,ITMP(5),IDUM
       CHARACTER*8 CDK(NMXDKS),CDKPRD(5,NMXDKS),CDUM
       JUNIT=ABS(IUNIT)
       INTEGER IUNIT,IOPT,IME,JUNIT,I,J,K,L,IDKY,ITMP(5),IDUM
       CHARACTER*8 CDK(NMXDKS),CDKPRD(5,NMXDKS),CDUM
       JUNIT=ABS(IUNIT)


@@ -39918,7 +39918,7 @@ C----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SETS INPUT PARAMETERS
 C----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     SETS INPUT PARAMETERS
 C----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION FAC,ANGLE
       INTEGER I,J,N,L
       CHARACTER*28 TITLE
       DOUBLE PRECISION FAC,ANGLE
       INTEGER I,J,N,L
       CHARACTER*28 TITLE


@@ -40533,7 +40533,7 @@ C----------------------------------------------------------------------
 C     Use the GUPI (Generic User Process Interface) run common block
 C     to initialise HERWIG -- Initialization for Les Houches interface
 C----------------------------------------------------------------------
 C     Use the GUPI (Generic User Process Interface) run common block
 C     to initialise HERWIG -- Initialization for Les Houches interface
 C----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER MAXPUP
       PARAMETER(MAXPUP=100)
       INTEGER IDBMUP,PDFGUP,PDFSUP,IDWTUP,NPRUP,LPRUP
       INTEGER MAXPUP
       PARAMETER(MAXPUP=100)
       INTEGER IDBMUP,PDFGUP,PDFSUP,IDWTUP,NPRUP,LPRUP


@@ -40633,7 +40633,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Subroutine to merge Higgs WW/ZZ decay modes for four body ME
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Subroutine to merge Higgs WW/ZZ decay modes for four body ME
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IH,I,NMODE,J,K
       LOGICAL REMOVE
       DOUBLE PRECISION BR
       INTEGER IH,I,NMODE,J,K
       LOGICAL REMOVE
       DOUBLE PRECISION BR


@@ -40743,7 +40743,7 @@ C     IOPT = 1 sets the weights for the different channels to their
 C              default values
 C     IOPT = 2 optimises the weights for the process selected
 C-----------------------------------------------------------------------
 C              default values
 C     IOPT = 2 optimises the weights for the process selected
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER I,IPRC,ICH,IOPT,ISTP,IWGT,IFER,IANT,IGAU,IQRK
       LOGICAL CALLED,TEV,LHC
       DOUBLE PRECISION CHNPST(IMAXCH,IMAXOP),D(IMAXOP),CHWGTS(IMAXCH),
       INTEGER I,IPRC,ICH,IOPT,ISTP,IWGT,IFER,IANT,IGAU,IQRK
       LOGICAL CALLED,TEV,LHC
       DOUBLE PRECISION CHNPST(IMAXCH,IMAXOP),D(IMAXOP),CHWGTS(IMAXCH),


@@ -41012,7 +41012,7 @@ C-----------------------------------------------------------------------
 C     Calculates the couplings for the SUSY decays for spin correlations
 C     and 3/4 body matrix elements
 C-----------------------------------------------------------------------
 C     Calculates the couplings for the SUSY decays for spin correlations
 C     and 3/4 body matrix elements
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWUALF,PRE,MCHAR(2),QIJPP(4,4),SIJPP(4,4),
      &     DIJ(2,2),QIJ(2,2),R(4,2),SIJ(2,2)
       INTEGER I,J,K,L,IH,IK,IL,IQ
       DOUBLE PRECISION HWUALF,PRE,MCHAR(2),QIJPP(4,4),SIJPP(4,4),
      &     DIJ(2,2),QIJ(2,2),R(4,2),SIJ(2,2)
       INTEGER I,J,K,L,IH,IK,IL,IQ


@@ -41232,7 +41232,7 @@ C-----------------------------------------------------------------------
 C     Initialise all the decay modes for three/four body MEs and spin
 C     correlations
 C-----------------------------------------------------------------------
 C     Initialise all the decay modes for three/four body MEs and spin
 C     correlations
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER I,J,K,NDKYST
 C--set the number of two and three body modes to zero
       N2MODE = 0
       INTEGER I,J,K,NDKYST
 C--set the number of two and three body modes to zero
       N2MODE = 0


@@ -41355,7 +41355,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Initialise the SUSY two body modes for spin correlations
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Initialise the SUSY two body modes for spin correlations
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER I,J,IL,IH,L,L1,IM,O(2),II,JJ,III,JJJ,KKK
       COMMON /HWSPNC/ SW,CW,TW,E,G,RT,ORT,MW,MZ,GS,AFN,AFG,AFC,OIJ,OIJP,
      &               OIJPP,HNN,HCC,HNC,HFF,HWW,HZZ,ZAB,HHB
       INTEGER I,J,IL,IH,L,L1,IM,O(2),II,JJ,III,JJJ,KKK
       COMMON /HWSPNC/ SW,CW,TW,E,G,RT,ORT,MW,MZ,GS,AFN,AFG,AFC,OIJ,OIJP,
      &               OIJPP,HNN,HCC,HNC,HFF,HWW,HZZ,ZAB,HHB


@@ -42569,7 +42569,7 @@ C-----------------------------------------------------------------------
 C     Initialise the top/SUSY three body decay modes
 C     gravitino and RPV modes added by Peter Richardson
 C-----------------------------------------------------------------------
 C     Initialise the top/SUSY three body decay modes
 C     gravitino and RPV modes added by Peter Richardson
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER I,J,K,L,L1,IL,IQ,IQ1,IQ2,IFR,SIFR,IH,IH1,IM,O(2),II,JJ,
      &     III,JJJ,KKK
       DOUBLE PRECISION SW,CW,TW,E,G,RT,ORT,MW,MZ,GS,AFN(2,12,2,4),
       INTEGER I,J,K,L,L1,IL,IQ,IQ1,IQ2,IFR,SIFR,IH,IH1,IM,O(2),II,JJ,
      &     III,JJJ,KKK
       DOUBLE PRECISION SW,CW,TW,E,G,RT,ORT,MW,MZ,GS,AFN(2,12,2,4),


@@ -44639,7 +44639,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Initialise the Higgs four body modes
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Initialise the Higgs four body modes
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER I,J,K,IL,IH,II,JJ
       DOUBLE PRECISION COL(2),SW,CW,TW,E,G,RT,ORT,MW,MZ,AFN(2,12,2,4),
      &     AFG(2,6,2),AFC(2,12,2,2),OIJ(2,4,2),OIJP(2,2,2),OIJPP(2,4,4),
       INTEGER I,J,K,IL,IH,II,JJ
       DOUBLE PRECISION COL(2),SW,CW,TW,E,G,RT,ORT,MW,MZ,AFN(2,12,2,4),
      &     AFG(2,6,2),AFC(2,12,2,2),OIJ(2,4,2),OIJP(2,2,2),OIJPP(2,4,4),


@@ -44741,7 +44741,7 @@ C-----------------------------------------------------------------------
 C  Reads in SUSY particle properties and decays,
 C  in format generated by ISAWIG
 C-----------------------------------------------------------------------
 C  Reads in SUSY particle properties and decays,
 C  in format generated by ISAWIG
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER I,J,K,IH,IHW,NSSP,NDEC,MDKYS
       DOUBLE PRECISION BETAH, WEINCOS,WEINSIN, MW,MZ, RMMAX
       DOUBLE PRECISION FTM,FTMUU(4),FTMDD(4),FTMTT(4),FTMBB(4),FTMU,FTMD
       INTEGER I,J,K,IH,IHW,NSSP,NDEC,MDKYS
       DOUBLE PRECISION BETAH, WEINCOS,WEINSIN, MW,MZ, RMMAX
       DOUBLE PRECISION FTM,FTMUU(4),FTMDD(4),FTMTT(4),FTMBB(4),FTMU,FTMD


@@ -44989,7 +44989,7 @@ C     IPROC = 1000,... ADDS SOFT UNDERLYING EVENT
 C           = 8000:  CREATES MINIMUM-BIAS EVENT
 C     SUPPRESSED BY ADDING 10000 TO IPROC
 C-----------------------------------------------------------------------
 C           = 8000:  CREATES MINIMUM-BIAS EVENT
 C     SUPPRESSED BY ADDING 10000 TO IPROC
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWREXP,ENFAC,TECM,SECM,SUMM,EMCL,BMP(5),BMR(3,3)
       INTEGER HWRINT,NETC,IBT,IDBT,ID1,ID2,ID3,KHEP,LHEP,NTRY,ICMS,
      & NPPBAR,MCHT,JCL,JD1,JD2,JD3,ICH,MODC,NCHT,INHEP(2),
       DOUBLE PRECISION HWREXP,ENFAC,TECM,SECM,SUMM,EMCL,BMP(5),BMR(3,3)
       INTEGER HWRINT,NETC,IBT,IDBT,ID1,ID2,ID3,KHEP,LHEP,NTRY,ICMS,
      & NPPBAR,MCHT,JCL,JD1,JD2,JD3,ICH,MODC,NCHT,INHEP(2),


@@ -45274,7 +45274,7 @@ C-----------------------------------------------------------------------
 C     GENERATES CYLINDRICAL PHASE SPACE USING THE METHOD OF JADACH
 C     RETURNS WITH NCL=0 IF UNSUCCESSFUL
 C-----------------------------------------------------------------------
 C     GENERATES CYLINDRICAL PHASE SPACE USING THE METHOD OF JADACH
 C     RETURNS WITH NCL=0 IF UNSUCCESSFUL
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWREXT,HWRUNG,HWUSQR,TECM,ESS,ALOGS,EPS,SUMX,
      & SUMY,PT,PX,PY,PT2,SUMPT2,SUMTM,XIMIN,XIMAX,YY,SUM1,SUM2,SUM3,
      & SUM4,EX,FY,DD,DYY,ZZ,E1,TM,SLOP,XI(NMXCL)
       DOUBLE PRECISION HWREXT,HWRUNG,HWUSQR,TECM,ESS,ALOGS,EPS,SUMX,
      & SUMY,PT,PX,PY,PT2,SUMPT2,SUMTM,XIMIN,XIMAX,YY,SUM1,SUM2,SUM3,
      & SUM4,EX,FY,DD,DYY,ZZ,E1,TM,SLOP,XI(NMXCL)


@@ -45418,7 +45418,7 @@ C     adjusts pointers,  sets CMMOM (for two-body mode) and resets RSTAB
 C     if necessary.  The branching ratios of any other IDKTMP decays are
 C     scaled by (1.-BRTMP)/(1.-BR_OLD)
 C-----------------------------------------------------------------------
 C     if necessary.  The branching ratios of any other IDKTMP decays are
 C     scaled by (1.-BRTMP)/(1.-BR_OLD)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWUPCM,BRTMP,SCALE,EPS
       INTEGER IDKTMP,IMETMP,IATMP,IBTMP,ICTMP,IDTMP,IETMP,IDKY,ITMP(5),
      & L,I,J,K,JPREV
       DOUBLE PRECISION HWUPCM,BRTMP,SCALE,EPS
       INTEGER IDKTMP,IMETMP,IATMP,IBTMP,ICTMP,IDTMP,IETMP,IDKY,ITMP(5),
      & L,I,J,K,JPREV


@@ -45574,7 +45574,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Chooses charged multiplicity NCHT at the p-pbar c.m. energy EPPBAR
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Chooses charged multiplicity NCHT at the p-pbar c.m. energy EPPBAR
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWMNBI,HWRGEN,EPPBAR,E0,ALOGS,RK,EK,AVN,SUM,R,
      & CUM(500)
       INTEGER NCHT,IMAX,I,N
       DOUBLE PRECISION HWMNBI,HWRGEN,EPPBAR,E0,ALOGS,RK,EK,AVN,SUM,R,
      & CUM(500)
       INTEGER NCHT,IMAX,I,N


@@ -45631,7 +45631,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C COMPUTES WEIGHT FOR MINIMUM-BIAS EVENT
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C COMPUTES WEIGHT FOR MINIMUM-BIAS EVENT
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION S,X,Y
       INTEGER IDB,IDT,IDBT
       IF (IERROR.NE.0) RETURN
       DOUBLE PRECISION S,X,Y
       INTEGER IDB,IDT,IDBT
       IF (IERROR.NE.0) RETURN


@@ -45693,7 +45693,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     subroutine for radiation in top decays
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     subroutine for radiation in top decays
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IHEP,KK,IPOS,NN,NHEP0,KK1,KK2,JMOH(NMXHEP)
       DOUBLE PRECISION HWDPWT
       EXTERNAL HWDPWT
       INTEGER IHEP,KK,IPOS,NN,NHEP0,KK1,KK2,JMOH(NMXHEP)
       DOUBLE PRECISION HWDPWT
       EXTERNAL HWDPWT


@@ -45744,7 +45744,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     subroutine for radiation in top production
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     subroutine for radiation in top production
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

 C--local variables
       INTEGER IMO(10),IFOUND,JMO(2),I,J,K,L,NSTART,NHEPX
 C--initialisation
 C--local variables
       INTEGER IMO(10),IFOUND,JMO(2),I,J,K,L,NSTART,NHEPX
 C--initialisation


@@ -45893,7 +45893,7 @@ C-----------------------------------------------------------------------
 C     Gaussian random number, mean A, standard deviation B.
 C     Generates uncorrelated pairs and throws one of them away.
 C-----------------------------------------------------------------------
 C     Gaussian random number, mean A, standard deviation B.
 C     Generates uncorrelated pairs and throws one of them away.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGAU,HWRGEN,A,B,X,TRASH
       INTEGER J
       EXTERNAL HWRGEN
       DOUBLE PRECISION HWRGAU,HWRGEN,A,B,X,TRASH
       INTEGER J
       EXTERNAL HWRGEN


@@ -46014,7 +46014,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Generates a random primary IP using a triple Gaussian distribution
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Generates a random primary IP using a triple Gaussian distribution
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGAU
       INTEGER I
       EXTERNAL HWRGAU
       DOUBLE PRECISION HWRGAU
       INTEGER I
       EXTERNAL HWRGAU


@@ -46101,7 +46101,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     DOES BRANCHING OF SPACELIKE PARTON KPAR
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     DOES BRANCHING OF SPACELIKE PARTON KPAR
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWBVMC,HWRGEN,HWRUNI,HWSTAB,HWUALF,HWUTAB,HWSGQQ,
      & HWSSUD,XLAST,QNOW,QLST,QP,QMIN,QLAM,QSAV,SMAX,SLST,SNOW,RN,SUDA,
      & SUDB,ZZ,ENOW,XI,PMOM,DIST(13),DMIN,X1,X2,REJFAC,OTHXI,OTHZ,QTMP,
       DOUBLE PRECISION HWBVMC,HWRGEN,HWRUNI,HWSTAB,HWUALF,HWUTAB,HWSGQQ,
      & HWSSUD,XLAST,QNOW,QLST,QP,QMIN,QLAM,QSAV,SMAX,SLST,SNOW,RN,SUDA,
      & SUDB,ZZ,ENOW,XI,PMOM,DIST(13),DMIN,X1,X2,REJFAC,OTHXI,OTHZ,QTMP,


@@ -46665,7 +46665,7 @@ C
 C     ID1.LT.0 ON RETURN MEANS NO PHASE SPACE
 C     ID1.EQ.0 ON RETURN FLAGS REJECTED BRANCHINGS
 C-----------------------------------------------------------------------
 C     ID1.LT.0 ON RETURN MEANS NO PHASE SPACE
 C     ID1.EQ.0 ON RETURN FLAGS REJECTED BRANCHINGS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWBVMC,HWRGEN,HWUALF,HWUAEM,QP,X,QQ,Z,WQG,WQV,
      & WQP,XQV,ZMIN,ZMAX,YMIN,YMAX,DELY,YY,PSUM,EZ,WQN,WR,ZR,WZ,ZZ,AZ,
      & PVAL,EY,DIST(13),PROB(13,100),PPHO
       DOUBLE PRECISION HWBVMC,HWRGEN,HWUALF,HWUAEM,QP,X,QQ,Z,WQG,WQV,
      & WQP,XQV,ZMIN,ZMAX,YMIN,YMAX,DELY,YY,PSUM,EZ,WQN,WR,ZR,WZ,ZZ,AZ,
      & PVAL,EY,DIST(13),PROB(13,100),PPHO


@@ -46874,7 +46874,7 @@ C   CERN_ROOT:[DOC]PDFLIB.TXT ON VXCERN
 C-----------------------------------------------------------------------
 C---BRW change 27/8/04: include Frixione's fix to reduce PDFSET calls
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C---BRW change 27/8/04: include Frixione's fix to reduce PDFSET calls
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWSGAM,X,SCALE,XOLD,QOLD,XMWN,QSCA,SS,SMIN,S,T,
      & TMIN,TMAX,VX,AA,VT,WT,UPV,DNV,SEA,STR,CHM,BTM,TOP,GLU,WX,XQSUM,
      & DMIN,TPMIN,TPMAX,DIST(13),G(2),Q0(5),QL(5),F(5),A(6,5),
       DOUBLE PRECISION HWSGAM,X,SCALE,XOLD,QOLD,XMWN,QSCA,SS,SMIN,S,T,
      & TMIN,TMAX,VX,AA,VT,WT,UPV,DNV,SEA,STR,CHM,BTM,TOP,GLU,WX,XQSUM,
      & DMIN,TPMIN,TPMAX,DIST(13),G(2),Q0(5),QL(5),F(5),A(6,5),


@@ -47262,7 +47262,7 @@ C---FIX TO CALL SCHULER-SJOSTRAND CODE
             IF(PARM(1).NE.PARMSAVE.OR.VAL(1).NE.VALSAVE)THEN
               PARMSAVE=PARM(1)
               VALSAVE=VAL(1)
             IF(PARM(1).NE.PARMSAVE.OR.VAL(1).NE.VALSAVE)THEN
               PARMSAVE=PARM(1)
               VALSAVE=VAL(1)

-              CALL PDFSET(PARM,VAL)
+              CALL PDFSET_HERWIG(PARM,VAL)

             ENDIF
             IF (X.LT.PDFXMN.AND.PDFWRX(IBEAM,1) .OR.
      &          X.GT.PDFXMX.AND.PDFWRX(IBEAM,2)) THEN
             ENDIF
             IF (X.LT.PDFXMN.AND.PDFWRX(IBEAM,1) .OR.
      &          X.GT.PDFXMX.AND.PDFWRX(IBEAM,2)) THEN


@@ -47362,7 +47362,7 @@ C---USE PDFLIB NUCLEON STRUCTURE FUNCTIONS
         IF(PARM(1).NE.PARMSAVE.OR.VAL(1).NE.VALSAVE)THEN
           PARMSAVE=PARM(1)
           VALSAVE=VAL(1)
         IF(PARM(1).NE.PARMSAVE.OR.VAL(1).NE.VALSAVE)THEN
           PARMSAVE=PARM(1)
           VALSAVE=VAL(1)

-          CALL PDFSET(PARM,VAL)
+          CALL PDFSET_HERWIG(PARM,VAL)

         ENDIF
         IF (X.LT.PDFXMN.AND.PDFWRX(IBEAM,1) .OR.
      &      X.GT.PDFXMX.AND.PDFWRX(IBEAM,2)) THEN
         ENDIF
         IF (X.LT.PDFXMN.AND.PDFWRX(IBEAM,1) .OR.
      &      X.GT.PDFXMX.AND.PDFWRX(IBEAM,2)) THEN


@@ -47667,7 +47667,7 @@ C-----------------------------------------------------------------------
 C     GENERATES X VALUES (IF GENEX)
 C     EVALUATES STRUCTURE FUNCTIONS AND ENFORCES CUTOFFS ON X
 C-----------------------------------------------------------------------
 C     GENERATES X VALUES (IF GENEX)
 C     EVALUATES STRUCTURE FUNCTIONS AND ENFORCES CUTOFFS ON X
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWBVMC,HWRUNI,X,QL
       INTEGER I,J
       LOGICAL GENEX
       DOUBLE PRECISION HWBVMC,HWRUNI,X,QL
       INTEGER I,J
       LOGICAL GENEX


@@ -47695,7 +47695,7 @@ C-----------------------------------------------------------------------
 C     CORRECTION TO GLUON STRUCTURE FUNCTION FOR BACKWARD EVOLUTION:
 C     G->Q-QBAR PART OF FORM FACTOR
 C-----------------------------------------------------------------------
 C     CORRECTION TO GLUON STRUCTURE FUNCTION FOR BACKWARD EVOLUTION:
 C     G->Q-QBAR PART OF FORM FACTOR
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWSGQQ,HWUALF,QSCA,GG
       EXTERNAL HWUALF
       GG=HWUALF(1,QSCA)**(-ONE/BETAF)
       DOUBLE PRECISION HWSGQQ,HWUALF,QSCA,GG
       EXTERNAL HWUALF
       GG=HWUALF(1,QSCA)**(-ONE/BETAF)


@@ -47718,7 +47718,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     MRST98 Leading order PDF's central and higher gluon + average
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     MRST98 Leading order PDF's central and higher gluon + average
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION X,Q,UPV,DNV,USEA,DSEA,STR,CHM,BOT,GLU,XMIN,XMAX,
      &     QSQMIN,QSQMAX,Q2,QQ(NQMRS),XXMRS(NXMRS),G(NPMRS),N0(NPMRS),
      &     XSAVE,Q2SAVE,XXX,A,B,FAC
       DOUBLE PRECISION X,Q,UPV,DNV,USEA,DSEA,STR,CHM,BOT,GLU,XMIN,XMAX,
      &     QSQMIN,QSQMAX,Q2,QQ(NQMRS),XXMRS(NXMRS),G(NPMRS),N0(NPMRS),
      &     XSAVE,Q2SAVE,XXX,A,B,FAC


@@ -47840,7 +47840,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     REPLACES SPACELIKE PARTONS BY SPECTATORS
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     REPLACES SPACELIKE PARTONS BY SPECTATORS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWUSQR,EMSQ,EMTR,EPAR,XPAR,QSQ,PCL(5)
       INTEGER KHEP,IP,JP,IDH,IDP,ISP,IDSPC,JHEP
       EXTERNAL HWUSQR
       DOUBLE PRECISION HWUSQR,EMSQ,EMTR,EPAR,XPAR,QSQ,PCL(5)
       INTEGER KHEP,IP,JP,IDH,IDP,ISP,IDSPC,JHEP
       EXTERNAL HWUSQR


@@ -47947,7 +47947,7 @@ CDECK  ID>, HWSSUD.
 C-----------------------------------------------------------------------
       FUNCTION HWSSUD(I)
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
       FUNCTION HWSSUD(I)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWSSUD,HWSGQQ,DMIN,QSCA,XLAST,DIST(13)
       INTEGER I,N0,IS,ID
       EXTERNAL HWSGQQ
       DOUBLE PRECISION HWSSUD,HWSGQQ,DMIN,QSCA,XLAST,DIST(13)
       INTEGER I,N0,IS,ID
       EXTERNAL HWSGQQ


@@ -48038,7 +48038,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     TRUE FOR VALENCE PARTON ID IN INCOMING HADRON INHAD
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     TRUE FOR VALENCE PARTON ID IN INCOMING HADRON INHAD
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER ID,IDHAD
       LOGICAL HWSVAL
       HWSVAL=.FALSE.
       INTEGER ID,IDHAD
       LOGICAL HWSVAL
       HWSVAL=.FALSE.


@@ -48069,7 +48069,7 @@ C     Running electromagnetic coupling constant.
 C     See R. Kleiss et al.: CERN yellow report 89-08, vol.3 p.129
 C     Hadronic component from: H. Burkhardt et al.: Z. Phys C43 (89) 497
 C-----------------------------------------------------------------------
 C     See R. Kleiss et al.: CERN yellow report 89-08, vol.3 p.129
 C     Hadronic component from: H. Burkhardt et al.: Z. Phys C43 (89) 497
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWUAEM,HWUAER,Q2,EPS,A1,B1,C1,A2,B2,C2,A3,B3,C3,
      & A4,B4,C4,AEMPI,EEL2,EMU2,ETAU2,ETOP2,REPIGG,X
       LOGICAL FIRST
       DOUBLE PRECISION HWUAEM,HWUAER,Q2,EPS,A1,B1,C1,A2,B2,C2,A3,B3,C3,
      & A4,B4,C4,AEMPI,EEL2,EMU2,ETAU2,ETOP2,REPIGG,X
       LOGICAL FIRST


@@ -48152,7 +48152,7 @@ C         .EQ.2  RATIO OF ABOVE TO ONE-LOOP
 C                WITH 5-FLAVOUR BETA, LAMBDA=QCDL3
 C         .EQ.3  ONE-LOOP WITH 5-FLAVOUR BETA, LAMBDA=QCDL3
 C-----------------------------------------------------------------------
 C                WITH 5-FLAVOUR BETA, LAMBDA=QCDL3
 C         .EQ.3  ONE-LOOP WITH 5-FLAVOUR BETA, LAMBDA=QCDL3
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWUALF,SCALE,KAFAC,B3,B4,B5,B6,C3,C4,C5,C6,C35,
      & C45,C65,D35,RHO,RAT,RLF,DRH,EPS
       INTEGER IOPT,ITN
       DOUBLE PRECISION HWUALF,SCALE,KAFAC,B3,B4,B5,B6,C3,C4,C5,C6,C35,
      & C45,C65,D35,RHO,RAT,RLF,DRH,EPS
       INTEGER IOPT,ITN


@@ -48244,7 +48244,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Returns the antiparticle of IPART; uses HERWIG numbering
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Returns the antiparticle of IPART; uses HERWIG numbering
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER HWUANT,IPART,IPDG,IANTI,OLDERR
       CHARACTER*8 CDUM
       OLDERR=IERROR
       INTEGER HWUANT,IPART,IPDG,IANTI,OLDERR
       CHARACTER*8 CDUM
       OLDERR=IERROR


@@ -48282,7 +48282,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Replaces all &'s in TXNAME by backslashes
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     Replaces all &'s in TXNAME by backslashes
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER I,J,L
       CHARACTER*1 Z
       Z=CHAR(92)
       INTEGER I,J,L
       CHARACTER*1 Z
       Z=CHAR(92)


@@ -48301,7 +48301,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     PRINTS OUT DATA ON PARTON SHOWER
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     PRINTS OUT DATA ON PARTON SHOWER
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER I,J
       IF (PRVTX) THEN
         WRITE(6,10) INHAD,XFACT
       INTEGER I,J
       IF (PRVTX) THEN
         WRITE(6,10) INHAD,XFACT


@@ -48333,7 +48333,7 @@ C-----------------------------------------------------------------------
 C     BOOST THE ENTIRE EVENT RECORD TO (IOPT=1) OR FROM (IOPT=0) ITS
 C     CENTRE-OF-MASS FRAME, WITH INCOMING HADRONS ON Z-AXIS
 C-----------------------------------------------------------------------
 C     BOOST THE ENTIRE EVENT RECORD TO (IOPT=1) OR FROM (IOPT=0) ITS
 C     CENTRE-OF-MASS FRAME, WITH INCOMING HADRONS ON Z-AXIS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION PBOOST(5),RBOOST(3,3)
       INTEGER IOPT,IHEP,BOOSTD,IHAD
       SAVE BOOSTD,PBOOST,RBOOST
       DOUBLE PRECISION PBOOST(5),RBOOST(3,3)
       INTEGER IOPT,IHEP,BOOSTD,IHAD
       SAVE BOOSTD,PBOOST,RBOOST


@@ -48386,7 +48386,7 @@ C     labels final fermion; type given as:
 C        I,J= 1- 6: d,u,s,c,b,t
 C           =11-16: e,nu_e,mu,nu_mu,tau,nu_tau
 C-----------------------------------------------------------------------
 C        I,J= 1- 6: d,u,s,c,b,t
 C           =11-16: e,nu_e,mu,nu_mu,tau,nu_tau
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION QSQ,CLF(7),POL1,POL2,QIF,VI,AI,VF,AF,PG,DQM,PMW,
      & DEN,XRE,XIM,XSQ,VI2,AI2,VF2,AF2,PG2,PG12,DQM2,PMW2,DEN2,XRE2,
      & XIM2,XSQ2,XRE12,XIM12
       DOUBLE PRECISION QSQ,CLF(7),POL1,POL2,QIF,VI,AI,VF,AF,PG,DQM,PMW,
      & DEN,XRE,XIM,XSQ,VI2,AI2,VF2,AF2,PG2,PG12,DQM2,PMW2,DEN2,XRE2,
      & XIM2,XSQ2,XRE12,XIM12


@@ -48501,7 +48501,7 @@ C       QORQQB(I) = .TRUE. if it is a quark or an antidiquark
 C       QBORQQ(I) = .TRUE. if it is an antiquark or a diquark
 C     And stores the particle decay tables: call HWUDPR to print them
 C-----------------------------------------------------------------------
 C       QBORQQ(I) = .TRUE. if it is an antiquark or a diquark
 C     And stores the particle decay tables: call HWUDPR to print them
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       COMMON/HWSEED/ISEED(2)
       INTEGER ISEED
       INTEGER NLAST,NNEXT,NLEFT,NREST,I,J,MMWIDE,MMLONG,MMHOFF,MMVOFF
       COMMON/HWSEED/ISEED(2)
       INTEGER ISEED
       INTEGER NLAST,NNEXT,NLEFT,NREST,I,J,MMWIDE,MMLONG,MMHOFF,MMVOFF


@@ -59117,7 +59117,7 @@ C
 C     Modified 16/01/01 by BRW to force particle on mass shell if
 C     p^2-m^2 < 10^-10 GeV^2 (rounding errors)
 C-----------------------------------------------------------------------
 C     Modified 16/01/01 by BRW to force particle on mass shell if
 C     p^2-m^2 < 10^-10 GeV^2 (rounding errors)
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,PMOM(4),DISP(4),PMOM2,SCALE,OFFSH
       INTEGER ID
       EXTERNAL HWRGEN
       DOUBLE PRECISION HWRGEN,PMOM(4),DISP(4),PMOM2,SCALE,OFFSH
       INTEGER ID
       EXTERNAL HWRGEN


@@ -59144,7 +59144,7 @@ C     Gives warnings if: a particle has no decay modes or antiparticle's
 C     modes are not the charge conjugates of the particles.
 C     (N.B. CP violation permits this).
 C-----------------------------------------------------------------------
 C     modes are not the charge conjugates of the particles.
 C     (N.B. CP violation permits this).
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWUPCM,HWUAEM,HWUALF,BRSUM,EPS,SCALE,
      & BRTMP(NMXDKS),FN,X,W,Q,FAC
       INTEGER HWUANT,I,IDKY,LAST,LTMP(NMXMOD),J,L,K,M,N,INDX(NMXMOD),
       DOUBLE PRECISION HWUPCM,HWUAEM,HWUALF,BRSUM,EPS,SCALE,
      & BRTMP(NMXDKS),FN,X,W,Q,FAC
       INTEGER HWUANT,I,IDKY,LAST,LTMP(NMXMOD),J,L,K,M,N,INDX(NMXMOD),


@@ -59462,7 +59462,7 @@ C                 Designed to be printed as landscape
 C     If (PRNWEB) HTML  to the file HW_decays/index.html
 C                                            /PART0000001.html etc.
 C-----------------------------------------------------------------------
 C     If (PRNWEB) HTML  to the file HW_decays/index.html
 C                                            /PART0000001.html etc.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER MMWIDE,MMLONG,MMHOFF,MMVOFF,IUNITT,IUNTW1,IUNTW2,I,NM,J,K,
      & L,M
       CHARACTER*1  Z
       INTEGER MMWIDE,MMLONG,MMHOFF,MMVOFF,IUNITT,IUNTW1,IUNTW2,I,NM,J,K,
      & L,M
       CHARACTER*1  Z


@@ -59800,7 +59800,7 @@ C     EDIT THE EVENT RECORD
 C     IF N>0 DELETE THE N ENTRIES IN IEDT FROM EVENT RECORD
 C     IF N<0 INSERT LINES AFTER THE -N ENTRIES IN IEDT
 C-----------------------------------------------------------------------
 C     IF N>0 DELETE THE N ENTRIES IN IEDT FROM EVENT RECORD
 C     IF N<0 INSERT LINES AFTER THE -N ENTRIES IN IEDT
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER N,IEDT(*),IMAP(0:NMXHEP),IHEP,I,J,I1,I2
       COMMON /HWUMAP/IMAP
 C---MOVE ENTRIES AND CALCULATE MAPPING OF POINTERS
       INTEGER N,IEDT(*),IMAP(0:NMXHEP),IHEP,I,J,I1,I2
       COMMON /HWUMAP/IMAP
 C---MOVE ENTRIES AND CALCULATE MAPPING OF POINTERS


@@ -59865,7 +59865,7 @@ C-----------------------------------------------------------------------
 C     Loads cross-section coefficients, for kinematically open channels,
 C     in llbar-->qqbar; lepton label IL=1-6: e,nu_e,mu,nu_mu,tau,nu_tau.
 C-----------------------------------------------------------------------
 C     Loads cross-section coefficients, for kinematically open channels,
 C     in llbar-->qqbar; lepton label IL=1-6: e,nu_e,mu,nu_mu,tau,nu_tau.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION Q2
       INTEGER IL,JL,IQ
       Q2=EMSCA**2
       DOUBLE PRECISION Q2
       INTEGER IL,JL,IQ
       Q2=EMSCA**2


@@ -59891,7 +59891,7 @@ C-----------------------------------------------------------------------
 C     MOVE A BLOCK OF ENTRIES IN THE EVENT RECORD
 C     N ENTRIES IN HEPEVT STARTING AT IFROM ARE MOVED TO AFTER ITO
 C-----------------------------------------------------------------------
 C     MOVE A BLOCK OF ENTRIES IN THE EVENT RECORD
 C     N ENTRIES IN HEPEVT STARTING AT IFROM ARE MOVED TO AFTER ITO
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER N,IFROM,ITO,IMAP(0:NMXHEP),LFROM,LTO,I,IEDT(NMXHEP),IHEP,
      $ JHEP,KHEP
       COMMON /HWUMAP/IMAP
       INTEGER N,IFROM,ITO,IMAP(0:NMXHEP),LFROM,LTO,I,IEDT(NMXHEP),IHEP,
      $ JHEP,KHEP
       COMMON /HWUMAP/IMAP


@@ -59940,7 +59940,7 @@ C                 Call HWUDPR to create particle property files in
 C                 the subdirectory HW_decays/
 C     ******* gives the event number 0000001 etc.
 C-----------------------------------------------------------------------
 C                 the subdirectory HW_decays/
 C     ******* gives the event number 0000001 etc.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER MMWIDE,MMLONG,MMHOFF,MMVOFF,I,IST,IS,ID,MS,J,K,IUNITW,
      & IUNITT
       CHARACTER*1 Z
       INTEGER MMWIDE,MMLONG,MMHOFF,MMVOFF,I,IST,IS,ID,MS,J,K,IUNITW,
      & IUNITT
       CHARACTER*1 Z


@@ -60373,7 +60373,7 @@ C-----------------------------------------------------------------------
 C     Subroutine to handle termination of HERWIG if reaches end of event
 C     file
 C-----------------------------------------------------------------------
 C     Subroutine to handle termination of HERWIG if reaches end of event
 C     file
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

 C--reset the number of events to the correct value
       NEVHEP = NEVHEP-1
 C--output information on the events
 C--reset the number of events to the correct value
       NEVHEP = NEVHEP-1
 C--output information on the events


@@ -60389,7 +60389,7 @@ C-----------------------------------------------------------------------
 C     FINALISES THE EVENT BY UNDOING THE LORENTZ BOOST IF THERE WAS ONE,
 C     CHECKING FOR ERRORS, AND PRINTING
 C-----------------------------------------------------------------------
 C     FINALISES THE EVENT BY UNDOING THE LORENTZ BOOST IF THERE WAS ONE,
 C     CHECKING FOR ERRORS, AND PRINTING
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IHEP
       LOGICAL CALLED
       COMMON/HWDBUG/CALLED
       INTEGER IHEP
       LOGICAL CALLED
       COMMON/HWDBUG/CALLED


@@ -60491,7 +60491,7 @@ C     IOPT= 1 GIVEN IPDG, RETURNS IWIG AND NWIG
 C     IOPT= 2 GIVEN IWIG, RETURNS IPDG AND NWIG
 C     IOPT= 3 GIVEN NWIG, RETURNS IPDG AND IWIG
 C-----------------------------------------------------------------------
 C     IOPT= 2 GIVEN IWIG, RETURNS IPDG AND NWIG
 C     IOPT= 3 GIVEN NWIG, RETURNS IPDG AND IWIG
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IOPT,IPDG,IWIG,I
       CHARACTER*8 NWIG
       IF (IOPT.EQ.1) THEN
       INTEGER IOPT,IPDG,IWIG,I
       CHARACTER*8 NWIG
       IF (IOPT.EQ.1) THEN


@@ -60549,7 +60549,7 @@ C-----------------------------------------------------------------------
 C     COMPUTES CONSTANTS AND LOOKUP TABLES
 C---BRW change 27/8/04: include Frixione's fix to reduce PDFSET calls
 C-----------------------------------------------------------------------
 C     COMPUTES CONSTANTS AND LOOKUP TABLES
 C---BRW change 27/8/04: include Frixione's fix to reduce PDFSET calls
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWBVMC,HWUALF,HWUPCM,XMIN,XMAX,XPOW,QR,DQKWT,
      & UQKWT,SQKWT,DIQWT,QMAX,PMAX,PTLIM,ETLIM,PGS,PTELM,X,QSCA,UPV,DNV,
      & USEA,DSEA,STR,CHM,BTM,TOP,GLU,VAL(20),CLMXPW,RCLPOW,TEST,RPM(2)
       DOUBLE PRECISION HWBVMC,HWUALF,HWUPCM,XMIN,XMAX,XPOW,QR,DQKWT,
      & UQKWT,SQKWT,DIQWT,QMAX,PMAX,PTLIM,ETLIM,PGS,PTELM,X,QSCA,UPV,DNV,
      & USEA,DSEA,STR,CHM,BTM,TOP,GLU,VAL(20),CLMXPW,RCLPOW,TEST,RPM(2)


@@ -61367,7 +61367,7 @@ C---FIX TO CALL SCHULER-SJOSTRAND CODE
               WRITE (6,'(10X,A)') 'THEIR WISHES.  SSph NO LONGER WORKS'
               STOP
             ELSE
               WRITE (6,'(10X,A)') 'THEIR WISHES.  SSph NO LONGER WORKS'
               STOP
             ELSE

-              CALL PDFSET(PARM,VAL)
+              CALL PDFSET_HERWIG(PARM,VAL)

               CALL STRUCTM(X,QSCA,UPV,DNV,USEA,DSEA,STR,CHM,BTM,TOP,GLU)
             ENDIF
           ENDIF
               CALL STRUCTM(X,QSCA,UPV,DNV,USEA,DSEA,STR,CHM,BTM,TOP,GLU)
             ENDIF
           ENDIF


@@ -61455,7 +61455,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     INITIALISES AN EVENT
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     INITIALISES AN EVENT
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWRGEN,HWRGET,DUMMY
       REAL TL
       LOGICAL CALLED,HWRLOG
       DOUBLE PRECISION HWRGEN,HWRGET,DUMMY
       REAL TL
       LOGICAL CALLED,HWRLOG


@@ -61677,7 +61677,7 @@ C-----------------------------------------------------------------------
 C     CHOOSES MASS ACCORDING TO BREIT-WIGNER DISTRIBUTION
 C--BRW fix 27/8/04: changed from mass to mass-squared BW formula
 C-----------------------------------------------------------------------
 C     CHOOSES MASS ACCORDING TO BREIT-WIGNER DISTRIBUTION
 C--BRW fix 27/8/04: changed from mass to mass-squared BW formula
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWUMBW,HWRGEN,WMX,TAU,GAM,T,TM
       INTEGER ID
 C--WMX IS MAX NUMBER OF WIDTHS FROM NOMINAL MASS
       DOUBLE PRECISION HWUMBW,HWRGEN,WMX,TAU,GAM,T,TM
       INTEGER ID
 C--WMX IS MAX NUMBER OF WIDTHS FROM NOMINAL MASS


@@ -61868,7 +61868,7 @@ C----------------------------------------------------------------------
 C     Prints contents of the GUPI (Generic User Process Interface)
 C     common block HEPEUP
 C----------------------------------------------------------------------
 C     Prints contents of the GUPI (Generic User Process Interface)
 C     common block HEPEUP
 C----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER MAXNUP
       PARAMETER (MAXNUP=500)
       INTEGER NUP,IDPRUP,IDUP,ISTUP,MOTHUP,ICOLUP
       INTEGER MAXNUP
       PARAMETER (MAXNUP=500)
       INTEGER NUP,IDPRUP,IDUP,ISTUP,MOTHUP,ICOLUP


@@ -61900,7 +61900,7 @@ C     ( rep. enhancement factor)^2  - used in cluster decays
 C     Finds iso-flavour hadrons and creates pointers for cluster decays.
 C     Sets CLDKWT(K) =(2J+1) spin weight normalizing largest value to 1.
 C-----------------------------------------------------------------------
 C     Finds iso-flavour hadrons and creates pointers for cluster decays.
 C     Sets CLDKWT(K) =(2J+1) spin weight normalizing largest value to 1.
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER NMXTMP
       PARAMETER (NMXTMP=20)
       DOUBLE PRECISION EPS,WTMX,REMMN,RWTMX,WTMP,RESTMP(91),WTMX2,
       INTEGER NMXTMP
       PARAMETER (NMXTMP=20)
       DOUBLE PRECISION EPS,WTMX,REMMN,RWTMX,WTMP,RESTMP(91),WTMX2,


@@ -62264,7 +62264,7 @@ C-----------------------------------------------------------------------
 C     RUNNING QUARK MASSES (MSBAR, 2-LOOP, 5 FLAVOUR, NO THRESHOLDS)
 C     ASSUMING RMASS(IQ) IS POLE MASS
 C-----------------------------------------------------------------------
 C     RUNNING QUARK MASSES (MSBAR, 2-LOOP, 5 FLAVOUR, NO THRESHOLDS)
 C     ASSUMING RMASS(IQ) IS POLE MASS
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       DOUBLE PRECISION HWUALF,SCALE,ALFAS,P0,C1,CC,MHAT(6),RQM(6)
       INTEGER IQ
       LOGICAL FIRST
       DOUBLE PRECISION HWUALF,SCALE,ALFAS,P0,C1,CC,MHAT(6),RQM(6)
       INTEGER IQ
       LOGICAL FIRST


@@ -62363,7 +62363,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C  Subroutine to output the contents of the spin common block
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C  Subroutine to output the contents of the spin common block
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER I
 C--write out the header
       WRITE(6,1000)
       INTEGER I
 C--write out the header
       WRITE(6,1000)


@@ -62394,7 +62394,7 @@ C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     MAKES PARTICLE TYPE 'NAME' STABLE
 C-----------------------------------------------------------------------
 C-----------------------------------------------------------------------
 C     MAKES PARTICLE TYPE 'NAME' STABLE
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER IPDG,IWIG
       CHARACTER*8 NAME
       CALL HWUIDT(3,IPDG,IWIG,NAME)
       INTEGER IPDG,IWIG
       CHARACTER*8 NAME
       CALL HWUIDT(3,IPDG,IWIG,NAME)


@@ -62597,7 +62597,7 @@ C                          300-399    FATAL, DUMP EVENT & TERMINATE RUN
 C                          400-499    FATAL, DUMP EVENT & STOP DEAD
 C                          500-       FATAL, STOP DEAD WITH NO DUMP
 C-----------------------------------------------------------------------
 C                          400-499    FATAL, DUMP EVENT & STOP DEAD
 C                          500-       FATAL, STOP DEAD WITH NO DUMP
 C-----------------------------------------------------------------------

-      INCLUDE 'HERWIG65.INC'
+      INCLUDE 'herwig65.inc'

       INTEGER ICODE
       CHARACTER*6 SUBRTN
       IF (ICODE.GE.0) WRITE (6,10) SUBRTN,ICODE
       INTEGER ICODE
       CHARACTER*6 SUBRTN
       IF (ICODE.GE.0) WRITE (6,10) SUBRTN,ICODE