compute_vorticity Subroutine

public subroutine compute_vorticity(equi, mesh_cano, boundaries, co, potv, nev, tiv, vortv, swpot, swdia)

Computes generalised vorticity \f[ \nabla\cdot\left[co\nabla_\perp\left(swpot\nabla_\perp\phi+swdia\tau\frac{\nabla_\perp p_i}{n}\right)\right] \f]

Arguments

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

Equilibrium

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

Mesh (canonical)

type(boundaries_braginskii_t), intent(in) :: boundaries

Boundary information for the BRAGINSKII model

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

Polarisation coefficient

real(kind=GP), dimension(mesh_cano%get_n_points()) :: potv

Values of electrostatic potential

real(kind=GP), dimension(mesh_cano%get_n_points()) :: nev

Values of density

real(kind=GP), dimension(mesh_cano%get_n_points()) :: tiv

Values of ion temperature

real(kind=GP), intent(out), dimension(mesh_cano%get_n_points()) :: vortv

Generalised values for vorticity on inner grid

real(kind=GP), intent(in), optional :: swpot

Switch in front of \nabla_\perp\phi term, default = 1

real(kind=GP), intent(in), optional :: swdia

Switch in front of \nabla_\perp p_i/n term, default = 1


Contents