In [2]:

    

from SimPEG import Mesh, np, EM
%pylab inline


    

    




Populating the interactive namespace from numpy and matplotlib

In [3]:

    

# problem = DC.Problem2D_CC(mesh)
cs = 2.5
npad = 6
hx = [(cs,npad, -1.3),(cs,160),(cs,npad, 1.3)]
hy = [(cs,npad, -1.3),(cs,20)]
mesh = Mesh.TensorMesh([hx, hy])
mesh = Mesh.TensorMesh([hx, hy],x0=[-mesh.hx[:6].sum()-0.25, -mesh.hy.sum()])


    

In [4]:

    

sigma0 = np.load("sigma0101NBound.npy")
sigma1= np.load("sigma0201NBound.npy")
sigma2 = np.load("sigma0301NBound.npy")
sigma3= np.load("sigma0401NBound.npy")
sigma4 = np.load("sigma0501NBound.npy")
sigma5= np.load("sigma0601NBound.npy")
sigma6 = np.load("sigma0701NBound.npy")
sigma7= np.load("sigma0801NBound.npy")
sigma8 = np.load("sigma0901NBound.npy")
sigma9=  np.load("sigma1001NBound.npy")
sigma10 = np.load("sigma1101NBound.npy")
sigma11= np.load("sigma1201NBound.npy")


    

In [5]:

    

drho2= (1./sigma1-1./sigma0) / (1./sigma0)*100.
drho3= (1./sigma2-1./sigma0) / (1./sigma0)*100.
drho4= (1./sigma3-1./sigma0) / (1./sigma0)*100.
drho5= (1./sigma4-1./sigma0) / (1./sigma0)*100.
drho6= (1./sigma5-1./sigma0) / (1./sigma0)*100.
drho7= (1./sigma6-1./sigma0) / (1./sigma0)*100.
drho8= (1./sigma7-1./sigma0) / (1./sigma0)*100.
drho9= (1./sigma8-1./sigma0) / (1./sigma0)*100.
drho10= (1./sigma9-1./sigma0) / (1./sigma0)*100.
drho11= (1./sigma10-1./sigma0) / (1./sigma0)*100.
drho12= (1./sigma11-1./sigma0) / (1./sigma0)*100.

dperc = (1./sigma1) / (1./sigma0)


    

In [6]:

    

fig, ax = plt.subplots(1,1, figsize = (10, 1.2))
dat = mesh.plotImage(1./sigma0, grid=False, ax=ax, pcolorOpts={"cmap":"jet"}, clim=(0, 200))
ax.set_ylim(-30, 0)
ax.set_xlim(-10, 290)
plt.colorbar(dat[0])


    

    
Out[6]:





<matplotlib.colorbar.Colorbar at 0x10e772290>






    

In [7]:

    

fig, ax = plt.subplots(1,1, figsize = (10, 1.2))
dat = mesh.plotImage(1./sigma1, grid=False, ax=ax, pcolorOpts={"cmap":"jet"}, clim=(0, 200))
ax.set_ylim(-30, 0)
ax.set_xlim(-10, 290)
plt.colorbar(dat[0])


    

    
Out[7]:





<matplotlib.colorbar.Colorbar at 0x10f6b7810>






    

In [8]:

    

fig, ax = plt.subplots(1,1, figsize = (10, 1.2))
dat = mesh.plotImage(drho2/100, grid=False, ax=ax, pcolorOpts={"cmap":"jet"}, clim=(-0.5, 0.5))
ax.set_ylim(-30, 0)
ax.set_xlim(-10, 290)
plt.colorbar(dat[0])


    

    
Out[8]:





<matplotlib.colorbar.Colorbar at 0x10fac7bd0>






    

In [9]:

    

fig, ax = plt.subplots(1,1, figsize = (10, 1.2))
dat = mesh.plotImage(drho2, grid=False, ax=ax, pcolorOpts={"cmap":"jet"}, clim=(-50, 50))
ax.set_ylim(-30, 0)
ax.set_xlim(-10, 290)
plt.colorbar(dat[0])


    

    
Out[9]:





<matplotlib.colorbar.Colorbar at 0x10fce26d0>






    

In [10]:

    

fig, ax = plt.subplots(1,1, figsize = (10, 1.2))
dat= EM.Static.Utils.plot2Ddata(mesh.gridCC, 1./sigma0, ax=ax, level=[80])
ax.set_ylim(-30, 0)
ax.set_xlim(-10, 290)
plt.colorbar(dat[1])


    

    
Out[10]:





<matplotlib.colorbar.Colorbar at 0x110005b50>






    

In [11]:

    

fig, ax = plt.subplots(1,1, figsize = (10, 1.2))
dat= EM.Static.Utils.plot2Ddata(mesh.gridCC, 1./sigma1, ax=ax, level=[80])
ax.set_ylim(-30, 0)
ax.set_xlim(-10, 290)
plt.colorbar(dat[1])


    

    
Out[11]:





<matplotlib.colorbar.Colorbar at 0x110245e50>






    

In [35]:

    

fig, ax = plt.subplots(1,1, figsize = (10, 1.2))
dat= EM.Static.Utils.plot2Ddata(mesh.gridCC, drho2, ax=ax, level=[-1],clim=(-50, 50))
ax.set_ylim(-30, 0)
ax.set_xlim(-10, 290)
plt.colorbar(dat[1])


    

    
Out[35]:





<matplotlib.colorbar.Colorbar at 0x116b20350>






    

In [13]:

    

fig, ax = plt.subplots(1,1, figsize = (10, 1.2))
dat= EM.Static.Utils.plot2Ddata(mesh.gridCC, drho2, ax=ax, level=[10],clim=(-32, 32))
ax.set_ylim(-30, 0)
ax.set_xlim(-10, 290)
plt.colorbar(dat[1])


    

    
Out[13]:





<matplotlib.colorbar.Colorbar at 0x111032910>






    

In [36]:

    

fig, ax = plt.subplots(1,1, figsize = (10, 1.2))
dat= EM.Static.Utils.plot2Ddata(mesh.gridCC, drho3, ax=ax, level=[-1],clim=(-50, 50))
ax.set_ylim(-30, 0)
ax.set_xlim(-10, 290)
plt.colorbar(dat[1])


    

    
Out[36]:





<matplotlib.colorbar.Colorbar at 0x116e0d6d0>






    

In [37]:

    

fig, ax = plt.subplots(1,1, figsize = (10, 1.2))
dat= EM.Static.Utils.plot2Ddata(mesh.gridCC, drho4, ax=ax, level=[-1],clim=(-50, 50))
ax.set_ylim(-30, 0)
ax.set_xlim(-10, 290)
plt.colorbar(dat[1])


    

    
Out[37]:





<matplotlib.colorbar.Colorbar at 0x1178a8cd0>






    

In [38]:

    

fig, ax = plt.subplots(1,1, figsize = (10, 1.2))
dat= EM.Static.Utils.plot2Ddata(mesh.gridCC, drho5, ax=ax, level=[-1],clim=(-50, 50))
ax.set_ylim(-30, 0)
ax.set_xlim(-10, 290)
plt.colorbar(dat[1])


    

    
Out[38]:





<matplotlib.colorbar.Colorbar at 0x117d8a890>






    

In [39]:

    

fig, ax = plt.subplots(1,1, figsize = (10, 1.2))
dat= EM.Static.Utils.plot2Ddata(mesh.gridCC, drho6, ax=ax, level=[-1],clim=(-50, 50))
ax.set_ylim(-30, 0)
ax.set_xlim(-10, 290)
plt.colorbar(dat[1])


    

    
Out[39]:





<matplotlib.colorbar.Colorbar at 0x117fd2b50>






    

In [40]:

    

fig, ax = plt.subplots(1,1, figsize = (10, 1.2))
dat= EM.Static.Utils.plot2Ddata(mesh.gridCC, drho7, ax=ax, level=[-1],clim=(-50, 50))
ax.set_ylim(-30, 0)
ax.set_xlim(-10, 290)
plt.colorbar(dat[1])


    

    
Out[40]:





<matplotlib.colorbar.Colorbar at 0x118b275d0>






    

In [41]:

    

fig, ax = plt.subplots(1,1, figsize = (10, 1.2))
dat= EM.Static.Utils.plot2Ddata(mesh.gridCC, drho8, ax=ax, level=[-1],clim=(-50, 50))
ax.set_ylim(-30, 0)
ax.set_xlim(-10, 290)
plt.colorbar(dat[1])


    

    
Out[41]:





<matplotlib.colorbar.Colorbar at 0x118c947d0>






    

In [42]:

    

fig, ax = plt.subplots(1,1, figsize = (10, 1.2))
dat= EM.Static.Utils.plot2Ddata(mesh.gridCC, drho9, ax=ax, level=[-1],clim=(-50, 50))
ax.set_ylim(-30, 0)
ax.set_xlim(-10, 290)
plt.colorbar(dat[1])


    

    
Out[42]:





<matplotlib.colorbar.Colorbar at 0x1198dec50>






    

In [43]:

    

fig, ax = plt.subplots(1,1, figsize = (10, 1.2))
dat= EM.Static.Utils.plot2Ddata(mesh.gridCC, drho10, ax=ax, level=[-1],clim=(-50, 50))
ax.set_ylim(-30, 0)
ax.set_xlim(-10, 290)
plt.colorbar(dat[1])


    

    
Out[43]:





<matplotlib.colorbar.Colorbar at 0x119c25150>






    

In [44]:

    

fig, ax = plt.subplots(1,1, figsize = (10, 1.2))
dat= EM.Static.Utils.plot2Ddata(mesh.gridCC, drho11, ax=ax, level=[-1],clim=(-50, 50))
ax.set_ylim(-30, 0)
ax.set_xlim(-10, 290)
plt.colorbar(dat[1])


    

    
Out[44]:





<matplotlib.colorbar.Colorbar at 0x119e545d0>






    

In [45]:

    

fig, ax = plt.subplots(1,1, figsize = (10, 1.2))
dat= EM.Static.Utils.plot2Ddata(mesh.gridCC, drho12, ax=ax, level=[-1],clim=(-50, 50))
ax.set_ylim(-30, 0)
ax.set_xlim(-10, 290)
plt.colorbar(dat[1])


    

    
Out[45]:





<matplotlib.colorbar.Colorbar at 0x11a1963d0>






    

In [ ]: