from em_examples import sphereElectrostatic_example as electrostatic_sphere

sig0 = 10.**-3.          # conductivity of the wholespace
sig1 = 10.**-1.         # conductivity of the conductive sphere
sig2 = 10.**-5.         # conductivity of the resistive sphere
R    = 50.          # radius of the sphere
E0   = 1.           # inducing field strength
n = 50             #level of discretisation
xr = np.linspace(-2.*R, 2.*R, n) # X-axis discretization
yr = xr.copy()      # Y-axis discretization
zr = np.r_[0]          # identical to saying `zr = np.array([0])`
XYZ = electrostatic_sphere.ndgrid(xr,yr,zr) # Space Definition
PlotOpt = 'Total'

xstart=-100.
ystart=-100.
xend=100.
yend=100.
nb_dipole=31
electrode_spacing=10

electrostatic_sphere.two_configurations_comparison(XYZ,sig0,sig1,sig2,R,R,E0,xstart,ystart,xend,yend,nb_dipole,electrode_spacing,PlotOpt)

plt.tight_layout()
plt.show()