Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=rp), | intent(in) | :: | v | |||
real(kind=rp), | intent(in) | :: | ne | ne is in m^-3 and below is converted to cm^-3 |
||
real(kind=rp), | intent(in) | :: | Te |
function CLogei(v,ne,Te)
REAL(rp), INTENT(IN) :: v
REAL(rp), INTENT(IN) :: ne
!! ne is in m^-3 and below is converted to cm^-3
REAL(rp), INTENT(IN) :: Te ! In Joules
REAL(rp) :: CLogei
REAL(rp) :: k=5._rp
REAL(rp) :: gam,p
gam=1/sqrt(1-v**2)
p=gam*v
if (cparams_ss%Clog_model.eq.'CONSTANT') then
CLogei = cparams_ss%Clog_const
else
CLogei = CLog0(ne,Te)+log(1+(2*p/VTe(Te))**k)/k
end if
end function CLogei