CO5BOLD can now 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: 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) 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).
real c_dust10) to control each dust formation scheme in detail. A parameter can be given as in
0.0in each case.
real c_dust01must be set to
1.0in 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.
1.0E-04means 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
qucarray for the metal is found in the input model.