smod_equil_gravito_mhd Submodule

This submodule defines an equilibrium in Cartesian geometry with a stratified equilibrium profile, giving rise to gravito-MHD waves. The geometry can be overridden using the parfile.

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

Note

Default values are given by

  • k2 =
  • k3 =
  • cte_p0 = 0.5 : used to set the pressure value.
  • alpha = 20 : used to constrain the density.
  • g = 0.5 : used to set the gravity constant.

and can all be changed in the parfile.


Uses


Variables

Type Visibility Attributes Name Initial
real(kind=dp), private :: B0

Functions

function rho0(x)

Arguments

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

Return Value real(kind=dp)

function drho0(x)

Arguments

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

Return Value real(kind=dp)

function T0()

Arguments

None

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 g0()

Arguments

None

Return Value real(kind=dp)


Module Procedures

module procedure /home/runner/work/legolas/legolas/docs/ford/module/smod_equil_gravito_mhd.html gravito_mhd_eq module subroutine gravito_mhd_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