smod_equil_harris_sheet Submodule

This submodule defines a resistive equilibrium with tearing modes created by a Harris sheet.

This equilibrium is taken from Shi et al. (2020), Oblique tearing mode instability: guide field and Hall effect. The Astrophysical Journal, 902:142. DOI.

Note

Default values are given by

  • k2 = 0.12
  • k3 = 0
  • cte_rho0 = 1
  • cte_T0 = 1
  • cte_B02 = 1
  • cte_B03 = 0 : guide field parameter.
  • alpha = 1 : used to set the width of the current sheet.
  • eq_bool = False : if True, an alternative force-free Harris sheet is used.

and can all be changed in the parfile.


Uses


Functions

function rho0()

Arguments

None

Return Value real(kind=dp)

function T0(x)

Arguments

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

Return Value real(kind=dp)

function dT0(x)

Arguments

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

Return Value real(kind=dp)

function B02(x)

Arguments

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

Return Value real(kind=dp)

function dB02(x)

Arguments

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

Return Value real(kind=dp)

function ddB02(x)

Arguments

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

Return Value real(kind=dp)

function B03(x)

Arguments

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

Return Value real(kind=dp)

function dB03(x)

Arguments

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

Return Value real(kind=dp)

function ddB03(x)

Arguments

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

Return Value real(kind=dp)

function B0(x)

Arguments

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

Return Value real(kind=dp)


Module Procedures

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

eq_bool >> if True, the alternative force-free Harris sheet is used

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