In [121]:

    

from SimPEG import Mesh, np, EM
%pylab inline


    

    




Populating the interactive namespace from numpy and matplotlib

In [122]:

    

# 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 [123]:

    

sigma0 = np.load("sigma0101NBound.npy")
sigma1= np.load("sigma0308NBound.npy")
sigma2 = np.load("sigma0915NBound.npy")
sigma3= np.load("sigma1215NBound.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 [124]:

    

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 [125]:

    

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


    

    
Out[125]:





<matplotlib.colorbar.Colorbar at 0x12a1943d0>






    

In [126]:

    

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


    

    
Out[126]:





<matplotlib.colorbar.Colorbar at 0x12a5d3d90>






    

In [127]:

    

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


    

    
Out[127]:





<matplotlib.colorbar.Colorbar at 0x12a904f90>






    

In [128]:

    

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


    

    
Out[128]:





<matplotlib.colorbar.Colorbar at 0x12ab271d0>






    

In [129]:

    

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[129]:





<matplotlib.colorbar.Colorbar at 0x12ae683d0>






    

In [130]:

    

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[130]:





<matplotlib.colorbar.Colorbar at 0x12b180550>






    

In [131]:

    

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


    

    
Out[131]:





<matplotlib.colorbar.Colorbar at 0x12b3716d0>






    

In [132]:

    

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


    

    
Out[132]:





<matplotlib.colorbar.Colorbar at 0x12b6e4750>






    

In [154]:

    

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


    

    
Out[154]:





(-10, 290)






    

In [153]:

    

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


    

    
Out[153]:





(-10, 290)






    

In [152]:

    

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


    

    
Out[152]:





(-10, 290)






    

In [149]:

    

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[149]:





(-10, 290)






    

In [151]:

    

plt.colorbar(dat[1])


    

    
Out[151]:





<matplotlib.colorbar.Colorbar at 0x130ae3cd0>






    

In [172]:

    

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[172]:





(-10, 290)






    

In [171]:

    

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


    

    
Out[171]:





(-10, 290)






    

In [173]:

    

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


    

    
Out[173]:





(-10, 290)






    

In [168]:

    

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


    

    
Out[168]:





(-10, 290)






    

In [169]:

    

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


    

    
Out[169]:





(-10, 290)






    

In [170]:

    

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


    

    
Out[170]:





(-10, 290)






    

In [139]:

    

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


    

    
Out[139]:





<matplotlib.colorbar.Colorbar at 0x12ecdb7d0>






    

In [140]:

    

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


    

    
Out[140]:





<matplotlib.colorbar.Colorbar at 0x12f1326d0>






    

In [141]:

    

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


    

    
Out[141]:





<matplotlib.colorbar.Colorbar at 0x12f215e10>






    

In [142]:

    

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


    

    
Out[142]:





<matplotlib.colorbar.Colorbar at 0x12f4e8f90>






    

In [143]:

    

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


    

    
Out[143]:





<matplotlib.colorbar.Colorbar at 0x12f83bf50>






    

In [144]:

    

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


    

    
Out[144]:





<matplotlib.colorbar.Colorbar at 0x12fba10d0>






    

In [145]:

    

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


    

    
Out[145]:





<matplotlib.colorbar.Colorbar at 0x12fddea50>






    

In [146]:

    

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


    

    
Out[146]:





<matplotlib.colorbar.Colorbar at 0x12ff35b10>






    

In [147]:

    

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


    

    
Out[147]:





<matplotlib.colorbar.Colorbar at 0x130477c90>






    

In [148]:

    

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


    

    
Out[148]:





<matplotlib.colorbar.Colorbar at 0x1307b0250>






    

In [ ]: