From 36f32540100e0551e6536626206c9182b4114d38 Mon Sep 17 00:00:00 2001
From: Dave Goodwin <dgg@caltech.edu>
Date: Tue, 7 Nov 2006 13:47:39 +0000
Subject: [PATCH] added method concentration(int k)

---
 Cantera/src/State.cpp | 10 ++++++++++
 Cantera/src/State.h   |  6 ++++++
 2 files changed, 16 insertions(+)

diff --git a/Cantera/src/State.cpp b/Cantera/src/State.cpp
index 9ec05ebe2..b74c65582 100644
--- a/Cantera/src/State.cpp
+++ b/Cantera/src/State.cpp
@@ -107,6 +107,16 @@ namespace Cantera {
         }
     }
 
+    doublereal State::concentration(int k) const {
+        if (k >= 0 && k < m_kk) {
+            return m_y[k] * m_dens * m_rmolwts[k] ;
+        }
+        else {
+            throw CanteraError("State:massFraction",
+                "illegal species index number");
+        }
+    }
+
     void State::setMassFractions(const doublereal* y) {
         doublereal norm = 0.0, sum = 0.0;
         int k;
diff --git a/Cantera/src/State.h b/Cantera/src/State.h
index 04089e0f9..a3d6c1ec8 100755
--- a/Cantera/src/State.h
+++ b/Cantera/src/State.h
@@ -156,6 +156,12 @@ namespace Cantera {
          */
 	void getConcentrations(doublereal* c) const;
 
+        /**
+         * Concentration of species k. If k is outside the valid
+         * range, an exception will be thrown.
+         */
+        doublereal concentration(int k) const;
+
 	/**
 	 * Set the concentrations to the specified values within the
 	 * phase.