CLogei_wu Function

  1. korc_collisions.f90
  2. korc_collisions
  3. CLogei_wu

public function CLogei_wu(params, ne, Te)

With units

Arguments

Type IntentOptional AttributesName
type(KORC_PARAMS), intent(in) :: params
real(kind=rp), intent(in) :: ne

ne is in m^-3 and below is converted to cm^-3
real(kind=rp), intent(in) :: Te

Return Value real(kind=rp)


Calls

proc~~clogei_wu~~CallsGraph proc~clogei_wu CLogei_wu proc~clog0_wu CLog0_wu proc~clogei_wu->proc~clog0_wu proc~vte_wu VTe_wu proc~clogei_wu->proc~vte_wu

Called by

proc~~clogei_wu~~CalledByGraph proc~clogei_wu CLogei_wu proc~load_params_ss load_params_ss proc~load_params_ss->proc~clogei_wu proc~initialize_collision_params initialize_collision_params proc~initialize_collision_params->proc~load_params_ss program~main main program~main->proc~initialize_collision_params

Contents

Source Code

CLogei_wu

Source Code

  function CLogei_wu(params,ne,Te)
    
    !! With units
    TYPE(KORC_PARAMS), INTENT(IN) 	:: params
    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_wu
    REAL(rp)  :: k=5._rp
    REAL(rp)  :: p

    p=sqrt(params%minimum_particle_g**2-1)

    if (cparams_ss%Clog_model.eq.'CONSTANT') then
       CLogei_wu = cparams_ss%Clog_const
    else
       CLogei_wu = CLog0_wu(ne,Te)+ &
            log(1+(2*p/(VTe_wu(Te)/C_C))**k)/k
    end if
  end function CLogei_wu