% f chart
clear all
Ta = [22.7 22.3 20.5 17.2 13.9 11.0 10.7 11.5 13.2 15.7 18.3 16.5]; %temp ambiente media mensual montevideo AMT
%Tm = [28.4 28.5 26.1 21.7 16.6 12.0 09.3 09.2 11.7 16.1 21.3 18.9]; %dato de ose (falso)
Tm = [23.1234 23.1755 22.3758 20.8911 19.1245 17.5520 16.6009 16.5290 17.3577 18.8604 20.6280 22.1848]; %resultados promedio mensual de este calculo
Ht = [21.5613 21.2250 19.0581 18.4067 12.7516 12.3633 12.1065 15.4710 17.1467 19.8355 19.0200 21.4903]*1e6; %irradiación media mensual sobre plano inclinado a 45º en J/m2
Ni = [31 28 31 30 31 30 31 31 30 31 30 31]; %días por mes
dt = 24*60*60*Ni; %segundos por mes
Xi = zeros(1,12);
Yi = Xi;
fi = Xi;

% parámetros:
MAI = 0.94; %colectores planos
FR = 0.9;
UL = 4.0/FR;
tauAn = 0.88;
eps = 0.8; %eficiencia intercambiador
cp = 3900; %de la mezcla de anticongelante por el circuito primario
Tw = 60; %temp mínima aceptable
CT = (11.6 + 1.18*Tw + 3.86*Tm - 2.32*Ta)./(100-Ta);

Vacs = 2800; %consumo diario de litros

Li = Vacs*1*4186*(45 - Tm).*Ni; %consumo por mes en J
tot=1; %cant casos
ff=zeros(1,tot);
gg=ff;

for i=1:tot

N = 20; %num colect

Vacu = 3000; %volumen acumulación

Ac = 1.9*N; %área acumulación
mp = N*100*1.03/3600; %lt/h * kg/lt * segundos/h
FpR = FR/(1 + (Ac*FR*UL/mp/cp)*(1/eps - 1));
CV = (Vacu/Ac/75)^(-0.25); %ajuste si no hay 75 litros por m2 de colector
Xi = FR*UL*FpR/FR*(100-Ta).*CV.*Ac.*CT.*dt./Li;
Yi = FR*tauAn*FpR/FR*MAI*Ht.*Ni*Ac./Li;
fi = 1.029*Yi - 0.065*Xi - 0.245*Yi.^2 + 0.0018*Xi.^2 + 0.0215*Yi.^3;

for j=1:12
    if fi(j)>1
        gg(i)=gg(i)+(fi(j)-1);
        fi(j)=1;
    end
end
ff(i) = sum(fi.*Li)/sum(Li);

end