Contains a dedicated type for the various settings of the ARPACK solvers. All variables that are used in different solver settings are defined, initialised and set in this module.
Module defining the assertion routine used by assert.fpp.
Contains data for a realistic solar atmosphere model taken from Avrett, E. H., & Loeser, R. (2008). Models of the solar chromosphere and transition region from SUMER and HRTS observations: formation of the extreme-ultraviolet spectrum of hydrogen, carbon, and oxygen. ApJS, 175(1), 229. link.
Contains types and routines to handle banded matrices. We use the same conventions as explained in the LAPACK guide http://www.netlib.org/lapack/lug/node124.html.
Contains types and routines to handle banded Hermitian matrices. We use the same conventions as explained in the LAPACK guide http://www.netlib.org/lapack/lug/node124.html.
This module contains various methods to check for small, NaN or negative values, equal values or inf values. Interfaces are provided for functionality with real and complex variables.
log10 temperature values from Melemma and Lundqvist (2002), solar metallicity
Parent module governing all equilibrium types and submodules. This module contains all equilibrium types and the initial declarations of the module subroutines. Every equilibrium submodule extends this module, implementing one of the module subroutines declared here. All "main" equilibrium configurations are set in the submodules. The ones that depend on "main" arrays, like radiative cooling, are set here through calls to their respective modules.
This submodule defines a simple, adiabatic homogeneous medium in Cartesian geometry. The geometry can be overridden using the parfile.
This submodule defines an equilibrium in cylindrical geometry with a constant axial current. The geometry can be overridden using the parfile.
This submodule defines a magnetic flux tube embedded in a uniform magnetic environment. In this case the flux tube is under coronal conditions where the subscript e denotes the outer region. More specifically the equilibrium is defined as . The geometry can be overridden in the parfile, and is cylindrical by default for .
This submodule defines a steady plane Couette flow in a Cartesian geometry with flow and viscosity.
This submodule defines an equilibrium in cylindrical geometry with an axial current profile (), modelling a solar coronal loop in which discrete Alfvén waves are present. The geometry can be overridden in the parfile.
This submodule defines flow driven instabilities in a Cartesian geometry. This equilibrium can not be called explicitly from the parfile, but rather acts as a "parent setup" for the Rayleigh-Taylor and Kelvin-Helmholtz submodules which use this specific kind of equilibrium but with different parameters. This submodule is called within its implicit children.
This submodule defines a Gold-Hoyle equilibrium in cylindrical geometry. This equilibrium configuration models a filament with a uniform twist such that all fieldlines perform an equal amount of turns around the cylinder axis. The geometry can be overridden in the parfile.
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 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 submodule defines a resistive equilibrium with tearing modes created by a Harris sheet.
This submodule defines an exponentially stratified medium in Cartesian geometry with a constant gravity term and magnetic shear. The geometry can be overridden in the parfile.
This submodule defines internal kink modes in force-free magnetic fields. The geometry is cylindrical with parabolic density and velocity profiles, The geometry can be overridden in the parfile.
This submodule defines an equilibrium in Cartesian geometry with a stratified equilibrium profile, representing a solar magnetic atmosphere. The equilibrium is isothermal with a constant magnetic field. The scale height is given by The geometry is fixed to Cartesian, boundaries can be overridden using the parfile.
This submodule defines an unperturbed magnetised jet model in cylindrical geometry, giving rise to Kelvin-Helmholtz and current-driven instabilities. The geometry is fixed for this problem; the cylinder wall is dependent on the equilibrium parameters and is given by 2rj.
This submodule defines Kelvin-Helmholtz instabilities in Cartesian geometry. This equilibrium is a specific case of the flow driven instabilities. No magnetic fields are considered in this case (pure HD).
This submodule defines an equilibrium containing magnetothermal instabilities in a cylindrical geometry. The geometry can be overridden in the parfile.
This submodule defines magneto-rotational instabilities in an accretion disk. Due to the special nature of this equilibrium x_start is hardcoded to one and can not be overridden in the parfile, the same goes for the geometry which is hardcoded to 'cylindrical'. The outer edge can be chosen freely. This equilibrium is chosen in such a way that the angular rotation is of order unity, implying Keplerian rotation. The thin-disk approximation is valid with small magnetic fields, but still large enough to yield magneto-rotational instabilities. Gravity is assumed to go like .
This submodule defines a magnetic flux tube embedded in a uniform magnetic environment. In this case the flux tube is under photospheric conditions where the subscript e denotes the outer region. More specifically the equilibrium is defined as . The geometry can be overridden in the parfile, and is cylindrical by default for .
This submodule defines a simple, homogeneous medium in Cartesian geometry with a constant resistivity value. The geometry can be overridden using the parfile.
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 submodule defines an equilibrium in Cartesian geometry with a flow profile and 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 submodule defines an inhomogeneous medium in Cartesian geometry with a constant resistivity value. Two (constant) density profiles are defined which are connected by a sine profile, representing the interface inbetween. This density profile allows for resonant absorption, the geometry can be overridden in the parfile.
This submodule defines a cylindrical equilibrium, resembling a rotating plasma cylinder. The geometry can be overridden using the parfile.
This submodule defines Rayleigh-Taylor instabilities in Cartesian geometry. This equilibrium is a specific case of the flow driven instabilities.
This submodule defines Rayleigh-Taylor and Kelvin-Helmholtz instabilities in Cartesian geometry. This equilibrium is a specific case of the flow driven instabilities.
This submodule defines Rayleigh-Taylor instabilities in rotating theta pinches. The straight cylinder approximation is used with a constant angular frequency. Density and pressure profiles decrease over the domain, with a uni-directional increasing magnetic field profile. Mode numbers correspond to HD Rayleigh-Taylor instabilities, while represent MHD RTIs. The geometry is hardcoded to 'cylindrical', the domain is forced to through division by x_end.
This submodule defines a cylindrical equilibrium in which a Suydam surface is present, such that this gives rises to Suydam cluster modes. The geometry can be overriden using the parfile.
This submodule defines a steady Taylor-Couette flow in a cylindrical geometry where a fluid is confined between two (rotating) coaxial cylinders (without a magnetic field).
This submodule defines a steady Taylor-Couette flow in a cylindrical geometry where a plasma is confined between two (rotating) coaxial cylinders with an azimuthal magnetic field and constant resistivity.
Submodule for user-defined equilibria. Look at the examples in the equilibria subdirectory or consult the website for more information.
Module containing all equilibrium parameters. All parameters used in the equilibrium submodules are defined here for convenience, including the wave numbers and . All of these values are NaN initially, such that variables that are not properly set propagate their value and are easy to spot in follow-up checks.
Module to explicitly handle exceptions. Depending on the application at hand we override what happens when an exception is raised, which is useful for testing purposes (no error stop if we expect something to fail). Loosely based on an example given in https://github.com/Goddard-Fortran-Ecosystem/pFUnit_demos/blob/main/Basic/src/throw.F90
Module defining convenient index retrieval functions on various arrays.
Module containing all grid-related things. Contains subroutines to create the base grid, Gaussian grid and scale factors. An integral of in can be approximated with where and are the weights and nodes of the Gaussian quadrature. The Gaussian grid is hence set up in every interval across the nodes as
Module containing Hall-related routines. Sets the Hall and electron inertia factors based on normalisations and specified profiles.
Module to handle parfile reading. Contains subroutines to retrieve the parfile based on the commandline arguments and to read the parfile, setting the global variables.
Module to inspect if certain conditions are fulfilled by doing additional sanity checks on the equilibrium configuration. For cylindrical geometries we check if is an integer and if the on-axis values obey regularity conditions. Equilibrium balance for both the Cartesian and cylindrical cases is checked.
Module responsible for integration of differential equations, useful when setting equilibria or integrating the equilibrium equation. Contains subroutines to numerically solve the following systems of differential equations: These are solved using a fifth-order Runge-Kutta method.
Module responsible for table interpolations and array lookups. Contains subroutines for table interpolations, numerical derivatives of arrays and lookup functions. Subroutines are loosely based on routines implemented in the MPI-AMRVAC code.
Module containing functions to solve linear systems .
Module that contains the implementation of nodes in the linked-list matrix representation.
Module that contains the implementation of a single row of nodes in the linked-list matrix representation.
Module that contains the datastructure for a linked-list matrix representation.
This module handles formatting of terminal-printed strings. Contains subroutines to colourise strings.
All physical constants used in the code are defined in this module. We include values both in SI units and in cgs units for convenience. All values are taken from the NRL Plasma Formulary.
Module containing radiative cooling-related routines. This module is responsible for initialising the radiative cooling variables and a correct handling of the cooling curves. If an interpolated cooling curve is selected this module calls the interpolation module to create one.
Module containing resistivity-related routines, calculates and sets the resistivity values based on the equilibrium configuration.
Module to set a realistic solar atmosphere, using tabulated density and temperature profiles (see mod_atmosphere_curves), in Cartesian geometries only.
Parent module for everything solver-related. Interfaces to the different submodules are defined here, and the solve_evp routine calls the correct solver based on parfile settings.
Main module for the Arnoldi-type solvers. Contains interfaces to the general Arnoldi procedures (general, shift-invert, etc.).
Module containing the implementation for the ARPACK general-type solver, that is, given the general eigenvalue problem find eigenvalues that satisfy a given criterion.
Module containing the implementation for the ARPACK shift-invert-type solver, that is, given the general eigenvalue problem choose a shift and solve the problem thereby finding eigenvalues of the shifted problem that satisfy a given criterion.
Submodule containing the implementation of the inverse iteration algorithm. TODO more docs
Submodule containing the implementation of the QR-cholesky algorithm. Using LAPACKS's zpbtrf and BLAS's zgbtrs, the original problem is written as a standard eigenvalue problem as where is positive definite and . Eventually a call to LAPACK's zgeev routine is done to obtain all eigenvalues and eigenvectors.
Submodule containing the implementation of the QR-invert algorithm. The original problem is written as a standard eigenvalue problem through . This is done using a LU decomposition via LAPACKS's zgbsv. Eventually a call to LAPACK's zgeev routine is done to obtain all eigenvalues and eigenvectors.
Submodule containing the implementation of the QZ-direct algorithm. We keep the general form of the eigenvalue problem and solve this directly by calling LAPACK's zggev3 routine.
Module that calculates the finite element basis functions. The different basis functions and their derivatives used throughout the code are calculated here. All routines defined in this module simply return the basis functions for a specific point r in the interval (rj_lo, rj_hi).
This module is responsible for calculating and setting the thermal conduction values based on the equilibrium configuration.
Module to provide timing facilities.
Contains various subroutines and functions to switch between linked-list matrix representations, banded matrix representations, and full array matrices.
Simple module, containing only version-related stuff. Versioning is done in a separate module to avoid cluttering the commit history of for example mod_global_variables or mod_output every time an update to the code is done. The Legolas version is added to the datfile as a string and has common MAJOR.MINOR.PATCH formatting. This means: