MPI-AMRVAC
3.1
The MPI - Adaptive Mesh Refinement - Versatile Advection Code (development version)
|
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, 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. |
Hydrodynamics physics module.
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 | ||
) |
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 | ||
) |
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) | 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 :
Definition at line 1191 of file mod_hd_phys.t.
subroutine, public mod_hd_phys::hd_check_params |
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 612 of file mod_hd_phys.t.
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 720 of file mod_hd_phys.t.
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 707 of file mod_hd_phys.t.
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 810 of file mod_hd_phys.t.
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 886 of file mod_hd_phys.t.
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 789 of file mod_hd_phys.t.
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 1000 of file mod_hd_phys.t.
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 | ||
) |
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 | ||
) |
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 | ||
) |
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 1013 of file mod_hd_phys.t.
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 480 of file mod_hd_phys.t.
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 | ||
) |
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 833 of file mod_hd_phys.t.
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 1080 of file mod_hd_phys.t.
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 1065 of file mod_hd_phys.t.
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 773 of file mod_hd_phys.t.
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 763 of file mod_hd_phys.t.
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 | ||
) |
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 1424 of file mod_hd_phys.t.
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 1410 of file mod_hd_phys.t.
subroutine, public mod_hd_phys::hd_phys_init |
Initialize the module.
Definition at line 155 of file mod_hd_phys.t.
subroutine mod_hd_phys::hd_physical_units |
Definition at line 559 of file mod_hd_phys.t.
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 | ||
) |
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 | ||
) |
subroutine mod_hd_phys::hd_te_images |
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 641 of file mod_hd_phys.t.
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 673 of file mod_hd_phys.t.
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 | ||
) |
subroutine mod_hd_phys::hd_write_info | ( | integer, intent(in) | fh | ) |
Write this module's parameters to a snapsoht.
Definition at line 137 of file mod_hd_phys.t.
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 492 of file mod_hd_phys.t.
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 1550 of file mod_hd_phys.t.
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 | ||
) |
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 | ||
) |
subroutine mod_hd_phys::tc_params_read_hd | ( | type(tc_fluid), intent(inout) | fl | ) |
Definition at line 460 of file mod_hd_phys.t.
integer, public, protected mod_hd_phys::e_ |
Index of the energy density (-1 if not present)
Definition at line 56 of file mod_hd_phys.t.
logical, public, protected mod_hd_phys::eq_state_units = .true. |
Definition at line 101 of file mod_hd_phys.t.
double precision, public, protected mod_hd_phys::h_ion_fr =1d0 |
Ionization fraction of H H_ion_fr = H+/(H+ + H)
Definition at line 87 of file mod_hd_phys.t.
double precision, public mod_hd_phys::hd_adiab = 1.0d0 |
The adiabatic constant.
Definition at line 71 of file mod_hd_phys.t.
logical, public, protected mod_hd_phys::hd_cak_force = .false. |
Whether CAK radiation line force is activated.
Definition at line 38 of file mod_hd_phys.t.
logical, public, protected mod_hd_phys::hd_dust = .false. |
Whether dust is added.
Definition at line 23 of file mod_hd_phys.t.
logical, public, protected mod_hd_phys::hd_energy = .true. |
Whether an energy equation is used.
Definition at line 11 of file mod_hd_phys.t.
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 81 of file mod_hd_phys.t.
double precision, public mod_hd_phys::hd_gamma = 5.d0/3.0d0 |
The adiabatic index.
Definition at line 68 of file mod_hd_phys.t.
logical, public, protected mod_hd_phys::hd_gravity = .false. |
Whether gravity is added.
Definition at line 29 of file mod_hd_phys.t.
integer, public, protected mod_hd_phys::hd_n_tracer = 0 |
Number of tracer species.
Definition at line 41 of file mod_hd_phys.t.
logical, public, protected mod_hd_phys::hd_partial_ionization = .false. |
Whether plasma is partially ionized.
Definition at line 44 of file mod_hd_phys.t.
logical, public, protected mod_hd_phys::hd_particles = .false. |
Whether particles module is added.
Definition at line 32 of file mod_hd_phys.t.
logical, public, protected mod_hd_phys::hd_radiative_cooling = .false. |
Whether radiative cooling is added.
Definition at line 19 of file mod_hd_phys.t.
logical, public, protected mod_hd_phys::hd_rotating_frame = .false. |
Whether rotating frame is activated.
Definition at line 35 of file mod_hd_phys.t.
logical, public, protected mod_hd_phys::hd_thermal_conduction = .false. |
Whether thermal conduction is added.
Definition at line 14 of file mod_hd_phys.t.
logical, public, protected mod_hd_phys::hd_trac = .false. |
Whether TRAC method is used.
Definition at line 77 of file mod_hd_phys.t.
integer, public, protected mod_hd_phys::hd_trac_type = 1 |
Definition at line 78 of file mod_hd_phys.t.
logical, public, protected mod_hd_phys::hd_viscosity = .false. |
Whether viscosity is added.
Definition at line 26 of file mod_hd_phys.t.
double precision, public, protected mod_hd_phys::he_abundance =0.1d0 |
Helium abundance over Hydrogen.
Definition at line 84 of file mod_hd_phys.t.
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 90 of file mod_hd_phys.t.
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 93 of file mod_hd_phys.t.
integer, dimension(:), allocatable, public, protected mod_hd_phys::mom |
Indices of the momentum density.
Definition at line 50 of file mod_hd_phys.t.
integer, public, protected mod_hd_phys::p_ |
Index of the gas pressure (-1 if not present) should equal e_.
Definition at line 59 of file mod_hd_phys.t.
type(rc_fluid), allocatable, public mod_hd_phys::rc_fl |
Definition at line 20 of file mod_hd_phys.t.
integer, public, protected mod_hd_phys::rho_ |
Index of the density (in the w array)
Definition at line 47 of file mod_hd_phys.t.
double precision, public, protected mod_hd_phys::rr =1d0 |
Definition at line 97 of file mod_hd_phys.t.
type(tc_fluid), allocatable, public mod_hd_phys::tc_fl |
Definition at line 15 of file mod_hd_phys.t.
integer, public, protected mod_hd_phys::tcoff_ |
Index of the cutoff temperature for the TRAC method.
Definition at line 65 of file mod_hd_phys.t.
integer, public, protected mod_hd_phys::te_ |
Indices of temperature.
Definition at line 62 of file mod_hd_phys.t.
type(te_fluid), allocatable, public mod_hd_phys::te_fl_hd |
Definition at line 16 of file mod_hd_phys.t.
integer, dimension(:), allocatable, public, protected mod_hd_phys::tracer |
Indices of the tracers.
Definition at line 53 of file mod_hd_phys.t.