Python ChemicalNetwork.add_species示例

示例#1

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)

示例#2

文件： run_ion_network.py 项目： cassidymwagner/dengo-merge

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)

示例#3

文件： write_ion_by_ion_network.py 项目： drreynolds/dengo-merge

    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))

示例#4

文件： run_ion_network.py 项目： drreynolds/dengo-merge

    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)