diff --git a/src/equil/ChemEquil.cpp b/src/equil/ChemEquil.cpp
index e3507a94a1..b42b491149 100644
--- a/src/equil/ChemEquil.cpp
+++ b/src/equil/ChemEquil.cpp
@@ -1027,10 +1027,8 @@ int ChemEquil::estimateEP_Brinkley(ThermoPhase& s, vector_fp& x,
             for (size_t m = 0; m < m_mm; m++) {
                 size_t kMSp = npos;
                 size_t kMSp2 = npos;
-                int nSpeciesWithElem = 0;
                 for (size_t k = 0; k < m_kk; k++) {
                     if (n_i_calc[k] > nCutoff && fabs(nAtoms(k,m)) > 0.001) {
-                        nSpeciesWithElem++;
                         if (kMSp != npos) {
                             kMSp2 = k;
                             double factor = fabs(nAtoms(kMSp,m) / nAtoms(kMSp2,m));
diff --git a/src/oneD/Sim1D.cpp b/src/oneD/Sim1D.cpp
index 07f27ade1c..0e1a4de72f 100644
--- a/src/oneD/Sim1D.cpp
+++ b/src/oneD/Sim1D.cpp
@@ -452,7 +452,6 @@ void Sim1D::solve(int loglevel, bool refine_grid)
     int new_points = 1;
     doublereal dt = m_tstep;
     m_nsteps = 0;
-    int soln_number = -1;
     finalize();
 
     while (new_points > 0) {
@@ -495,7 +494,6 @@ void Sim1D::solve(int loglevel, bool refine_grid)
                                  "After successful Newton solve");
                 }
                 ok = true;
-                soln_number++;
             } else {
                 debuglog("    failure. \n", loglevel);
                 if (loglevel > 6) {