In [340]:

    

from SimPEG import Mesh, np, EM
%pylab inline


    

    




Populating the interactive namespace from numpy and matplotlib

In [341]:

    

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

    

sigma0 = np.load("sigma0101Oct.npy")
sigma1=  np.load("sigma0721Oct.npy")
sigma2 = np.load("sigma0722Oct.npy")
sigma3=  np.load("sigma0723Oct.npy")
sigma4 = np.load("sigma0724Oct.npy")
sigma5=  np.load("sigma0725Oct.npy")
sigma6 = np.load("sigma0726Oct.npy")
sigma7=  np.load("sigma0727Oct.npy")
sigma8 = np.load("sigma0728Oct.npy")
sigma9 = np.load("sigma0810Oct.npy")
sigma10= np.load("sigma0820Oct.npy")
sigma11= np.load("sigma0830Oct.npy")
sigma12= np.load("sigma0915Oct.npy")
sigma13= np.load("sigma1030Oct.npy")
sigma14= np.load("sigma1130Oct.npy")
sigma15= np.load("sigma1218Oct.npy")
sigma16= np.load("sigma1228Oct.npy")


    

In [343]:

    

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

    

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





<matplotlib.colorbar.Colorbar at 0x15974d3d0>






    

In [345]:

    

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





<matplotlib.colorbar.Colorbar at 0x159a82250>






    

In [346]:

    

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





<matplotlib.colorbar.Colorbar at 0x159c8bd10>






    

In [347]:

    

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





<matplotlib.colorbar.Colorbar at 0x159ea2e90>






    

In [348]:

    

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





<matplotlib.colorbar.Colorbar at 0x15a053f90>






    

In [349]:

    

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





<matplotlib.colorbar.Colorbar at 0x15a49bed0>






    

In [350]:

    

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])
#print drho2


    

    
Out[350]:





<matplotlib.colorbar.Colorbar at 0x15b809b10>






    

In [351]:

    

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





<matplotlib.colorbar.Colorbar at 0x15ba37690>






    

In [352]:

    

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





<matplotlib.colorbar.Colorbar at 0x15bcee590>






    

In [353]:

    

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





<matplotlib.colorbar.Colorbar at 0x15ca0ae10>






    

In [354]:

    

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





<matplotlib.colorbar.Colorbar at 0x15cc36e50>






    

In [355]:

    

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





<matplotlib.colorbar.Colorbar at 0x15d07a690>






    

In [356]:

    

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





<matplotlib.colorbar.Colorbar at 0x15d2ab850>






    

In [357]:

    

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





<matplotlib.colorbar.Colorbar at 0x15d4e57d0>






    

In [358]:

    

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





<matplotlib.colorbar.Colorbar at 0x15e037410>






    

In [359]:

    

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





<matplotlib.colorbar.Colorbar at 0x15e49e690>






    

In [360]:

    

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





<matplotlib.colorbar.Colorbar at 0x15e6dc690>






    

In [361]:

    

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


    

    
Out[361]:





<matplotlib.colorbar.Colorbar at 0x15ea3ec90>






    

In [362]:

    

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


    

    
Out[362]:





<matplotlib.colorbar.Colorbar at 0x15ed9cfd0>






    

In [363]:

    

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


    

    
Out[363]:





<matplotlib.colorbar.Colorbar at 0x15f0d6b10>






    

In [364]:

    

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





<matplotlib.colorbar.Colorbar at 0x15f304c50>






    

In [365]:

    

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





<matplotlib.colorbar.Colorbar at 0x15f5412d0>






    

In [ ]: