#species = ['A', 'B', 'C', 'D']
species = ['A2', 'B2', 'C2', 'D2']#, 'B2', 'C10A10']
enthalpy = [0.0, 0.0, 1.5e6/2, 0.0]
S_minf = [300.0, 0.1, 0.1, 0.1, 0.7]#, 0.3, 0.5]
S_pinf = [900.0, 0.0, 0.0, 0.0, 1.0]#, 0.0, 0.0]
c = Chemistry(species, enthalpy)
# reaction = {'equation':'A = B', 'A': 1e9, 'b': 1.0, 'Ta':14000.0, 'Q': 1.5e6}
# c.add_reaction(reaction)
# reaction = {'equation':'D = C', 'A': 1e9, 'b': 1.0, 'Ta':14000.0, 'Q': 1.5e6}
# c.add_reaction(reaction)
# reaction = {'equation':'C = D', 'A': 1e9, 'b': 1.0, 'Ta':14000.0, 'Q': -1.5e6}
# c.add_reaction(reaction)
# reaction = {'equation':'B = A', 'A': 1e9, 'b': 1.0, 'Ta':14000.0, 'Q': -1.5e6}
# c.add_reaction(reaction)
#reaction = {'equation':'A2 = C2', 'A': 1e9, 'b': 1.0, 'Ta':14000.0}
#c.add_reaction(reaction)
reaction = {'equation':'C2 = A2', 'A': 1e9, 'b': 1.0, 'Ta':14000.0}
# reaction = {'equation':'H2 + 0.5*O2 = H2O', 'A': 1e9, 'b': 1.0, 'Ta':14000.0, 'Q': -1.5e6}
c.add_reaction(reaction)
#reaction = {'equation':'2*C = A + B', 'A': 1e9, 'b': 1.0, 'Ta':14000.0, 'Q': -1.5e6}
#c.add_reaction(reaction)
flame = FlameBase(c, S_minf, S_pinf)
flame.dt = 1e-3
flame.solve()
