Resistivity Transform Graph
4 layered earth model

All calculated with python (Matplotlib, numpy)

Back Again !!!

After a long time I am back again. Due to master stuff at University I was away for a long time period but now I feel free so that I decided to add a few things. Like signal contribution section used in Geoelectrics 

which was one of my homeworks. Below figure is computed for dipole - dipole geometry.

I give the code below for the ones who want to try and see alive.

'''

Created on Mar 7, 2011

@author: Burak Nebioglu

'''

from pylab import *

from scipy import mgrid

import math as m

z=0

Sigma=10000000000000

I=5000000000

a=0.4

b=0.4

c=0.4

def imsave(filename, X, **kwargs):

    """ Homebrewed imsave to have nice colors... """

    figsize=(array(X.shape)/10.0)[::-1]

    rcParams.update({'figure.figsize':figsize})

    fig = figure(figsize=figsize,dpi=200)

    axes([0,0,1,1]) # Make the plot occupy the whole canvas

    axis('off')

    fig.set_size_inches(figsize)

    colorbar(imshow(X,origin='lower', **kwargs))

    savefig(filename, facecolor='white', edgecolor='white', dpi=200)

    close(fig)

x,y=mgrid[-2.5:2.5:0.1,-1:0:0.1]

denk1=(x*(x-a)+y**2+z**2)/((x**2+y**2+z**2)**(3/2.)*((x-a)**2+y**2+z**2)**(3/2.))

print denk1

denk2=((x-a)*(x-a-b-c)+y**2+z**2)/(((x-a)**2+y**2+z**2)**(3/2.)*((x-a-b-c)**2+y**2+z**2)**(3/2.))

print denk2

denk3=(x*(x-a-b)+y**2+z**2)/((x**2+y**2+z**2)**(3/2.)*((x-a-b)**2+y**2+z**2)**(3/2.))

print denk3

denk4=((x-a-b)*(x-a-b-c)+y**2+z**2)/(((x-a-b)**2+y**2+z**2)**(3/2.)*((x-a-b-c)**2+y**2+z**2)**(3/2.))

print denk4

'''

za=((Sigma*I)/(2*m.pi))*((x*(x-a)+y**2+z**2)/((x**2+y**2+z**2)**3/2.*((x-a)**2+y**2+z**2)**3/2.)-(((x-a)*(x-a-b-c)+y**2+z**2)/(((x-a)**2+y**2+z**2)**3/2.*((x-a-b-c)**2+y**2+z**2)**3/2.)))

'''

Z=((Sigma*I)/(2*m.pi))*(denk1-denk2-denk3+denk4) # Create the data to be plotted

imsave('last.png', transpose(Z), cmap=cm.jet )