smod_equil_gold_hoyle – Legolas v. 2.1.1

This submodule defines a Gold-Hoyle equilibrium in cylindrical geometry. This equilibrium configuration models a filament with a uniform twist such that all fieldlines perform an equal amount of turns around the cylinder axis. The geometry can be overridden in the parfile.

This equilibrium is taken from Van der Linden, R. A. M., Goossens, M. (1991). "The thermal continuum in coronal loops: instability criteria and the influence of perpendicular thermal conduction." Solar physics, 134.2, 247-273.

Note

Default values are given by

k2 = 1
k3 = 1
cte_T0 = 1 : used to set the temperature.
cte_rho0 = 1 : used to set the density
alpha = 20 : used in the magnetic field components
cooling_curve = 'rosner'
parallel thermal conduction, no perpendicular conduction

and can all be changed in the parfile.

Note

To reproduce the three profiles of the original paper you can supply one of the following normalisations:

cold plasma: unit_magneticfield = 10 Gauss, unit_density = 1.6726e-15 g/cm3, unit_length = 1.0e10 cm, corresponding to a temperature of 5.7e4 K.
cool plasma: unit_magneticfield = 22.5 Gauss, unit_density = 1.6726e-15 g/cm3, unit_length = 1.0e10 cm, corresponding to a temperature of 2.9e5 K (default).
hot plasma: unit_magneticfield = 67 Gauss, unit_density = 1.6726e-15 g/cm3, unit_length = 1.0e11 cm, corresponding to a temperature of 2.6e6 K.

Uses

- mod_equilibrium_params
- Ancestors:
- mod_equilibrium

Functions

rho0 T0 B02 dB02 B03 dB03

Module Procedures

gold_hoyle_eq

Functions

function rho0()

Arguments

None

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)				::	r

Return Value real(kind=dp)

function dB02(r)

Arguments

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

Return Value real(kind=dp)

function B03(r)

Arguments

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

Return Value real(kind=dp)

function dB03(r)

Arguments

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

Return Value real(kind=dp)

Module Procedures

module procedure gold_hoyle_eq module module subroutine gold_hoyle_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_gold_hoyle Submodule

Contents

Functions

Module Procedures

Note

Note

Uses

Ancestors:

Contents

Functions

Module Procedures

Functions

function rho0()

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)

function B03(r)

Arguments

Return Value real(kind=dp)

function dB03(r)

Arguments

Return Value real(kind=dp)

Module Procedures

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

Arguments

smod_equil_gold_hoyle Submodule

Contents

Functions

Module Procedures

Note

Note

Uses

Ancestors:

Contents

Functions

Module Procedures

Functions

function rho0()

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)

function B03(r)

Arguments

Return Value real(kind=dp)

function dB03(r)

Arguments

Return Value real(kind=dp)

Module Procedures

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

Arguments

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