MPI-AMRVAC  3.1
The MPI - Adaptive Mesh Refinement - Versatile Advection Code (development version)
Data Types | Modules | Functions/Subroutines | Variables
mod_twofl_phys.t File Reference

Go to the source code of this file.

Data Types

interface  mod_twofl_phys::mask_subroutine2
 

Modules

module  mod_twofl_phys
 Magneto-hydrodynamics module.
 

Functions/Subroutines

subroutine mod_twofl_phys::twofl_init_hyper (files)
 
subroutine mod_twofl_phys::twofl_write_info (fh)
 Write this module's parameters to a snapsoht. More...
 
subroutine, public mod_twofl_phys::twofl_phys_init ()
 
subroutine mod_twofl_phys::twofl_te_images
 
subroutine mod_twofl_phys::twofl_sts_set_source_tc_c_mhd (ixIL, ixOL, w, x, wres, fix_conserve_at_step, my_dt, igrid, nflux)
 
subroutine mod_twofl_phys::twofl_sts_set_source_tc_c_hd (ixIL, ixOL, w, x, wres, fix_conserve_at_step, my_dt, igrid, nflux)
 
double precision function mod_twofl_phys::twofl_get_tc_dt_mhd_c (w, ixIL, ixOL, dxD, x)
 
double precision function mod_twofl_phys::twofl_get_tc_dt_hd_c (w, ixIL, ixOL, dxD, x)
 
subroutine mod_twofl_phys::twofl_tc_handle_small_e_c (w, x, ixIL, ixOL, step)
 
subroutine mod_twofl_phys::twofl_sts_set_source_tc_n_hd (ixIL, ixOL, w, x, wres, fix_conserve_at_step, my_dt, igrid, nflux)
 
subroutine mod_twofl_phys::twofl_tc_handle_small_e_n (w, x, ixIL, ixOL, step)
 
double precision function mod_twofl_phys::twofl_get_tc_dt_hd_n (w, ixIL, ixOL, dxD, x)
 
subroutine mod_twofl_phys::tc_n_params_read_hd (fl)
 
subroutine mod_twofl_phys::rc_params_read_n (fl)
 
subroutine mod_twofl_phys::tc_c_params_read_mhd (fl)
 
subroutine mod_twofl_phys::tc_c_params_read_hd (fl)
 
subroutine mod_twofl_phys::rc_params_read_c (fl)
 
subroutine mod_twofl_phys::set_equi_vars_grid_faces (igrid, x, ixIL, ixOL)
 sets the equilibrium variables More...
 
subroutine mod_twofl_phys::set_equi_vars_grid (igrid)
 sets the equilibrium variables More...
 
double precision function, dimension(ixo^s, 1:nwc) mod_twofl_phys::convert_vars_splitting (ixIL, ixOL, w, x, nwc)
 
subroutine mod_twofl_phys::grav_params_read (files)
 copied from mod_gravity More...
 
subroutine mod_twofl_phys::associate_dump_hyper ()
 
subroutine mod_twofl_phys::twofl_check_params
 
subroutine mod_twofl_phys::twofl_physical_units ()
 
subroutine mod_twofl_phys::twofl_check_w (primitive, ixIL, ixOL, w, flag)
 
subroutine, public mod_twofl_phys::twofl_to_conserved (ixIL, ixOL, w, x)
 Transform primitive variables into conservative ones. More...
 
subroutine, public mod_twofl_phys::twofl_to_primitive (ixIL, ixOL, w, x)
 Transform conservative variables into primitive ones. More...
 
subroutine mod_twofl_phys::twofl_ei_to_e_c (ixIL, ixOL, w, x)
 Transform internal energy to total energy. More...
 
subroutine mod_twofl_phys::twofl_e_to_ei_c (ixIL, ixOL, w, x)
 Transform total energy to internal energy. More...
 
subroutine mod_twofl_phys::twofl_ei_to_e_n (ixIL, ixOL, w, x)
 
subroutine mod_twofl_phys::twofl_e_to_ei_n (ixIL, ixOL, w, x)
 Transform total energy to internal energy. More...
 
subroutine mod_twofl_phys::twofl_handle_small_values (primitive, w, x, ixIL, ixOL, subname)
 
subroutine mod_twofl_phys::twofl_get_cmax (w, x, ixIL, ixOL, idim, cmax)
 Calculate cmax_idim=csound+abs(v_idim) within ixO^L. More...
 
subroutine mod_twofl_phys::twofl_get_a2max (w, x, ixIL, ixOL, a2max)
 
subroutine mod_twofl_phys::twofl_get_tcutoff_n (ixIL, ixOL, w, x, tco_local, Tmax_local)
 get adaptive cutoff temperature for TRAC (Johnston 2019 ApJL, 873, L22) More...
 
subroutine mod_twofl_phys::twofl_get_tcutoff_c (ixIL, ixOL, w, x, Tco_local, Tmax_local)
 get adaptive cutoff temperature for TRAC (Johnston 2019 ApJL, 873, L22) More...
 
subroutine mod_twofl_phys::twofl_get_h_speed_one (wprim, x, ixIL, ixOL, idim, Hspeed)
 get H speed for H-correction to fix the carbuncle problem at grid-aligned shock front More...
 
subroutine mod_twofl_phys::twofl_get_h_speed_species (wprim, x, ixIL, ixOL, idim, Hspeed)
 get H speed for H-correction to fix the carbuncle problem at grid-aligned shock front More...
 
subroutine mod_twofl_phys::twofl_get_cbounds_one (wLC, wRC, wLp, wRp, x, ixIL, ixOL, idim, Hspeed, cmax, cmin)
 Estimating bounds for the minimum and maximum signal velocities. More...
 
subroutine mod_twofl_phys::twofl_get_csound_prim_c (w, x, ixIL, ixOL, idim, csound)
 Calculate fast magnetosonic wave speed. More...
 
subroutine mod_twofl_phys::twofl_get_csound_prim_n (w, x, ixIL, ixOL, idim, csound)
 Calculate fast magnetosonic wave speed. More...
 
subroutine mod_twofl_phys::twofl_get_cbounds_species (wLC, wRC, wLp, wRp, x, ixIL, ixOL, idim, Hspeed, cmax, cmin)
 Estimating bounds for the minimum and maximum signal velocities. More...
 
subroutine mod_twofl_phys::twofl_get_ct_velocity (vcts, wLp, wRp, ixIL, ixOL, idim, cmax, cmin)
 prepare velocities for ct methods More...
 
subroutine mod_twofl_phys::twofl_get_csound_c_idim (w, x, ixIL, ixOL, idim, csound)
 
subroutine mod_twofl_phys::twofl_get_csound_prim (w, x, ixIL, ixOL, idim, csound)
 Calculate fast magnetosonic wave speed when cbounds_species=false. More...
 
subroutine twofl_get_csound2_primitive (w, x, ixIL, ixOL, csound2)
 
subroutine mod_twofl_phys::twofl_get_csound2 (w, x, ixIL, ixOL, csound2)
 
subroutine mod_twofl_phys::twofl_get_csound2_adiab (w, x, ixIL, ixOL, csound2)
 
subroutine mod_twofl_phys::twofl_get_csound (w, x, ixIL, ixOL, idim, csound)
 
subroutine mod_twofl_phys::twofl_get_csound2_from_pthermal (w, x, ixIL, ixOL, pth_c, pth_n, csound2)
 
subroutine mod_twofl_phys::twofl_get_csound_n (w, x, ixIL, ixOL, csound)
 
subroutine mod_twofl_phys::twofl_get_temperature_from_eint_n (w, x, ixIL, ixOL, res)
 separate routines so that it is faster Calculate temperature=p/rho when in e_ the internal energy is stored More...
 
subroutine mod_twofl_phys::twofl_get_temperature_from_eint_n_with_equi (w, x, ixIL, ixOL, res)
 
subroutine mod_twofl_phys::twofl_get_temperature_n_equi (w, x, ixIL, ixOL, res)
 
subroutine mod_twofl_phys::twofl_get_rho_n_equi (w, x, ixIL, ixOL, res)
 
subroutine mod_twofl_phys::twofl_get_pe_n_equi (w, x, ixIL, ixOL, res)
 
subroutine mod_twofl_phys::twofl_get_temperature_from_etot_n (w, x, ixIL, ixOL, res)
 Calculate temperature=p/rho when in e_ the total energy is stored this does not check the values of twofl_energy and twofl_internal_e, twofl_energy = .true. and twofl_internal_e = .false. also check small_values is avoided. More...
 
subroutine mod_twofl_phys::twofl_get_temperature_from_etot_n_with_equi (w, x, ixIL, ixOL, res)
 
subroutine mod_twofl_phys::twofl_get_temperature_from_eint_c (w, x, ixIL, ixOL, res)
 separate routines so that it is faster Calculate temperature=p/rho when in e_ the internal energy is stored More...
 
subroutine mod_twofl_phys::twofl_get_temperature_from_eint_c_with_equi (w, x, ixIL, ixOL, res)
 
subroutine mod_twofl_phys::twofl_get_temperature_c_equi (w, x, ixIL, ixOL, res)
 
subroutine mod_twofl_phys::twofl_get_rho_c_equi (w, x, ixIL, ixOL, res)
 
subroutine mod_twofl_phys::twofl_get_pe_c_equi (w, x, ixIL, ixOL, res)
 
subroutine mod_twofl_phys::twofl_get_temperature_from_etot_c (w, x, ixIL, ixOL, res)
 Calculate temperature=p/rho when in e_ the total energy is stored this does not check the values of twofl_energy and twofl_internal_e, twofl_energy = .true. and twofl_internal_e = .false. also check small_values is avoided. More...
 
subroutine mod_twofl_phys::twofl_get_temperature_from_eki_c (w, x, ixIL, ixOL, res)
 
subroutine mod_twofl_phys::twofl_get_temperature_from_etot_c_with_equi (w, x, ixIL, ixOL, res)
 
subroutine mod_twofl_phys::twofl_get_temperature_from_eki_c_with_equi (w, x, ixIL, ixOL, res)
 
subroutine mod_twofl_phys::twofl_get_csound2_adiab_n (w, x, ixIL, ixOL, csound2)
 
subroutine mod_twofl_phys::twofl_get_csound2_n_from_conserved (w, x, ixIL, ixOL, csound2)
 
subroutine mod_twofl_phys::twofl_get_csound2_n_from_primitive (w, x, ixIL, ixOL, csound2)
 
subroutine mod_twofl_phys::twofl_get_csound2_adiab_c (w, x, ixIL, ixOL, csound2)
 
subroutine, public mod_twofl_phys::twofl_get_csound2_c_from_conserved (w, x, ixIL, ixOL, csound2)
 
subroutine mod_twofl_phys::twofl_get_flux (wC, w, x, ixIL, ixOL, idim, f)
 Calculate fluxes within ixO^L. More...
 
subroutine mod_twofl_phys::twofl_add_source (qdt, dtfactor, ixIL, ixOL, wCT, wCTprim, w, x, qsourcesplit, active)
 w[iws]=w[iws]+qdt*S[iws,wCT] where S is the source based on wCT within ixO More...
 
subroutine mod_twofl_phys::add_pe_n0_divv (qdt, ixIL, ixOL, wCT, w, x)
 
subroutine mod_twofl_phys::add_pe_c0_divv (qdt, ixIL, ixOL, wCT, w, x)
 
subroutine mod_twofl_phys::add_geom_pdivv (qdt, ixIL, ixOL, v, p, w, x, ind)
 
subroutine mod_twofl_phys::get_lorentz (ixIL, ixOL, w, JxB)
 Compute the Lorentz force (JxB) More...
 
subroutine mod_twofl_phys::add_source_lorentz_work (qdt, ixIL, ixOL, w, wCT, x)
 
subroutine mod_twofl_phys::twofl_get_v_n (w, x, ixIL, ixOL, v)
 Calculate v_n vector. More...
 
subroutine, public mod_twofl_phys::get_rhon_tot (w, x, ixIL, ixOL, rhon)
 
subroutine, public mod_twofl_phys::twofl_get_pthermal_n (w, x, ixIL, ixOL, pth)
 
subroutine mod_twofl_phys::twofl_get_pthermal_n_primitive (w, x, ixIL, ixOL, pth)
 
subroutine, public mod_twofl_phys::twofl_get_v_n_idim (w, x, ixIL, ixOL, idim, v)
 Calculate v component. More...
 
subroutine mod_twofl_phys::internal_energy_add_source_n (qdt, ixIL, ixOL, wCT, w, x)
 
subroutine mod_twofl_phys::twofl_get_v_c (w, x, ixIL, ixOL, v)
 Calculate v_c vector. More...
 
subroutine, public mod_twofl_phys::get_rhoc_tot (w, x, ixIL, ixOL, rhoc)
 
subroutine, public mod_twofl_phys::twofl_get_pthermal_c (w, x, ixIL, ixOL, pth)
 
subroutine mod_twofl_phys::twofl_get_pthermal_c_primitive (w, x, ixIL, ixOL, pth)
 
subroutine, public mod_twofl_phys::twofl_get_v_c_idim (w, x, ixIL, ixOL, idim, v)
 Calculate v_c component. More...
 
subroutine mod_twofl_phys::internal_energy_add_source_c (qdt, ixIL, ixOL, wCT, w, x, ie)
 
subroutine mod_twofl_phys::twofl_handle_small_ei_c (w, x, ixIL, ixOL, ie, subname)
 handle small or negative internal energy More...
 
subroutine mod_twofl_phys::twofl_handle_small_ei_n (w, x, ixIL, ixOL, ie, subname)
 handle small or negative internal energy More...
 
subroutine mod_twofl_phys::add_source_b0split (qdt, ixIL, ixOL, wCT, w, x)
 Source terms after split off time-independent magnetic field. More...
 
subroutine mod_twofl_phys::add_source_res1 (qdt, ixIL, ixOL, wCT, w, x)
 Add resistive source to w within ixO Uses 3 point stencil (1 neighbour) in each direction, non-conservative. If the fourthorder precompiler flag is set, uses fourth order central difference for the laplacian. Then the stencil is 5 (2 neighbours). More...
 
subroutine mod_twofl_phys::add_source_res2 (qdt, ixIL, ixOL, wCT, w, x)
 Add resistive source to w within ixO Uses 5 point stencil (2 neighbours) in each direction, conservative. More...
 
subroutine mod_twofl_phys::add_source_hyperres (qdt, ixIL, ixOL, wCT, w, x)
 Add Hyper-resistive source to w within ixO Uses 9 point stencil (4 neighbours) in each direction. More...
 
subroutine mod_twofl_phys::add_source_glm (qdt, ixIL, ixOL, wCT, w, x)
 
subroutine mod_twofl_phys::add_source_powel (qdt, ixIL, ixOL, wCT, w, x)
 Add divB related sources to w within ixO corresponding to Powel. More...
 
subroutine mod_twofl_phys::add_source_janhunen (qdt, ixIL, ixOL, wCT, w, x)
 
subroutine mod_twofl_phys::add_source_linde (qdt, ixIL, ixOL, wCT, w, x)
 
subroutine, public mod_twofl_phys::get_normalized_divb (w, ixIL, ixOL, divb)
 get dimensionless div B = |divB| * volume / area / |B| More...
 
subroutine, public mod_twofl_phys::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. More...
 
subroutine mod_twofl_phys::gravity_add_source (qdt, ixIL, ixOL, wCT, w, x, energy, qsourcesplit, active)
 w[iw]=w[iw]+qdt*S[wCT,qtC,x] where S is the source based on wCT within ixO More...
 
subroutine mod_twofl_phys::gravity_get_dt (w, ixIL, ixOL, dtnew, dxD, x)
 
subroutine mod_twofl_phys::twofl_get_dt (w, ixIL, ixOL, dtnew, dxD, x)
 If resistivity is not zero, check diffusion time limit for dt. More...
 
pure logical function mod_twofl_phys::has_collisions ()
 
subroutine mod_twofl_phys::twofl_add_source_geom (qdt, dtfactor, ixIL, ixOL, wCT, w, x)
 
subroutine twofl_get_p_c_total (w, x, ixIL, ixOL, p)
 
subroutine mod_twofl_phys::twofl_get_temp_c_pert_from_etot (w, x, ixIL, ixOL, res)
 
double precision function, dimension(ixo^s) mod_twofl_phys::twofl_mag_en_all (w, ixIL, ixOL)
 Compute 2 times total magnetic energy. More...
 
double precision function, dimension(ixo^s) mod_twofl_phys::twofl_mag_i_all (w, ixIL, ixOL, idir)
 Compute full magnetic field by direction. More...
 
double precision function, dimension(ixo^s) mod_twofl_phys::twofl_mag_en (w, ixIL, ixOL)
 Compute evolving magnetic energy. More...
 
double precision function, dimension(ixo^s) mod_twofl_phys::twofl_kin_en_n (w, ixIL, ixOL)
 compute kinetic energy of neutrals More...
 
subroutine mod_twofl_phys::twofl_get_temp_n_pert_from_etot (w, x, ixIL, ixOL, res)
 
double precision function, dimension(ixo^s) mod_twofl_phys::twofl_kin_en_c (w, ixIL, ixOL)
 compute kinetic energy of charges w are conserved variables More...
 
subroutine mod_twofl_phys::twofl_getv_hall (w, x, ixIL, ixOL, vHall)
 
subroutine mod_twofl_phys::twofl_modify_wlr (ixIL, ixOL, qt, wLC, wRC, wLp, wRp, s, idir)
 
subroutine mod_twofl_phys::twofl_boundary_adjust (igrid, psb)
 
subroutine mod_twofl_phys::fixdivb_boundary (ixGL, ixOL, w, x, iB)
 
subroutine, public mod_twofl_phys::twofl_clean_divb_multigrid (qdt, qt, active)
 
subroutine mod_twofl_phys::twofl_update_faces (ixIL, ixOL, qt, qdt, wprim, fC, fE, sCT, s, vcts)
 
subroutine mod_twofl_phys::update_faces_average (ixIL, ixOL, qt, qdt, fC, fE, sCT, s)
 get electric field though averaging neighors to update faces in CT More...
 
subroutine mod_twofl_phys::update_faces_contact (ixIL, ixOL, qt, qdt, wp, fC, fE, sCT, s, vcts)
 update faces using UCT contact mode by Gardiner and Stone 2005 JCP 205, 509 More...
 
subroutine mod_twofl_phys::update_faces_hll (ixIL, ixOL, qt, qdt, fE, sCT, s, vcts)
 update faces More...
 
subroutine mod_twofl_phys::get_resistive_electric_field (ixIL, ixOL, sCT, s, jce)
 calculate eta J at cell edges More...
 
subroutine, public mod_twofl_phys::twofl_face_to_center (ixOL, s)
 calculate cell-center values from face-center values More...
 
subroutine, public mod_twofl_phys::b_from_vector_potential (ixIsL, ixIL, ixOL, ws, x)
 calculate magnetic field from vector potential More...
 
subroutine mod_twofl_phys::hyperdiffusivity_get_dt (w, ixIL, ixOL, dtnew, dxD, x)
 
subroutine mod_twofl_phys::add_source_hyperdiffusive (qdt, ixIL, ixOL, w, wCT, x)
 
subroutine add_density_hyper_source (index_rho)
 
double precision function, dimension(ixo^s, 1:nwc) mod_twofl_phys::dump_hyperdiffusivity_coef_y (ixIL, ixOL, w, x, nwc)
 
double precision function, dimension(ixo^s, 1:nwc) mod_twofl_phys::dump_hyperdiffusivity_coef_z (ixIL, ixOL, w, x, nwc)
 
double precision function, dimension(ixop^s, 1:nw) mod_twofl_phys::dump_hyperdiffusivity_coef_dim (ixIL, ixOPL, w, x, ii)
 
double precision function, dimension(ixo^s, 1:nwc) mod_twofl_phys::dump_coll_terms (ixIL, ixOL, w, x, nwc)
 
subroutine, public mod_twofl_phys::get_gamma_ion_rec (ixIL, ixOL, w, x, gamma_rec, gamma_ion)
 
subroutine, public mod_twofl_phys::get_alpha_coll (ixIL, ixOL, w, x, alpha)
 
subroutine mod_twofl_phys::get_alpha_coll_plasma (ixIL, ixOL, w, x, alpha)
 
subroutine mod_twofl_phys::calc_mult_factor1 (ixIL, ixOL, step_dt, JJ, res)
 
subroutine mod_twofl_phys::twofl_implicit_coll_terms_update (dtfactor, qdt, qtC, psb, psa)
 Implicit solve of psb=psa+dtfactor*dt*F_im(psb) More...
 
subroutine mod_twofl_phys::twofl_evaluate_implicit (qtC, psa)
 inplace update of psa==>F_im(psa) More...
 
subroutine mod_twofl_phys::coll_terms (ixIL, ixOL, w, x)
 
subroutine mod_twofl_phys::twofl_explicit_coll_terms_update (qdt, ixIL, ixOL, w, wCT, x)
 
subroutine mod_twofl_phys::rfactor_c (w, x, ixIL, ixOL, Rfactor)
 

Variables

integer, parameter, public mod_twofl_phys::eq_energy_tot =2
 
integer, parameter, public mod_twofl_phys::eq_energy_int =1
 
integer, parameter, public mod_twofl_phys::eq_energy_none =0
 
integer, parameter, public mod_twofl_phys::eq_energy_ki =3
 
integer, public, protected mod_twofl_phys::twofl_eq_energy = EQ_ENERGY_TOT
 
logical, public, protected mod_twofl_phys::twofl_hyperdiffusivity = .false.
 Whether hyperdiffusivity is used. More...
 
logical, public, protected mod_twofl_phys::twofl_dump_hyperdiffusivity_coef = .false.
 
double precision, dimension(:), allocatable, public, protected mod_twofl_phys::c_shk
 
double precision, dimension(:), allocatable, public, protected mod_twofl_phys::c_hyp
 
logical, public, protected mod_twofl_phys::twofl_thermal_conduction_c = .false.
 Whether thermal conduction is used. More...
 
logical, public, protected mod_twofl_phys::twofl_radiative_cooling_c = .false.
 Whether radiative cooling is added. More...
 
type(rc_fluid), allocatable, public mod_twofl_phys::rc_fl_c
 
logical, public, protected mod_twofl_phys::twofl_viscosity = .false.
 Whether viscosity is added. More...
 
logical, public, protected mod_twofl_phys::twofl_gravity = .false.
 Whether gravity is added: common flag for charges and neutrals. More...
 
logical, public, protected mod_twofl_phys::twofl_dump_full_vars = .false.
 whether dump full variables (when splitting is used) in a separate dat file More...
 
logical, public, protected mod_twofl_phys::twofl_hall = .false.
 Whether Hall-MHD is used. More...
 
type(tc_fluid), allocatable, public mod_twofl_phys::tc_fl_c
 
type(te_fluid), allocatable, public mod_twofl_phys::te_fl_c
 
logical, public, protected mod_twofl_phys::twofl_thermal_conduction_n = .false.
 
logical, public, protected mod_twofl_phys::twofl_radiative_cooling_n = .false.
 
logical, public, protected mod_twofl_phys::twofl_trac = .false.
 Whether TRAC method is used. More...
 
logical, public, protected mod_twofl_phys::twofl_glm = .false.
 Whether GLM-MHD is used. More...
 
integer, public, protected mod_twofl_phys::twofl_trac_type =1
 Which TRAC method is used
More...
 
double precision, public, protected mod_twofl_phys::twofl_trac_mask = 0.d0
 Height of the mask used in the TRAC method. More...
 
logical, public, protected mod_twofl_phys::source_split_divb = .false.
 Whether divB cleaning sources are added splitting from fluid solver. More...
 
double precision, public mod_twofl_phys::twofl_glm_alpha = 0.5d0
 GLM-MHD parameter: ratio of the diffusive and advective time scales for div b taking values within [0, 1]. More...
 
logical, public, protected mod_twofl_phys::twofl_4th_order = .false.
 MHD fourth order. More...
 
integer, public mod_twofl_phys::rho_c_
 Index of the density (in the w array) More...
 
integer, dimension(:), allocatable, public mod_twofl_phys::mom_c
 Indices of the momentum density. More...
 
integer, public mod_twofl_phys::e_c_ =-1
 Index of the energy density (-1 if not present) More...
 
integer, public mod_twofl_phys::tcoff_c_
 Index of the cutoff temperature for the TRAC method. More...
 
integer, public mod_twofl_phys::tweight_c_
 
integer, public, protected mod_twofl_phys::psi_
 Indices of the GLM psi. More...
 
logical, public mod_twofl_phys::has_equi_rho_c0 = .false.
 equi vars flags More...
 
logical, public mod_twofl_phys::has_equi_pe_c0 = .false.
 
integer, public mod_twofl_phys::equi_rho_c0_ = -1
 equi vars indices in the stateequi_vars array More...
 
integer, public mod_twofl_phys::equi_pe_c0_ = -1
 
logical, public mod_twofl_phys::twofl_equi_thermal_c = .false.
 
logical, public mod_twofl_phys::twofl_equi_thermal = .false.
 
integer, public mod_twofl_phys::rho_n_
 
integer, dimension(:), allocatable, public mod_twofl_phys::mom_n
 
integer, public mod_twofl_phys::e_n_
 
integer, public mod_twofl_phys::tcoff_n_
 
integer, public mod_twofl_phys::tweight_n_
 
logical, public mod_twofl_phys::has_equi_rho_n0 = .false.
 
logical, public mod_twofl_phys::has_equi_pe_n0 = .false.
 
integer, public mod_twofl_phys::equi_rho_n0_ = -1
 
integer, public mod_twofl_phys::equi_pe_n0_ = -1
 
double precision, public mod_twofl_phys::twofl_alpha_coll = 0d0
 collisional alpha More...
 
logical, public mod_twofl_phys::twofl_alpha_coll_constant = .true.
 
logical, public mod_twofl_phys::twofl_coll_inc_te = .true.
 whether include thermal exchange collisional terms More...
 
logical, public mod_twofl_phys::twofl_coll_inc_ionrec = .false.
 whether include ionization/recombination inelastic collisional terms More...
 
logical, public mod_twofl_phys::twofl_equi_thermal_n = .false.
 
double precision, public mod_twofl_phys::dtcollpar = -1d0
 
logical, public, protected mod_twofl_phys::twofl_dump_coll_terms = .false.
 whether dump collisional terms in a separte dat file More...
 
double precision, public, protected mod_twofl_phys::he_abundance = 0d0
 
double precision, public, protected mod_twofl_phys::rc = 2d0
 
double precision, public, protected mod_twofl_phys::rn = 1d0
 
double precision, public mod_twofl_phys::twofl_gamma = 5.d0/3.0d0
 The adiabatic index. More...
 
double precision, public mod_twofl_phys::twofl_adiab = 1.0d0
 The adiabatic constant. More...
 
double precision, public mod_twofl_phys::twofl_eta = 0.0d0
 The MHD resistivity. More...
 
double precision, public mod_twofl_phys::twofl_eta_hyper = 0.0d0
 The MHD hyper-resistivity. More...
 
double precision, public mod_twofl_phys::twofl_etah = 0.0d0
 The MHD Hall coefficient. More...
 
character(len=std_len), public, protected mod_twofl_phys::typedivbfix = 'linde'
 Method type to clean divergence of B. More...
 
character(len=std_len), public, protected mod_twofl_phys::type_ct = 'uct_contact'
 Method type of constrained transport. More...
 
logical, public, protected mod_twofl_phys::twofl_divb_4thorder = .false.
 Whether divB is computed with a fourth order approximation. More...
 
logical, public mod_twofl_phys::clean_initial_divb = .false.
 clean initial divB More...
 
logical, public mod_twofl_phys::divbwave = .true.
 Add divB wave in Roe solver. More...
 
logical, dimension(2 *^nd), public, protected mod_twofl_phys::boundary_divbfix =.true.
 To control divB=0 fix for boundary. More...
 
integer, dimension(2 *^nd), public, protected mod_twofl_phys::boundary_divbfix_skip =0
 To skip * layer of ghost cells during divB=0 fix for boundary. More...
 
logical, public, protected mod_twofl_phys::b0field_forcefree =.true.
 B0 field is force-free. More...
 
procedure(implicit_mult_factor_subroutine), pointer mod_twofl_phys::calc_mult_factor => null()
 
procedure(mask_subroutine), pointer, public mod_twofl_phys::usr_mask_alpha => null()
 
procedure(mask_subroutine2), pointer, public mod_twofl_phys::usr_mask_gamma_ion_rec => null()
 

Function/Subroutine Documentation

◆ add_density_hyper_source()

subroutine add_source_hyperdiffusive::add_density_hyper_source ( integer, intent(in)  index_rho)

Definition at line 6308 of file mod_twofl_phys.t.

Here is the call graph for this function:

◆ twofl_get_csound2_primitive()

subroutine twofl_get_csound_prim::twofl_get_csound2_primitive ( 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 
)

Definition at line 2632 of file mod_twofl_phys.t.

Here is the call graph for this function:

◆ twofl_get_p_c_total()

subroutine twofl_add_source_geom::twofl_get_p_c_total ( 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)  p 
)

Definition at line 4758 of file mod_twofl_phys.t.

Here is the call graph for this function: