mod_hd_phys Module Reference

Hydrodynamics physics module. More...

Functions/Subroutines
subroutine	hd_write_info (fh)
	Write this module's parameters to a snapsoht. More...
subroutine, public	hd_phys_init ()
	Initialize the module. More...
subroutine	hd_te_images
subroutine	hd_sts_set_source_tc_hd (ixIL, ixOL, w, x, wres, fix_conserve_at_step, my_dt, igrid, nflux)
double precision function	hd_get_tc_dt_hd (w, ixIL, ixOL, dxD, x)
subroutine	hd_tc_handle_small_e (w, x, ixIL, ixOL, step)
subroutine	tc_params_read_hd (fl)
subroutine	hd_get_rho (w, x, ixIL, ixOL, rho)
subroutine	rc_params_read (fl)
subroutine, public	hd_check_params
subroutine	hd_physical_units
subroutine, public	hd_check_w (primitive, ixIL, ixOL, w, flag)
	Returns logical argument flag where values are ok. More...
subroutine, public	hd_to_conserved (ixIL, ixOL, w, x)
	Transform primitive variables into conservative ones. More...
subroutine, public	hd_to_primitive (ixIL, ixOL, w, x)
	Transform conservative variables into primitive ones. More...
subroutine	hd_ei_to_e (ixIL, ixOL, w, x)
	Transform internal energy to total energy. More...
subroutine	hd_e_to_ei (ixIL, ixOL, w, x)
	Transform total energy to internal energy. More...
subroutine	e_to_rhos (ixIL, ixOL, w, x)
subroutine	rhos_to_e (ixIL, ixOL, w, x)
subroutine	hd_get_v_idim (w, x, ixIL, ixOL, idim, v)
	Calculate v_i = m_i / rho within ixO^L. More...
subroutine	hd_get_v (w, x, ixIL, ixOL, v)
	Calculate velocity vector v_i = m_i / rho within ixO^L. More...
subroutine	hd_get_cmax (w, x, ixIL, ixOL, idim, cmax)
	Calculate cmax_idim = csound + abs(v_idim) within ixO^L. More...
subroutine	hd_get_a2max (w, x, ixIL, ixOL, a2max)
subroutine	hd_get_tcutoff (ixIL, ixOL, w, x, tco_local, Tmax_local)
	get adaptive cutoff temperature for TRAC (Johnston 2019 ApJL, 873, L22) More...
subroutine	hd_get_cbounds (wLC, wRC, wLp, wRp, x, ixIL, ixOL, idim, Hspeed, cmax, cmin)
	Calculate cmax_idim = csound + abs(v_idim) within ixO^L. More...
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. More...
subroutine, public	hd_get_pthermal (w, x, ixIL, ixOL, pth)
	Calculate thermal pressure=(gamma-1)*(e-0.5*m**2/rho) within ixO^L. More...
subroutine	hd_get_temperature_from_etot (w, x, ixIL, ixOL, res)
	Calculate temperature=p/rho when in e_ the total energy is stored. More...
subroutine	hd_get_temperature_from_eint (w, x, ixIL, ixOL, res)
	Calculate temperature=p/rho when in e_ the internal energy is stored. More...
subroutine	hd_get_flux_cons (w, x, ixIL, ixOL, idim, f)
subroutine	hd_get_flux (wC, w, x, ixIL, ixOL, idim, f)
subroutine	hd_add_source_geom (qdt, dtfactor, ixIL, ixOL, wCT, wprim, w, x)
	Add geometrical source terms to w. More...
subroutine	hd_add_source (qdt, dtfactor, ixIL, ixOL, wCT, wCTprim, w, x, qsourcesplit, active)
subroutine	hd_get_dt (w, ixIL, ixOL, dtnew, dxD, x)
double precision function, dimension(ixo^s), public	hd_kin_en (w, ixIL, ixOL, inv_rho)
double precision function, dimension(ixo^s)	hd_inv_rho (w, ixIL, ixOL)
subroutine	hd_handle_small_values (primitive, w, x, ixIL, ixOL, subname)
subroutine	rfactor_from_temperature_ionization (w, x, ixIL, ixOL, Rfactor)
subroutine	rfactor_from_constant_ionization (w, x, ixIL, ixOL, Rfactor)
subroutine	hd_update_temperature (ixIL, ixOL, wCT, w, x)

Variables
logical, public, protected	hd_energy = .true.
	Whether an energy equation is used. More...
logical, public, protected	hd_thermal_conduction = .false.
	Whether thermal conduction is added. More...
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. More...
type(rc_fluid), allocatable, public	rc_fl
logical, public, protected	hd_dust = .false.
	Whether dust is added. More...
logical, public, protected	hd_viscosity = .false.
	Whether viscosity is added. More...
logical, public, protected	hd_gravity = .false.
	Whether gravity is added. More...
logical, public, protected	hd_particles = .false.
	Whether particles module is added. More...
logical, public, protected	hd_rotating_frame = .false.
	Whether rotating frame is activated. More...
logical, public, protected	hd_cak_force = .false.
	Whether CAK radiation line force is activated. More...
integer, public, protected	hd_n_tracer = 0
	Number of tracer species. More...
logical, public, protected	hd_partial_ionization = .false.
	Whether plasma is partially ionized. More...
integer, public, protected	rho_
	Index of the density (in the w array) More...
integer, dimension(:), allocatable, public, protected	mom
	Indices of the momentum density. More...
integer, dimension(:), allocatable, public, protected	tracer
	Indices of the tracers. More...
integer, public, protected	e_
	Index of the energy density (-1 if not present) More...
integer, public, protected	p_
	Index of the gas pressure (-1 if not present) should equal e_. More...
integer, public, protected	te_
	Indices of temperature. More...
integer, public, protected	tcoff_
	Index of the cutoff temperature for the TRAC method. More...
double precision, public	hd_gamma = 5.d0/3.0d0
	The adiabatic index. More...
double precision, public	hd_adiab = 1.0d0
	The adiabatic constant. More...
logical, public, protected	hd_trac = .false.
	Whether TRAC method is used. More...
integer, public, protected	hd_trac_type = 1
logical, public, protected	hd_force_diagonal = .false.
	Allows overruling default corner filling (for debug mode, since otherwise corner primitives fail) More...
double precision, public, protected	he_abundance =0.1d0
	Helium abundance over Hydrogen. More...
double precision, public, protected	h_ion_fr =1d0
	Ionization fraction of H H_ion_fr = H+/(H+ + H) More...
double precision, public, protected	he_ion_fr =1d0
	Ionization fraction of He He_ion_fr = (He2+ + He+)/(He2+ + He+ + He) More...
double precision, public, protected	he_ion_fr2 =1d0
	Ratio of number He2+ / number He+ + He2+ He_ion_fr2 = He2+/(He2+ + He+) More...
double precision, public, protected	rr =1d0
logical, public, protected	eq_state_units = .true.

Detailed Description

Hydrodynamics physics module.

Function/Subroutine Documentation

e_to_rhos()

subroutine mod_hd_phys::e_to_rhos	(	integer, intent(in)	ixI,
		integer, intent(in)	L,
		integer, intent(in)	ixO,
			L,
		double precision, dimension(ixi^s, nw)	w,
		double precision, dimension(ixi^s, 1:ndim), intent(in)	x
	)

Definition at line 733 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_add_source()

subroutine mod_hd_phys::hd_add_source	(	double precision, intent(in)	qdt,
		double precision, intent(in)	dtfactor,
		integer, intent(in)	ixI,
		integer, intent(in)	L,
		integer, intent(in)	ixO,
			L,
		double precision, dimension(ixi^s, 1:nw), intent(in)	wCT,
		double precision, dimension(ixi^s,1:nw), intent(in)	wCTprim,
		double precision, dimension(ixi^s, 1:nw), intent(inout)	w,
		double precision, dimension(ixi^s, 1:ndim), intent(in)	x,
		logical, intent(in)	qsourcesplit,
		logical, intent(inout)	active
	)

Definition at line 1331 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_add_source_geom()

subroutine mod_hd_phys::hd_add_source_geom	(	double precision, intent(in)	qdt,
		double precision, intent(in)	dtfactor,
		integer, intent(in)	ixI,
		integer, intent(in)	L,
		integer, intent(in)	ixO,
			L,
		double precision, dimension(ixi^s, 1:nw), intent(inout)	wCT,
		double precision, dimension(ixi^s,1:nw), intent(inout)	wprim,
		double precision, dimension(ixi^s, 1:nw), intent(inout)	w,
		double precision, dimension(ixi^s, 1:ndim), intent(in)	x
	)

Add geometrical source terms to w.

Notice that the expressions of the geometrical terms depend only on ndir, not ndim. Eg, they are the same in 2.5D and in 3D, for any geometry.

Ileyk : to do :

• address the source term for the dust in case (coordinate == spherical)

Definition at line 1190 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_check_params()

subroutine, public mod_hd_phys::hd_check_params

Definition at line 537 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_check_w()

subroutine, public mod_hd_phys::hd_check_w	(	logical, intent(in)	primitive,
		integer, intent(in)	ixI,
		integer, intent(in)	L,
		integer, intent(in)	ixO,
			L,
		double precision, dimension(ixi^s, nw), intent(in)	w,
		logical, dimension(ixi^s,1:nw), intent(inout)	flag
	)

Returns logical argument flag where values are ok.

Definition at line 613 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_e_to_ei()

subroutine mod_hd_phys::hd_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 721 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_ei_to_e()

subroutine mod_hd_phys::hd_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 708 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_get_a2max()

subroutine mod_hd_phys::hd_get_a2max	(	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(ndim), intent(inout)	a2max
	)

4th order

Definition at line 811 of file mod_hd_phys.t.

hd_get_cbounds()

subroutine mod_hd_phys::hd_get_cbounds	(	double precision, dimension(ixi^s, nw), intent(in)	wLC,
		double precision, dimension(ixi^s, nw), intent(in)	wRC,
		double precision, dimension(ixi^s, nw), intent(in)	wLp,
		double precision, dimension(ixi^s, nw), intent(in)	wRp,
		double precision, dimension(ixi^s, 1:ndim), intent(in)	x,
		integer, intent(in)	ixI,
		integer, intent(in)	L,
		integer, intent(in)	ixO,
			L,
		integer, intent(in)	idim,
		double precision, dimension(ixi^s,1:number_species), intent(in)	Hspeed,
		double precision, dimension(ixi^s,1:number_species), intent(inout)	cmax,
		double precision, dimension(ixi^s,1:number_species), intent(inout), optional	cmin
	)

Calculate cmax_idim = csound + abs(v_idim) within ixO^L.

Definition at line 887 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_get_cmax()

subroutine mod_hd_phys::hd_get_cmax	(	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,
		integer, intent(in)	idim,
		double precision, dimension(ixi^s), intent(inout)	cmax
	)

Calculate cmax_idim = csound + abs(v_idim) within ixO^L.

Definition at line 790 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_get_csound2()

subroutine, public mod_hd_phys::hd_get_csound2	(	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), intent(out)	csound2
	)

Calculate the square of the thermal sound speed csound2 within ixO^L. csound2=gamma*p/rho.

Definition at line 1001 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_get_dt()

subroutine mod_hd_phys::hd_get_dt	(	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, intent(inout)	dtnew,
		double precision, intent(in)	dx,
		double precision, intent(in)	D,
		double precision, dimension(ixi^s, 1:^nd), intent(in)	x
	)

Definition at line 1394 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_get_flux()

subroutine mod_hd_phys::hd_get_flux	(	double precision, dimension(ixi^s, 1:nw), intent(in)	wC,
		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,
		integer, intent(in)	idim,
		double precision, dimension(ixi^s, nwflux), intent(out)	f
	)

Definition at line 1132 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_get_flux_cons()

subroutine mod_hd_phys::hd_get_flux_cons	(	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,
		integer, intent(in)	idim,
		double precision, dimension(ixi^s, nwflux), intent(out)	f
	)

Definition at line 1093 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_get_pthermal()

subroutine, public mod_hd_phys::hd_get_pthermal	(	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
	)

Calculate thermal pressure=(gamma-1)*(e-0.5*m**2/rho) within ixO^L.

Definition at line 1014 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_get_rho()

subroutine mod_hd_phys::hd_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 481 of file mod_hd_phys.t.

hd_get_tc_dt_hd()

double precision function mod_hd_phys::hd_get_tc_dt_hd	(	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, intent(in)	dx,
		double precision, intent(in)	D,
		double precision, dimension(ixi^s,1:ndim), intent(in)	x
	)

Definition at line 412 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_get_tcutoff()

subroutine mod_hd_phys::hd_get_tcutoff	(	integer, intent(in)	ixI,
		integer, intent(in)	L,
		integer, intent(in)	ixO,
			L,
		double precision, dimension(ixi^s,1:nw), intent(inout)	w,
		double precision, dimension(ixi^s,1:ndim), intent(in)	x,
		double precision, intent(out)	tco_local,
		double precision, intent(out)	Tmax_local
	)

get adaptive cutoff temperature for TRAC (Johnston 2019 ApJL, 873, L22)

iijima et al. 2021, LTRAC method

Definition at line 834 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_get_temperature_from_eint()

subroutine mod_hd_phys::hd_get_temperature_from_eint	(	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)	res
	)

Calculate temperature=p/rho when in e_ the internal energy is stored.

Definition at line 1079 of file mod_hd_phys.t.

hd_get_temperature_from_etot()

subroutine mod_hd_phys::hd_get_temperature_from_etot	(	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)	res
	)

Calculate temperature=p/rho when in e_ the total energy is stored.

Definition at line 1064 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_get_v()

subroutine mod_hd_phys::hd_get_v	(	double precision, dimension(ixi^s,nw), intent(in)	w,
		double precision, dimension(ixi^s,1:^nd), intent(in)	x,
		integer, intent(in)	ixI,
		integer, intent(in)	L,
		integer, intent(in)	ixO,
			L,
		double precision, dimension(ixi^s,1:ndir), intent(out)	v
	)

Calculate velocity vector v_i = m_i / rho within ixO^L.

Definition at line 774 of file mod_hd_phys.t.

hd_get_v_idim()

subroutine mod_hd_phys::hd_get_v_idim	(	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,
		integer, intent(in)	idim,
		double precision, dimension(ixi^s), intent(out)	v
	)

Calculate v_i = m_i / rho within ixO^L.

Definition at line 764 of file mod_hd_phys.t.

hd_handle_small_values()

subroutine mod_hd_phys::hd_handle_small_values	(	logical, intent(in)	primitive,
		double precision, dimension(ixi^s,1:nw), intent(inout)	w,
		double precision, dimension(ixi^s,1:ndim), intent(in)	x,
		integer, intent(in)	ixI,
		integer, intent(in)	L,
		integer, intent(in)	ixO,
			L,
		character(len=*), intent(in)	subname
	)

Definition at line 1455 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_inv_rho()

double precision function, dimension(ixo^s) mod_hd_phys::hd_inv_rho	(	double precision, dimension(ixi^s, nw), intent(in)	w,
		integer, intent(in)	ixI,
		integer, intent(in)	L,
		integer, intent(in)	ixO,
			L
	)

Definition at line 1445 of file mod_hd_phys.t.

hd_kin_en()

double precision function, dimension(ixo^s), public mod_hd_phys::hd_kin_en	(	double precision, dimension(ixi^s, nw), intent(in)	w,
		integer, intent(in)	ixI,
		integer, intent(in)	L,
		integer, intent(in)	ixO,
			L,
		double precision, dimension(ixo^s), intent(in), optional	inv_rho
	)

Definition at line 1431 of file mod_hd_phys.t.

hd_phys_init()

subroutine, public mod_hd_phys::hd_phys_init

Initialize the module.

Definition at line 156 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_physical_units()

subroutine mod_hd_phys::hd_physical_units

Definition at line 560 of file mod_hd_phys.t.

hd_sts_set_source_tc_hd()

subroutine mod_hd_phys::hd_sts_set_source_tc_hd	(	integer, intent(in)	ixI,
		integer, intent(in)	L,
		integer, intent(in)	ixO,
			L,
		double precision, dimension(ixi^s,1:nw), intent(inout)	w,
		double precision, dimension(ixi^s,1:ndim), intent(in)	x,
		double precision, dimension(ixi^s,1:nw), intent(inout)	wres,
		logical, intent(in)	fix_conserve_at_step,
		double precision, intent(in)	my_dt,
		integer, intent(in)	igrid,
		integer, intent(in)	nflux
	)

Definition at line 400 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_tc_handle_small_e()

subroutine mod_hd_phys::hd_tc_handle_small_e	(	double precision, dimension(ixi^s,1:nw), intent(inout)	w,
		double precision, dimension(ixi^s,1:ndim), intent(in)	x,
		integer, intent(in)	ixI,
		integer, intent(in)	L,
		integer, intent(in)	ixO,
			L,
		integer, intent(in)	step
	)

Definition at line 427 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_te_images()

subroutine mod_hd_phys::hd_te_images

Definition at line 380 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_to_conserved()

subroutine, public mod_hd_phys::hd_to_conserved	(	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 primitive variables into conservative ones.

Definition at line 642 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_to_primitive()

subroutine, public mod_hd_phys::hd_to_primitive	(	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 conservative variables into primitive ones.

Definition at line 674 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_update_temperature()

subroutine mod_hd_phys::hd_update_temperature	(	integer, intent(in)	ixI,
		integer, intent(in)	L,
		integer, intent(in)	ixO,
			L,
		double precision, dimension(ixi^s,1:nw), intent(in)	wCT,
		double precision, dimension(ixi^s,1:nw), intent(inout)	w,
		double precision, dimension(ixi^s,1:ndim), intent(in)	x
	)

Definition at line 1580 of file mod_hd_phys.t.

Here is the call graph for this function:

hd_write_info()

subroutine mod_hd_phys::hd_write_info

(

integer, intent(in)

fh

)

Write this module's parameters to a snapsoht.

Definition at line 138 of file mod_hd_phys.t.

rc_params_read()

subroutine mod_hd_phys::rc_params_read

(

type(rc_fluid), intent(inout)

fl

)

Name of cooling curve

Name of cooling method

Fixed temperature not lower than tlow

Lower limit of temperature

Add cooling source in a split way (.true.) or un-split way (.false.)

Definition at line 493 of file mod_hd_phys.t.

rfactor_from_constant_ionization()

subroutine mod_hd_phys::rfactor_from_constant_ionization	(	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)	Rfactor
	)

Definition at line 1569 of file mod_hd_phys.t.

rfactor_from_temperature_ionization()

subroutine mod_hd_phys::rfactor_from_temperature_ionization	(	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)	Rfactor
	)

Definition at line 1553 of file mod_hd_phys.t.

Here is the call graph for this function:

rhos_to_e()

subroutine mod_hd_phys::rhos_to_e	(	integer, intent(in)	ixI,
		integer, intent(in)	L,
		integer, intent(in)	ixO,
			L,
		double precision, dimension(ixi^s, nw)	w,
		double precision, dimension(ixi^s, 1:ndim), intent(in)	x
	)

Definition at line 748 of file mod_hd_phys.t.

Here is the call graph for this function:

tc_params_read_hd()

subroutine mod_hd_phys::tc_params_read_hd

(

type(tc_fluid), intent(inout)

fl

)

Definition at line 462 of file mod_hd_phys.t.

Variable Documentation

e_

integer, public, protected mod_hd_phys::e_

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

Definition at line 57 of file mod_hd_phys.t.

eq_state_units

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

Definition at line 102 of file mod_hd_phys.t.

h_ion_fr

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

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

Definition at line 88 of file mod_hd_phys.t.

hd_adiab

double precision, public mod_hd_phys::hd_adiab = 1.0d0

The adiabatic constant.

Definition at line 72 of file mod_hd_phys.t.

hd_cak_force

logical, public, protected mod_hd_phys::hd_cak_force = .false.

Whether CAK radiation line force is activated.

Definition at line 39 of file mod_hd_phys.t.

hd_dust

logical, public, protected mod_hd_phys::hd_dust = .false.

Whether dust is added.

Definition at line 24 of file mod_hd_phys.t.

hd_energy

logical, public, protected mod_hd_phys::hd_energy = .true.

Whether an energy equation is used.

Definition at line 12 of file mod_hd_phys.t.

hd_force_diagonal

logical, public, protected mod_hd_phys::hd_force_diagonal = .false.

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

Definition at line 82 of file mod_hd_phys.t.

hd_gamma

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

The adiabatic index.

Definition at line 69 of file mod_hd_phys.t.

hd_gravity

logical, public, protected mod_hd_phys::hd_gravity = .false.

Whether gravity is added.

Definition at line 30 of file mod_hd_phys.t.

hd_n_tracer

integer, public, protected mod_hd_phys::hd_n_tracer = 0

Number of tracer species.

Definition at line 42 of file mod_hd_phys.t.

hd_partial_ionization

logical, public, protected mod_hd_phys::hd_partial_ionization = .false.

Whether plasma is partially ionized.

Definition at line 45 of file mod_hd_phys.t.

hd_particles

logical, public, protected mod_hd_phys::hd_particles = .false.

Whether particles module is added.

Definition at line 33 of file mod_hd_phys.t.

hd_radiative_cooling

logical, public, protected mod_hd_phys::hd_radiative_cooling = .false.

Whether radiative cooling is added.

Definition at line 20 of file mod_hd_phys.t.

hd_rotating_frame

logical, public, protected mod_hd_phys::hd_rotating_frame = .false.

Whether rotating frame is activated.

Definition at line 36 of file mod_hd_phys.t.

hd_thermal_conduction

logical, public, protected mod_hd_phys::hd_thermal_conduction = .false.

Whether thermal conduction is added.

Definition at line 15 of file mod_hd_phys.t.

hd_trac

logical, public, protected mod_hd_phys::hd_trac = .false.

Whether TRAC method is used.

Definition at line 78 of file mod_hd_phys.t.

hd_trac_type

integer, public, protected mod_hd_phys::hd_trac_type = 1

Definition at line 79 of file mod_hd_phys.t.

hd_viscosity

logical, public, protected mod_hd_phys::hd_viscosity = .false.

Whether viscosity is added.

Definition at line 27 of file mod_hd_phys.t.

he_abundance

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

Helium abundance over Hydrogen.

Definition at line 85 of file mod_hd_phys.t.

he_ion_fr

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

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

Definition at line 91 of file mod_hd_phys.t.

he_ion_fr2

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

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

Definition at line 94 of file mod_hd_phys.t.

mom

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

Indices of the momentum density.

Definition at line 51 of file mod_hd_phys.t.

p_

integer, public, protected mod_hd_phys::p_

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

Definition at line 60 of file mod_hd_phys.t.

rc_fl

type(rc_fluid), allocatable, public mod_hd_phys::rc_fl

Definition at line 21 of file mod_hd_phys.t.

rho_

integer, public, protected mod_hd_phys::rho_

Index of the density (in the w array)

Definition at line 48 of file mod_hd_phys.t.

rr

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

Definition at line 98 of file mod_hd_phys.t.

tc_fl

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

Definition at line 16 of file mod_hd_phys.t.

tcoff_

integer, public, protected mod_hd_phys::tcoff_

Index of the cutoff temperature for the TRAC method.

Definition at line 66 of file mod_hd_phys.t.

te_

integer, public, protected mod_hd_phys::te_

Indices of temperature.

Definition at line 63 of file mod_hd_phys.t.

te_fl_hd

type(te_fluid), allocatable, public mod_hd_phys::te_fl_hd

Definition at line 17 of file mod_hd_phys.t.

tracer

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

Indices of the tracers.

Definition at line 54 of file mod_hd_phys.t.