Metal dust radiation

[genchiaki]

Metal/dust chemistry/cooling is included into grackle by pull request #78.
This modification make enzo compatible with it.
Each new feature can be switched on/off with corresponding parameters:

MultiSpecies = 4
switch to solve chemistry for 3 extra primordial species (D-, HD+, and HeH+)
MetalChemistry
switch to solve metal chemistry
GrainGrowth
switch to solve the growth of dust grains
MultiMetals
This switch enables to consider multiple metal sources (e.g. supernovae with different metal yields).
Now this version includes

  0. metal/dust abundances of local ISM (Pollack et al. 1994)
  1-4. Pop III normal core-collapse supernovae (Nozawa et al. 2007)
    with progenitor masses 13, 20, 25 and 30 Msun
  5-8. Pop III faint supernovae (Marassi et al. 2014)
    with progenitor masses 13, 50 and 80 Msun
  9-10. Pop III pair-instability supernovae (Nozawa et al. 2007)
    with progenitor masses 170 and 200 Msun
  11. simple dust model (only include silicate and graphite; Yajima et al. 2019)

MetalAbundances
If multi_metals = 0, you can select a single metal/dust yield from the list above.
DustSpecies
You can consider multiple grain species

  1. enstatite + amorphous carbon
  2. + metallic silicon + metallic iron + forsterite + magnetite + silica + magnesia + troilite + alumina
  3. + water ice + volatile organics + refractory organics

DustTemperatureMulti
switch to consider different dust temperature for each species if dust_species > 0
DustSublimation
switch to include dust sublimation

UseDustDensityField
enables to use dust density fields
VariableGamma
enables to consider non-uniform adiabatic indices in hydro solvers

Also photochemistry of metal species is included
RadiativeTransferHDIDiss
photodissociation of HD molecules
RadiativeTransferMetalIon
photoionization of C and O atoms
RadiativeTransferMetalDiss
photodissociation of CO, OH and H2O molecules

[jwise77]

During a discussion in the Enzo/Enzo-E workshop, we decided to address this PR after #160 because they might conflict because they add new metal fields.