diff --git a/meshLibrary/Make/options b/meshLibrary/Make/options
index fa4782054ab3ac8b358ca728064c580c0356f36f..5c336ce378a32641adb227b87fc3bb1958fe572a 100644
--- a/meshLibrary/Make/options
+++ b/meshLibrary/Make/options
@@ -14,10 +14,6 @@ ifeq (OpenFOAM,$(findstring OpenFOAM,$(WM_PROJECT)))
     endif
 endif
 
-ifeq (v, $(findstring v, $(WM_PROJECT_VERSION)))
-    CFMESH_MACROS += -DOpenCFDSpecific
-endif
-
 ifeq ($(WM_PROJECT), foam)
     VER := $(shell expr `echo $(WM_PROJECT_VERSION)` \>= 3.2)
     ifeq ($(VER), 1)
diff --git a/meshLibrary/utilities/helperClasses/geometry/quadricFitting/quadricFittingI.H b/meshLibrary/utilities/helperClasses/geometry/quadricFitting/quadricFittingI.H
index 98f156841668f80a400549caf2b0472d6f477dd8..e2dbf5f7675d5ce57aa29f22e789df580e067b78 100644
--- a/meshLibrary/utilities/helperClasses/geometry/quadricFitting/quadricFittingI.H
+++ b/meshLibrary/utilities/helperClasses/geometry/quadricFitting/quadricFittingI.H
@@ -52,11 +52,7 @@ void quadricFitting::calculateNormalVector()
 
     //- estimate the normal as the eigenvector associated
     //- to the smallest eigenvalue
-    # ifdef OpenCFDSpecific
     normal_ = eigenVectors(mat, ev).x();
-    # else
-    normal_ = eigenVector(mat, ev[0]);
-    # endif
 }
 
 void quadricFitting::calculateCoordinateSystem()
diff --git a/meshLibrary/utilities/smoothers/geometry/meshOptimizer/tetMeshOptimisation/tetMeshOptimisation.C b/meshLibrary/utilities/smoothers/geometry/meshOptimizer/tetMeshOptimisation/tetMeshOptimisation.C
index 2f07cabc4b8f9c653e31c370859972e30c04fe2a..600d68b108c749c2e99833962a85d16e7b0e76f0 100644
--- a/meshLibrary/utilities/smoothers/geometry/meshOptimizer/tetMeshOptimisation/tetMeshOptimisation.C
+++ b/meshLibrary/utilities/smoothers/geometry/meshOptimizer/tetMeshOptimisation/tetMeshOptimisation.C
@@ -495,11 +495,7 @@ void tetMeshOptimisation::optimiseBoundaryVolumeOptimizer
                         if( mag(ev[2]) > (mag(ev[1]) + mag(ev[0])) )
                         {
                             //- ordinary surface vertex
-                            # ifdef OpenCFDSpecific
                             vector normal = eigenVectors(nt, ev).z();
-                            # else
-                            vector normal = eigenVector(nt, ev[2]);
-                            # endif
 
                             normal /= (mag(normal)+VSMALL);
                             disp -= (disp & normal) * normal;
@@ -507,19 +503,11 @@ void tetMeshOptimisation::optimiseBoundaryVolumeOptimizer
                         else if( mag(ev[1]) > 0.5 * (mag(ev[2]) + mag(ev[0])) )
                         {
                             //- this vertex is on an edge
-                            # ifdef OpenCFDSpecific
                             vector normal1 = eigenVectors(nt, ev).y();
-                            # else
-                            vector normal1 = eigenVector(nt, ev[1]);
-                            # endif
 
                             normal1 /= (mag(normal1)+VSMALL);
 
-                            # ifdef OpenCFDSpecific
                             vector normal2 = eigenVectors(nt, ev).z();
-                            # else
-                            vector normal2 = eigenVector(nt, ev[2]);
-                            # endif
 
                             normal2 /= (mag(normal2)+VSMALL);
 
diff --git a/meshLibrary/utilities/triSurfaceTools/triSurface2DCheck/triSurface2DCheck.C b/meshLibrary/utilities/triSurfaceTools/triSurface2DCheck/triSurface2DCheck.C
index 051efadf36874c832e3b77f75fdcec8f8c39201d..00c31ae9341aff8313cf4281a701125f9558ba1e 100644
--- a/meshLibrary/utilities/triSurfaceTools/triSurface2DCheck/triSurface2DCheck.C
+++ b/meshLibrary/utilities/triSurfaceTools/triSurface2DCheck/triSurface2DCheck.C
@@ -86,19 +86,11 @@ bool triSurface2DCheck::is2DSurface() const
 
     //- calculate the plane normal as a cross prduct of the two
     //- eigenVectors spanning the plane
-    # ifdef OpenCFDSpecific
     const vector n
     (
         eigenVectors(covarianceMatrix_, eigenVal).y() ^
         eigenVectors(covarianceMatrix_, eigenVal).z()
     );
-    # else
-    const vector n
-    (
-        eigenVector(covarianceMatrix_, eigenVal[1]) ^
-        eigenVector(covarianceMatrix_, eigenVal[2])
-    );
-    # endif
 
     //- check if the plane is in the x-y plane of the coordinate system
     if( mag(n.x()) > SMALL || mag(n.y()) > SMALL )