CO5BOLD can handle a number of additional density arrays. They can be used to describe e.g. the mass density of dust distribution moments or number densities of molecules. These species are properly advected with the gas density. There is also already a small number of dust/molecule formation models available. These models have to be improved in the future and the influence on the radiation field (opacities, radiation pressure on dust) has to be taken into account.
none, None: No handling of any dust/molecule density at all.
nosource: Skip source term step for dust/molecules entirely, but do the transport.
dust_simple_01: Simple and unrealistic 'dust' formation model (only for testing of
the numerics).
co_component01_01: Simple CO formation (from Matthias Steffen)
with one component only but realistic time scales.
dust_k3mon_01: Simple C-rich dust formation (routines from Susanne Höfner)
with one component only but realistic time scales.
dust_k3mon_02: Simple C-rich dust formation (routines from Susanne Höfner)
with two components for dust density and free carbon density.
dust_k3mon_03: Simple forsterite dust formation (based on routines from Susanne Höfner)
with two components for dust density and free "forsterite monomer" density.
dust_bins_01: Multi-size-bin forsterite dust formation
(based on Rossow's equations, see Rossow, 1978)
with one bin for the monomers and several bins for the different grain sizes.
dust_moment04_c2: C-rich dust chemistry, 4 moments (routines from Susanne Höfner).
chemreacnet: chemical reaction networks (routines from Sven Wedemeyer-Böhm and Inga Kamp).
hiontd: time-dependent hydrogen ionization (routines from Jorrit Leenaarts and Sven Wedemeyer-Böhm).
PP)
or less (e.g. PPmimo) aggressive.
It can be chosen e.g. with
reconstruction=vanLeer, dustreconstruction=Superbee.
reconstruction=FRweno, dustreconstruction=FRmono.
reconstruction=FRweno, dustreconstruction=PP.
reconstruction=HBweno, dustreconstruction=PP.
reconstruction=HBweno, dustreconstruction=PPmimo.
real c_dust01 to c_dust10)
to control each dust formation scheme in detail. A parameter can be given as in
0.0 in each case.
real c_dust01 must be set to 1.0 in order to
activate advection of particle densities for the CHEM (chemreacnet) and the
HION (hiontd) module. The default value is 0.0, i.e. advection is switched off.
![$\epsilon_\mathrm{M} = [M] / [H]$](img215.png){kind=small}
). A value of 1.0E-04 means there
are 
hydrogen atoms for every metal atom. The metal number density
for each grid cell is then derived via
(assuming a pure
hydrogen gas).
This parameter is ignored when a quc array for the metal is found in the
input model.