Calculates the Energy flux by neoclassical particle transport
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
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 |
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