In [1]:

    

from SimPEG import Mesh, np, EM
%pylab inline


    

    




Populating the interactive namespace from numpy and matplotlib

In [2]:

    

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

    

sigma0 = np.load("sigma0101NBound.npy")
sigma1=  np.load("sigma0721NBound.npy")
sigma2 = np.load("sigma0722NBound.npy")
sigma3=  np.load("sigma0723NBound.npy")
sigma4 = np.load("sigma0724NBound.npy")
sigma5=  np.load("sigma0725NBound.npy")
sigma6 = np.load("sigma0726NBound.npy")
sigma7=  np.load("sigma0727NBound.npy")
sigma8 = np.load("sigma0730NBound.npy")
sigma9 = np.load("sigma0801NBound.npy")
sigma10= np.load("sigma0805NBound.npy")
sigma11= np.load("sigma0810NBound.npy")
sigma12= np.load("sigma0815NBound.npy")
sigma13= np.load("sigma0820NBound.npy")
sigma14= np.load("sigma0825NBound.npy")
sigma15= np.load("sigma0901NBound.npy")
sigma16= np.load("sigma0915NBound.npy")


    

In [4]:

    

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.
drho13= (1./sigma12-1./sigma0) / (1./sigma0)*100.
drho14= (1./sigma13-1./sigma0) / (1./sigma0)*100.
drho15= (1./sigma14-1./sigma0) / (1./sigma0)*100.
drho16= (1./sigma15-1./sigma0) / (1./sigma0)*100.
drho17= (1./sigma16-1./sigma0) / (1./sigma0)*100.

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


    

In [5]:

    

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





<matplotlib.colorbar.Colorbar at 0x10c51c390>






    

In [6]:

    

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





<matplotlib.colorbar.Colorbar at 0x10c7cf790>






    

In [7]:

    

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





<matplotlib.colorbar.Colorbar at 0x10f9c0b50>






    

In [8]:

    

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





<matplotlib.colorbar.Colorbar at 0x10fbdb550>






    

In [9]:

    

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





<matplotlib.colorbar.Colorbar at 0x10fff48d0>






    

In [10]:

    

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





<matplotlib.colorbar.Colorbar at 0x110490690>






    

In [11]:

    

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


    

    
Out[11]:





<matplotlib.colorbar.Colorbar at 0x1106cd790>






    

In [12]:

    

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


    

    
Out[12]:





<matplotlib.colorbar.Colorbar at 0x1111593d0>






    

In [13]:

    

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


    

    
Out[13]:





<matplotlib.colorbar.Colorbar at 0x111526250>






    

In [14]:

    

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


    

    
Out[14]:





<matplotlib.colorbar.Colorbar at 0x11209b710>






    

In [15]:

    

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


    

    
Out[15]:





<matplotlib.colorbar.Colorbar at 0x1123d31d0>






    

In [16]:

    

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


    

    
Out[16]:





<matplotlib.colorbar.Colorbar at 0x11260bd10>






    

In [17]:

    

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


    

    
Out[17]:





<matplotlib.colorbar.Colorbar at 0x112838e50>






    

In [18]:

    

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


    

    
Out[18]:





<matplotlib.colorbar.Colorbar at 0x1129b8b90>






    

In [19]:

    

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


    

    
Out[19]:





<matplotlib.colorbar.Colorbar at 0x112ed7b50>






    

In [20]:

    

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


    

    
Out[20]:





<matplotlib.colorbar.Colorbar at 0x113a1c790>






    

In [21]:

    

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


    

    
Out[21]:





<matplotlib.colorbar.Colorbar at 0x113c4ec50>






    

In [22]:

    

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


    

    
Out[22]:





<matplotlib.colorbar.Colorbar at 0x113fc56d0>






    

In [23]:

    

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


    

    
Out[23]:





<matplotlib.colorbar.Colorbar at 0x114972190>






    

In [24]:

    

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


    

    
Out[24]:





<matplotlib.colorbar.Colorbar at 0x114e2e750>






    

In [25]:

    

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


    

    
Out[25]:





<matplotlib.colorbar.Colorbar at 0x11519e3d0>






    

In [26]:

    

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


    

    
Out[26]:





<matplotlib.colorbar.Colorbar at 0x1154db350>






    

In [ ]: