diff --git a/simpegEM/Utils/Ana/FEM.py b/simpegEM/Analytics/FEM.py
similarity index 100%
rename from simpegEM/Utils/Ana/FEM.py
rename to simpegEM/Analytics/FEM.py
diff --git a/simpegEM/Utils/Ana/TEM.py b/simpegEM/Analytics/TEM.py
similarity index 100%
rename from simpegEM/Utils/Ana/TEM.py
rename to simpegEM/Analytics/TEM.py
diff --git a/simpegEM/Utils/Ana/__init__.py b/simpegEM/Analytics/__init__.py
similarity index 100%
rename from simpegEM/Utils/Ana/__init__.py
rename to simpegEM/Analytics/__init__.py
diff --git a/simpegEM/FDEM/FDEM.py b/simpegEM/FDEM/FDEM.py
index 3381a6ab..d137060c 100644
--- a/simpegEM/FDEM/FDEM.py
+++ b/simpegEM/FDEM/FDEM.py
@@ -1,7 +1,7 @@
 from SimPEG import Survey, Problem, Utils, np, sp, Solver as SimpegSolver
 from scipy.constants import mu_0
 from SurveyFDEM import SurveyFDEM, FieldsFDEM
-from simpegEM.Utils import Sources
+from simpegEM import Sources
 from simpegEM.Base import BaseEMProblem
 
 def omega(freq):
diff --git a/simpegEM/Utils/Sources/CircularLoop.py b/simpegEM/Sources/CircularLoop.py
similarity index 100%
rename from simpegEM/Utils/Sources/CircularLoop.py
rename to simpegEM/Sources/CircularLoop.py
diff --git a/simpegEM/Utils/Sources/__init__.py b/simpegEM/Sources/__init__.py
similarity index 100%
rename from simpegEM/Utils/Sources/__init__.py
rename to simpegEM/Sources/__init__.py
diff --git a/simpegEM/Utils/Sources/magneticDipole.py b/simpegEM/Sources/magneticDipole.py
similarity index 100%
rename from simpegEM/Utils/Sources/magneticDipole.py
rename to simpegEM/Sources/magneticDipole.py
diff --git a/simpegEM/TDEM/BaseTDEM.py b/simpegEM/TDEM/BaseTDEM.py
index d12e776e..7626d944 100644
--- a/simpegEM/TDEM/BaseTDEM.py
+++ b/simpegEM/TDEM/BaseTDEM.py
@@ -1,6 +1,6 @@
 from SimPEG import Solver, Problem
 from SimPEG.Problem import BaseTimeProblem
-from simpegEM.Utils import Sources
+from simpegEM import Sources
 from scipy.constants import mu_0
 from SimPEG.Utils import sdiag, mkvc
 from SimPEG import Utils, Mesh
diff --git a/simpegEM/TDEM/SurveyTDEM.py b/simpegEM/TDEM/SurveyTDEM.py
index 629658de..2b0b2ce5 100644
--- a/simpegEM/TDEM/SurveyTDEM.py
+++ b/simpegEM/TDEM/SurveyTDEM.py
@@ -1,6 +1,6 @@
 from SimPEG import Utils, Survey, np
 from SimPEG.Survey import BaseSurvey
-from simpegEM.Utils import Sources
+from simpegEM import Sources
 from BaseTDEM import FieldsTDEM
 
 
@@ -106,7 +106,7 @@ class TxTDEM(Survey.BaseTx):
         else:
             raise Exception('Unknown mesh for CircularLoop')
 
-        return {"b": mesh.edgeCurl*MVP}        
+        return {"b": mesh.edgeCurl*MVP}
 
     def getJs(self, mesh, time):
         return None
diff --git a/simpegEM/Tests/test_FDEM_analytics.py b/simpegEM/Tests/test_FDEM_analytics.py
index 2928c50c..732b51ce 100644
--- a/simpegEM/Tests/test_FDEM_analytics.py
+++ b/simpegEM/Tests/test_FDEM_analytics.py
@@ -58,7 +58,7 @@ class FDEM_analyticTests(unittest.TestCase):
 
         P = self.mesh.getInterpolationMat(XYZ, 'Fz')
 
-        an = EM.Utils.Ana.FEM.hzAnalyticDipoleF(x, self.Tx0.freq, self.sig)
+        an = EM.Analytics.FEM.hzAnalyticDipoleF(x, self.Tx0.freq, self.sig)
 
         diff = np.log10(np.abs(P*np.imag(u[self.Tx0, 'b']) - mu_0*np.imag(an)))
 
diff --git a/simpegEM/Tests/test_TDEM_forward_Analytic.py b/simpegEM/Tests/test_TDEM_forward_Analytic.py
index 8400dcf2..4557c5cf 100644
--- a/simpegEM/Tests/test_TDEM_forward_Analytic.py
+++ b/simpegEM/Tests/test_TDEM_forward_Analytic.py
@@ -41,7 +41,7 @@ def halfSpaceProblemAnaDiff(meshType, sig_half=1e-2, rxOffset=50., bounds=[1e-5,
     sigma = np.log(sigma[active])
     prb.pair(survey)
 
-    bz_ana = mu_0*EM.Utils.Ana.hzAnalyticDipoleT(rx.locs[0][0]+1e-3, rx.times, sig_half)
+    bz_ana = mu_0*EM.Analytics.hzAnalyticDipoleT(rx.locs[0][0]+1e-3, rx.times, sig_half)
 
     bz_calc = survey.dpred(sigma)
 
diff --git a/simpegEM/Utils/__init__.py b/simpegEM/Utils/__init__.py
deleted file mode 100644
index c9e5e23b..00000000
--- a/simpegEM/Utils/__init__.py
+++ /dev/null
@@ -1,3 +0,0 @@
-import Sources
-import Ana
-# import Solver
\ No newline at end of file
diff --git a/simpegEM/__init__.py b/simpegEM/__init__.py
index 92c6ae92..90a1dcf7 100644
--- a/simpegEM/__init__.py
+++ b/simpegEM/__init__.py
@@ -1,5 +1,6 @@
 # from EM import *
-import Utils
 import TDEM
 import FDEM
-import Base
\ No newline at end of file
+import Base
+import Sources
+import Analytics