program main

implicit none

    ! -------- Valores de entrada -------- "
    Real :: G_tot = 578                          !Radiación total incidente"
    Real :: G_b = 374                            !Radiación directa"
    
    Real :: Tem_a = 282                          !25 + 273.15, Temperatura ambiente"
    Real :: Tem_sky = 282                        !25 + 273.15, Temperatura del cielo"
    Real :: Tem_f_in = 282                       !70 + 273.15, Temperatura de entrada del fluido"
    Real :: Tem_pv = 355.45                      !83.75 + 273.15 Temperatura de entrada del fluido"
    
    ! -------- Geometría del PTC y SRC-PVT -------- "
    Real :: D = 0.06        !Mod"                !Diámetro de la tubería por donde circula el HTF"
    Real :: A_ap = 60       !Mod"                !Área de apertura del concentrador"
    Real :: A_pv = 6        !Mod"                !Área del PV"
    Real :: A_abs = 3       !Mod"                !Área del Absorbedor"
    Real :: A_conc = 94     !Mod"                !Área del concentrador"
    Real :: CR_pv = 10      !Mod"                !Relación de concentración del concentrador"
    
    ! -------- Propiedades físicas del SRC-PVT y Concentrador -------- "
    Real :: alpha_pv = 0.97                      !Admitancia del PV"
    Real :: alpha_abs = 0.9                      !Admitancia del Absorbedor"
    Real :: alpha_conc =0.03                     !Admitancia del Concentrador"
    Real :: eta_opt = 0.83                       !Eficiencia óptica"
    !"eta_pv = 0.2998"                    "Eficiencia del PV"
    Real :: IAM_th = 1                           !Ángulo de incicencia térmica"
    Real :: IAM_elec = 1                         !Ángulo de incicencia eléctrica"
    
    Real :: epsilon_pv = 0.2                     !Emisividad del PV"
    Real :: epsilon_abs = 0.2                    !Emisividad del Absorbedor"
    Real :: epsilon_conc = 0.3                   !Emisividad del Concentrador"
    Real :: sigma = 5.67E-8                 !Constante de Stefan-Boltzmann"
    Real :: P = 1.01325                          !Presión atmosferica en bar"
    
    ! -------- Propiedades del HTF -------- "
    Real :: m_dot_f = 0.1   !Mod"                !Flujo másico del agua"
    Real :: C_p_f = 4190
    !"C_p_f = cp(Water,T=Tem_f_in,P=P) "    "Calor específico del agua a 70 °C - C_p_f = 4190"
    
    ! -------- Coeficientes y Parámetros de las resistencias térmicas conductivas -------- "
    Real :: t_pv = 0.002                         !Grosor del PV"
    Real :: k_pv = 50                            !Coeficiente de conductividad del PV"
    Real :: t_abs = 0.002                        !Grosor del Absorbedor"
    Real :: k_abs = 205                          !Coeficiente de conductividad del Absorbedor"
    Real :: t_sub = 0.05                         !Grosor del Sustrato"
    Real :: k_sub = 250                          !Coeficiente de conductividad del Sustrato"

    ! -------- Coeficientes y Parámetros de las resistencias térmicas convectivas -------- "
    Real :: h_f = 417.4                          !Coeficiente de convección del HTF"
    Real :: h_pv = 33.49                         !Coeficiente de convección del PV"
    Real :: h_abs = 47.36                        !Coeficiente de convección del Absorbedor"
    Real :: h_conc = 1.579                       !Coeficiente de convección del Concentrador"
    
Real :: R_cond_pv, R_cond_abs, R_cond_sub

Real :: R_conv_pv, R_conv_abs

Real :: u_agua, NUT, E, eta_pv

    ! -------- Cálculo de resitencias térmicas conductivas -------- "
!(23)" 
R_cond_pv = t_pv / (k_pv * A_pv)

!(19)" 
R_cond_abs = t_abs / (k_abs * A_abs)

!(20)"
R_cond_sub = t_sub / (k_sub * A_pv)

    ! -------- Cálculo de resitencias térmicas convectivas -------- "
!(9)" 
R_conv_pv = 1 / (h_pv * A_pv)

!(14)"
R_conv_abs = 1 / (h_abs * A_abs)

    ! -------- Cálculo de la velocidad del HTF (agua) -------- "
u_agua = 0.03537

    ! -------- Cálculo del número NUT y Efectividad de la transferencia de calor -------- "
NUT = 5.557

E = 0.9961

    ! -------- Cálculo de la eficiencia del PV en función de la relación de concentración y temperatura del PV -------- "
eta_pv = 0.3222

print *,'Sigma', sigma
print *,'C_p_f', C_p_f
print *,''

print *,'Resistencias conductivas'
print *,'Cond_pv', R_cond_pv
print *,'Cond_abs', R_cond_abs
print *,'Cond_sub', R_cond_sub
print *,''

print *,'Resistencias convectivas'
print *,'Conv_pv', R_conv_pv
print *,'Conv_abs', R_conv_abs
print *,''

print *,'Cálculos adicionales'
print *,'u_agua', u_agua
print *,'NUT', NUT
print *,'E', E
print *,'eta_pv', eta_pv
print *,''

stop


end program main