diff --git a/include/cantera/zeroD/ReactorNet.h b/include/cantera/zeroD/ReactorNet.h
index 5a3d8b86e..2ba1b7fed 100644
--- a/include/cantera/zeroD/ReactorNet.h
+++ b/include/cantera/zeroD/ReactorNet.h
@@ -218,15 +218,6 @@ public:
     }
 
 protected:
-    void connect(size_t i, size_t j) {
-        m_connect[j*m_reactors.size() + i] = 1;
-        m_connect[i*m_reactors.size() + j] = 1;
-    }
-
-    bool connected(size_t i, size_t j) {
-        return (m_connect[m_reactors.size()*i + j] == 1);
-    }
-
     /**
      * Initialize the reactor network. Called automatically the first time
      * advance or step is called.
@@ -264,7 +255,6 @@ protected:
     //! output.
     std::vector<size_t> m_sensIndex;
 
-    vector_int m_connect;
     vector_fp m_ydot;
 
     std::vector<bool> m_iown;
diff --git a/src/zeroD/ReactorNet.cpp b/src/zeroD/ReactorNet.cpp
index 769dd2b86..22342ac96 100644
--- a/src/zeroD/ReactorNet.cpp
+++ b/src/zeroD/ReactorNet.cpp
@@ -82,43 +82,7 @@ void ReactorNet::initialize()
         }
     }
 
-    m_connect.resize(m_reactors.size()*m_reactors.size(), 0);
     m_ydot.resize(m_nv,0.0);
-    size_t i, j, nin, nout, nw;
-    ReactorBase* r, *rj;
-    for (i = 0; i < m_reactors.size(); i++) {
-        r = m_reactors[i];
-        for (j = 0; j < m_reactors.size(); j++) {
-            if (i == j) {
-                connect(i,j);
-            } else {
-                rj = m_reactors[j];
-                nin = rj->nInlets();
-                for (n = 0; n < nin; n++) {
-                    if (&rj->inlet(n).out() == r) {
-                        connect(i,j);
-                    }
-                }
-                nout = rj->nOutlets();
-                for (n = 0; n < nout; n++) {
-                    if (&rj->outlet(n).in() == r) {
-                        connect(i,j);
-                    }
-                }
-                nw = rj->nWalls();
-                for (n = 0; n < nw; n++) {
-                    if (&rj->wall(n).left() == rj
-                            && &rj->wall(n).right() == r) {
-                        connect(i,j);
-                    } else if (&rj->wall(n).left() == r
-                               && &rj->wall(n).right() == rj) {
-                        connect(i,j);
-                    }
-                }
-            }
-        }
-    }
-
     m_atol.resize(neq());
     fill(m_atol.begin(), m_atol.end(), m_atols);
     m_integ->setTolerances(m_rtol, neq(), DATA_PTR(m_atol));