import matplotlib.pyplot as plt
import numpy as np

# Datos
parametros = ['Longitud de\nla tubería', 'Ancho de\nla cara', 'Diámetro de\nla tubería', 'Flujo másico\ndel HTF']
ef_inicial = [0.1804, 0.1946, 0.1813, 0.1607]
ef_final = [0.1786, 0.1795, 0.1585, 0.18]
diferencias = [0.0018, 0.0151, 0.0228, 0.0193]

diferencias_porcentuales = [(dif / ini) * 100 for dif, ini in zip(diferencias, ef_inicial)]

# Gráfico de barras para comparar eficiencias iniciales y finales
x = np.arange(len(parametros))  # Posiciones de las barras
width = 0.35  # Ancho de las barras

fig, ax = plt.subplots(1, 2, figsize=(20, 8))  # 1 fila, 2 columnas

# Gráfico 1: Eficiencias iniciales y finales
ax[0].bar(x - width/2, ef_inicial, width, label='Eficiencia inicial', color='blue', alpha=0.7)
ax[0].bar(x + width/2, ef_final, width, label='Eficiencia final', color='green', alpha=0.7)
ax[0].set_xticks(x)
ax[0].set_xticklabels(parametros, fontsize=20)
ax[0].set_ylabel('Eficiencia eléctrica', fontsize=20)
ax[0].set_title('', fontsize=20)
ax[0].legend(loc='upper left', fontsize=20)
ax[0].grid(axis='y', linestyle='--', alpha=0.7)
ax[0].tick_params(axis='y', labelsize=20)

# Remarcar los bordes del gráfico 1
for spine in ax[0].spines.values():
    spine.set_edgecolor('black')
    spine.set_linewidth(2)

# Ajustar límites del eje Y para que no cubra la gráfica
ax[0].set_ylim(0, max(max(ef_inicial), max(ef_final)) * 1.2)

# Gráfico 2: Diferencias porcentuales
ax[1].bar(x, diferencias_porcentuales, width, color='orange', alpha=0.7)  # Usar el mismo ancho de barras
ax[1].set_xticks(x)
ax[1].set_xticklabels(parametros, fontsize=20)
ax[1].set_ylabel('\nVariación de la eficiencia (%)', fontsize=20)
ax[1].set_title('', fontsize=20)
ax[1].grid(axis='y', linestyle='--', alpha=0.7)
ax[1].tick_params(axis='y', labelsize=20)
ax[1].tick_params(axis='x', labelsize=20)

# Remarcar los bordes del gráfico 2
for spine in ax[1].spines.values():
    spine.set_edgecolor('black')
    spine.set_linewidth(2)

# Ajustar diseño
plt.tight_layout()
plt.show()