# -*- coding: utf-8 -*-
"""
Editor de Spyder

Este es un archivo temporal
"""
#!/usr/bin/env python

"""
Cálculo del campo magnético generado por dos bobinas sobre el eje de simetría.
"""

import numpy as np
from scipy import *
#%% Parámetros

R = 4.9                 ;"radio de las bobinas (en cm)" 
d = 3.0                 ;"Separación entre las bobinas (en cm)"
L1 = 11.3               ;"Longitud de la bobina 1 (en cm)"
n1 = 1000.0             ;"Espiras por unidad de longitud para la bobina 1 (por m)"
L2 = 11.3               ;"Longitud de la bobina 2 (en cm)"
n2 = 1000.0             ;"Espiras por unidad de longitud para la bobina 2 (por m)"
#mu0 = 4*3.14*10**(-7)   ;"mu0 en Tesla*m/A"  
mu0 = 4*3.14*10**(-3)   ;"mu0 en Gauss*m/A"

#%% Función que calcula el campo magnético (ver Anexo-2B.pdf)
 
def B(I,z):
    "Campo B sobre la sonda Hall (ubicada en la posición z) en función de I (corriente) "
    # Convierto a I en un array
    I=conj(conj(I))
    
    # Campo de la bobina 1
    B1=0.5*I*n1*mu0*((L1 - z)/np.sqrt((z - L1)**2 + R**2) + z/np.sqrt((z)**2 + R**2))
   
    # Campo de la bobina 2
    B2=0.5*I*n2*mu0*((L1 + L2 + d - z)/np.sqrt((z - L1 - L2 - d)**2 + R**2) - (L1 + d - z)/np.sqrt((z - L1 - d)**2 + R**2))
    
    # Campo total
    return B1 + B2

#%% Mis mediciones:

# Defino el vector corriente (unidades: A)
I=(0.32,0.3,0.25,0.2,0.15,0.1,0.05,0.01)

# Posición en el eje z  (unidades: cm)
z=12

#%% Cálculo la intensidad del campo para cada corriente en Gauss
print("campo B (Gauss):")

print(B(I,z))