diff --git a/SimPEG/__init__.py b/SimPEG/__init__.py
index ea9faf33..e8946f88 100644
--- a/SimPEG/__init__.py
+++ b/SimPEG/__init__.py
@@ -1,5 +1,4 @@
-from TensorMesh import TensorMesh
-from LogicallyOrthogonalMesh import LogicallyOrthogonalMesh
+import mesh
 import utils
 import inverse
 from Solver import Solver
diff --git a/SimPEG/forward/DCProblem/DCProblem.py b/SimPEG/forward/DCProblem/DCProblem.py
index b8115598..1df8897b 100644
--- a/SimPEG/forward/DCProblem/DCProblem.py
+++ b/SimPEG/forward/DCProblem/DCProblem.py
@@ -1,4 +1,4 @@
-from SimPEG import TensorMesh
+from SimPEG.mesh import TensorMesh
 from SimPEG.forward import Problem, SyntheticProblem
 from SimPEG.tests import checkDerivative
 from SimPEG.utils import ModelBuilder, sdiag
diff --git a/SimPEG/BaseMesh.py b/SimPEG/mesh/BaseMesh.py
similarity index 99%
rename from SimPEG/BaseMesh.py
rename to SimPEG/mesh/BaseMesh.py
index c57828e9..29c43dd1 100644
--- a/SimPEG/BaseMesh.py
+++ b/SimPEG/mesh/BaseMesh.py
@@ -1,5 +1,5 @@
 import numpy as np
-from utils import mkvc
+from SimPEG.utils import mkvc
 
 
 class BaseMesh(object):
diff --git a/SimPEG/mesh/Cyl1DMesh.py b/SimPEG/mesh/Cyl1DMesh.py
new file mode 100644
index 00000000..5fc4deb6
--- /dev/null
+++ b/SimPEG/mesh/Cyl1DMesh.py
@@ -0,0 +1,330 @@
+import numpy as np
+import scipy.sparse as sp
+from scipy.constants import pi
+from SimPEG.utils import mkvc, ndgrid, sdiag
+
+class Cyl1DMesh(object):
+    """
+    Cyl1DMesh is a mesh class for cylindrically symmetric 1D problems  
+    """ 
+
+    _meshType = 'CYL1D'
+
+    def __init__(self, h, z0=None):
+        assert len(h) == 2, "len(h) must equal 2"
+        if z0 is not None:
+            assert z0.size == 1, "z0.size must equal 1"
+
+        for i, h_i in enumerate(h):
+            assert type(h_i) == np.ndarray, ("h[%i] is not a numpy array." % i)
+            assert len(h_i.shape) == 1, ("h[%i] must be a 1D numpy array." % i)
+
+        # Ensure h contains 1D vectors
+        self._h = [mkvc(x) for x in h]
+
+        if z0 is None:
+            z0 = 0
+        self._z0 = z0
+
+    ####################################################
+    # Mesh properties
+    ####################################################
+
+    def h():
+        doc = "list containing the width of each cell"
+        def fget(self):
+            return self._h
+        return locals()
+    h = property(**h())
+
+    def z0():
+        doc = "The z-origin"
+        def fget(self):
+            return self._z0
+        return locals()
+    z0 = property(**z0())
+
+    def hr():
+        doc = "Width of the cells in the r direction"
+        def fget(self):
+            return self._h[0]
+        return locals()
+    hr = property(**hr())
+
+    def hz():
+        doc = "Width of the cells in the z direction"
+        def fget(self):
+            return self._h[1]
+        return locals()
+    hz = property(**hz())
+
+    ####################################################
+    # Counting
+    ####################################################
+
+    def nCr():
+        doc = "Number of cells in the radial direction"
+        fget = lambda self: self.hr.size
+        return locals()
+    nCr = property(**nCr())
+
+    def nCz():
+        doc = "Number of cells in the z direction"
+        fget = lambda self: self.hz.size
+        return locals()
+    nCz = property(**nCz())
+
+    def nC():
+        doc = "Total number of cells"
+        fget = lambda self: self.nCr * self.nCz
+        return locals()
+    nC = property(**nC())
+
+    def nNr():
+        doc = "Number of nodes in the radial direction"
+        fget = lambda self: self.hr.size
+        return locals()
+    nNr = property(**nNr())
+
+    def nNz():
+        doc = "Number of nodes in the radial direction"
+        fget = lambda self: self.hz.size + 1
+        return locals()
+    nNz = property(**nNz())
+
+    def nN():
+        doc = "Total number of nodes"
+        fget = lambda self: self.nNr * self.nNz
+        return locals()
+    nN = property(**nN())
+
+    def nFr():
+        doc = "Number of r faces"
+        fget = lambda self: self.nNr * self.nCz
+        return locals()
+    nFr = property(**nFr())
+
+    def nFz():
+        doc = "Number of z faces"
+        fget = lambda self: self.nNz * self.nCr
+        return locals()
+    nFz = property(**nFz())
+
+    def nF():
+        doc = "Total number of faces"
+        fget = lambda self: self.nFr + self.nFz
+        return locals()
+    nF = property(**nF())
+
+    def nE():
+        doc = "Number of edges"
+        fget = lambda self: self.nN
+        return locals()
+    nE = property(**nE())
+
+    ####################################################
+    # Vectors & Grids
+    ####################################################
+
+    def vectorNr():
+        doc = "Nodal grid vector (1D) in the r direction"
+        fget = lambda self: self.hr.cumsum()
+        return locals()
+    vectorNr = property(**vectorNr())
+
+    def vectorNz():
+        doc = "Nodal grid vector (1D) in the z direction"
+        fget = lambda self: np.r_[0, self.hz.cumsum()] + self._z0
+        return locals()
+    vectorNz = property(**vectorNz())
+
+    def vectorCCr():
+        doc = "Cell centered grid vector (1D) in the r direction"
+        fget = lambda self: np.r_[0, self.hr.cumsum()[1:] - self.hr[1:]/2]
+        return locals()
+    vectorCCr = property(**vectorCCr())
+
+    def vectorCCz():
+        doc = "Cell centered grid vector (1D) in the z direction"
+        fget = lambda self: self.hz.cumsum() - self.hz/2 + self._z0 
+        return locals()
+    vectorCCz = property(**vectorCCz())
+
+    def gridCC():
+        doc = "Cell-centered grid"
+        def fget(self):
+            if self._gridCC is None:
+                self._gridCC = ndgrid([self.vectorCCr, self.vectorCCz])
+            return self._gridCC
+        return locals()
+    _gridCC = None
+    gridCC = property(**gridCC())
+
+    def gridN():
+        doc = "Nodal grid"
+        def fget(self):
+            if self._gridN is None:
+                self._gridN = ndgrid([self.vectorNr, self.vectorNz])
+            return self._gridN
+        return locals()
+    _gridN = None
+    gridN = property(**gridN())
+
+    def gridFr():
+        doc = "r face grid"
+        def fget(self):
+            if self._gridFr is None:
+                self._gridFr = ndgrid([self.vectorNr, self.vectorCCz])
+            return self._gridFr
+        return locals()
+    _gridFr = None    
+    gridFr = property(**gridFr())
+
+    def gridFz():
+        doc = "z face grid"
+        def fget(self):
+            if self._gridFz is None:
+                self._gridFz = ndgrid([self.vectorCCr, self.vectorNz])
+            return self._gridFz
+        return locals()
+    _gridFz = None    
+    gridFz = property(**gridFz())
+
+    ####################################################
+    # Geometries
+    ####################################################
+
+    def edge():
+        doc = "Edge lengths"
+        def fget(self):
+            if self._edge is None:
+                self._edge = 2*pi*self.gridN[:,0]
+            return self._edge
+        return locals()
+    _edge = None
+    edge = property(**edge())
+
+    def area():
+        doc = "Face areas"
+        def fget(self):
+            if self._area is None:
+                areaR = np.kron(self.hz, 2*pi*self.vectorNr)
+                areaZ = np.kron(np.ones_like(self.vectorNz),pi*(self.vectorNr**2 - np.r_[0, self.vectorNr[:-1]]**2))
+                self._area = np.r_[areaR, areaZ]
+            return self._area
+        return locals()
+    _area = None
+    area = property(**area())
+
+    def vol():
+        doc = "Volume of each cell"
+        def fget(self):
+            if self._vol is None:
+                az = pi*(self.vectorNr**2 - np.r_[0, self.vectorNr[:-1]]**2)
+                self._vol = np.kron(self.hz,az)
+            return self._vol
+        return locals()
+    _vol = None
+    vol = property(**vol())
+
+    ####################################################
+    # Operators
+    ####################################################
+
+    def edgeCurl():
+        doc = "The edgeCurl property."
+        def fget(self):
+            if self._edgeCurl is None:
+                #1D Difference matricies
+                dr = sp.spdiags((np.ones((self.nCr+1, 1))*[-1, 1]).T, [-1,0], self.nCr, self.nCr, format="csr")
+                dz = sp.spdiags((np.ones((self.nCz+1, 1))*[-1, 1]).T, [0,1], self.nCz, self.nCz+1, format="csr")
+
+                #2D Difference matricies
+                Dr = sp.kron(sp.eye(self.nNz), dr)
+                Dz = -sp.kron(dz, sp.eye(self.nCr))  #Not sure about this negative
+
+                #Edge curl operator
+                self._edgeCurl = sp.diags(1/self.area,0)*sp.vstack((Dz, Dr))*sp.diags(self.edge,0)
+            return self._edgeCurl
+        return locals()
+    _edgeCurl = None
+    edgeCurl = property(**edgeCurl())
+
+    def aveE2CC():
+        doc = "Averaging operator from cell edges to cell centres"
+        def fget(self):
+            if self._aveE2CC is None:
+                az = sp.spdiags(0.5*np.ones((2, self.nNz)), [-1,0], self.nNz, self.nCz, format='csr')
+                ar = sp.spdiags(0.5*np.ones((2, self.nCr)), [0, 1], self.nCr, self.nCr, format='csr')
+                ar[0,0] = 1
+                self._aveE2CC = sp.kron(az, ar).T
+            return self._aveE2CC
+        return locals()
+    _aveE2CC = None
+    aveE2CC = property(**aveE2CC())
+
+    def aveF2CC():
+        doc = "Averaging operator from cell faces to cell centres"
+        def fget(self):
+            if self._aveF2CC is None:
+                az = sp.spdiags(0.5*np.ones((2, self.nNz)), [-1,0], self.nNz, self.nCz, format='csr')
+                ar = sp.spdiags(0.5*np.ones((2, self.nCr)), [0, 1], self.nCr, self.nCr, format='csr')
+                ar[0,0] = 1
+                Afr = sp.kron(sp.eye(self.nCz),ar)
+                Afz = sp.kron(az,sp.eye(self.nCr))
+                self._aveF2CC = sp.vstack((Afr,Afz)).T
+            return self._aveF2CC
+        return locals()
+    _aveF2CC = None
+    aveF2CC = property(**aveF2CC())
+
+    ####################################################
+    # Methods
+    ####################################################
+
+
+    def getMass(self, materialProp=None, loc='e'):
+        """ Produces mass matricies.
+
+        :param str loc: Average to location: 'e'-edges, 'f'-faces
+        :param None,float,numpy.ndarray materialProp: property to be averaged (see below)
+        :rtype: scipy.sparse.csr.csr_matrix
+        :return: M, the mass matrix
+
+        materialProp can be::
+
+            None            -> takes materialProp = 1 (default)
+            float           -> a constant value for entire domain
+            numpy.ndarray   -> if materialProp.size == self.nC
+                                    3D property model
+                               if materialProp.size = self.nCz
+                                    1D (layered eath) property model
+        """
+        if materialProp is None:
+            materialProp = np.ones(self.nC)
+        elif type(materialProp) is float:
+            materialProp = np.ones(self.nC)*materialProp
+        elif materialProp.shape == (self.nCz,):
+            materialProp = materialProp.repeat(self.nCr)
+        materialProp = mkvc(materialProp)
+        assert materialProp.shape == (self.nC,), "materialProp incorrect shape"
+
+        if loc=='e':
+            Av = self.aveE2CC
+        elif loc=='f':
+            Av = self.aveF2CC
+        else:
+            raise ValueError('Invalid loc')
+
+        diag = Av.T * (self.vol * mkvc(materialProp))
+
+        return sdiag(diag)
+
+    def getEdgeMass(self, materialProp=None):
+        """mass matrix for products of edge functions w'*M(materialProp)*e"""
+        return self.getMass(loc='e', materialProp=materialProp)
+
+    def getFaceMass(self, materialProp=None):
+        """mass matrix for products of face functions w'*M(materialProp)*f"""
+        return self.getMass(loc='f', materialProp=materialProp)
+
diff --git a/SimPEG/DiffOperators.py b/SimPEG/mesh/DiffOperators.py
similarity index 100%
rename from SimPEG/DiffOperators.py
rename to SimPEG/mesh/DiffOperators.py
diff --git a/SimPEG/InnerProducts.py b/SimPEG/mesh/InnerProducts.py
similarity index 100%
rename from SimPEG/InnerProducts.py
rename to SimPEG/mesh/InnerProducts.py
diff --git a/SimPEG/LogicallyOrthogonalMesh.py b/SimPEG/mesh/LogicallyOrthogonalMesh.py
similarity index 99%
rename from SimPEG/LogicallyOrthogonalMesh.py
rename to SimPEG/mesh/LogicallyOrthogonalMesh.py
index 5d2b7ffa..7b1e1bca 100644
--- a/SimPEG/LogicallyOrthogonalMesh.py
+++ b/SimPEG/mesh/LogicallyOrthogonalMesh.py
@@ -3,7 +3,7 @@ from BaseMesh import BaseMesh
 from DiffOperators import DiffOperators
 from InnerProducts import InnerProducts
 from LomView import LomView
-from utils import mkvc, ndgrid, volTetra, indexCube, faceInfo
+from SimPEG.utils import mkvc, ndgrid, volTetra, indexCube, faceInfo
 
 # Some helper functions.
 length2D = lambda x: (x[:, 0]**2 + x[:, 1]**2)**0.5
diff --git a/SimPEG/LomView.py b/SimPEG/mesh/LomView.py
similarity index 99%
rename from SimPEG/LomView.py
rename to SimPEG/mesh/LomView.py
index 4b4f36a3..2a9b242b 100644
--- a/SimPEG/LomView.py
+++ b/SimPEG/mesh/LomView.py
@@ -2,7 +2,7 @@ import numpy as np
 import matplotlib.pyplot as plt
 import matplotlib
 from mpl_toolkits.mplot3d import Axes3D
-from utils import mkvc
+from SimPEG.utils import mkvc
 
 
 class LomView(object):
diff --git a/SimPEG/TensorMesh.py b/SimPEG/mesh/TensorMesh.py
similarity index 96%
rename from SimPEG/TensorMesh.py
rename to SimPEG/mesh/TensorMesh.py
index 2015838c..a3329c14 100644
--- a/SimPEG/TensorMesh.py
+++ b/SimPEG/mesh/TensorMesh.py
@@ -3,7 +3,7 @@ from BaseMesh import BaseMesh
 from TensorView import TensorView
 from DiffOperators import DiffOperators
 from InnerProducts import InnerProducts
-from utils import ndgrid, mkvc
+from SimPEG.utils import ndgrid, mkvc
 
 
 class TensorMesh(BaseMesh, TensorView, DiffOperators, InnerProducts):
diff --git a/SimPEG/TensorView.py b/SimPEG/mesh/TensorView.py
similarity index 99%
rename from SimPEG/TensorView.py
rename to SimPEG/mesh/TensorView.py
index ba4f1b63..687a520d 100644
--- a/SimPEG/TensorView.py
+++ b/SimPEG/mesh/TensorView.py
@@ -2,7 +2,7 @@ import numpy as np
 import matplotlib.pyplot as plt
 import matplotlib
 from mpl_toolkits.mplot3d import Axes3D
-from utils import mkvc
+from SimPEG.utils import mkvc
 
 
 class TensorView(object):
diff --git a/SimPEG/mesh/__init__.py b/SimPEG/mesh/__init__.py
new file mode 100644
index 00000000..3b8e1eef
--- /dev/null
+++ b/SimPEG/mesh/__init__.py
@@ -0,0 +1,8 @@
+from Cyl1DMesh import Cyl1DMesh
+from TensorMesh import TensorMesh
+from LogicallyOrthogonalMesh import LogicallyOrthogonalMesh
+from BaseMesh import BaseMesh
+from TensorView import TensorView
+from LomView import LomView
+from InnerProducts import InnerProducts
+from DiffOperators import DiffOperators
diff --git a/SimPEG/tests/TestUtils.py b/SimPEG/tests/TestUtils.py
index 79340da7..9b2158c4 100644
--- a/SimPEG/tests/TestUtils.py
+++ b/SimPEG/tests/TestUtils.py
@@ -2,7 +2,8 @@ import numpy as np
 import matplotlib.pyplot as plt
 from pylab import norm
 from SimPEG.utils import mkvc
-from SimPEG import TensorMesh, utils, LogicallyOrthogonalMesh
+from SimPEG import utils
+from SimPEG.mesh import TensorMesh, LogicallyOrthogonalMesh
 import numpy as np
 import unittest
 
diff --git a/SimPEG/tests/test_LogicallyOrthogonalMesh.py b/SimPEG/tests/test_LogicallyOrthogonalMesh.py
index 523b18d2..7f656aae 100644
--- a/SimPEG/tests/test_LogicallyOrthogonalMesh.py
+++ b/SimPEG/tests/test_LogicallyOrthogonalMesh.py
@@ -1,10 +1,7 @@
 import numpy as np
 import unittest
-import sys
-sys.path.append('../')
-from TensorMesh import TensorMesh
-from LogicallyOrthogonalMesh import LogicallyOrthogonalMesh
-from utils import ndgrid
+from SimPEG.mesh import TensorMesh, LogicallyOrthogonalMesh
+from SimPEG.utils import ndgrid
 
 
 class BasicLOMTests(unittest.TestCase):
diff --git a/SimPEG/tests/test_basemesh.py b/SimPEG/tests/test_basemesh.py
index 60373011..5762de57 100644
--- a/SimPEG/tests/test_basemesh.py
+++ b/SimPEG/tests/test_basemesh.py
@@ -1,7 +1,6 @@
 import unittest
 import sys
-sys.path.append('../')
-from BaseMesh import BaseMesh
+from SimPEG.mesh import BaseMesh
 import numpy as np
 
 
diff --git a/SimPEG/tests/test_forward_DCproblem.py b/SimPEG/tests/test_forward_DCproblem.py
index b5e6d844..6ffa4650 100644
--- a/SimPEG/tests/test_forward_DCproblem.py
+++ b/SimPEG/tests/test_forward_DCproblem.py
@@ -1,6 +1,6 @@
 import numpy as np
 import unittest
-from SimPEG import TensorMesh
+from SimPEG.mesh import TensorMesh
 from SimPEG.utils import ModelBuilder, sdiag
 from SimPEG.forward import Problem, SyntheticProblem
 from SimPEG.forward.DCProblem import DCProblem, DCutils
diff --git a/SimPEG/tests/test_forward_problem.py b/SimPEG/tests/test_forward_problem.py
index b7180dca..d913a697 100644
--- a/SimPEG/tests/test_forward_problem.py
+++ b/SimPEG/tests/test_forward_problem.py
@@ -1,6 +1,6 @@
 import numpy as np
 import unittest
-from SimPEG import TensorMesh
+from SimPEG.mesh import TensorMesh
 from SimPEG.forward import Problem
 from TestUtils import checkDerivative
 from scipy.sparse.linalg import dsolve
diff --git a/SimPEG/tests/test_tensorMesh.py b/SimPEG/tests/test_tensorMesh.py
index 300e6634..9544dab6 100644
--- a/SimPEG/tests/test_tensorMesh.py
+++ b/SimPEG/tests/test_tensorMesh.py
@@ -1,6 +1,6 @@
 import numpy as np
 import unittest
-from SimPEG import TensorMesh
+from SimPEG.mesh import TensorMesh
 from TestUtils import OrderTest
 from scipy.sparse.linalg import dsolve
 
diff --git a/SimPEG/utils/ModelBuilder.py b/SimPEG/utils/ModelBuilder.py
index 170da1bf..d0512680 100644
--- a/SimPEG/utils/ModelBuilder.py
+++ b/SimPEG/utils/ModelBuilder.py
@@ -15,6 +15,7 @@ def getIndecesBlock(p0,p1,ccMesh):
         ccMesh represents the cell-centered mesh
 
         The points p0 and p1 must live in the the same dimensional space as the mesh.
+
     """
 
     # Validation: p0 and p1 live in the same dimensional space
@@ -131,7 +132,7 @@ def scalarConductivity(ccMesh,pFunction):
 
 if __name__ == '__main__':
 
-    from SimPEG import TensorMesh
+    from SimPEG.mesh import TensorMesh
     from matplotlib import pyplot as plt
 
     # Define the mesh
diff --git a/docs/api_BaseMesh.rst b/docs/api_BaseMesh.rst
index 66249953..e5f81299 100644
--- a/docs/api_BaseMesh.rst
+++ b/docs/api_BaseMesh.rst
@@ -3,6 +3,6 @@
 Base Mesh
 *********
 
-.. automodule:: SimPEG.BaseMesh
+.. automodule:: SimPEG.mesh.BaseMesh
     :members:
     :undoc-members:
diff --git a/docs/api_Cyl1DMesh.rst b/docs/api_Cyl1DMesh.rst
new file mode 100644
index 00000000..f0109b0d
--- /dev/null
+++ b/docs/api_Cyl1DMesh.rst
@@ -0,0 +1,8 @@
+.. _api_Cyl1DMesh:
+
+Cylindrical 1D Mesh
+*******************
+
+.. automodule:: SimPEG.mesh.Cyl1DMesh
+    :members:
+    :undoc-members:
\ No newline at end of file
diff --git a/docs/api_DiffOperators.rst b/docs/api_DiffOperators.rst
index 57763a33..022fe564 100644
--- a/docs/api_DiffOperators.rst
+++ b/docs/api_DiffOperators.rst
@@ -3,6 +3,6 @@
 Differential Operators
 **********************
 
-.. automodule:: SimPEG.DiffOperators
+.. automodule:: SimPEG.mesh.DiffOperators
     :members:
     :undoc-members:
diff --git a/docs/api_InnerProducts.rst b/docs/api_InnerProducts.rst
index 6dfd2a77..90642807 100644
--- a/docs/api_InnerProducts.rst
+++ b/docs/api_InnerProducts.rst
@@ -3,6 +3,6 @@
 Inner Products
 **************
 
-.. automodule:: SimPEG.InnerProducts
+.. automodule:: SimPEG.mesh.InnerProducts
     :members:
     :undoc-members:
diff --git a/docs/api_LOMView.rst b/docs/api_LOMView.rst
index 61630c26..cfbfa6fa 100644
--- a/docs/api_LOMView.rst
+++ b/docs/api_LOMView.rst
@@ -3,6 +3,6 @@
 LOM View
 ********
 
-.. automodule:: SimPEG.LomView
+.. automodule:: SimPEG.mesh.LomView
     :members:
     :undoc-members:
diff --git a/docs/api_LogicallyOrthogonalMesh.rst b/docs/api_LogicallyOrthogonalMesh.rst
index 3b8296d7..af9e73d0 100644
--- a/docs/api_LogicallyOrthogonalMesh.rst
+++ b/docs/api_LogicallyOrthogonalMesh.rst
@@ -3,6 +3,6 @@
 Logically Orthogonal Mesh
 *************************
 
-.. automodule:: SimPEG.LogicallyOrthogonalMesh
+.. automodule:: SimPEG.mesh.LogicallyOrthogonalMesh
     :members:
     :undoc-members:
diff --git a/docs/api_TensorMesh.rst b/docs/api_TensorMesh.rst
index 9eb346cb..36fb925d 100644
--- a/docs/api_TensorMesh.rst
+++ b/docs/api_TensorMesh.rst
@@ -3,6 +3,6 @@
 Tensor Mesh
 ***********
 
-.. automodule:: SimPEG.TensorMesh
+.. automodule:: SimPEG.mesh.TensorMesh
     :members:
     :undoc-members:
diff --git a/docs/api_TensorView.rst b/docs/api_TensorView.rst
index 6ae80276..c6b5a9b8 100644
--- a/docs/api_TensorView.rst
+++ b/docs/api_TensorView.rst
@@ -3,6 +3,6 @@
 Tensor View
 ***********
 
-.. automodule:: SimPEG.TensorView
+.. automodule:: SimPEG.mesh.TensorView
     :members:
     :undoc-members:
diff --git a/docs/examples/mesh/plot_TensorMesh.py b/docs/examples/mesh/plot_TensorMesh.py
index 94c627f8..e773ed13 100644
--- a/docs/examples/mesh/plot_TensorMesh.py
+++ b/docs/examples/mesh/plot_TensorMesh.py
@@ -1,6 +1,6 @@
 import numpy as np
 import matplotlib.pyplot as plt
-from SimPEG import TensorMesh
+from SimPEG.mesh import TensorMesh
 
 pad = 7
 padfactor = 1.4
diff --git a/docs/examples/mesh/plot_grid_2D.py b/docs/examples/mesh/plot_grid_2D.py
index 258fbdfc..4c6d5da5 100644
--- a/docs/examples/mesh/plot_grid_2D.py
+++ b/docs/examples/mesh/plot_grid_2D.py
@@ -1,6 +1,6 @@
 import numpy as np
 import matplotlib.pyplot as plt
-from SimPEG import TensorMesh
+from SimPEG.mesh import TensorMesh
 
 h1 = np.linspace(.1,.5,3)
 h2 = np.linspace(.1,.5,5)
diff --git a/docs/examples/mesh/plot_grid_3D.py b/docs/examples/mesh/plot_grid_3D.py
index caf8dbd2..6278802c 100644
--- a/docs/examples/mesh/plot_grid_3D.py
+++ b/docs/examples/mesh/plot_grid_3D.py
@@ -1,6 +1,6 @@
 import numpy as np
 import matplotlib.pyplot as plt
-from SimPEG import TensorMesh
+from SimPEG.mesh import TensorMesh
 
 h1 = np.linspace(.1,.5,3)
 h2 = np.linspace(.1,.5,5)
diff --git a/docs/examples/mesh/plot_image_2D.py b/docs/examples/mesh/plot_image_2D.py
index c7bcb0d5..9fef3f1b 100644
--- a/docs/examples/mesh/plot_image_2D.py
+++ b/docs/examples/mesh/plot_image_2D.py
@@ -1,6 +1,6 @@
 import numpy as np
 import matplotlib.pyplot as plt
-from SimPEG import TensorMesh
+from SimPEG.mesh import TensorMesh
 
 n = 20
 h = np.ones(n)/n
diff --git a/docs/examples/mesh/plot_image_3D.py b/docs/examples/mesh/plot_image_3D.py
index d00756a1..c67b0d84 100644
--- a/docs/examples/mesh/plot_image_3D.py
+++ b/docs/examples/mesh/plot_image_3D.py
@@ -1,6 +1,6 @@
 import numpy as np
 import matplotlib.pyplot as plt
-from SimPEG import TensorMesh
+from SimPEG.mesh import TensorMesh
 
 n = 20
 h = np.ones(n)/n
diff --git a/docs/examples/mesh/plot_nodes.py b/docs/examples/mesh/plot_nodes.py
index 5d2c912b..ef11b351 100644
--- a/docs/examples/mesh/plot_nodes.py
+++ b/docs/examples/mesh/plot_nodes.py
@@ -1,6 +1,6 @@
 import numpy as np
 import matplotlib.pyplot as plt
-from SimPEG import TensorMesh
+from SimPEG.mesh import TensorMesh
 
 x0 = np.zeros(2)
 h1 = np.linspace(.1,.5,3)
diff --git a/docs/index.rst b/docs/index.rst
index 4c01c4cd..692363a2 100644
--- a/docs/index.rst
+++ b/docs/index.rst
@@ -26,6 +26,7 @@ Meshing & Operators
    api_TensorMesh
    api_TensorView
    api_LogicallyOrthogonalMesh
+   api_Cyl1DMesh
    api_LOMView
    api_DiffOperators
    api_InnerProducts