smod_equil_coronal_flux_tube Submodule

This submodule defines a magnetic flux tube embedded in a uniform magnetic environment. In this case the flux tube is under coronal conditions where the subscript e denotes the outer region. More specifically the equilibrium is defined as . The geometry can be overridden in the parfile, and is cylindrical by default for .

This equilibrium is taken from chapter 6, fig. 6.7 in Roberts, Bernard (2019). MHD Waves in the Solar Atmosphere. Cambridge University Press. DOI.

Note

For best results, it is recommended to enable mesh accumulation.

Note

Default values are given by

  • k2 = 0
  • k3 = 2
  • cte_rho0 = 1 : density value for the inner tube.
  • cte_p0 = 1 : pressure value for the inner tube.
  • r0 = 1 : radius of the inner tube.

and can all be changed in the parfile.



Variables

Type Visibility Attributes Name Initial
real(kind=dp), private :: rho_e
real(kind=dp), private :: p_e
real(kind=dp), private :: B_0
real(kind=dp), private :: B_e

Functions

function rho0(r)

Arguments

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

Return Value real(kind=dp)

function T0(r)

Arguments

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

Return Value real(kind=dp)

function p0(r)

Arguments

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

Return Value real(kind=dp)

function B03(r)

Arguments

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

Return Value real(kind=dp)


Module Procedures

module procedure /home/runner/work/legolas/legolas/docs/ford/module/smod_equil_coronal_flux_tube.html coronal_flux_tube_eq module subroutine coronal_flux_tube_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