# print Levels of O III
print(oiii_elj_data['j_v'])
# print Energy Levels (cm-1) of O III
print(oiii_elj_data['ej'])
# get citations for Energy Levels (Ej) Reference L7288
citation = atomneb.get_elj_reference_citation(atom_elj_file, 'L7288')
# print citations for Energy Levels (Ej) Reference L7288
print(citation)
atom = 'o'
ion = 'iii'
reference = 'SSB14'
# read Collision Strengths (Omegaij) of O III from Reference SSB14
oiii_omij_data = atomneb.read_omij(atom_omij_file,
atom,
ion,
reference=reference)
# print Level 1 of Collision Strengths (Omegaij) of O III
print(oiii_omij_data['level1'])
# print Level 2 of Collision Strengths (Omegaij) of O III
print(oiii_omij_data['level2'])
# print Strength[1] of Collision Strengths (Omegaij) of O III
print(oiii_omij_data['strength'][0])
atom = 'o'
ion = 'iii'
# list all Collision Strengths (Omegaij) data for O III
list_oiii_omij_data = atomneb.search_omij(atom_omij_file, atom, ion)
# print all Collision Strengths (Omegaij) of O III
print(list_oiii_omij_data)
data_dir = os.path.join('atomic-data', 'chianti70')
data_rc_dir = os.path.join('atomic-data-rc')
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)
# read Energy Levels (Ej) of O III upto level number 6 oiii_elj_data = atomneb.read_elj(atom_elj_file, atom, ion, level_num=6) # print Levels of O III print(oiii_elj_data['j_v']) # print Energy Levels (cm-1) of O III print(oiii_elj_data['ej']) # get citations for Energy Levels (Ej) Reference o_iii_elj citation = atomneb.get_elj_reference_citation(atom_elj_file, 'o_iii_elj') # print citations for Energy Levels (Ej) Reference o_iii_elj print(citation) atom = 'o' ion = 'iii' # read Collision Strengths (Omegaij) of O III oiii_omij_data = atomneb.read_omij(atom_omij_file, atom, ion) # print Level 1 of Collision Strengths (Omegaij) of O III print(oiii_omij_data['level1']) # print Level 2 of Collision Strengths (Omegaij) of O III print(oiii_omij_data['level2']) # print Strength[1] of Collision Strengths (Omegaij) of O III print(oiii_omij_data['strength'][0]) atom = 'o' ion = 'iii' # list all Collision Strengths (Omegaij) data for O III list_oiii_omij_data = atomneb.search_omij(atom_omij_file, atom, ion) # print all Collision Strengths (Omegaij) of O III print(list_oiii_omij_data) atom = 'o'
