xfmt = md.DateFormatter("%H:%M:%S")
plt.figure(1)
imshow_UTC(year, month, day, tiempo_central, cociente, escalas, "inferno", 3)
# plot_datetime(year, month, day,t_MVA, B_MVA, 'cyan', '-', 1, 0.5) #no sé por qué se superpone mal, tiene mal los tiempos.
for tt in timestamps:
plt.axvline(x=tt, color="g") # plotea los tiempos t1t2t3t4
plt.title(
f"Heatmap del cociente de lambdas en distintas escalas temporales \n para el día {day}-{month}-{year}"
)
plt.xlabel("Tiempo en el que está centrado (hh:mm:ss)")
plt.ylabel("Radio (s) \n |B| (nT)")
plt.figure(2)
plot_datetime(year, month, day, t_cut, B_cut, "red", "-", 1, 1)
plt.ylabel("|B| (nT)")
plt.xlabel("Tiempo UTC (hh:mm:ss)")
fig = plt.figure(
3, constrained_layout=True
) # Lo bueno de esta forma es que puedo hacer que solo algunos compartan eje
fig.subplots_adjust(top=0.93,
bottom=0.07,
left=0.05,
right=0.95,
hspace=0.005,
wspace=0.15)
fig.set_size_inches(15,
10) # con este tamaño ocupa toda la pantalla de la laptop
import numpy as np
import matplotlib.pyplot as plt
from funciones import fechas
from funciones_plot import plot_datetime
from importar_datos import importar_mag_1s
"""
Plotea sólo los datos de MAG de baja resolución
"""
np.set_printoptions(precision=4)
year, month, day, doy = fechas()
mag, t, B, posicion = importar_mag_1s(year, month, day, 0.1, 24)
plt.figure()
plot_datetime(year, month, day, t, np.linalg.norm(B, axis=1))
plt.xlabel("t (UTC)")
plt.ylabel("|B|")
plt.ylim(ymin=0, ymax=70)
plt.title("MAG lowres hdec")
plt.show(block=False)
inicio = donde(t, ti)
fin = donde(t, tf)
t1, t2, t3, t4 = importar_t1t2t3t4()
Bu = np.zeros((180, 4))
for i in range(180):
paso = t1 - 180 + 0.0001 * i
inicio_up = donde(t, paso)
fin_up = donde(t, t1)
B_upstream = np.mean(B[inicio_up:fin_up, :], axis=0) # nT
Bu[i, 0:3] = B_upstream
Bu[i, 3] = np.linalg.norm(B_upstream)
plt.figure(1)
plot_datetime(year, month, day, t[fin_up - 180 : fin_up], Bu[:, 0])
plot_datetime(year, month, day, t[fin_up - 180 : fin_up], Bu[:, 1], colour="C1")
plot_datetime(year, month, day, t[fin_up - 180 : fin_up], Bu[:, 2], colour="C2")
plt.xlabel("t (UTC)")
plt.ylabel("componentes Bup")
plt.figure(2)
plot_datetime(year, month, day, t[fin_up - 180 : fin_up], Bu[:, 3])
plt.xlabel("t (UTC)")
plt.ylabel("|B_up|")
Bd = np.zeros((200, 4))
for i in range(200):
paso = t4 + 0.0001 * i
inicio_down = donde(t, t4)
year, month, day, doy = fechas()
ti, tf = tiempos()
mag, t, B, posicion = importar_mag(year, month, day, ti, tf)
inicio = donde(t, ti)
fin = donde(t, tf)
if any(posicion[inicio:fin, 2]) > 0:
print("Norte")
B_norm = np.linalg.norm(B, axis=1)
B_cut = B_norm[inicio:fin]
t_cut = t[inicio:fin]
plt.figure()
plt.plot(t_cut, B_cut)
plt.xlabel("t (hdec)")
plt.ylabel("|B|")
plt.title("MAG hires hdec")
plt.figure()
plot_datetime(year, month, day, t_cut, B_cut)
plt.xlabel("t (UTC)")
plt.ylabel("|B|")
plt.title("MAG hires UTC")
plt.show(block=False)
n = int(ti * 32 * 3600)
mag, t, B, posicion = importar_mag(year, month, day, ti, tf)
inicio = donde(t, ti)
fin = donde(t, tf)
B_norm = np.linalg.norm(B, axis=1)
inicio_MVA = donde(t, tiempo_central[0])
fin_MVA = donde(t, tiempo_central[-1])
B_MVA = B_norm[inicio_MVA:fin_MVA]
t_MVA = t[inicio_MVA:fin_MVA]
plt.figure()
imshow_UTC(year, month, day, tiempo_central, cociente, escalas, "inferno")
# plot_datetime(year, month, day,t_MVA, B_MVA, 'cyan', '-', 1, 0.5) #no sé por qué se superpone mal, tiene mal los tiempos.
for tt in timestamps[1:3]:
plt.axvline(x=tt, color="g") # plotea los tiempos t2 y t3
plt.title("Heatmap del cociente de lambdas en distintas escalas temporales")
plt.xlabel("Tiempo en el que está centrado (hh:mm:ss)")
plt.ylabel("Radio (s) \n |B| (nT)")
plt.figure()
plot_datetime(year, month, day, t, B_norm, "red", "-", 1, 1)
plt.ylabel("|B| (nT)")
plt.xlabel("Tiempo UTC (hh:mm:ss)")
"""
pensar criterio para la escala maxima
"""