MCSDIF    MCSDIF    MCSDIF    MCSDIF    MCSDIF    MCSDIF    MCSDIF    
          **********************************************
                    **************************
                              ******

SUBROUTINE MCSDIF (P, T, KDIM, RMCWRK, DJK)
Returns the binary diffusion coefficients given pressure and
temperature.

INPUT
P         - Real scalar, pressure.
               cgs units, dynes/cm**2
T         - Real scalar, temperature.
               cgs units, K
KDIM      - Integer scalar, actual first dimension of DJK(KDIM,KK).
RMCWRK(*) - Real workspace array; dimension at least LENRMC.

OUTPUT
DJK(*,*)  - Real matrix, binary diffusion coefficients;
            dimension at least KK, the total species count, for
            both the first and second dimensions.
               cgs units, cm**2/s
            CJK(J,K) is the diffusion coefficient of species J
            in species K.