CLogee_wu Function

  1. korc_collisions.f90
  2. korc_collisions
  3. CLogee_wu

public function CLogee_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~~clogee_wu~~CallsGraph proc~clogee_wu CLogee_wu proc~clog0_wu CLog0_wu proc~clogee_wu->proc~clog0_wu proc~vte_wu VTe_wu proc~clogee_wu->proc~vte_wu

Called by

proc~~clogee_wu~~CalledByGraph proc~clogee_wu CLogee_wu proc~load_params_ss load_params_ss proc~load_params_ss->proc~clogee_wu proc~gammac_wu Gammac_wu proc~load_params_ss->proc~gammac_wu proc~gammac_wu->proc~clogee_wu proc~initialize_collision_params initialize_collision_params proc~initialize_collision_params->proc~load_params_ss proc~load_params_ms load_params_ms proc~initialize_collision_params->proc~load_params_ms proc~load_params_ms->proc~gammac_wu program~main main program~main->proc~initialize_collision_params

Contents

Source Code

CLogee_wu

Source Code

  function CLogee_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) 				:: CLogee_wu
    REAL(rp)  :: k=5._rp

    if (cparams_ss%Clog_model.eq.'HESSLOW') then
       CLogee_wu = CLog0_wu(ne,Te)+ &
            log(1+(2*(params%minimum_particle_g-1)/ &
            (VTe_wu(Te)/C_C)**2)**(k/2._rp))/k

    else if (cparams_ss%Clog_model.eq.'CONSTANT') then
       CLogee_wu = cparams_ss%Clog_const

    else if (cparams_ss%Clog_model.eq.'MCDEVITT') then
       CLogee_wu = CLog0_wu(ne,Te)+ &
            log(1+(2*(params%minimum_particle_g-1)/ &
            (VTe_wu(Te)/C_C)**2)**(k/2._rp))/k
    end if
       
  end function CLogee_wu