radiation_ray

Module

Description

$Id$

Radiation (solves transfer equation along rays).

The direction of the ray is given by the vector (lrad,mrad,nrad), and the parameters radx0,rady0,radz0 gives the maximum number of steps of the direction vector in the corresponding direction.

This module currently does not work for fully periodic domains. The z-direction has to be non-periodic

Note that Qrad is the heating rate (Q=I-S) and not the cooling rate used in Heinemann et al. 2006. Furthermore, Frad is the radiative flux multiplied by kappa*rho, not actually the radiative flux itself.

TODO: Calculate weights properly

** AUTOMATIC CPARAM.INC GENERATION ************************ Declare (for generation of cparam.inc) the number of f array variables and auxiliary variables added by this module

CPARAM logical, parameter :: lradiation = .true., lradiation_ray = .true.

MVAR CONTRIBUTION 0 MAUX CONTRIBUTION 2

Quick access

Types:

qbound

Variables:

a_sb, ampl_pert, c_gam, dff_new, flim, idiag_egas_max, idiag_egas_rms, idiag_erad_max, idiag_erad_rms, idiag_fmax, idiag_frms, inflow, initrad, k_b_radiation, kappa_es_radiation, opas, pertee, qbound

Routines:

bc_ee_inflow_x(), bc_ee_outflow_x(), calc_diagnostics_radiation(), calc_pencils_radiation(), de_dt(), dradiation_dt(), flux_limiter(), get_slices_radiation(), init_equil(), init_rad(), initialize_radiation(), pencil_criteria_radiation(), pencil_interdep_radiation(), radiative_cooling(), radiative_pressure(), radtransfer(), read_radiation_init_pars(), read_radiation_run_pars(), register_radiation(), rprint_radiation(), write_radiation_init_pars(), write_radiation_run_pars()

Needed modules

Types

  • type  radiation/qbound
    Type fields:

    • % qpoint [real] :: Qpt

    • % set [logical]

    • % val [real]

Variables

  • radiation/qbound [private]