neoclassics_calc_Gamma_flux Function

public function neoclassics_calc_Gamma_flux(densities, temperatures, dr_densities, dr_temperatures)

Calculates the Energy flux by neoclassical particle transport

Arguments

Type IntentOptional AttributesName
real(kind=dp), dimension(4):: densities
real(kind=dp), dimension(4):: temperatures
real(kind=dp), dimension(4):: dr_densities
real(kind=dp), dimension(4):: dr_temperatures

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


Contents


Source Code

    function neoclassics_calc_Gamma_flux(densities, temperatures, dr_densities, dr_temperatures)
        !! Calculates the Energy flux by neoclassical particle transport
        !
        ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

        real(dp),dimension(4) :: neoclassics_calc_Gamma_flux, densities, dr_densities, z, temperatures, dr_temperatures


        z = (/-1.0,1.0,1.0,2.0/)

        neoclassics_calc_Gamma_flux = - densities * D111 * ((dr_densities/densities - z * Er/temperatures)+ &
                        (D112/D111-3.0/2.0) * dr_temperatures/temperatures )

    end function neoclassics_calc_Gamma_flux