parallel_divergence_nonaligned Subroutine

public subroutine parallel_divergence_nonaligned(comm, equi, mesh, map, uflx, pardivu, order)

Computes parallel divergence from parallel flux defined at 3D staggered positions Result (pardivu) is defined on full grid Based on symmetric formulation of Guenter's scheme (2nd and 4th order available)

Arguments

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

MPI communicator

class(equilibrium_t), intent(inout) :: equi

Equilibrium

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

Mesh within poloidal plane

type(parallel_map_t), intent(in) :: map

Parallel map

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

Parallel flux

real(kind=GP), intent(out), dimension(mesh%get_n_points_inner()) :: pardivu

Parallel divergence of uflx

integer, intent(in), optional :: order

Order of scheme default: 2


Contents