SKICOV    SKICOV    SKICOV    SKICOV    SKICOV    SKICOV    SKICOV    
          **********************************************
                    **************************
                              ******

SUBROUTINE SKICOV (IR, NDIM, ISKWRK, RSKWRK, NCOVI, KCOVI, CPARI)
Returns the coverage species index numbers and their coverage
parameters for reaction IR.

INPUT
IR        - Integer scalar, surface reaction index.
NDIM      - Integer scalar, first dimension of array CPAR, the
            coverage parameters; NDIM must be at least NSCOV,
            the total number of coverage parameters.
ISKWRK(*) - Integer workspace array; dimension at least LENISK.
RSKWRK(*) - Real    workspace array; dimension at least LENRSK.

OUTPUT
NCOVI     - Integer scalar, total number of species that modify the
            rate of reaction IR through coverage dependence.
KCOVI(*)  - Integer array, species indices for the NCOVI species that
            modify the rate of a coverage dependence reaction;
            dimension at least KKTOT, the total species count.
CPARI(*,*)- Real matrix, coverage parameters for the coverage species
            of reaction IR;
            dimension at least NSCOV for the first, the number of
            coverage parameters required, and at least KKTOT for the
            second, the total species count.