#
#  This is the NASA thermodynamic database, which is available for download
#  from http://www.galcit.caltech.edu/EDL/public/thermo.html. The original 
#  sources are
#
#  S. Gordon and B.J. McBride, "Computer Program for Calculation of Complex 
#  Chemical Equilibrium Composition, Rocket Performance, Incident and 
#  Reflected Shocks and Chapman-Jouguet Detonations", NASA Report SP-273, 1971.
#
#  B.J. McBride, S. Gordon, and M.A. Reno, "Coefficients for Calculating 
#  Thermodynamic and Transport Properties of Individual Species", NASA 
#  Report TM-4513, October 1993.
#
#-----------------------------------------------------------------------------

ideal_gas(name = "air_below6000K",
      elements = " N O E ",
      species = """ N2 O2 N O NO NO+ Electron N+ O+ N2+ O2+ """,
      reactions = "all",
      initial_state = state(temperature = 298.15,
                            pressure = OneAtm,
                            mole_fractions = 'N2:0.78847, O2:0.21153')    )

species(name = "Electron",
    atoms = " E:1 ",
    thermo = (
       NASA( [  200.00,  1000.00], [  2.500000000E+00,   0.000000000E+00, 
                0.000000000E+00,   0.000000000E+00,   0.000000000E+00,
               -7.453750000E+02,  -1.172469020E+01] ),
       NASA( [ 1000.00,  6000.00], [  2.500000000E+00,   0.000000000E+00, 
                0.000000000E+00,   0.000000000E+00,   0.000000000E+00,
               -7.453750000E+02,  -1.172469020E+01] )
             ),
    note = "gas L10/92"
       )

species(name = "N",
    atoms = " N:1 ",
    thermo = (
       NASA( [  200.00,  1000.00], [  2.500000000E+00,   0.000000000E+00, 
                0.000000000E+00,   0.000000000E+00,   0.000000000E+00,
                5.610463780E+04,   4.193909320E+00] ),
       NASA( [ 1000.00,  6000.00], [  2.415942930E+00,   1.748906000E-04, 
               -1.190236670E-07,   3.022623870E-11,  -2.036097900E-15,
                5.613377480E+04,   4.649609860E+00] )
             ),
    note = "L 6/88"
       )

species(name = "N+",
    atoms = " N:1  E:-1 ",
    thermo = (
       NASA( [  200.00,  1000.00], [  2.802694450E+00,  -1.447589110E-03, 
                2.771183800E-06,  -2.401873520E-09,   7.808399310E-13,
                2.255752440E+05,   3.578778350E+00] ),
       NASA( [ 1000.00,  6000.00], [  2.511129670E+00,   3.464417510E-06, 
               -1.594269380E-08,   7.248656630E-12,  -6.445014260E-16,
                2.256243400E+05,   4.927676610E+00] )
             ),
    note = "L 7/88"
       )

species(name = "NO",
    atoms = " N:1  O:1 ",
    thermo = (
       NASA( [  200.00,  1000.00], [  4.218598960E+00,  -4.639881240E-03, 
                1.104430490E-05,  -9.340555070E-09,   2.805548740E-12,
                9.845099640E+03,   2.280610010E+00] ),
       NASA( [ 1000.00,  6000.00], [  3.260712340E+00,   1.191011350E-03, 
               -4.291226460E-07,   6.944814630E-11,  -4.032956810E-15,
                9.921431320E+03,   6.369005180E+00] )
             ),
    note = "TPIS89"
       )

species(name = "NO+",
    atoms = " N:1  O:1  E:-1 ",
    thermo = (
       NASA( [  200.00,  1000.00], [  3.693012310E+00,  -1.342291580E-03, 
                2.673433950E-06,  -1.026093080E-09,  -6.956104920E-14,
                1.181030550E+05,   3.091266980E+00] ),
       NASA( [ 1000.00,  6000.00], [  2.945877020E+00,   1.403252600E-03, 
               -4.955031960E-07,   7.959489730E-11,  -4.720766680E-15,
                1.182443400E+05,   6.706446410E+00] )
             ),
    note = "TPIS89"
       )

species(name = "N2",
    atoms = " N:2 ",
    thermo = (
       NASA( [  200.00,  1000.00], [  3.531005280E+00,  -1.236609870E-04, 
               -5.029994370E-07,   2.435306120E-09,  -1.408812350E-12,
               -1.046976280E+03,   2.967474680E+00] ),
       NASA( [ 1000.00,  6000.00], [  2.952576260E+00,   1.396900570E-03, 
               -4.926316910E-07,   7.860103670E-11,  -4.607553210E-15,
               -9.239486450E+02,   5.871892520E+00] )
             ),
    note = "TPIS78"
       )

species(name = "N2+",
    atoms = " N:2  E:-1 ",
    thermo = (
       NASA( [  200.00,  1000.00], [  3.775407110E+00,  -2.064591570E-03, 
                4.757523010E-06,  -3.156642280E-09,   6.705099730E-13,
                1.804811150E+05,   2.693221860E+00] ),
       NASA( [ 1000.00,  6000.00], [  3.586613630E+00,   2.530719490E-04, 
                1.847782140E-07,  -4.552572230E-11,   3.268180290E-15,
                1.803909940E+05,   3.095841500E+00] )
             ),
    note = "TPIS89"
       )

species(name = "O",
    atoms = " O:1 ",
    thermo = (
       NASA( [  200.00,  1000.00], [  3.168267100E+00,  -3.279318840E-03, 
                6.643063960E-06,  -6.128066240E-09,   2.112659710E-12,
                2.912225920E+04,   2.051933460E+00] ),
       NASA( [ 1000.00,  6000.00], [  2.543636970E+00,  -2.731624860E-05, 
               -4.190295200E-09,   4.954818450E-12,  -4.795536940E-16,
                2.922601200E+04,   4.922294570E+00] )
             ),
    note = "L 1/90"
       )

species(name = "O+",
    atoms = " O:1  E:-1 ",
    thermo = (
       NASA( [  200.00,  1000.00], [  2.500000000E+00,   0.000000000E+00, 
                0.000000000E+00,   0.000000000E+00,   0.000000000E+00,
                1.879352910E+05,   4.393376890E+00] ),
       NASA( [ 1000.00,  6000.00], [  2.487733170E+00,   2.176600160E-05, 
               -1.089558060E-08,   1.259092120E-12,   1.373167200E-16,
                1.879399650E+05,   4.461340910E+00] )
             ),
    note = "L 1/90"
       )

species(name = "O2",
    atoms = " O:2 ",
    thermo = (
       NASA( [  200.00,  1000.00], [  3.782456360E+00,  -2.996734150E-03, 
                9.847302000E-06,  -9.681295080E-09,   3.243728360E-12,
               -1.063943560E+03,   3.657675730E+00] ),
       NASA( [ 1000.00,  6000.00], [  3.660960830E+00,   6.563655230E-04, 
               -1.411494850E-07,   2.057976580E-11,  -1.299132480E-15,
               -1.215977250E+03,   3.415361840E+00] )
             ),
    note = "TPIS89"
       )

species(name = "O2+",
    atoms = " O:2  E:-1 ",
    thermo = (
       NASA( [  200.00,  1000.00], [  4.610171670E+00,  -6.359519520E-03, 
                1.424256240E-05,  -1.209979230E-08,   3.709568780E-12,
                1.397422290E+05,  -2.013268740E-01] ),
       NASA( [ 1000.00,  6000.00], [  3.316759220E+00,   1.115222440E-03, 
               -3.834925560E-07,   5.727846870E-11,  -2.776483810E-15,
                1.398768230E+05,   5.447264760E+00] )
             ),
    note = "TPIS89"
       )