|
| subroutine | mod_thermal_conduction::tc_init_params (phys_gamma) |
| |
| subroutine, public | mod_thermal_conduction::tc_get_mhd_params (fl, read_mhd_params) |
| | Init TC coefficients: MHD case.
|
| |
| subroutine, public | mod_thermal_conduction::tc_get_hd_params (fl, read_hd_params) |
| | Init TC coefficients: HD case.
|
| |
| double precision function, public | mod_thermal_conduction::get_tc_dt_mhd (w, ixil, ixol, dxd, x, fl) |
| | Get the explicit timestep for the TC (mhd implementation) Note: for multi-D MHD (1D MHD will use HD fall-back)
|
| |
| subroutine, public | mod_thermal_conduction::sts_set_source_tc_mhd (ixil, ixol, w, x, wres, fix_conserve_at_step, my_dt, igrid, nflux, fl) |
| | anisotropic thermal conduction with slope limited symmetric scheme Sharma 2007 Journal of Computational Physics 227, 123
|
| |
| subroutine | mod_thermal_conduction::set_source_tc_mhd (ixil, ixol, w, x, fl, qvec, rho, te, alpha) |
| |
| subroutine | mod_thermal_conduction::set_source_tc_mhd_geo (ixil, ixol, w, x, fl, qvec, rho, te, alpha) |
| |
| double precision function, dimension(ixi^s) | mod_thermal_conduction::slope_limiter (f, ixil, ixol, idims, pm, tc_slope_limiter) |
| |
| double precision function, public | mod_thermal_conduction::get_tc_dt_hd (w, ixil, ixol, dxd, x, fl) |
| | Get the explicit timestep for the TC (hd implementation) Note: also used in 1D MHD (or for neutrals in twofl)
|
| |
| subroutine, public | mod_thermal_conduction::sts_set_source_tc_hd (ixil, ixol, w, x, wres, fix_conserve_at_step, my_dt, igrid, nflux, fl) |
| |
| subroutine | mod_thermal_conduction::set_source_tc_hd (ixil, ixol, w, x, fl, qvec, rho, te) |
| |
| subroutine | mod_thermal_conduction::set_source_tc_hd_geo (ixil, ixol, w, x, fl, qvec, rho, te) |
| |
| subroutine, public | mod_thermal_conduction::tc_patch_negative_eint (w, x, ixil, ixol, te, ie, fl) |
| | Patch cells where e_int <= 0 by neighbor-averaging the temperature AND repairing the conserved internal energy w(:, ie) in place. Called after get_temperature_from_eint in the TC source routines when tc_patch_eint is .true. During STS RKL2 Chebyshev substeps the polynomial can overshoot e_int to negative values at low-density coronal cells. A negative T would produce NaN via sqrt(T) in the Spitzer conductivity. This routine replaces those cells' T with the average of valid (e_int > 0) neighbors and writes the matching e_int back to w(ix, ie) via the EoS inverse helper eint_nH_from_T, so the next STS substep does not re-read negative w and amplify it.
|
| |
