neoclassics_calc_nu_star Function

  1. neoclassics_module.f90
  2. neoclassics_module
  3. neoclassics_calc_nu_star

public function neoclassics_calc_nu_star() result(nu_star)

Uses

Calculates the normalized collision frequency

Arguments

None

Return Value real(kind=dp), dimension(4,no_roots)


Contents

Source Code

neoclassics_calc_nu_star

Source Code

    function neoclassics_calc_nu_star() result(nu_star)
        !! Calculates the normalized collision frequency
        !
        ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

        use const_and_precisions, only: e_,c_,me_,mp_
        use physics_variables, only: rmajor

        real(dp), dimension(no_roots) ::  K
        real(dp), dimension(4,no_roots) :: v,nu_star,KK
        real(dp), dimension(4) :: mass

        K = roots

        KK(1,:) = K * temperatures(1)
        KK(2,:) = K * temperatures(2)
        KK(3,:) = K * temperatures(3)
        KK(4,:) = K * temperatures(4)

        mass = (/me_,mp_*2.0d0,mp_*3.0d0,mp_*4.0d0/)

        v(1,:) = c_ * sqrt(1.0d0-(KK(1,:)/(mass(1) * c_**2)+1)**(-1))
        v(2,:) = c_ * sqrt(1.0d0-(KK(2,:)/(mass(2) * c_**2)+1)**(-1))
        v(3,:) = c_ * sqrt(1.0d0-(KK(3,:)/(mass(3) * c_**2)+1)**(-1))
        v(4,:) = c_ * sqrt(1.0d0-(KK(4,:)/(mass(4) * c_**2)+1)**(-1))

        nu_star = rmajor * nu/(iota*v)

    end function neoclassics_calc_nu_star