smod_equil_internal_kink_instability – Legolas v. 2.1.1

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.

Note

Default values are given by

k2 = 1
k3 = $0.16\alpha$
cte_rho0 = 1 : used as prefactor in setting the density.
cte_v03 = 1 : used as prefactor in setting the z-component of velocity.
cte_p0 = 3 : used to set the pressure.
alpha = 5 / x_end : used in the Bessel functions.

and can all be changed in the parfile.

Uses

- mod_equilibrium_params
- Ancestors:
- mod_equilibrium

Variables

Functions

rho0 drho0 T0 dT0 v03 dv03 B02 dB02 B03 dB03 J0 J1 J2 dJ0 dJ1

Module Procedures

internal_kink_eq

Variables

Type	Visibility	Attributes		Name		Initial
real(kind=dp),	private		::	a0

Functions

function rho0(r)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	r

Return Value real(kind=dp)

function drho0(r)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	r

Return Value real(kind=dp)

function T0(r)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	r

Return Value real(kind=dp)

function dT0(r)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	r

Return Value real(kind=dp)

function v03(r)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	r

Return Value real(kind=dp)

function dv03(r)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	r

Return Value real(kind=dp)

function B02(r)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	r

Return Value real(kind=dp)

function dB02(r)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	r

Return Value real(kind=dp)

function B03(r)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	r

Return Value real(kind=dp)

function dB03(r)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	r

Return Value real(kind=dp)

function J0(r)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	r

Return Value real(kind=dp)

function J1(r)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	r

Return Value real(kind=dp)

function J2(r)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	r

Return Value real(kind=dp)

function dJ0(r)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	r

Return Value real(kind=dp)

function dJ1(r)

Arguments

Type	Intent	Optional	Attributes		Name
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

Type	Intent		Name
type(settings_t),	intent(inout)	::	settings
type(grid_t),	intent(inout)	::	grid
type(background_t),	intent(inout)	::	background
type(physics_t),	intent(inout)	::	physics

smod_equil_internal_kink_instability Submodule

Contents

Variables

Functions

Module Procedures

Note

Uses

Ancestors:

Contents

Variables

Functions

Module Procedures

Variables

Functions

function rho0(r)

Arguments

Return Value real(kind=dp)

function drho0(r)

Arguments

Return Value real(kind=dp)

function T0(r)

Arguments

Return Value real(kind=dp)

function dT0(r)

Arguments

Return Value real(kind=dp)

function v03(r)

Arguments

Return Value real(kind=dp)

function dv03(r)

Arguments

Return Value real(kind=dp)

function B02(r)

Arguments

Return Value real(kind=dp)

function dB02(r)

Arguments

Return Value real(kind=dp)

function B03(r)

Arguments

Return Value real(kind=dp)

function dB03(r)

Arguments

Return Value real(kind=dp)

function J0(r)

Arguments

Return Value real(kind=dp)

function J1(r)

Arguments

Return Value real(kind=dp)

function J2(r)

Arguments

Return Value real(kind=dp)

function dJ0(r)

Arguments

Return Value real(kind=dp)

function dJ1(r)

Arguments

Return Value real(kind=dp)

Module Procedures

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

Arguments

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