diff --git a/include/cantera/oneD/StFlow.h b/include/cantera/oneD/StFlow.h
index 5cd8a6002..0ad50f605 100644
--- a/include/cantera/oneD/StFlow.h
+++ b/include/cantera/oneD/StFlow.h
@@ -131,7 +131,7 @@ public:
      */
     void setMassFraction(size_t j, size_t k, doublereal y) {
         m_fixedy(k,j) = y;
-        m_do_species[k] = true; // false;
+        m_do_species[k] = true;
     }
 
     //! The fixed temperature value at point j.
@@ -469,21 +469,10 @@ protected:
     //             member data
     //---------------------------------------------------------
 
-    // inlet
-    doublereal m_inlet_u;
-    doublereal m_inlet_V;
-    doublereal m_inlet_T;
-    doublereal m_rho_inlet;
-    vector_fp m_yin;
-
-    // surface
-    doublereal m_surface_T;
-
     doublereal m_press;        // pressure
 
     // grid parameters
     vector_fp m_dz;
-    //vector_fp m_z;
 
     // mixture thermo properties
     vector_fp m_rho;
@@ -492,7 +481,6 @@ protected:
     // species thermo properties
     vector_fp m_wt;
     vector_fp m_cp;
-    vector_fp m_enth;
 
     // transport properties
     vector_fp m_visc;
@@ -504,7 +492,6 @@ protected:
 
     // production rates
     Array2D m_wdot;
-    vector_fp m_surfdot;
 
     size_t m_nsp;
 
@@ -514,8 +501,6 @@ protected:
 
     MultiJac* m_jac;
 
-    bool m_ok;
-
     // boundary emissivities for the radiation calculations
     doublereal m_epsilon_left;
     doublereal m_epsilon_right;
@@ -537,10 +522,6 @@ protected:
     // vector which contains the values of the radiative heat loss
     vector_fp m_qdotRadiation;
 
-    // solution estimate
-    //vector_fp m_zest;
-    //Array2D   m_yest;
-
     // fixed T and Y values
     Array2D   m_fixedy;
     vector_fp m_fixedtemp;
diff --git a/src/oneD/StFlow.cpp b/src/oneD/StFlow.cpp
index 2a4871c50..384944626 100644
--- a/src/oneD/StFlow.cpp
+++ b/src/oneD/StFlow.cpp
@@ -18,17 +18,12 @@ namespace Cantera
 
 StFlow::StFlow(IdealGasPhase* ph, size_t nsp, size_t points) :
     Domain1D(nsp+4, points),
-    m_inlet_u(0.0),
-    m_inlet_V(0.0),
-    m_inlet_T(-1.0),
-    m_surface_T(-1.0),
     m_press(-1.0),
     m_nsp(nsp),
     m_thermo(0),
     m_kin(0),
     m_trans(0),
     m_jac(0),
-    m_ok(false),
     m_epsilon_left(0.0),
     m_epsilon_right(0.0),
     m_do_soret(false),
@@ -69,7 +64,6 @@ StFlow::StFlow(IdealGasPhase* ph, size_t nsp, size_t points) :
     m_multidiff.resize(m_nsp*m_nsp*m_points);
     m_flux.resize(m_nsp,m_points);
     m_wdot.resize(m_nsp,m_points, 0.0);
-    m_surfdot.resize(m_nsp, 0.0);
     m_ybar.resize(m_nsp);
     m_qdotRadiation.resize(m_points, 0.0);
 
@@ -116,7 +110,6 @@ void StFlow::resize(size_t ncomponents, size_t points)
     m_rho.resize(m_points, 0.0);
     m_wtm.resize(m_points, 0.0);
     m_cp.resize(m_points, 0.0);
-    m_enth.resize(m_points, 0.0);
     m_visc.resize(m_points, 0.0);
     m_tcon.resize(m_points, 0.0);