smod_equil_gravito_acoustic Submodule

This submodule defines an equilibrium in Cartesian geometry with a stratified equilibrium profile, giving rise to gravito-acoustic waves. No magnetic fields are included, such that this treats the hydrodynamic regime. The geometry can be overridden using the parfile.

This equilibrium is taken from section 7.2.3, p. 242 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 = 1 : used to set the pressure value.
  • alpha = 20.42 : used to constrain the density.
  • g = 0.5 : used to set the gravity constant.

and can all be changed in the parfile.


Uses


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

Arguments

None

Return Value real(kind=dp)


Module Procedures

module procedure /home/runner/work/legolas/legolas/docs/ford/module/smod_equil_gravito_acoustic.html gravito_acoustic_eq module subroutine gravito_acoustic_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