NCELLS = 1
density = 1.0
#* 1.67e-24
temperature = np.logspace(2, 8, NCELLS)
temperature[:] = 5e6
X = 1e-8
ion_by_ion = ChemicalNetwork(write_intermediate = False,
stop_time = 3.1557e13)
HI = AtomicSpecies("H", 0.0)
HII = AtomicSpecies("H", 1.0)
HeI = AtomicSpecies("He", 0.0)
HeII = AtomicSpecies("He", 1.0)
de = species_registry['de']
ion_by_ion.add_species(de)
ion_by_ion.add_species(HI)
ion_by_ion.add_species(HII)
ion_by_ion.add_species(HeI)
ion_by_ion.add_species(HeII)
#ion_by_ion.add_species('O_1')
#ion_by_ion.add_species('O_2')
#ion_by_ion.add_species('O_3')
#ion_by_ion.add_species('O_4')
#ion_by_ion.add_species('O_5')
#ion_by_ion.add_species('O_6')
for atom in ["H", "He", "O"]: #"C", "N", "O", "Ne", "Si"]:
s, c, r = setup_ionization(atom, photo_background="HM12")
ion_by_ion.add_collection(s, c, r)
NCELLS = 1
density = 1.0
#* 1.67e-24
temperature = np.logspace(2, 8, NCELLS)
temperature[:] = 5e6
X = 1e-8
ion_by_ion = ChemicalNetwork(write_intermediate=False, stop_time=3.1557e13)
HI = AtomicSpecies("H", 0.0)
HII = AtomicSpecies("H", 1.0)
HeI = AtomicSpecies("He", 0.0)
HeII = AtomicSpecies("He", 1.0)
de = species_registry['de']
ion_by_ion.add_species(de)
ion_by_ion.add_species(HI)
ion_by_ion.add_species(HII)
ion_by_ion.add_species(HeI)
ion_by_ion.add_species(HeII)
ion_by_ion.add_species('O_1')
ion_by_ion.add_species('O_2')
ion_by_ion.add_species('O_3')
ion_by_ion.add_species('O_4')
ion_by_ion.add_species('O_5')
ion_by_ion.add_species('O_6')
for atom in ["H", "He", "O"]: #"C", "N", "O", "Ne", "Si"]:
s, c, r = setup_ionization(atom, photo_background="HM12")
ion_by_ion.add_collection(s, c, r)
reaction_registry, \
cooling_registry
from dengo.primordial_rates import setup_primordial
import dengo.primordial_cooling
from dengo.ion_by_ion import setup_ionization
from dengo.chemistry_constants import tiny, kboltz, mh
import numpy as np
NCELLS = 1
density = 1e-3 * 1.67e-24
temperature = np.logspace(4, 6.7, NCELLS)
temperature[:] = 5e6
X = 1e-3
ion_by_ion = ChemicalNetwork(write_intermediate=False, stop_time=3.1557e13)
ion_by_ion.add_species("de")
#for atom in ["O", "C", "Si", "Mg", "N", "S", "He", "Ne", "H"]:
for atom in ["H", "O", "He"]:
s, c, r = setup_ionization(atom)
ion_by_ion.add_collection(s, c, r)
#ion_by_ion.add_collection(s, [], r)
# ion_by_ion.add_cooling('compton')
#s, c, r = setup_primordial()
#ion_by_ion.add_collection(s, c, r)
# This defines the temperature range for the rate tables
ion_by_ion.init_temperature((1e0, 1e12))
ChemicalNetwork, \
reaction_registry, \
cooling_registry
from dengo.ion_by_ion import setup_ionization
from dengo.chemistry_constants import tiny, kboltz, mh
import numpy as np
NCELLS = 1
density = 1.0
#* 1.67e-24
temperature = np.logspace(2, 8, NCELLS)
temperature[:] = 5e6
X = 1e-8
ion_by_ion = ChemicalNetwork(write_intermediate=False, stop_time=3.1557e13)
ion_by_ion.add_species("de")
ion_by_ion.add_species('H_1')
ion_by_ion.add_species('H_2')
ion_by_ion.add_species('He_1')
ion_by_ion.add_species('He_2')
#ion_by_ion.add_species('O_1')
#ion_by_ion.add_species('O_2')
#ion_by_ion.add_species('O_3')
#ion_by_ion.add_species('O_4')
#ion_by_ion.add_species('O_5')
#ion_by_ion.add_species('O_6')
for atom in ["H", "He", "O"]: #"C", "N", "O", "Ne", "Si"]:
s, c, r = setup_ionization(atom, photo_background="hm12")
ion_by_ion.add_collection(s, c, r)