def test_create_ChargeStates_dictionary(): elements = ['H', 'He', 'Fe'] AtomicData = read_atomic_data(elements) CSD = create_ChargeStates_dictionary(elements, 2.5879e8, AtomicData) for element in elements: print() print(element+':') print(CSD[element])
def test_create_ChargeStates_dictionary(): elements = ['H', 'He', 'Fe'] AtomicData = read_atomic_data(elements) CSD = create_ChargeStates_dictionary(elements, 2.5879e8, AtomicData) for element in elements: print() print(element + ':') print(CSD[element])
#------------------------------------------------------------------------------ # read AtomicData #------------------------------------------------------------------------------ AtomicData = sunnei.read_atomic_data(screen_output=True) # Define the plasma parameters elements = ['He', 'O', 'Mg'] #natom = AtomicData[element]['nstates'] - 1 te_ini = 4.0e+4 te0 = 2.0e+6 ne0 = 1.0e+7 te_list = [te0, te0] ne_list = [ne0, ne0 * 0.9] # Create a dictionary to save results ChargeStateDic = sunnei.create_ChargeStates_dictionary(elements) ChargeStateDic_pre = sunnei.create_ChargeStates_dictionary(elements) ChargeStateDic_next = sunnei.create_ChargeStates_dictionary(elements) #ChargeStateTime = [ChargeStateDic] #------------------------------------------------------------------------------ # Tests #------------------------------------------------------------------------------ # Initial conditions, set plasma temperature, density and dt # Start from any ionization states, e.g., Te = 4.0d4 K, time = 0 ind = sunnei.core.func_index_te(te_ini, AtomicData['temperatures']) print('************************************') print('TEST STARTING: elements =', elements)
#------------------------------------------------------------------------------ # read AtomicData #------------------------------------------------------------------------------ AtomicData = sunnei.read_atomic_data(screen_output=True) # Define the plasma parameters elements = ['He', 'O', 'Mg'] #natom = AtomicData[element]['nstates'] - 1 te_ini = 4.0e+4 te0 = 2.0e+6 ne0 = 1.0e+7 te_list = [te0, te0] ne_list = [ne0, ne0*0.9] # Create a dictionary to save results ChargeStateDic = sunnei.create_ChargeStates_dictionary(elements) ChargeStateDic_pre = sunnei.create_ChargeStates_dictionary(elements) ChargeStateDic_next = sunnei.create_ChargeStates_dictionary(elements) #ChargeStateTime = [ChargeStateDic] #------------------------------------------------------------------------------ # Tests #------------------------------------------------------------------------------ # Initial conditions, set plasma temperature, density and dt # Start from any ionization states, e.g., Te = 4.0d4 K, time = 0 ind = sunnei.core.func_index_te(te_ini, AtomicData['temperatures']) print('************************************') print('TEST STARTING: elements =', elements)