compute_diffusive_flux Subroutine

public subroutine compute_diffusive_flux(equi, mesh, dphi, polygon, conn, u, diff_flux, co)

Computes the diffusive flux through a perpendicular boundary polygon, i.e. \int_S v \nabla u dS, where nabla is gradient within plane Sign convention is that flux out of inner domain is positive

Arguments

Type IntentOptional Attributes Name
class(equilibrium_t), intent(in) :: equi

Equilibrium

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

Mesh

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

Toroidal extent of surface

type(polygon_in_mesh_t), intent(in) :: polygon

Polygon

type(csrmat_t), intent(in) :: conn

Connectivity matrix

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

Field

real(kind=GP), intent(out), dimension(polygon%get_np_total()) :: diff_flux

Diffusive flux through boundary polygon points (its associated area)

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

Coefficient


Contents