def __init__(self, species, completeLevels, tempRange, cutoffMethod, numberDensity=0): self.species = species self.completeLevels = completeLevels self.tempRange = tempRange self.numberDensity = numberDensity if type(cutoffMethod) in [float, int]: self.cutoffMethod = 'fixedLowering' self.cutoffEnergy = cutoffMethod else: self.cutoffMethod = cutoffMethod self.constants = util.Constants() self.epsilon = [] self.J = [] for config in self.completeLevels: for level in config.term.levels: if level.energy is not None: self.epsilon.append(level.energy * self.constants.Cm_1ToJoules) self.J.append(level.J)
def __init__(self, species, configuration, level): self.species = species self.configuration = configuration self.level = level self.constants = util.Constants() self.energy = self.hydrogenicApprox()
def __init__(self, species): self.name = species self.constants = util.Constants() if '2' in self.name: self.readDiatomProperties() else: self.readAtomProperties()
def __init__(self, species, completeLevels, tempRange, cutoffEnergy): self.species = species self.completeLevels = completeLevels self.tempRange = tempRange self.cutoffEnergy = cutoffEnergy self.constants = util.Constants() self.epsilon = [] self.J = [] for config in self.completeLevels: for level in config.term.levels: self.epsilon.append(level.energy * self.constants.Cm_1ToJoules) self.J.append(level.J)
def __init__(self, nMax, species, NIST, theory, allSpecies): self.nMax = nMax self.species = species self.NIST = NIST self.theory = theory self.allSpecies = allSpecies self.constants = util.Constants() self.ions = [] for ionSpecies in self.allSpecies: speciesObj = ionSpecies[0] if speciesObj.noElectrons == self.species.noElectrons and speciesObj.name != self.species.name: self.ions.append(ionSpecies) for ionSpecies in self.allSpecies: speciesObj = ionSpecies[0] if speciesObj.atomicNumber == self.species.atomicNumber and speciesObj.noElectrons == self.species.noElectrons - 1: self.speciesIon = ionSpecies break
def __init__(self, nMax, species, NIST, theory, allSpecies): self.nMax = nMax self.species = species self.NIST = NIST self.theory = theory self.constants = util.Constants() self.speciesIon = None self.comprableIons = [] for ionSpecies in allSpecies: speciesObj = ionSpecies[0] if speciesObj.noElectrons == self.species.noElectrons and speciesObj.name != self.species.name: self.comprableIons.append(ionSpecies) elif speciesObj.atomicNumber == self.species.atomicNumber and speciesObj.noElectrons == self.species.noElectrons - 1: self.speciesIon = ionSpecies if len(self.comprableIons) == 0: print('missing ion data for ', self.species.name)
def __init__(self, nMax, species, NIST, theory, allSpecies): self.nMax = nMax self.species = species self.NIST = NIST self.theory = theory self.constants = util.Constants() self.speciesIon = None with open('hi.dat', 'w') as outFile: for config in self.NIST: outFile.write(config.ID + '\n') self.comprableIons = [] for ionSpecies in allSpecies: speciesObj = ionSpecies[0] if speciesObj.noElectrons == self.species.noElectrons and speciesObj.name != self.species.name: self.comprableIons.append(ionSpecies) elif speciesObj.atomicNumber == self.species.atomicNumber and speciesObj.noElectrons == self.species.noElectrons - 1: self.speciesIon = ionSpecies if len(self.comprableIons) == 0: print('missing ion data for ', self.species.name)
23.70, 22.74, 21.26, 21.04, 20.91, 20.83, 20.94, 21.80, 23.29, 24.35, 25.74, 28.15, 31.89, 37.19, 44.04, 52.19, 61.11, 70.00, 84.37, 88.54, 90.36, 89.94, 87.66, 84.00, 79.46, 69.37, 51.52, 35.92, 30.6, 26.30 ]] Gordon1999 = [[ 100, 200, 500, 700, 1000, 2000, 3000, 5000, 7000, 10000, 15000, 20000 ], [ 23.703, 22.734, 21.257, 21.040, 20.915, 20.827, 20.937, 21.799, 22.708, 23.150, 23.868, 24.722 ]] constants = util.Constants() #oi = Species('oi', 1000) #oii = Species('oii', 1000) #ni = Species('ni', 1000) #nii = Species('nii', 1000) #ci = Species('ci', 1000) #oi.calcCpRange(range(1000, 51000, 1000)) #oii.calcCpRange(range(1000, 51000, 1000)) #ni.calcCpRange(range(1000, 51000, 1000)) #nii.calcCpRange(range(1000, 51000, 1000)) #he = Species('he', 1000) #for level in he.states:
def __init__(self, species): self.name = species self.constants = util.Constants() self.readProperties()
def __init__(self, species): self.constants = util.Constants() self.species = species