smod_equil_internal_kink_instability Submodule

This submodule defines internal kink modes in force-free magnetic fields. The geometry is cylindrical with parabolic density and velocity profiles, The geometry can be overridden in the parfile.

This equilibrium is taken from section III.B in Goedbloed, J. P. "The Spectral Web of stationary plasma equilibria. II. Internal modes." Physics of Plasmas 25.3 (2018): 032110.



Contents


Variables

TypeVisibilityAttributesNameInitial
real(kind=dp), private :: a0

Functions

function rho0(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function drho0(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function T0(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function dT0(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function v03(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function dv03(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function B02(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function dB02(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function B03(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function dB03(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function J0(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function J1(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function J2(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function dJ0(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)

function dJ1(r)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: r

Return Value real(kind=dp)


Module Procedures

module procedure internal_kink_eq module module subroutine internal_kink_eq(settings, grid, background, physics)

Arguments

TypeIntentOptionalAttributesName
type(settings_t), intent(inout) :: settings
type(grid_t), intent(inout) :: grid
type(background_t), intent(inout) :: background
type(physics_t), intent(inout) :: physics