SUBROUTINE RRK(T,RMS,PS1,PS2,TP,ANS)
      implicit none
C     ****
C     ****      EVALUATES  ANS = K (DIFFUSION COEFFICIENT)
C     ****
      include "params.h"
      include "cons.h"
      real :: rmn4s,rmn2d,rmno,brn2d,cee
      COMMON/MASS/RMN4S,RMN2D,RMNO,BRN2D,CEE
!
! Args:
      real,intent(in) :: T(ZIMXP,ZKMXP),RMS(ZIMXP,ZKMXP),
     |  PS1(ZIMXP,ZKMXP),PS2(ZIMXP,ZKMXP),TP(ZIMXP,ZKMXP)
      real,intent(out) :: ANS(ZIMXP,ZKMXP)
!
! Local:
      integer :: i,k
C     ****
C     ****      SET UP EXPS ARRAY
C     ****
      DO 1 K=1,KMAX
        DO 2 I=1,LEN1
          ANS(I,K) = EXPS(K)
    2   CONTINUE
    1 CONTINUE

      write(6,"('rrk: cee=',e12.4)") cee

      DO 3 I=1,LEN2
        ANS(I,1) = 1.42E17*C(84)*T(I,1)/(C(81)*ANS(I,1)*.5*(RMS(I,1)+
     A   RMS(I,2))*(PS2(I,1)/RMASS(2)*SQRT(TP(I,1))*(1.-0.064*
     B   ALOG10(TP(I,1)))**2*CEE+18.6*(1.-PS1(I,1)-PS2(I,1))/RMASS(3)
     C   +18.1*PS1(I,1)/RMASS(1)))
    3 CONTINUE
      RETURN
      END
C