import species
import matplotlib.pyplot as plt
############### SORT LEVELS ##################
constants = util.Constants()
nMax = 30
speciesList = ['O', 'O+', 'F+', 'Ne++']
allSpecies = []
for speciesStr in speciesList:
speciesObj = species.Species(speciesStr)
allSpecies.append([speciesObj, spectra.readNISTSpectra(speciesObj)])
O = species.Species('O')
use = O
NIST = spectra.readNISTSpectra(use)
theory = spectra.calculateExpectedStates(use, nMax)
calcEnergy = spectra.CalcEnergy(nMax, use, NIST, theory, allSpecies)
completeLevels, calculatedLevels = calcEnergy.populateTheory()
calculatedLevels = spectra.sortSpectra(calculatedLevels, 'calculated')
completeLevels = spectra.sortSpectra(completeLevels, 'complete')
NISTSorted = spectra.sortSpectra(NIST, 'NIST')
import spectra
import thermo
import util
import species
import matplotlib.pyplot as plt
constants = util.Constants()
nMax = 10
O = species.Species('O')
#He = species.Species('He')
use = O
NIST = spectra.readNISTSpectra(use)
theory = spectra.calculateExpectedStates(use, nMax)
completeLevels = spectra.calculateEnergies(nMax, use, NIST, theory)
tempRange = range(300, 50000, 50)
asdf = thermo.Thermo(O, NIST, tempRange, 250)
CpRange = asdf.getCpRange()
#print(CpRange)
#print(len(CpRange), len(tempRange))
Capitelli2005oi = [[
100, 200, 500, 700, 1000, 2000, 3000, 5000, 10000, 12000, 13000, 14000,
15000, 16000, 17000, 18000, 19000, 20000, 22000, 23000, 24000, 25000,
26000, 27000, 28000, 30000, 34000, 40000, 44000, 50000