diff --git a/opm/simulators/wells/WellHelpers.cpp b/opm/simulators/wells/WellHelpers.cpp
index 41cfc5f60..676ded12f 100644
--- a/opm/simulators/wells/WellHelpers.cpp
+++ b/opm/simulators/wells/WellHelpers.cpp
@@ -25,8 +25,6 @@
 
 #include <opm/common/OpmLog/OpmLog.hpp>
 
-#include <opm/grid/cpgrid/GridHelpers.hpp>
-
 #include <opm/simulators/wells/ParallelWellInfo.hpp>
 
 #include <vector>
@@ -159,47 +157,6 @@ void sumDistributedWellEntries(Dune::DynamicMatrix<Scalar>& mat,
     std::copy(pos, allEntries.end(), &(vec[0]));
 }
 
-template <int dim, class C2F, class FC>
-std::array<double, dim>
-getCubeDim(const C2F& c2f,
-           FC         begin_face_centroids,
-           int        cell)
-{
-    std::array<std::vector<double>, dim> X;
-    {
-        const std::vector<double>::size_type
-            nf = std::distance(c2f[cell].begin(),
-                               c2f[cell].end  ());
-
-        for (int d = 0; d < dim; ++d) {
-            X[d].reserve(nf);
-        }
-    }
-
-    for (const auto& f : c2f[cell]) {
-        using Opm::UgGridHelpers::increment;
-        using Opm::UgGridHelpers::getCoordinate;
-
-        const FC& fc = increment(begin_face_centroids, f, dim);
-
-        for (int d = 0; d < dim; ++d) {
-            X[d].push_back(getCoordinate(fc, d));
-        }
-    }
-
-    std::array<double, dim> cube;
-    for (int d = 0; d < dim; ++d) {
-        using VI = std::vector<double>::iterator;
-        using PVI = std::pair<VI,VI>;
-
-        const PVI m = std::minmax_element(X[d].begin(), X[d].end());
-
-        cube[d] = *m.second - *m.first;
-    }
-
-    return cube;
-}
-
 template <class DenseMatrix>
 DenseMatrix transposeDenseDynMatrix(const DenseMatrix& M)
 {
diff --git a/opm/simulators/wells/WellHelpers.hpp b/opm/simulators/wells/WellHelpers.hpp
index 0f661d050..4ba8ab2c7 100644
--- a/opm/simulators/wells/WellHelpers.hpp
+++ b/opm/simulators/wells/WellHelpers.hpp
@@ -78,13 +78,6 @@ void sumDistributedWellEntries(Dune::DynamicMatrix<Scalar>& mat,
                                const Comm& comm);
 
 
-
-template <int dim, class C2F, class FC>
-std::array<double, dim>
-getCubeDim(const C2F& c2f,
-           FC         begin_face_centroids,
-           int        cell);
-
 // explicit transpose of a dense matrix due to compilation problems
 // used for calculating quasiimpes well weights
 template <class DenseMatrix>