smod_equil_discrete_alfven Submodule

This submodule defines an equilibrium in cylindrical geometry with an axial current profile ($\nu = 2$), modelling a solar coronal loop in which discrete Alfvén waves are present. The geometry can be overridden in the parfile.

This equilibrium is taken from Keppens, Rony, Ronald AM Van Der Linden, and Marcel Goossens. "Non-adiabatic discrete Alfven waves in coronal loops and prominences.", Solar physics 144.2 (1993): 267-281.

Contents

Variables

x_end

Functions

rho0 drho0 T0 dT0 B02 dB02 B03 p0 dp0

Module Procedures

discrete_alfven_eq

Variables

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

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 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()

Arguments

None

Return Value real(kind=dp)

function p0(r)

Arguments

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

Return Value real(kind=dp)

function dp0(r)

Arguments

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

Return Value real(kind=dp)

Module Procedures

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

Arguments

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