smod_equil_resistive_tearing Submodule

This submodule defines an equilibrium in Cartesian geometry with a constant resistivity value. Parameters are taken in such a way as to allow for resistive tearing modes. The geometry can be overridden using the parfile.

This equilibrium is taken from section 14.3, p. 551 in Goedbloed, H., Keppens, R., & Poedts, S. (2019). Magnetohydrodynamics of Laboratory and Astrophysical Plasmas. Cambridge University Press. DOI.

Contents

Functions

rho0 T0 B02 dB02 ddB02 B03 dB03 ddB03 B0

Module Procedures

resistive_tearing_modes_eq

Functions

function rho0()

Arguments

None

Return Value real(kind=dp)

function T0(x)

Arguments

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

Return Value real(kind=dp)

function B02(x)

Arguments

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

Return Value real(kind=dp)

function dB02(x)

Arguments

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

Return Value real(kind=dp)

function ddB02(x)

Arguments

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

Return Value real(kind=dp)

function B03(x)

Arguments

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

Return Value real(kind=dp)

function dB03(x)

Arguments

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

Return Value real(kind=dp)

function ddB03(x)

Arguments

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

Return Value real(kind=dp)

function B0(x)

Arguments

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

Return Value real(kind=dp)

Module Procedures

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

Sets the equilibrium

Arguments

Type	Intent	Optional	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