CF_FIO Function

private function CF_FIO(params, v)

Arguments

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

Return Value real(kind=rp)


Calls

proc~~cf_fio~~CallsGraph proc~cf_fio CF_FIO proc~psi psi proc~cf_fio->proc~psi proc~h_j h_j proc~cf_fio->proc~h_j proc~clogee CLogee proc~cf_fio->proc~clogee proc~clog0 CLog0 proc~clogee->proc~clog0 proc~vte VTe proc~clogee->proc~vte

Called by

proc~~cf_fio~~CalledByGraph proc~cf_fio CF_FIO proc~nu_s_fio nu_S_FIO proc~nu_s_fio->proc~cf_fio proc~define_collisions_time_step define_collisions_time_step proc~define_collisions_time_step->proc~nu_s_fio program~main main program~main->proc~define_collisions_time_step

Contents

Source Code


Source Code

  function CF_FIO(params,v)
    TYPE(KORC_PARAMS), INTENT(IN) 	:: params
    REAL(rp), INTENT(IN) 	:: v
    REAL(rp) 				:: CF_FIO
    REAL(rp) 				:: CF_temp
    REAL(rp) 				:: x
    INTEGER :: i
    REAL(rp)  :: k=5._rp

    x = v/cparams_ss%VTe
    CF_FIO  = cparams_ss%Gammac*psi(x)/cparams_ss%Te

    if (params%bound_electron_model.eq.'HESSLOW') then
       CF_temp=CF_FIO
       do i=1,1
          CF_temp=CF_temp+CF_FIO*cparams_ms%nz(i)/cparams_ms%ne* &
               (cparams_ms%Zo(i)-cparams_ms%Zj(i))/ &
               CLogee(v,cparams_ss%ne,cparams_ss%Te)* &
               (log(1+h_j(i,v)**k)/k-v**2) 
       end do
       CF_FIO=CF_temp
       
    end if
    
  end function CF_FIO