logLei Module Function

pure module function logLei(Te, ne, Z, removeDisc) result(res)

Calculate Coulomb logarithm for electron-ion collisions (NRL Formulary 2013 page 34 equation b). Assumes Te > TiZm_e/m_i

Arguments

Type IntentOptional Attributes Name
real(kind=rk), intent(in) :: Te

Electron temperature in eV

real(kind=rk), intent(in) :: ne

Electron density in m^{-3}

real(kind=rk), intent(in) :: Z

Ion charge

logical, intent(in), optional :: removeDisc

Remove the discontinuity at 10eV by moving the branch switch to Te = Z 2 * e2 eV. Defaults to .false.

Return Value real(kind=rk)