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])
Exemple #2
0
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])
Exemple #3
0
#------------------------------------------------------------------------------
# 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)