smod_equil_magnetothermal_instabilities – Legolas v. 2.1.1

This submodule defines an equilibrium containing magnetothermal instabilities in a cylindrical geometry. The geometry can be overridden in the parfile.

This equilibrium is taken from Van der Linden, R. A. M., Goossens, M., & Hood, A. W. (1992). The relevance of the ballooning approximation for magnetic, thermal, and coalesced magnetothermal instabilities. Solar physics, 140(2), 317-342.

Note

Default values are given by

k2 = 0
k3 = 1
cte_T0 = 1 : used to set the temperature.
cooling_curve = 'rosner'
parallel thermal conduction, no perpendicular conduction

and normalisations given by

unit_temperature = 2.6e6 K
unit_magneticfield = 10 Gauss
unit_length = 1e8 cm

and can all be changed in the parfile.

Note

The default setup handles CASE I in the original paper. For CASE II, you can set k2 = 10. To reproduce the thermal continuum plot, set unit_length = 2.44e9 cm in CASE I.

Uses

- mod_equilibrium_params
- Ancestors:
- mod_equilibrium

Functions

p0 rho0 drho0 T0 B02 dB02

Module Procedures

magnetothermal_instability_eq

Functions

function p0(r)

Arguments

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

Return Value real(kind=dp)

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

Arguments

None

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)

Module Procedures

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

Sets the equilibrium

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_magnetothermal_instabilities Submodule

Contents

Functions

Module Procedures

Note

Note

Uses

Ancestors:

Contents

Functions

Module Procedures

Functions

function p0(r)

Arguments

Return Value real(kind=dp)

function rho0(r)

Arguments

Return Value real(kind=dp)

function drho0(r)

Arguments

Return Value real(kind=dp)

function T0()

Arguments

Return Value real(kind=dp)

function B02(r)

Arguments

Return Value real(kind=dp)

function dB02(r)

Arguments

Return Value real(kind=dp)

Module Procedures

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

Arguments

smod_equil_magnetothermal_instabilities Submodule

Contents

Functions

Module Procedures

Note

Note

Uses

Ancestors:

Contents

Functions

Module Procedures

Functions

function p0(r)

Arguments

Return Value real(kind=dp)

function rho0(r)

Arguments

Return Value real(kind=dp)

function drho0(r)

Arguments

Return Value real(kind=dp)

function T0()

Arguments

Return Value real(kind=dp)

function B02(r)

Arguments

Return Value real(kind=dp)

function dB02(r)

Arguments

Return Value real(kind=dp)

Module Procedures

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

Arguments

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