mod_physics_vars_srhd Module


Uses

  • module~~mod_physics_vars_srhd~~UsesGraph module~mod_physics_vars_srhd mod_physics_vars_srhd module~mod_global_parameters mod_global_parameters module~mod_physics_vars_srhd->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

Helium abundance over Hydrogen

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

Index of the energy density

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

derived values from adiabatic index

double precision, public :: gamma_to_gamma_1

derived values from adiabatic index

double precision, public :: inv_gamma_1

derived values from adiabatic index

integer, public :: lfac_

Index of the Lorentz factor

integer, public, allocatable :: mom(:)

Indices of the momentum density

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

Index of the gas pressure 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 = "srhd"

String describing the physics type of the simulation

integer, public :: rho_

Index of the density (in the w array) as primitive or conserved

logical, public :: srhd_eos = .false.

Whether synge eos is used

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

The adiabatic index

integer, public :: srhd_n_tracer = 0

Number of tracer species

logical, public :: srhd_particles = .false.

Whether particles module is added

double precision, public :: srhd_small_density = 0.0d0

Smallest rest-mass density allowed

double precision, public :: srhd_small_pressure = 0.0d0

Smallest gas pressure allowed. Both the lower bracket of the pressure root-find in con2prim and the floor fix_prim_state applies to a reconstructed face state. Zero, the default, leaves both unfloored.

logical, public :: srhd_source_usr = .false.

switch for source user

integer, public :: tracer(1)

Indices of the tracers

integer, public :: xi_

Index of the inertia