calculate_fluxes_fullfluxsurface Subroutine

public subroutine calculate_fluxes_fullfluxsurface(dens, upar, pres, pflux_vs_s, vflux_vs_s, qflux_vs_s, pflux_kxkyzts, vflux_kxkyzts, qflux_kxkyzts)

Calculate the total particle, momentum and heat fluxes (pflux_vs_s, vflux_vs_s, qflux_vs_s) and the contributions from a given (kx,ky,z) location (pflux_kxkyz, vflux_kxkyz, qflux_kxkyz) inputs are the particle density (dens), parallel flow (upar) and pressure (pres)

assume a single flux annulus

Set all fluxes to zero just in case Obtain the y-component of the electric field that appears as a factor in the flux expression due to the radial component of the ExB velocity

Calculate Jacobian for spacial integral

Arguments

Type IntentOptional Attributes Name
complex, intent(in), dimension(:, :, -nzgrid:, :) :: dens
complex, intent(in), dimension(:, :, -nzgrid:, :) :: upar
complex, intent(in), dimension(:, :, -nzgrid:, :) :: pres
real, intent(out), dimension(:) :: pflux_vs_s
real, intent(out), dimension(:) :: vflux_vs_s
real, intent(out), dimension(:) :: qflux_vs_s
real, intent(out), dimension(:, :, -nzgrid:, :, :) :: pflux_kxkyzts
real, intent(out), dimension(:, :, -nzgrid:, :, :) :: vflux_kxkyzts
real, intent(out), dimension(:, :, -nzgrid:, :, :) :: qflux_kxkyzts