In this part of the parameter file the radiation transport module has to be selected. Depending on this selection a couple of additional parameters have to be specified. Table 12 gives a list of the parameters and the modules they apply to. The standard routines are now in the MSrad module for local models and the SHORTrad module for global ``StarinaBox'' models. The LHDrad module is not maintained very much anymore.

character radscheme
:
None
: Skip radiation transport entirely.
LHDrad
:
(old ``supergiant module'')
It uses long characteristics and is restricted to an equidistant grid and
open boundaries at all surfaces.
Note that the switch Drhd_r01=1
has to be set during compilation
(see Sect. 3.6).
MSrad
:
(``solar module'')
It uses long characteristics. The lateral boundaries have to be periodic.
Top and bottom can be closed or open.
Note that the switch Drhd_r02=1
has to be set during compilation
(see Sect. 3.6).
SHORTrad
:
(new ``supergiant module'')
It uses short characteristics and is restricted to an equidistant grid and
open boundaries at all surfaces.
Note that the switch Drhd_r03=1
has to be set during compilation
(see Sect. 3.6).
integer n_radminiter
:
1
.
In the implicit case it is set to a higher value (typically 5
).
integer n_raditer
:
n_raditer
iterations will (probably) needed.
The parameter can be set e.g. with
1
.
For stars with shorter radiative time scales values around 10
may be considered.
All three radiation transport modules understand this parameter.
integer n_radmaxiter
:
n_radmaxiter
will either be set to a value somewhat larger
than the recommended number of iterations (n_raditer
)
or to 0
which disables the check for too many iterations completely.
This can be safely allowed in many cases and has the advantage that there is
no need to save the initial model before calling the radiation transport module,
which saves some memory.
To disable the iteration of the radiation transport substep set
n_radminiter
=n_raditer
=n_radmaxiter
=1
.
All three radiation transport modules understand this parameter.
character radraybase
:
LHDrad
or SHORTrad
the orientation of the base axis system
can be selected e.g. with
radscheme=SHORTrad
are:
unity
: (default) During all time steps and radiative substeps
the direction of the rays stays the same.
random
: At each time step (and radiative substep) a new base axis system is
chosen at random
randomgroup
: At each new time step a new base axis system is
chosen at random. It is kept for all radiative substeps.
alternate
: Alternate between unity and random orientation matrix.
drift01
: Slow drift, algorithm 1.
drift02
: Slow drift, algorithm 2.
radraystar
) it is advisable to vary the directions of the rays
(by choosing radraybase
=random
or randomgroup
) to cover
the entire sphere at least over a longer time.
radscheme=MSrad
are:
lobatto
:
dblgaus
:
character radraystar
:
LHDrad
or SHORTrad
the list of ray directions (i.e. the number of rays and their coordinates)
relative to the base axis system can be specified with e.g.
x1
: (N=1) one single ray along x1 axis (not enough to specify fluxes in all directions)
x2
: (N=1) one single ray along x2 axis (not enough to specify fluxes in all directions)
x3
: (N=1) one single ray along x3 axis (not enough to specify fluxes in all directions)
oktaeder
: (N=3, default) octahedron
tetraeder
: (N=4) tetrahedron
cube
: (N=4)
ikosaeder
: (N=6) icosahedron
dodekaeder
: (N=10) dodecahedron
list01
, list01(3)
: Choose ray systems from a list (oktahedrons, tetrahedrons).
If character radraybase
is set to unity
the rays
will only be aligned to the axes or diagonals and thus avoid the
timeconsuming interpolation step of the shortcharacteristics method.
integer n_radtheta
:
MSrad
module the ray directions have to specified in a different way:
The number of ray sets in theta direction can be chosen with e.g.
integer n_radphi
:
MSrad
module
the number of ray sets in phi direction can be set e.g. with
integer n_radsubray
:
MSrad
module
the number of rays per cell (with the same direction)
can be specified e.g. with
integer n_radthickpoint
:
MSrad
module the lower part of the model can be computed in diffusion approximation.
The number of points in diffusion approximation can be set with e.g.
0
means that the diffusion approximation is not used
in any part of the model.
integer n_radtaurefine
:
LHDrad
and the MSrad
module the number of points on the rays
can be finer than the number of points in the basic numerical grid.
The refinement can be set e.g. with
integer n_radband
:
1
: Use the grey part of the table
2
: Use the other (possibly nongrey, frequencydependent) part of the table
3
: Use a continuum band plus an infrared band with CO opacity, calculated with
CO equilibrium density
4
: Use a continuum band plus an infrared band with CO opacity, calculated with
CO density resulting from timedependent chemistry
MSrad
module so far can handle nongrey tables.
integer n_radrsyslevel
:
integer n_radoutput
:
real c_radtcool
:
N_radband
0) of MSrad.
It can be activated with the parameter
real c_raddcool
:
N_radband
0) of MSrad.
It can be activated with the parameter
real c_radscool
:
N_radband
0) of MSrad.
It can be activated with the parameter
real c_radtinci
:
real c_raddinci
:
real c_radimplicitmu
:
LHDrad
and the SHORTrad
module support
implicit radiation transport.
It can be activated with the parameter
0.0
: Fully explicit radiation transport
(possible with all modules)
0.0
C 1.0
: Partly implicit radiation transport
0.5
: Radiation transport timecentered
1.0
: Fully implicit radiation transport
real c_raditereps
:
LHDrad
module only)
the requested convergence accuracy of the iteration can be set e.g. with
real c_raditerstep
:
LHDrad
module only)
the step size of the iteration can be restricted with e.g.
0.0
1.0
: Restricted step size
1.0
: No restriction, standard step size
1.0
: Extra large steps
real c_radtvisdtau
:
LHDrad
module
the limit in delta optical depth (rho*kappa*dx) below which the
``radiative temperature viscosity'' (=temperature smoothing) is to be applied
can be set with e.g.
c_radtvisdtau
0.0
.
real c_radtvis
:
LHDrad
module
the amount of the ``radiative temperature viscosity'' (=temperature smoothing)
can be specified e.g. with
c_radtvis
0.0
).
real c_tminlimit
: