# -*- coding: utf-8 -*-
"""
Created on Sat Sep 12 14:25:03 2020

@author: Fabrice
"""
import numpy as np
import matplotlib.pyplot as plt

##perméabilité magnétique
mu0 = 4*np.pi*1e-7

#intensité du courant en ampères
I=100

#grille cartésienne (attention, très près du fil, car B diminue fortement !)
X = np.linspace(-5e-3 , 5e-3 , 1000)
Y = X.copy()
X , Y = np.meshgrid(X,Y)

#intensité du champ B
B = mu0*I/(2*np.pi*np.sqrt(X**2+Y**2))

#graphique
fig = plt.figure (figsize = (8 , 8))
ax = fig.add_subplot(111)

#contours=lignes de champ en tesla
levels = (0.005, 0.01 , 0.02 , 0.03 , 0.04)

#dessine les contours
cp = ax.contour(X , Y, B, levels = levels, colors = 'k' , linewidths = 2)
ax.clabel(cp , fontsize = 10 , colors = 'b')

#dessine la position du fil au centre
ax.plot(0 , 0, marker = '*')

#derniers paramètres du graphique et dessin final
ax.axis('equal')
plt.title("Lignes de champ magnétique, intensité en tesla, d'un fil rectiligne infini, I = " +str(I)+' A', fontsize=10)
ax.set_xlabel('x(m)')
ax.set_ylabel('y(m)')
plt.grid()
plt.show()