c
c###################################################################
c
      subroutine jrrec2(pt, pcf, kncf, pfit, kermsg)
c
c     this is a subroutine to calculate a reaction rate coefficient
c     in cm[3]/s as function of electron tempreature in ev for
c     radiative recombination for specific transitions to excited
c     states of neutral helium.
c
c     pcf is the coefficient data array, where
c
c     pcf(1)  =   threshold energy (ev)
c
c     pcf(2)  =   coefficient a (nl)
c
c     pt = tepreature (ev)
c
c     kermsg = blank if no errors
c
c     pfit = reaction rate coefficient in cm[3]/s
c
c     written by j. j. smith , iaea atomic and molecular data unit
c
c------------------------------------------------------------------------
c
      double precision pt, pcf, pfit
      dimension pcf(2)
      character*(*) kermsg
c
      kermsg = ' '
      eth=pcf(1)
      a=pcf(2)
c
      beta = eth/pt
c
      pfit = a * 1.0e-14 * (beta**1.5) * exp(beta) * exint(beta)
      return
c
      end