Hydrodynamics physics module. More...

Functions/Subroutines
subroutine, public	hd_phys_init ()
	Initialize the module.

subroutine, public	hd_check_params

subroutine, public	hd_check_w (primitive, ixil, ixol, w, flag)
	Returns logical argument flag where values are ok.

subroutine, public	hd_to_conserved (ixil, ixol, w, x)
	Transform primitive variables into conservative ones.

subroutine, public	hd_to_primitive (ixil, ixol, w, x)
	Transform conservative variables into primitive ones.

subroutine, public	hd_get_csound2 (w, x, ixil, ixol, csound2)
	Calculate the square of the thermal sound speed csound2 within ixO^L. csound2=gamma*p/rho.

subroutine, public	hd_get_pthermal (w, x, ixil, ixol, pth)
	Calculate thermal pressure=(gamma-1)(e-0.5m**2/rho) within ixO^L.

double precision function, dimension(ixo^s), public	hd_kin_en (w, ixil, ixol, inv_rho)

Variables
logical, public, protected	hd_energy = .true.
	Whether an energy equation is used.

logical, public, protected	hd_thermal_conduction = .false.
	Whether thermal conduction is added.

type(tc_fluid), allocatable, public	tc_fl

type(te_fluid), allocatable, public	te_fl_hd

logical, public, protected	hd_radiative_cooling = .false.
	Whether radiative cooling is added.

type(rc_fluid), allocatable, public	rc_fl

logical, public, protected	hd_dust = .false.
	Whether dust is added.

logical, public, protected	hd_viscosity = .false.
	Whether viscosity is added.

logical, public, protected	hd_gravity = .false.
	Whether gravity is added.

logical, public, protected	hd_particles = .false.
	Whether particles module is added.

logical, public, protected	hd_rotating_frame = .false.
	Whether rotating frame is activated.

logical, public, protected	hd_cak_force = .false.
	Whether CAK radiation line force is activated.

integer, public, protected	hd_n_tracer = 0
	Number of tracer species.

logical, public, protected	hd_partial_ionization = .false.
	Whether plasma is partially ionized.

integer, public, protected	rho_
	Index of the density (in the w array)

integer, dimension(:), allocatable, public, protected	mom
	Indices of the momentum density.

integer, dimension(:), allocatable, public, protected	tracer
	Indices of the tracers.

integer, public, protected	e_
	Index of the energy density (-1 if not present)

integer, public, protected	p_
	Index of the gas pressure (-1 if not present) should equal e_.

integer, public, protected	te_
	Indices of temperature.

integer, public, protected	tcoff_
	Index of the cutoff temperature for the TRAC method.

double precision, public	hd_gamma = 5.d0/3.0d0
	The adiabatic index.

double precision, public	hd_adiab = 1.0d0
	The adiabatic constant.

logical, public, protected	hd_trac = .false.
	Whether TRAC method is used.

integer, public, protected	hd_trac_type = 1

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

logical, public, protected	eq_state_units = .true.

Detailed Description

Hydrodynamics physics module.