mod_rmhd_phys Module Reference

Radiation-magneto-hydrodynamics module. More...

Functions/Subroutines
subroutine, public	rmhd_phys_init ()
subroutine, public	rmhd_set_mg_bounds
	Set the boundaries for the diffusion of E.
subroutine, public	rmhd_ei_to_e (ixil, ixol, w, x)
	Transform internal energy to total energy.
subroutine, public	rmhd_e_to_ei (ixil, ixol, w, x)
	Transform total energy to internal energy.
subroutine, public	rmhd_get_v (w, x, ixil, ixol, v)
	Calculate v vector.
subroutine, public	rmhd_get_pradiation (w, x, ixil, ixol, prad, nth)
	Calculate radiation pressure within ixO^L.
subroutine, public	rmhd_get_pthermal_plus_pradiation (w, x, ixil, ixol, pth_plus_prad)
	Calculates the sum of the gas pressure and the max Prad tensor element.
subroutine, public	rmhd_get_tgas (w, x, ixil, ixol, tgas)
	Calculates gas temperature.
subroutine, public	rmhd_get_trad (w, x, ixil, ixol, trad)
	Calculates radiation temperature.
subroutine, public	rmhd_get_rho (w, x, ixil, ixol, rho)
subroutine, public	get_normalized_divb (w, ixil, ixol, divb)
	get dimensionless div B = |divB| * volume / area / |B|
subroutine, public	get_current (w, ixil, ixol, idirmin, current)
	Calculate idirmin and the idirmin:3 components of the common current array make sure that dxlevel(^D) is set correctly.
double precision function, dimension(ixo^s), public	rmhd_mag_en_all (w, ixil, ixol)
	Compute 2 times total magnetic energy.
subroutine, public	rmhd_clean_divb_multigrid (qdt, qt, active)
subroutine, public	rmhd_face_to_center (ixol, s)
	calculate cell-center values from face-center values
subroutine, public	b_from_vector_potential (ixisl, ixil, ixol, ws, x)
	calculate magnetic field from vector potential

Variables
double precision, public	rmhd_gamma = 5.d0/3.0d0
	The adiabatic index.
double precision, public	rmhd_adiab = 1.0d0
	The adiabatic constant.
double precision, public	rmhd_eta = 0.0d0
	The MHD resistivity.
double precision, public	rmhd_eta_hyper = 0.0d0
	The MHD hyper-resistivity.
double precision, public	rmhd_etah = 0.0d0
	Hall resistivity.
double precision, public, protected	small_r_e = 0.d0
	The smallest allowed radiation energy.
double precision, public, protected	rmhd_trac_mask = 0.d0
	Height of the mask used in the TRAC method.
double precision, public	rmhd_glm_alpha = 0.5d0
	GLM-MHD parameter: ratio of the diffusive and advective time scales for div b taking values within [0, 1].
double precision, public	hypertc_kappa
	The thermal conductivity kappa in hyperbolic thermal conduction.
double precision, public, protected	he_abundance =0.1d0
	Helium abundance over Hydrogen.
double precision, public, protected	h_ion_fr =1d0
	Ionization fraction of H H_ion_fr = H+/(H+ + H)
double precision, public, protected	he_ion_fr =1d0
	Ionization fraction of He He_ion_fr = (He2+ + He+)/(He2+ + He+ + He)
double precision, public, protected	he_ion_fr2 =1d0
	Ratio of number He2+ / number He+ + He2+ He_ion_fr2 = He2+/(He2+ + He+)
double precision, public, protected	rr =1d0
integer, public	equi_rho0_ = -1
	equi vars indices in the stateequi_vars array
integer, public	equi_pe0_ = -1
integer, public, protected	rmhd_n_tracer = 0
	Number of tracer species.
integer, public, protected	rho_
	Index of the density (in the w array)
integer, dimension(:), allocatable, public, protected	mom
	Indices of the momentum density.
integer, public, protected	c
	Indices of the momentum density for the form of better vectorization.
integer, public, protected	m
integer, public, protected	c_
integer, public, protected	e_
	Index of the energy density (-1 if not present)
integer, public, protected	r_e
	Index of the radiation energy.
integer, public, protected	b
integer, public, protected	p_
	Index of the gas pressure (-1 if not present) should equal e_.
integer, public, protected	q_
	Index of the heat flux q.
integer, public, protected	psi_
	Indices of the GLM psi.
integer, public, protected	te_
	Indices of temperature.
integer, public, protected	tcoff_
	Index of the cutoff temperature for the TRAC method.
integer, public, protected	tweight_
integer, dimension(:), allocatable, public, protected	tracer
	Indices of the tracers.
integer, dimension(2 *^nd), public, protected	boundary_divbfix_skip =0
	To skip * layer of ghost cells during divB=0 fix for boundary.
logical, public, protected	rmhd_energy = .true.
	Whether an energy equation is used.
logical, public, protected	rmhd_thermal_conduction = .false.
	Whether thermal conduction is used.
logical, public, protected	rmhd_hyperbolic_thermal_conduction = .false.
	Whether thermal conduction is used.
logical, public, protected	rmhd_viscosity = .false.
	Whether viscosity is added.
logical, public, protected	rmhd_gravity = .false.
	Whether gravity is added.
logical, public, protected	rmhd_particles = .false.
	Whether particles module is added.
logical, public, protected	rmhd_glm = .false.
	Whether GLM-MHD is used to control div B.
logical, public, protected	rmhd_glm_extended = .true.
	Whether extended GLM-MHD is used with additional sources.
logical, public, protected	rmhd_trac = .false.
	Whether TRAC method is used.
integer, public, protected	rmhd_trac_type =1
	Which TRAC method is used.
integer, public, protected	rmhd_trac_finegrid =4
	Distance between two adjacent traced magnetic field lines (in finest cell size)
logical, public, protected	source_split_divb = .false.
	Whether divB cleaning sources are added splitting from fluid solver.
logical, public, protected	rmhd_partial_ionization = .false.
	Whether plasma is partially ionized.
logical, public, protected	rmhd_cak_force = .false.
	Whether CAK radiation line force is activated.
logical, public, protected	rmhd_4th_order = .false.
	MHD fourth order.
logical, public	has_equi_rho0 = .false.
	whether split off equilibrium density
logical, public	has_equi_pe0 = .false.
	whether split off equilibrium thermal pressure
logical, public	rmhd_equi_thermal = .false.
logical, public, protected	rmhd_dump_full_vars = .false.
	whether dump full variables (when splitting is used) in a separate dat file
integer, public, protected	rmhd_divb_nth = 1
	Whether divB is computed with a fourth order approximation.
logical, public	divbwave = .true.
	Add divB wave in Roe solver.
logical, public	clean_initial_divb = .false.
	clean initial divB
character(len=8), public	rmhd_radiation_formalism = 'fld'
	Formalism to treat radiation.
character(len=8), public	rmhd_pressure = 'Trad'
	In the case of no rmhd_energy, how to compute pressure.
logical, public, protected	rmhd_radiation_force = .true.
	Treat radiation fld_Rad_force.
logical, public, protected	rmhd_energy_interact = .true.
	Treat radiation-gas energy interaction.
logical, public, protected	rmhd_radiation_diffusion = .true.
	Treat radiation energy diffusion.
logical, public, protected	rmhd_radiation_advection = .true.
	Treat radiation advection.
double precision, public	kbmpmua4
	kb/(m_p mu)* 1/a_rad**4,
logical, public, protected	eq_state_units = .true.
logical, dimension(2 *^nd), public, protected	boundary_divbfix =.true.
	To control divB=0 fix for boundary.
logical, public, protected	b0field_forcefree =.true.
	B0 field is force-free.
logical, public	partial_energy = .false.
	Whether an internal or hydrodynamic energy equation is used.
character(len=std_len), public, protected	typedivbfix = 'linde'
	Method type to clean divergence of B.
character(len=std_len), public, protected	type_ct = 'uct_contact'
	Method type of constrained transport.
type(tc_fluid), allocatable, public	tc_fl
	type of fluid for thermal conduction
type(te_fluid), allocatable, public	te_fl_rmhd
	type of fluid for thermal emission synthesis
procedure(sub_convert), pointer, public	rmhd_to_primitive => null()
procedure(sub_convert), pointer, public	rmhd_to_conserved => null()
procedure(sub_get_pthermal), pointer, public	rmhd_get_pthermal => null()
procedure(sub_get_pthermal), pointer, public	rmhd_get_temperature => null()

Detailed Description

Radiation-magneto-hydrodynamics module.

Function/Subroutine Documentation

b_from_vector_potential()

subroutine, public mod_rmhd_phys::b_from_vector_potential	(	integer, intent(in)	ixis,
		integer, intent(in)	l,
		integer, intent(in)	ixi,
			l,
		integer, intent(in)	ixo,
			l,
		double precision, dimension(ixis^s,1:nws), intent(inout)	ws,
		double precision, dimension(ixi^s,1:ndim), intent(in)	x
	)

calculate magnetic field from vector potential

Definition at line 4894 of file mod_rmhd_phys.t.

Here is the call graph for this function:

get_current()

subroutine, public mod_rmhd_phys::get_current	(	double precision, dimension(ixi^s,1:nw), intent(in)	w,
		integer, intent(in)	ixi,
		integer, intent(in)	l,
		integer, intent(in)	ixo,
			l,
		integer, intent(out)	idirmin,
		double precision, dimension(ixi^s,7-2*ndir:3)	current
	)

Calculate idirmin and the idirmin:3 components of the common current array make sure that dxlevel(^D) is set correctly.

Definition at line 3303 of file mod_rmhd_phys.t.

Here is the call graph for this function:

get_normalized_divb()

subroutine, public mod_rmhd_phys::get_normalized_divb	(	double precision, dimension(ixi^s,1:nw), intent(in)	w,
		integer, intent(in)	ixi,
		integer, intent(in)	l,
		integer, intent(in)	ixo,
			l,
		double precision, dimension(ixi^s)	divb
	)

get dimensionless div B = |divB| * volume / area / |B|

Definition at line 3273 of file mod_rmhd_phys.t.

Here is the call graph for this function:

rmhd_clean_divb_multigrid()

subroutine, public mod_rmhd_phys::rmhd_clean_divb_multigrid	(	double precision, intent(in)	qdt,
		double precision, intent(in)	qt,
		logical, intent(inout)	active
	)

Parameters

[in]	qdt	Current time step
[in]	qt	Current time
[in,out]	active	Output if the source is active

Definition at line 4133 of file mod_rmhd_phys.t.

Here is the call graph for this function:

rmhd_e_to_ei()

subroutine, public mod_rmhd_phys::rmhd_e_to_ei	(	integer, intent(in)	ixi,
		integer, intent(in)	l,
		integer, intent(in)	ixo,
			l,
		double precision, dimension(ixi^s, nw), intent(inout)	w,
		double precision, dimension(ixi^s, 1:ndim), intent(in)	x
	)

Transform total energy to internal energy.

Definition at line 1257 of file mod_rmhd_phys.t.

rmhd_ei_to_e()

subroutine, public mod_rmhd_phys::rmhd_ei_to_e	(	integer, intent(in)	ixi,
		integer, intent(in)	l,
		integer, intent(in)	ixo,
			l,
		double precision, dimension(ixi^s, nw), intent(inout)	w,
		double precision, dimension(ixi^s, 1:ndim), intent(in)	x
	)

Transform internal energy to total energy.

Definition at line 1230 of file mod_rmhd_phys.t.

rmhd_face_to_center()

subroutine, public mod_rmhd_phys::rmhd_face_to_center	(	integer, intent(in)	ixo,
		integer, intent(in)	l,
		type(state)	s
	)

calculate cell-center values from face-center values

Definition at line 4832 of file mod_rmhd_phys.t.

rmhd_get_pradiation()

subroutine, public mod_rmhd_phys::rmhd_get_pradiation	(	double precision, dimension(ixi^s, 1:nw), intent(in)	w,
		double precision, dimension(ixi^s, 1:ndim), intent(in)	x,
		integer, intent(in)	ixi,
		integer, intent(in)	l,
		integer, intent(in)	ixo,
			l,
		double precision, dimension(ixo^s, 1:ndim, 1:ndim), intent(out)	prad,
		integer, intent(in)	nth
	)

Calculate radiation pressure within ixO^L.

Definition at line 2178 of file mod_rmhd_phys.t.

Here is the call graph for this function:

rmhd_get_pthermal_plus_pradiation()

subroutine, public mod_rmhd_phys::rmhd_get_pthermal_plus_pradiation	(	double precision, dimension(ixi^s, 1:nw), intent(in)	w,
		double precision, dimension(ixi^s, 1:ndim), intent(in)	x,
		integer, intent(in)	ixi,
		integer, intent(in)	l,
		integer, intent(in)	ixo,
			l,
		double precision, dimension(ixi^s), intent(out)	pth_plus_prad
	)

Calculates the sum of the gas pressure and the max Prad tensor element.

filter cmax

Definition at line 2198 of file mod_rmhd_phys.t.

Here is the call graph for this function:

rmhd_get_rho()

subroutine, public mod_rmhd_phys::rmhd_get_rho	(	double precision, dimension(ixi^s,1:nw), intent(in)	w,
		double precision, dimension(ixi^s,1:ndim), intent(in)	x,
		integer, intent(in)	ixi,
		integer, intent(in)	l,
		integer, intent(in)	ixo,
			l,
		double precision, dimension(ixi^s), intent(out)	rho
	)

Definition at line 2762 of file mod_rmhd_phys.t.

Here is the call graph for this function:

rmhd_get_tgas()

subroutine, public mod_rmhd_phys::rmhd_get_tgas	(	double precision, dimension(ixi^s, 1:nw), intent(in)	w,
		double precision, dimension(ixi^s, 1:ndim), intent(in)	x,
		integer, intent(in)	ixi,
		integer, intent(in)	l,
		integer, intent(in)	ixo,
			l,
		double precision, dimension(ixi^s), intent(out)	tgas
	)

Calculates gas temperature.

Definition at line 2247 of file mod_rmhd_phys.t.

rmhd_get_trad()

subroutine, public mod_rmhd_phys::rmhd_get_trad	(	double precision, dimension(ixi^s, 1:nw), intent(in)	w,
		double precision, dimension(ixi^s, 1:ndim), intent(in)	x,
		integer, intent(in)	ixi,
		integer, intent(in)	l,
		integer, intent(in)	ixo,
			l,
		double precision, dimension(ixi^s), intent(out)	trad
	)

Calculates radiation temperature.

Definition at line 2260 of file mod_rmhd_phys.t.

Here is the call graph for this function:

rmhd_get_v()

subroutine, public mod_rmhd_phys::rmhd_get_v	(	double precision, dimension(ixi^s,nw), intent(in)	w,
		double precision, dimension(ixi^s,1:ndim), intent(in)	x,
		integer, intent(in)	ixi,
		integer, intent(in)	l,
		integer, intent(in)	ixo,
			l,
		double precision, dimension(ixi^s,ndir), intent(out)	v
	)

Calculate v vector.

Definition at line 1384 of file mod_rmhd_phys.t.

Here is the call graph for this function:

rmhd_mag_en_all()

double precision function, dimension(ixo^s), public mod_rmhd_phys::rmhd_mag_en_all	(	double precision, dimension(ixi^s, nw), intent(in)	w,
		integer, intent(in)	ixi,
		integer, intent(in)	l,
		integer, intent(in)	ixo,
			l
	)

Compute 2 times total magnetic energy.

Definition at line 3657 of file mod_rmhd_phys.t.

rmhd_phys_init()

subroutine, public mod_rmhd_phys::rmhd_phys_init

set radiation energy

Initiate radiation-closure module

Definition at line 277 of file mod_rmhd_phys.t.

Here is the call graph for this function:

rmhd_set_mg_bounds()

subroutine, public mod_rmhd_phys::rmhd_set_mg_bounds

Set the boundaries for the diffusion of E.

Definition at line 922 of file mod_rmhd_phys.t.

Variable Documentation

b

integer, public, protected mod_rmhd_phys::b

Definition at line 74 of file mod_rmhd_phys.t.

b0field_forcefree

logical, public, protected mod_rmhd_phys::b0field_forcefree =.true.

B0 field is force-free.

Definition at line 176 of file mod_rmhd_phys.t.

boundary_divbfix

logical, dimension(2*^nd), public, protected mod_rmhd_phys::boundary_divbfix =.true.

To control divB=0 fix for boundary.

Definition at line 174 of file mod_rmhd_phys.t.

boundary_divbfix_skip

integer, dimension(2*^nd), public, protected mod_rmhd_phys::boundary_divbfix_skip =0

To skip * layer of ghost cells during divB=0 fix for boundary.

Definition at line 93 of file mod_rmhd_phys.t.

c

integer public protected mod_rmhd_phys::c

Indices of the momentum density for the form of better vectorization.

Indices of the magnetic field for the form of better vectorization.

Definition at line 68 of file mod_rmhd_phys.t.

c_

integer public protected mod_rmhd_phys::c_

Definition at line 68 of file mod_rmhd_phys.t.

clean_initial_divb

logical, public mod_rmhd_phys::clean_initial_divb = .false.

clean initial divB

Definition at line 149 of file mod_rmhd_phys.t.

divbwave

logical, public mod_rmhd_phys::divbwave = .true.

Add divB wave in Roe solver.

Definition at line 147 of file mod_rmhd_phys.t.

e_

integer, public, protected mod_rmhd_phys::e_

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

Definition at line 70 of file mod_rmhd_phys.t.

eq_state_units

logical, public, protected mod_rmhd_phys::eq_state_units = .true.

Definition at line 172 of file mod_rmhd_phys.t.

equi_pe0_

integer, public mod_rmhd_phys::equi_pe0_ = -1

Definition at line 60 of file mod_rmhd_phys.t.

equi_rho0_

integer, public mod_rmhd_phys::equi_rho0_ = -1

equi vars indices in the stateequi_vars array

Definition at line 59 of file mod_rmhd_phys.t.

h_ion_fr

double precision, public, protected mod_rmhd_phys::h_ion_fr =1d0

Ionization fraction of H H_ion_fr = H+/(H+ + H)

Definition at line 43 of file mod_rmhd_phys.t.

has_equi_pe0

logical, public mod_rmhd_phys::has_equi_pe0 = .false.

whether split off equilibrium thermal pressure

Definition at line 138 of file mod_rmhd_phys.t.

has_equi_rho0

logical, public mod_rmhd_phys::has_equi_rho0 = .false.

whether split off equilibrium density

Definition at line 136 of file mod_rmhd_phys.t.

he_abundance

double precision, public, protected mod_rmhd_phys::he_abundance =0.1d0

Helium abundance over Hydrogen.

Definition at line 40 of file mod_rmhd_phys.t.

he_ion_fr

double precision, public, protected mod_rmhd_phys::he_ion_fr =1d0

Ionization fraction of He He_ion_fr = (He2+ + He+)/(He2+ + He+ + He)

Definition at line 46 of file mod_rmhd_phys.t.

he_ion_fr2

double precision, public, protected mod_rmhd_phys::he_ion_fr2 =1d0

Ratio of number He2+ / number He+ + He2+ He_ion_fr2 = He2+/(He2+ + He+)

Definition at line 49 of file mod_rmhd_phys.t.

hypertc_kappa

double precision, public mod_rmhd_phys::hypertc_kappa

The thermal conductivity kappa in hyperbolic thermal conduction.

Definition at line 36 of file mod_rmhd_phys.t.

kbmpmua4

double precision, public mod_rmhd_phys::kbmpmua4

kb/(m_p mu)* 1/a_rad**4,

Definition at line 166 of file mod_rmhd_phys.t.

m

integer, public, protected mod_rmhd_phys::m

Definition at line 68 of file mod_rmhd_phys.t.

mom

integer, dimension(:), allocatable, public, protected mod_rmhd_phys::mom

Indices of the momentum density.

Definition at line 66 of file mod_rmhd_phys.t.

p_

integer, public, protected mod_rmhd_phys::p_

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

Definition at line 76 of file mod_rmhd_phys.t.

partial_energy

logical, public mod_rmhd_phys::partial_energy = .false.

Whether an internal or hydrodynamic energy equation is used.

Definition at line 180 of file mod_rmhd_phys.t.

psi_

integer, public, protected mod_rmhd_phys::psi_

Indices of the GLM psi.

Definition at line 80 of file mod_rmhd_phys.t.

q_

integer, public, protected mod_rmhd_phys::q_

Index of the heat flux q.

Definition at line 78 of file mod_rmhd_phys.t.

r_e

integer, public, protected mod_rmhd_phys::r_e

Index of the radiation energy.

Definition at line 72 of file mod_rmhd_phys.t.

rho_

integer, public, protected mod_rmhd_phys::rho_

Index of the density (in the w array)

Definition at line 64 of file mod_rmhd_phys.t.

rmhd_4th_order

logical, public, protected mod_rmhd_phys::rmhd_4th_order = .false.

MHD fourth order.

Definition at line 134 of file mod_rmhd_phys.t.

rmhd_adiab

double precision, public mod_rmhd_phys::rmhd_adiab = 1.0d0

The adiabatic constant.

Definition at line 19 of file mod_rmhd_phys.t.

rmhd_cak_force

logical, public, protected mod_rmhd_phys::rmhd_cak_force = .false.

Whether CAK radiation line force is activated.

Definition at line 132 of file mod_rmhd_phys.t.

rmhd_divb_nth

integer, public, protected mod_rmhd_phys::rmhd_divb_nth = 1

Whether divB is computed with a fourth order approximation.

Definition at line 143 of file mod_rmhd_phys.t.

rmhd_dump_full_vars

logical, public, protected mod_rmhd_phys::rmhd_dump_full_vars = .false.

whether dump full variables (when splitting is used) in a separate dat file

Definition at line 141 of file mod_rmhd_phys.t.

rmhd_energy

logical, public, protected mod_rmhd_phys::rmhd_energy = .true.

Whether an energy equation is used.

Definition at line 106 of file mod_rmhd_phys.t.

rmhd_energy_interact

logical, public, protected mod_rmhd_phys::rmhd_energy_interact = .true.

Treat radiation-gas energy interaction.

Definition at line 157 of file mod_rmhd_phys.t.

rmhd_equi_thermal

logical, public mod_rmhd_phys::rmhd_equi_thermal = .false.

Definition at line 139 of file mod_rmhd_phys.t.

rmhd_eta

double precision, public mod_rmhd_phys::rmhd_eta = 0.0d0

The MHD resistivity.

Definition at line 21 of file mod_rmhd_phys.t.

rmhd_eta_hyper

double precision, public mod_rmhd_phys::rmhd_eta_hyper = 0.0d0

The MHD hyper-resistivity.

Definition at line 23 of file mod_rmhd_phys.t.

rmhd_etah

double precision, public mod_rmhd_phys::rmhd_etah = 0.0d0

Hall resistivity.

Definition at line 25 of file mod_rmhd_phys.t.

rmhd_gamma

double precision, public mod_rmhd_phys::rmhd_gamma = 5.d0/3.0d0

The adiabatic index.

Definition at line 17 of file mod_rmhd_phys.t.

rmhd_get_pthermal

procedure(sub_get_pthermal), pointer, public mod_rmhd_phys::rmhd_get_pthermal => null()

Definition at line 199 of file mod_rmhd_phys.t.

rmhd_get_temperature

procedure(sub_get_pthermal), pointer, public mod_rmhd_phys::rmhd_get_temperature => null()

Definition at line 201 of file mod_rmhd_phys.t.

rmhd_glm

logical, public, protected mod_rmhd_phys::rmhd_glm = .false.

Whether GLM-MHD is used to control div B.

Definition at line 118 of file mod_rmhd_phys.t.

rmhd_glm_alpha

double precision, public mod_rmhd_phys::rmhd_glm_alpha = 0.5d0

GLM-MHD parameter: ratio of the diffusive and advective time scales for div b taking values within [0, 1].

Definition at line 34 of file mod_rmhd_phys.t.

rmhd_glm_extended

logical, public, protected mod_rmhd_phys::rmhd_glm_extended = .true.

Whether extended GLM-MHD is used with additional sources.

Definition at line 120 of file mod_rmhd_phys.t.

rmhd_gravity

logical, public, protected mod_rmhd_phys::rmhd_gravity = .false.

Whether gravity is added.

Definition at line 114 of file mod_rmhd_phys.t.

rmhd_hyperbolic_thermal_conduction

logical, public, protected mod_rmhd_phys::rmhd_hyperbolic_thermal_conduction = .false.

Whether thermal conduction is used.

Definition at line 110 of file mod_rmhd_phys.t.

rmhd_n_tracer

integer, public, protected mod_rmhd_phys::rmhd_n_tracer = 0

Number of tracer species.

Definition at line 62 of file mod_rmhd_phys.t.

rmhd_partial_ionization

logical, public, protected mod_rmhd_phys::rmhd_partial_ionization = .false.

Whether plasma is partially ionized.

Definition at line 130 of file mod_rmhd_phys.t.

rmhd_particles

logical, public, protected mod_rmhd_phys::rmhd_particles = .false.

Whether particles module is added.

Definition at line 116 of file mod_rmhd_phys.t.

rmhd_pressure

character(len=8), public mod_rmhd_phys::rmhd_pressure = 'Trad'

In the case of no rmhd_energy, how to compute pressure.

Definition at line 153 of file mod_rmhd_phys.t.

rmhd_radiation_advection

logical, public, protected mod_rmhd_phys::rmhd_radiation_advection = .true.

Treat radiation advection.

Definition at line 161 of file mod_rmhd_phys.t.

rmhd_radiation_diffusion

logical, public, protected mod_rmhd_phys::rmhd_radiation_diffusion = .true.

Treat radiation energy diffusion.

Definition at line 159 of file mod_rmhd_phys.t.

rmhd_radiation_force

logical, public, protected mod_rmhd_phys::rmhd_radiation_force = .true.

Treat radiation fld_Rad_force.

Definition at line 155 of file mod_rmhd_phys.t.

rmhd_radiation_formalism

character(len=8), public mod_rmhd_phys::rmhd_radiation_formalism = 'fld'

Formalism to treat radiation.

Definition at line 151 of file mod_rmhd_phys.t.

rmhd_thermal_conduction

logical, public, protected mod_rmhd_phys::rmhd_thermal_conduction = .false.

Whether thermal conduction is used.

Definition at line 108 of file mod_rmhd_phys.t.

rmhd_to_conserved

procedure(sub_convert), pointer, public mod_rmhd_phys::rmhd_to_conserved => null()

Definition at line 197 of file mod_rmhd_phys.t.

rmhd_to_primitive

procedure(sub_convert), pointer, public mod_rmhd_phys::rmhd_to_primitive => null()

Definition at line 196 of file mod_rmhd_phys.t.

rmhd_trac

logical, public, protected mod_rmhd_phys::rmhd_trac = .false.

Whether TRAC method is used.

Definition at line 122 of file mod_rmhd_phys.t.

rmhd_trac_finegrid

integer, public, protected mod_rmhd_phys::rmhd_trac_finegrid =4

Distance between two adjacent traced magnetic field lines (in finest cell size)

Definition at line 126 of file mod_rmhd_phys.t.

rmhd_trac_mask

double precision, public, protected mod_rmhd_phys::rmhd_trac_mask = 0.d0

Height of the mask used in the TRAC method.

Definition at line 31 of file mod_rmhd_phys.t.

rmhd_trac_type

integer, public, protected mod_rmhd_phys::rmhd_trac_type =1

Which TRAC method is used.

Definition at line 124 of file mod_rmhd_phys.t.

rmhd_viscosity

logical, public, protected mod_rmhd_phys::rmhd_viscosity = .false.

Whether viscosity is added.

Definition at line 112 of file mod_rmhd_phys.t.

rr

double precision, public, protected mod_rmhd_phys::rr =1d0

Definition at line 53 of file mod_rmhd_phys.t.

small_r_e

double precision, public, protected mod_rmhd_phys::small_r_e = 0.d0

The smallest allowed radiation energy.

Definition at line 29 of file mod_rmhd_phys.t.

source_split_divb

logical, public, protected mod_rmhd_phys::source_split_divb = .false.

Whether divB cleaning sources are added splitting from fluid solver.

Definition at line 128 of file mod_rmhd_phys.t.

tc_fl

type(tc_fluid), allocatable, public mod_rmhd_phys::tc_fl

type of fluid for thermal conduction

Definition at line 192 of file mod_rmhd_phys.t.

tcoff_

integer, public, protected mod_rmhd_phys::tcoff_

Index of the cutoff temperature for the TRAC method.

Definition at line 84 of file mod_rmhd_phys.t.

te_

integer, public, protected mod_rmhd_phys::te_

Indices of temperature.

Definition at line 82 of file mod_rmhd_phys.t.

te_fl_rmhd

type(te_fluid), allocatable, public mod_rmhd_phys::te_fl_rmhd

type of fluid for thermal emission synthesis

Definition at line 194 of file mod_rmhd_phys.t.

tracer

integer, dimension(:), allocatable, public, protected mod_rmhd_phys::tracer

Indices of the tracers.

Definition at line 87 of file mod_rmhd_phys.t.

tweight_

integer, public, protected mod_rmhd_phys::tweight_

Definition at line 85 of file mod_rmhd_phys.t.

type_ct

character(len=std_len), public, protected mod_rmhd_phys::type_ct = 'uct_contact'

Method type of constrained transport.

Definition at line 188 of file mod_rmhd_phys.t.

typedivbfix

character(len=std_len), public, protected mod_rmhd_phys::typedivbfix = 'linde'

Method type to clean divergence of B.

Definition at line 186 of file mod_rmhd_phys.t.