print(list_oiii_omij_reference)

atom = 'o'
ion = 'iii'
reference = 'SSB14'
# get citations for Collision Strengths (Omegaij) of O III with reference SSB14
citation = atomneb.get_omij_reference_citation(atom_omij_file, atom, ion,
                                               reference)
# print citations for Collision Strengths (Omegaij) of O III with reference SSB14
print(citation)

atom = 'o'
ion = 'iii'
reference = 'FFT04'
# read Transition Probabilities (Aij) of O III with reference FFT04
oiii_aij_data = atomneb.read_aij(atom_aij_file, atom, ion, reference)
# print Transition Probabilities (Aij) of O III with reference FFT04
print(oiii_aij_data['aij'][0])

atom = 'o'
ion = 'iii'
reference = 'FFT04'
# get citations for Transition Probabilities (Aij) of O III with reference FFT04
citation = atomneb.get_aij_reference_citation(atom_aij_file, atom, ion,
                                              reference)
# print citations for Transition Probabilities (Aij) of O III with reference FFT04
print(citation)

atom = 'o'
ion = 'iii'
# list all Transition Probabilities (Aij) data for O III
Exemplo n.º 2
0
atom_elj_file = os.path.join(base_dir, data_dir, 'AtomElj.fits')
atom_omij_file = os.path.join(base_dir, data_dir, 'AtomOmij.fits')
atom_aij_file = os.path.join(base_dir, data_dir, 'AtomAij.fits')
atom_rc_sh95_file = os.path.join(base_dir, data_rc_dir, 'rc_SH95.fits')

atom = 'h'
ion = 'ii'  # H I Rec
hi_rc_data = atomneb.read_aeff_sh95(atom_rc_sh95_file, atom, ion)

atom = 'o'
ion = 'iii'  # [O III]
o_iii_elj = atomneb.read_elj(atom_elj_file, atom, ion,
                             level_num=5)  # read Energy Levels (Ej)
o_iii_omij = atomneb.read_omij(atom_omij_file, atom,
                               ion)  # read Collision Strengths (Omegaij)
o_iii_aij = atomneb.read_aij(atom_aij_file, atom,
                             ion)  # read Transition Probabilities (Aij)

levels5007 = '3,4/'
temperature = np.float64(10000.0, )
density = np.float64(5000.0)
iobs5007 = np.float64(1200.0)
abb5007 = np.float64(0.0)

emis = pyequib.calc_emissivity(temperature=temperature,
                               density=density,
                               atomic_levels=levels5007,
                               elj_data=o_iii_elj,
                               omij_data=o_iii_omij,
                               aij_data=o_iii_aij)
print('Emissivity(O III 5007):', emis)
Exemplo n.º 3
0
# print all Collision Strengths (Omegaij) References for O III
print(list_oiii_omij_reference)

atom = 'o'
ion = 'iii'
reference = 'CHI52'
# get citations for Collision Strengths (Omegaij) of O III with reference CHI52
citation = atomneb.get_omij_reference_citation(atom_omij_file, atom, ion, reference)
# print citations for Collision Strengths (Omegaij) of O III with reference CHI52
print(citation)

atom = 'o'
ion = 'iii'
reference = 'CHI52'
# read Transition Probabilities (Aij) of O III with reference CHI52
oiii_aij_data = atomneb.read_aij(atom_aij_file, atom, ion)
# print Transition Probabilities (Aij) of O III with reference CHI52
print(oiii_aij_data['aij'][0])

atom = 'o'
ion = 'iii'
reference = 'CHI52'
# get citations for Transition Probabilities (Aij) of O III with reference CHI52
citation = atomneb.get_aij_reference_citation(atom_aij_file, atom, ion, reference)
# print citations for Transition Probabilities (Aij) of O III with reference CHI52
print(citation)

atom = 'o'
ion = 'iii'
# list all Transition Probabilities (Aij) data for O III
list_oiii_aij_data = atomneb.search_aij(atom_aij_file, atom, ion)