drift_surface_integral Subroutine

public subroutine drift_surface_integral(comm, equi, mesh, polars, dphi, u, drift_flux, v)

Computes surface integral of flux in drift form, i.e. \int dS v/B^2 B\times\nabla u via integration along flux surface, i.e. = \int v/B du/d\theta J dphi d\theta

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: comm

MPI-communicator

class(equilibrium_t) :: equi

EQuilibrium

type(mesh_cart_t), intent(in) :: mesh

Mesh within poloidal plane

type(polars_t), intent(in) :: polars

Polar grid and matrices

real(kind=GP), intent(in) :: dphi

Toroidal grid distance

real(kind=GP), intent(in), dimension(mesh%get_n_points()) :: u

Field that gardient acts upon

real(kind=GP), intent(out), dimension(polars%grid%get_nrho()) :: drift_flux

Drift flux

real(kind=GP), intent(in), optional, dimension(mesh%get_n_points()) :: v

Secondary field, default values = 1


Contents