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