Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
type(KORC_PARAMS), | intent(in) | :: | params | |||
real(kind=rp), | intent(in) | :: | v | |||
real(kind=rp), | intent(in) | :: | ne | |||
real(kind=rp), | intent(in) | :: | Te | |||
real(kind=rp), | intent(in), | DIMENSION(cparams_ms%num_impurity_species) | :: | nimp | ||
real(kind=rp), | intent(in) | :: | Zeff |
function CB_ei_SD_FIO(params,v,ne,Te,nimp,Zeff)
TYPE(KORC_PARAMS), INTENT(IN) :: params
REAL(rp), INTENT(IN) :: v
REAL(rp), INTENT(IN) :: ne
REAL(rp), INTENT(IN) :: Te
REAL(rp), INTENT(IN) :: Zeff
REAL(rp), DIMENSION(cparams_ms%num_impurity_species), INTENT(IN) :: nimp
REAL(rp) :: CB_ei_SD_FIO
REAL(rp) :: CB_ei_temp
REAL(rp) :: x
INTEGER :: i
x = v/VTe(Te)
CB_ei_SD_FIO = (0.5_rp*Gammacee(v,ne,Te)/v)* &
(Zeff*CLogei(v,ne,Te)/CLogee(v,ne,Te))
if (params%bound_electron_model.eq.'HESSLOW') then
CB_ei_temp=CB_ei_SD_FIO
do i=1,cparams_ms%num_impurity_species
CB_ei_temp=CB_ei_temp+CB_ei_SD_FIO*nimp(i)/(ne* &
Zeff*CLogei(v,ne,Te))*g_j(i,v)
end do
CB_ei_SD_FIO=CB_ei_temp
end if
end function CB_ei_SD_FIO