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 |
function CF_SD_FIO(params,v,ne,Te,nimp)
TYPE(KORC_PARAMS), INTENT(IN) :: params
REAL(rp), INTENT(IN) :: v
REAL(rp), INTENT(IN) :: ne
REAL(rp), INTENT(IN) :: Te
REAL(rp),DIMENSION(cparams_ms%num_impurity_species), INTENT(IN) :: nimp
REAL(rp) :: CF_SD_FIO
REAL(rp) :: CF_temp
REAL(rp) :: x
INTEGER :: i
REAL(rp) :: k=5._rp
x = v/VTe(Te)
CF_SD_FIO = Gammacee(v,ne,Te)*psi(x)/Te
if (params%bound_electron_model.eq.'HESSLOW') then
CF_temp=CF_SD_FIO
do i=1,cparams_ms%num_impurity_species
CF_temp=CF_temp+CF_SD_FIO*nimp(i)/ne* &
(cparams_ms%Zo(i)-cparams_ms%Zj(i))/ &
CLogee(v,ne,Te)*(log(1+h_j(i,v)**k)/k-v**2)
end do
CF_SD_FIO=CF_temp
end if
end function CF_SD_FIO