mod_physics_vars_hd Module


Uses

  • module~~mod_physics_vars_hd~~UsesGraph module~mod_physics_vars_hd mod_physics_vars_hd module~mod_global_parameters mod_global_parameters module~mod_physics_vars_hd->module~mod_global_parameters module~mod_basic_types mod_basic_types module~mod_global_parameters->module~mod_basic_types module~mod_connectivity mod_connectivity module~mod_global_parameters->module~mod_connectivity module~mod_constants mod_constants module~mod_global_parameters->module~mod_constants module~mod_physicaldata mod_physicaldata module~mod_global_parameters->module~mod_physicaldata module~mod_variables mod_variables module~mod_global_parameters->module~mod_variables mpi mpi module~mod_global_parameters->mpi module~mod_variables->module~mod_basic_types

Variables

Type Visibility Attributes Name Initial
double precision, public :: He_abundance = 0.1d0

The helium abundance

integer, public, parameter :: dp = kind(0.0d0)
integer, public :: e_

Index of the energy density (-1 if not present)

integer, public, parameter :: flux_default = 0

Indicates a normal flux

integer, public, parameter :: flux_tvdlf = 1

Indicates the flux should be treated with tvdlf

integer, public, allocatable :: flux_type(:,:)

Array per direction per variable, which can be used to specify that certain fluxes have to be treated differently

double precision, public :: gamma_1
double precision, public :: hd_adiab = 1.0d0

The adiabatic constant

logical, public :: hd_energy = .true.

Whether an energy equation is used

logical, public :: hd_force_diagonal = .true.

Allows overruling default corner filling (for debug mode, since otherwise corner primitives fail)

double precision, public :: hd_gamma = 5.d0/3.0d0

The adiabatic index

logical, public :: hd_gravity = .false.

Whether to use gravity

integer, public :: hd_n_tracer = 0

Number of tracer species

logical, public :: hd_partial_ionization = .false.

Whether plasma is partially ionized

logical, public :: hd_particles = .false.

Whether particles module is added

logical, public :: hd_radiative_cooling = .false.

switch for radiative cooling

double precision, public :: inv_gamma_1
integer, public, allocatable :: mom(:)

Indices of the momentum density

integer, public, parameter :: nw_flux = 2+ndim+1
integer, public, parameter :: nw_phys = 2+ndim+1
integer, public :: p_

Index of the gas pressure (-1 if not present) should equal e_

logical, public :: phys_energy = .false.

Solve energy equation or not

logical, public :: phys_equi_pe = .false.

if equilibrium pressure is splitted

double precision, public :: phys_gamma = 5.d0/3.d0
logical, public :: phys_internal_e = .false.

Solve internal energy instead of total energy

logical, public :: phys_partial_ionization = .false.

Solve partially ionized one-fluid plasma

logical, public :: phys_req_diagonal = .true.

Whether the physics routines require diagonal ghost cells, for example for computing a curl.

logical, public :: phys_total_energy = .false.

Solve total energy equation or not

integer, public :: phys_wider_stencil = 0

To use wider stencils in flux calculations. A value of 1 will extend it by one cell in both directions, in any dimension

character(len=name_len), public :: physics_type = "hd"

String describing the physics type of the simulation

integer, public :: rho_

Index of the density (in the w array)

integer, public :: tracer(1)

Indices of the tracers