diff --git a/SimPEG/Tests/test_TreeMesh.py b/SimPEG/Tests/test_TreeMesh.py
index 2d5367c1..cbeeaf6a 100644
--- a/SimPEG/Tests/test_TreeMesh.py
+++ b/SimPEG/Tests/test_TreeMesh.py
@@ -249,7 +249,6 @@ class TestOcTreeObjects(unittest.TestCase):
         self.assertTrue(q('c4fXp') is q('c2fXm').children[1,0])
         self.assertTrue(q('c4fXp').parent is q('c2fXm'))
 
-        #TEST EDGES!
 
     def test_gridCC(self):
         x = np.r_[0.25,0.75]
@@ -290,30 +289,55 @@ class TestOcTreeObjects(unittest.TestCase):
         z = np.r_[0.5,0.5,0.5,0.5]
         self.assertTrue(np.linalg.norm((np.c_[x,y,z]-self.M.gridFy).flatten()) == 0)
 
+        x = np.r_[0.125,0.375,0.75,0.125,0.375,0.125,0.375,0.75,0.125,0.375,0.125,0.375,0.125,0.375]
+        y = np.r_[0,0,0,0.5,0.5,1,1,1,0,0,0.5,0.5,1,1]
+        z = np.r_[0.25,0.25,0.5,0.25,0.25,0.25,0.25,0.5,0.75,0.75,0.75,0.75,0.75,0.75]
+        self.assertTrue(np.linalg.norm((np.c_[x,y,z]-self.Mr.gridFy).flatten()) == 0)
+
     def test_gridFz(self):
         x = np.r_[0.25,0.75,0.25,0.75]
         y = np.r_[0.5,0.5,0.5,0.5]
         z = np.r_[0,0,1.,1.]
         self.assertTrue(np.linalg.norm((np.c_[x,y,z]-self.M.gridFz).flatten()) == 0)
 
+        x = np.r_[0.125,0.375,0.75,0.125,0.375,0.125,0.375,0.125,0.375,0.125,0.375,0.75,0.125,0.375]
+        y = np.r_[0.25,0.25,0.5,0.75,0.75,0.25,0.25,0.75,0.75,0.25,0.25,0.5,0.75,0.75]
+        z = np.r_[0,0,0,0,0,0.5,0.5,0.5,0.5,1,1,1,1,1]
+        self.assertTrue(np.linalg.norm((np.c_[x,y,z]-self.Mr.gridFz).flatten()) == 0)
+
+
     def test_gridEx(self):
         x = np.r_[0.25,0.75,0.25,0.75,0.25,0.75,0.25,0.75]
         y = np.r_[0,0,1.,1.,0,0,1.,1.]
         z = np.r_[0,0,0,0,1.,1.,1.,1.]
         self.assertTrue(np.linalg.norm((np.c_[x,y,z]-self.M.gridEx).flatten()) == 0)
 
+        x = np.r_[0.125,0.375,0.75,0.125,0.375,0.125,0.375,0.75,0.125,0.375,0.125,0.375,0.125,0.375,0.125,0.375,0.75,0.125,0.375,0.125,0.375,0.75]
+        y = np.r_[0,0,0,0.5,0.5,1,1,1,0,0,0.5,0.5,1,1,0,0,0,0.5,0.5,1,1,1]
+        z = np.r_[0,0,0,0,0,0,0,0,0.5,0.5,0.5,0.5,0.5,0.5,1,1,1,1,1,1,1,1]
+        self.assertTrue(np.linalg.norm((np.c_[x,y,z]-self.Mr.gridEx).flatten()) == 0)
+
     def test_gridEy(self):
         x = np.r_[0,0.5,1,0,0.5,1]
         y = np.r_[0.5,0.5,0.5,0.5,0.5,0.5]
         z = np.r_[0,0,0,1.,1.,1.]
         self.assertTrue(np.linalg.norm((np.c_[x,y,z]-self.M.gridEy).flatten()) == 0)
 
+        x = np.r_[0,0.25,0.5,1,0,0.25,0.5,0,0.25,0.5,0,0.25,0.5,0,0.25,0.5,1,0,0.25,0.5]
+        y = np.r_[0.25,0.25,0.25,0.5,0.75,0.75,0.75,0.25,0.25,0.25,0.75,0.75,0.75,0.25,0.25,0.25,0.5,0.75,0.75,0.75]
+        z = np.r_[0,0,0,0,0,0,0,0.5,0.5,0.5,0.5,0.5,0.5,1,1,1,1,1,1,1]
+        self.assertTrue(np.linalg.norm((np.c_[x,y,z]-self.Mr.gridEy).flatten()) == 0)
+
     def test_gridEz(self):
         x = np.r_[0,0.5,1,0,0.5,1]
         y = np.r_[0,0,0,1.,1.,1.]
         z = np.r_[0.5,0.5,0.5,0.5,0.5,0.5]
         self.assertTrue(np.linalg.norm((np.c_[x,y,z]-self.M.gridEz).flatten()) == 0)
 
+        x = np.r_[0,0.25,0.5,1,0  ,0.25,0.5,0,0.25,0.5,1,0,0.25,0.5,0  ,0.25,0.5,0  ,0.25,0.5]
+        y = np.r_[0,0   ,0  ,0,0.5,0.5 ,0.5,1,1   ,1  ,1,0,0   ,0  ,0.5,0.5 ,0.5,1  ,1   ,1  ]
+        z = np.r_[0.25,0.25,0.25,0.5,0.25,0.25,0.25,0.25,0.25,0.25,0.5,0.75,0.75,0.75,0.75,0.75,0.75,0.75,0.75,0.75]
+        self.assertTrue(np.linalg.norm((np.c_[x,y,z]-self.Mr.gridEz).flatten()) == 0)
 
 
 class TestQuadTreeObjects(unittest.TestCase):