mod_rt_sc_slab.t File Reference

Go to the source code of this file.

Modules
module	mod_rt_sc_slab
	Public face of the short-characteristics (SC) radiative transfer and its coupling to the (M)HD.

Functions/Subroutines
subroutine, public	mod_rt_sc_slab::rt_sc_activate ()
	Switch the RT on. The physics module calls this from its initialisation when hd_radiative_transfer or mhd_radiative_transfer is set. It reads &rt_list, reserves the intensity slots of the wavefront and registers the output variables Qrad and Frad, with Iout as well for the replicated solver and the validation mode, which fill it. rt_sc_setup follows once the units and the domain are known.
subroutine, public	mod_rt_sc_slab::rt_sc_pin_refine (igrid, level, my_refine, my_coarsen)
	Hold every block overlapping an RT band on rt_sc_amr_level. Called from forcedrefine_grid after usr_refine_grid, so it overrides both the estimator and the user.
subroutine, public	mod_rt_sc_slab::rt_sc_selftest ()
	Self-test of the real sweep on an analytic slab. J must approach B where the slab is optically thick, stay finite in the tenuous top, where the gray FLD coefficient D = c lambda/(kappa rho) would diverge, and stay bounded by the source everywhere. The slab is horizontally infinite, so the test sets thermal side inflows, which are exact for it, a thermal bottom and an open top, and restores the run's boundaries afterwards.
subroutine, public	mod_rt_sc_slab::rt_sc_emergent_tau1 ()
	Horizontal mean of Fy on the tau = rt_sc_tau_emergent surface, the photosphere, into rt_sc_Ftau1. The flux controller needs this value, since Frad is not conserved above tau = 1 and the flux at the band top reads larger. Only the first band is searched, since it is the one that holds the photosphere, and tau comes from the band top as in the KBA heat blend.
subroutine, public	mod_rt_sc_slab::rt_sc_add_source (qdt, ixil, ixol, wct, w, x, thin_share, e_internal)
	Unsplit heating source of one block and stage, e += qdt Q, with Q read from i_qrad, which is zero outside the bands. The heating is limited so that one stage changes the gas energy by at most rt_sc_cfl of itself.
pure elemental double precision function, public	mod_rt_sc_slab::rt_sc_thin_weight (te_code, y_code, p_code)
	Share of the radiative losses taken by the optically thin treatment, with SC taking the rest, 1 - w_thin. The thin cooling and the user SC source both call this, so they share one definition. Inside a band w_thin is the product of a temperature gate and a pressure gate, so a cell goes to SC when it is cool or dense, as in the opacity switch. The pressure gate max(0, 1 - (p/p_TR)^2) of MURaM's Qthin is exactly zero in the photosphere. Outside every band, or without RT, w_thin = 1, because there is no SC there to take over. The step to one at the band edge marks where SC stops, and the gas there is coronal, with the temperature gate already close to one.