cjl For Herwig version 6.202, May-2001
      subroutine herwiginterface
c----------------------------------------------------------------------
c   program initialization
c----------------------------------------------------------------------
c
c  HERWIG common blocks (`external` loads in  block data)
c-----------------------------------------------------------------------
cc#include  "herwig.inc"
#include "HERWIG62.INC"
      external hwudat
c-----------------------------------------------------------------------
c  local variables
c-----------------------------------------------------------------------
      integer       status, luherc, ierr, i, mm
cjl
      integer ii,jj,j1,j2,j3
      INTEGER IPDG,IWIG
      CHARACTER*8 NAME
c------------------------------------------------------------------------------ 
c  set defaults
c------------------------------------------------------------------------------
      call hwigin                       ! set HERWIG defaults
c
      maxpr         = 0                 ! Turn off screen data output
      lwsud         = 0
      lwevt         = 0
      iprint        = 1
      PART1         = 'P       '
      PART2         = 'PBAR    '
      IPROC         = 1500
      maxev         = 100
      maxer         = 100
      pbeam1        = 900.0
      pbeam2        = 900.0
      lhwout        = 6
cjl sudord=2 for v 6.1
      sudord        = 2
cjl set default tlout to large negative, so timing-out feature is inoperative;
cjl   on systems where herwig's call timel works, user can reset tlout to 10
cjl   (or whatever) in talkto if desired (see Cern KERNLIB short writup Z007).
      tlout         = -10000000.0
cjl turn off mixing in Herwig - do it in QQ
      mixing = .false.
cjl set stable all particles with b or c quarks - let QQ do them
      do ii=1,nres
        jj=abs(iflav(ii))
        j1=mod(jj,10)
        jj=jj/10
        j2=mod(jj,10)
        jj=jj/10
        j3=mod(jj,10)
        if(j1.eq.4.or.j2.eq.4.or.j3.eq.4) then
cjl  copied from hwusta.f
          CALL HWUIDT(3,IPDG,IWIG,rname(ii))
          RSTAB(IWIG)=.TRUE.
        else if(j1.eq.5.or.j2.eq.5.or.j3.eq.5) then
          CALL HWUIDT(3,IPDG,IWIG,rname(ii))
          RSTAB(IWIG)=.TRUE.
        end if
      end do
c
      WRITE (LHWOUT,10)
   10 FORMAT(/10X,
     &    'HERWIG: ALL PARTICLE TYPE WITH A B QUARK SET STABLE, USE QQ')
      WRITE (LHWOUT,12)
 12   FORMAT(/10X,
     &    'HERWIG: ALL PARTICLE TYPE WITH A C QUARK SET STABLE, USE QQ')
cjl set default pltcut at 1.D-11, so default is to NOT decay K_short and
c     weak-decaying strange baryons in Herwig - do it in geant; can change
c     in talk-to
      pltcut=1.D-11
cjl set default pdf to MRS set G (replace Owens set 1.1 from 1991)
      modpdf(1)=41
      modpdf(2)=41
      autpdf(1)='MRS'
      autpdf(2)='MRS'
c
      return
      end

C
c
      Subroutine cdfreadsusy(filnam,iunit)
      character*(*)  filnam
      integer status
      integer sigtlu, sifrlu
      integer iunit
c-----------------------------------------------------------------------
c  open the SUSY particle properties file
c-----------------------------------------------------------------------

      status        = sigtlu(iunit)
      write(*,*) ' '
      write (*,6660)
 6660 format (' SUSY filename is ')
      write (*,*) filnam
      write (*,6664) iunit
 6664 format ('      and unit number is ',i6)
C
      open(unit=iunit,file=filnam,status='old',err=101
     &     ,FORM='FORMATTED',ACCESS='SEQUENTIAL')
      goto 102
c-----------------------------------------------------------------------
c  error: can't open SUSY file
c-----------------------------------------------------------------------
101   continue
      write(*,*) ' '
      write(*,*) ' ********************************************* '
      write(*,*) ' '
      write(*,*) 
     $     ' Error in opening SUSY file for READING'
      write(*,*)
     $     ' Program will abort.'
      write(*,*) ' '
      write(*,*) ' ********************************************* '
      write(*,*) ' '
102   continue
      end
c
c-----------------------------------------------------------------------
c  dummy hwanal, called from HerwigModule.cc, can be useful
      subroutine hwanal(ihwcod)
c  input ihwcod = 1, set to zero if want this event skipped
cc#include  "herwig.inc"
#include "HERWIG62.INC" 
      integer ihwcod
c
c
c skip if ierror ne 0
      if(ierror.ne.0) return
c
cc  do whatever here
cc  if want to skip this event, say  ihwcod = 0
c
      end
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c  dummy needed
      SUBROUTINE HWAEND
C     USER ROUTINE FOR TERMINAL CALCULATIONS, HISTOGRAM OUTPUT, ETC
C----------------------------------------------------------------------
      END