### Fit Function: JRREC1

 Fit Function \begin{align*} \end{align*} Comments

#### Fortran

Arguments
namedescriptionunitstype(s)
pt requested photon energy eV real, dimension(:)
pcf coefficient data array real, dimension(9)
kncf number of coefficients in the data array integer
pfit cross section cm2 real, dimension(:)
kermsg error message character
Return values
namedescriptionunitstype(s)
pfit cross section cm2 real, dimension(:)
Code
c
c###################################################################
c
subroutine jrrec1(pt, pcf, kncf, pfit, kermsg)
c
c     this is a subroutine to calculate a reaction rate coefficient
c     in cm/s as function of electron tempreature in ev for
c
c     pcf is the coefficient data array, where
c
c     pcf(1) =   the threshold energy for the transition
c
c     pcf(2)  =   coefficient a(nl)
c
c     pcf(3)  =   coefficient chi(nl)
c
c     pt = tempreature (ev)
c
c     kermsg = blank if no errors
c
c     pfit = reaction rate coefficient in cm/s
c
c     written by j. j. smith , iaea atomic and molecular data unit
c
c------------------------------------------------------------------------
c
double precision pt, pcf, pfit
real n
dimension pcf(3)
character*(*) kermsg
data ry/13.58/

eth=pcf(1)
a=pcf(2)
chi=pcf(3)
beta=eth/pt
c
pfit = a * 1.0e-14 * sqrt(eth/ry) * (beta**1.5) / (beta + chi)
return
c
end

#### Python

Arguments
namedescriptionunitstype(s)
pt requested photon energy eV float, np.ndarray
pcf coefficient data array float, np.ndarray
kncf number of coefficients in the data array int
pfit cross section cm2 float, np.ndarray
kermsg error message str
Return values
namedescriptionunitstype(s)
pfit cross section cm2 float, np.ndarray
Code