smod_equil_harris_sheet Submodule

  1. smod_equil_harris_sheet.f08
  2. smod_equil_harris_sheet

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.


Uses


Contents

Functions

Module Procedures


Functions

function rho0()

Arguments

None

Return Value real(kind=dp)

function T0(x)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: x

Return Value real(kind=dp)

function dT0(x)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: x

Return Value real(kind=dp)

function B02(x)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: x

Return Value real(kind=dp)

function dB02(x)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: x

Return Value real(kind=dp)

function ddB02(x)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: x

Return Value real(kind=dp)

function B03(x)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: x

Return Value real(kind=dp)

function dB03(x)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: x

Return Value real(kind=dp)

function ddB03(x)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: x

Return Value real(kind=dp)

function B0(x)

Arguments

TypeIntentOptionalAttributesName
real(kind=dp), intent(in) :: x

Return Value real(kind=dp)


Module Procedures

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

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

Arguments

TypeIntentOptionalAttributesName
type(settings_t), intent(inout) :: settings
type(grid_t), intent(inout) :: grid
type(background_t), intent(inout) :: background
type(physics_t), intent(inout) :: physics