zonal_neumann_core Subroutine

public subroutine zonal_neumann_core(comm_handler, mesh, equi, polars, u, v)

Sets Zonal Neumann boundary conditions on u (or v), u_b = , after local Neumann b.c. d/dr u_b_N = 0 had been applied to u. u is flux-surface averaged and the core value is used as a Dirichlet b.c.. For data handling reasons, in- and output can be separated with variable v = .

Arguments

Type IntentOptional Attributes Name
type(comm_handler_t), intent(in) :: comm_handler

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

Mesh
class(equilibrium_t), intent(in) :: equi

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

Polar grid
real(kind=GP), intent(inout), dimension(mesh%get_n_points()) :: u

Variable for which boundary conditions are computed and set
real(kind=GP), intent(inout), optional, dimension(mesh%get_n_points()) :: v

Variable for which boundary conditions are set instead of u

Contents