smod_equil_gold_hoyle Submodule

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


Functions

function rho0()

Arguments

None

Return Value real(kind=dp)

function T0()

Arguments

None

Return Value real(kind=dp)

function B02(r)

Arguments

Type IntentOptional Attributes Name
real(kind=dp) :: r

Return Value real(kind=dp)

function dB02(r)

Arguments

Type IntentOptional Attributes Name
real(kind=dp) :: r

Return Value real(kind=dp)

function B03(r)

Arguments

Type IntentOptional Attributes Name
real(kind=dp) :: r

Return Value real(kind=dp)

function dB03(r)

Arguments

Type IntentOptional Attributes Name
real(kind=dp) :: r

Return Value real(kind=dp)


Module Procedures

module procedure /home/runner/work/legolas/legolas/docs/ford/module/smod_equil_gold_hoyle.html gold_hoyle_eq module subroutine gold_hoyle_eq(settings, grid, background, physics)

Arguments

Type IntentOptional 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