In [ ]:

    

# latest version at ../0.4


    

In [ ]:

    

from pylab import *
import copy


    

In [ ]:

    

x,y,z = genfromtxt('f1.dat').T
N = int(sqrt(x.shape[0]))


    

In [ ]:

    

x = x.reshape(N, N)
y = y.reshape(N, N)
z = z.reshape(N, N)
mark1 = copy.deepcopy(z)
mark2 = copy.deepcopy(z)
cell1 = copy.deepcopy(z)


    

In [ ]:

    

for i in range (0, N-1):
    z_point = z[i] * z[i+1]
    for j in range (0, N):
        if (z_point[j] < 0) :
            mark1[i][j] = 1.0
        else:
            mark1[i][j] = 0.0


    

In [ ]:

    

z = z.T
x = x.T
y = y.T


    

In [ ]:

    

for i in range (0, N-1):
    z_point = z[i] * z[i+1]
    for j in range (0, N):
        if (z_point[j] < 0) :
            mark2[j][i] = 1.0
        else:
            mark2[j][i] = 0.0


    

In [ ]:

    

print mark1[0][226] + mark1[0][227] + mark2[0][226] + mark2[1][226]


    

In [ ]:

    

for i in range (0, N-1):
    for j in range (0, N-):
        cell1[i][j] = mark1[i][j] + mark1[i][j+1] + mark2[i][j] + mark2[i+1][j]


    

In [ ]:

    

x,y,z = genfromtxt('f2.dat').T
N = int(sqrt(x.shape[0]))


    

In [ ]:

    

x = x.reshape(N, N)
y = y.reshape(N, N)
z = z.reshape(N, N)
mark1 = copy.deepcopy(z)
mark2 = copy.deepcopy(z)
cell2 = copy.deepcopy(z)


    

In [ ]:

    

for i in range (0, N-1):
    z_point = z[i] * z[i+1]
    for j in range (0, N-1):
        if (z_point[j] < 0) :
            mark2[j][i] = 1.0
        else:
            mark2[j][i] = 0.0


    

In [ ]:

    

z = z.T
x = x.T
y = y.T


    

In [ ]:

    

for i in range (0, N-1):
    z_point = z[i] * z[i+1]
    for j in range (0, N):
        if (z_point[j] < 0) :
            mark2[j][i] = 1.0
        else:
            mark2[j][i] = 0.0


    

In [ ]:

    

for i in range (0, N-1):
    for j in range (0, N-1):
        cell2[i][j] = mark1[i][j] + mark1[i][j+1] + mark2[i][j] + mark2[i+1][j]


    

In [ ]:

    

cell = cell1*cell2


    

In [ ]:

    

for i in range (0, N-1):
    for j in range (0, N-1):
        if (cell[i][j] != 0):
            print i, j