%SCRIPTS PARA EL ANÁLISIS DE OLAS ÉPOCA SECA SEGÚN EL ESTADO DE MAREA
%(FLUJO/REFLUJO)
%FASE DE FLUJO
% Para histogramas de HS
dato=olas(t_f,8);
d=0.4;   %es la distancia para ubicar el texto sobre las barras

hs_data=hist(dato,0:2:14);
[fv,cv]=size(hs_data);
tot_hs=sum(hs_data(fv,:));
for j=1:cv
    hs_data(fv+1,j)=hs_data(fv,j)*100/tot_hs;
end
hs_data=hs_data';
hs_data(:,2)=roundn(hs_data(:,2),-1);

str1=num2str(hs_data(1,2));
str2=num2str(hs_data(2,2));
str3=num2str(hs_data(3,2));
str4=num2str(hs_data(4,2));
str5=num2str(hs_data(5,2));
str6=num2str(hs_data(6,2));
str7=num2str(hs_data(7,2));
str8=num2str(hs_data(8,2));

figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
hist(dato,0:2:14)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',0:2:14,'Fontsize',11,'FontWeight','bold')
title({'Oleaje - Época Seca';'Flujo'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Altura Significativa (m)','FontName','Arial','FontSize',12,'FontWeight','bold')

text(-0.5,hs_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(1.5,hs_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.5,hs_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(5.5,hs_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(7.5,hs_data(5,1)+d,str5,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(9.5,hs_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(11.5,hs_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(13.5,hs_data(8,1)+d,str8,'Fontsize',12,'Color','magenta','FontWeight','bold')

gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')

print('-dpng','-r400','HS_f.png')

clear j d fv cv tot_hs dato str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 hs_data


% Para histogramas de HMAX
dato=olas(t_f,9);
d=0.4;   %es la distancia para ubicar el texto sobre las barras

hmax_data=hist(dato,0:3:18);
[fv,cv]=size(hmax_data);
tot_hmax=sum(hmax_data(fv,:));
for j=1:cv
    hmax_data(fv+1,j)=hmax_data(fv,j)*100/tot_hmax;
end
hmax_data=hmax_data';
hmax_data(:,2)=roundn(hmax_data(:,2),-1);

str1=num2str(hmax_data(1,2));
str2=num2str(hmax_data(2,2));
str3=num2str(hmax_data(3,2));
str4=num2str(hmax_data(4,2));
str5=num2str(hmax_data(5,2));
str6=num2str(hmax_data(6,2));
str7=num2str(hmax_data(7,2));
%str8=num2str(hmax_data(8,2));
%str9=num2str(hmax_data(9,2));
%str10=num2str(hmax_data(10,2));
%str11=num2str(hmax_data(11,2));

figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
hist(dato,0:3:18)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',0:3:18,'Fontsize',11,'FontWeight','bold')
title({'Oleaje - Época Seca';'Flujo'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Altura Máxima (m)','FontName','Arial','FontSize',12,'FontWeight','bold')

text(-1,hmax_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(2,hmax_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(5,hmax_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(8,hmax_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(11,hmax_data(5,1)+d,str5,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(14,hmax_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(17,hmax_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')

gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')

print('-dpng','-r400','Hmax_f.png')

clear ans j d fv cv tot_hmax dato str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 hmax_data


% Para histogramas de TS
dato=olas(t_f,10);
d=1;   %es la distancia para ubicar el texto sobre las barras

Ts_data=hist(dato,3:1:8);
[fv,cv]=size(Ts_data);
tot_Ts=sum(Ts_data(fv,:));
for j=1:cv
    Ts_data(fv+1,j)=Ts_data(fv,j)*100/tot_Ts;
end
Ts_data=Ts_data';
Ts_data(:,2)=roundn(Ts_data(:,2),-1);

str1=num2str(Ts_data(1,2));
str2=num2str(Ts_data(2,2));
str3=num2str(Ts_data(3,2));
str4=num2str(Ts_data(4,2));
str5=num2str(Ts_data(5,2));
str6=num2str(Ts_data(6,2));
%str7=num2str(Ts_data(7,2));
%str8=num2str(Ts_data(8,2));

figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
hist(dato,3:1:8)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',3:1:8,'Fontsize',11,'FontWeight','bold')
title({'Oleaje - Época Seca','Flujo'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Periodo Significativo (s)','FontName','Arial','FontSize',12,'FontWeight','bold')

text(2.7,Ts_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.7,Ts_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(4.7,Ts_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(5.7,Ts_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(6.7,Ts_data(5,1)+d,str5,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(7.7,Ts_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')

gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')

print('-dpng','-r400','TS_f.png')

clear j d ans fv cv tot_Ts dato str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 Ts_data


% Para histogramas de TMAX
dato=olas(t_f,11);
d=2;   %es la distancia para ubicar el texto sobre las barras

Tmax_data=hist(dato,3:1:7);
[fv,cv]=size(Tmax_data);
tot_Tmax=sum(Tmax_data(fv,:));
for j=1:cv
    Tmax_data(fv+1,j)=Tmax_data(fv,j)*100/tot_Tmax;
end
Tmax_data=Tmax_data';
Tmax_data(:,2)=roundn(Tmax_data(:,2),-1);

str1=num2str(Tmax_data(1,2));
str2=num2str(Tmax_data(2,2));
str3=num2str(Tmax_data(3,2));
str4=num2str(Tmax_data(4,2));
str5=num2str(Tmax_data(5,2));

figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
hist(dato,3:1:7)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',3:1:7,'Fontsize',11,'FontWeight','bold')
title({'Oleaje - Época Seca','Flujo'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Periodo Máximo (s)','FontName','Arial','FontSize',12,'FontWeight','bold')

text(2.7,Tmax_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.7,Tmax_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(4.7,Tmax_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(5.7,Tmax_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(6.7,Tmax_data(5,1)+d,str5,'Fontsize',12,'Color','magenta','FontWeight','bold')

gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')

print('-dpng','-r400','Tmax_f.png')

clear d ans j fv cv tot_Tmax dato str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 Tmax_data

%DIRECCIONES
dato=olas(t_f,12);
d=0.4;

f_dir=my_angle(dato);
dir_data=hist(f_dir,1:1:8);
[fd,cd]=size(dir_data);
tot_d=sum(dir_data(fd,:));
for j=1:cd
    dir_data(fd+1,j)=dir_data(fd,j)*100/tot_d;
end
dir_data=dir_data';
dir_data(:,2)=roundn(dir_data(:,2),-1);

str1=num2str(dir_data(1,2));
str2=num2str(dir_data(2,2));
str3=num2str(dir_data(3,2));
str4=num2str(dir_data(4,2));
str5=num2str(dir_data(5,2));
str6=num2str(dir_data(6,2));
str7=num2str(dir_data(7,2));
str8=num2str(dir_data(8,2));

figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
hist(f_dir,1:1:8)
grid on
set(gca,'XTickLabel',{'N';'NE';'E';'SE';'S';'SO';'O';'NO'},'FontSize',11,'FontWeight','bold')
title({'Oleaje - Época Seca';'Flujo'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Dirección','FontName','Arial','FontSize',12,'FontWeight','bold')

text(0.7,dir_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(1.7,dir_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(2.7,dir_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.7,dir_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(4.7,dir_data(5,1)+d,str5,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(5.7,dir_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(6.7,dir_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(7.7,dir_data(8,1)+d,str8,'Fontsize',12,'Color','magenta','FontWeight','bold')

gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')

print('-dpng','-r400','olas_dir_f.png')
f_dir_f=f_dir;
clear d ans j fd cd tot_d Est str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 dir_data dato

mf_dir_f=find(f_dir_f==2 | f_dir_f==3);
mf_hs_f=olas(mf_dir_f,8);
mf_hmax_f=olas(mf_dir_f,9);
mf_ts_f=olas(mf_dir_f,10);
mf_tmax_f=olas(mf_dir_f,11);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% FLUJO
%%%%% Para histogramas de HS
dato=mf_hs_f;
d=0.5;   %es la distancia para ubicar el texto sobre las barras

hs_data=hist(dato,0:2:14);
[fv,cv]=size(hs_data);
tot_hs=sum(hs_data(fv,:));
for j=1:cv
    hs_data(fv+1,j)=hs_data(fv,j)*100/tot_hs;
end
hs_data=hs_data';
hs_data(:,2)=roundn(hs_data(:,2),-1);

str1=num2str(hs_data(1,2));
str2=num2str(hs_data(2,2));
str3=num2str(hs_data(3,2));
str4=num2str(hs_data(4,2));
str5=num2str(hs_data(5,2));
str6=num2str(hs_data(6,2));
str7=num2str(hs_data(7,2));
str8=num2str(hs_data(8,2));

figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
hist(dato,0:2:14)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',0:2:14,'Fontsize',11,'FontWeight','bold')
title({'Oleaje - Época Seca';'Flujo (Dir NE-E)'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Altura Significativa (m)','FontName','Arial','FontSize',12,'FontWeight','bold')

text(-0.5,hs_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(1.5,hs_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.5,hs_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(5.5,hs_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(7.5,hs_data(5,1)+d,str5,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(9.5,hs_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(11.5,hs_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(13.5,hs_data(8,1)+d,str8,'Fontsize',12,'Color','magenta','FontWeight','bold')

gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')

print('-dpng','-r400','mf_HS_f.png')

clear ans j d fv cv tot_hs dato str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 hs_data


% Para histogramas de HMAX
dato=mf_hmax_f;
d=0.5;   %es la distancia para ubicar el texto sobre las barras

hmax_data=hist(dato,0:3:18);
[fv,cv]=size(hmax_data);
tot_hmax=sum(hmax_data(fv,:));
for j=1:cv
    hmax_data(fv+1,j)=hmax_data(fv,j)*100/tot_hmax;
end
hmax_data=hmax_data';
hmax_data(:,2)=roundn(hmax_data(:,2),-1);

str1=num2str(hmax_data(1,2));
str2=num2str(hmax_data(2,2));
str3=num2str(hmax_data(3,2));
str4=num2str(hmax_data(4,2));
str5=num2str(hmax_data(5,2));
str6=num2str(hmax_data(6,2));
str7=num2str(hmax_data(7,2));

figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
hist(dato,0:3:18)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',0:3:18,'Fontsize',11,'FontWeight','bold')
title({'Oleaje - Época Seca';'Flujo (Dir NE-E)'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Altura Máxima (m)','FontName','Arial','FontSize',12,'FontWeight','bold')

text(-1,hmax_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(2,hmax_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(5,hmax_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(8,hmax_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(11,hmax_data(5,1)+d,str5,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(14,hmax_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(17,hmax_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')

gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')

print('-dpng','-r400','mf_Hmax_f.png')

clear ans j d fv cv tot_hmax dato str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 hmax_data

% Para histogramas de TS
dato=mf_ts_f;
d=0.3;   %es la distancia para ubicar el texto sobre las barras

Ts_data=hist(dato,3:0.1:3.6);
[fv,cv]=size(Ts_data);
tot_Ts=sum(Ts_data(fv,:));
for j=1:cv
    Ts_data(fv+1,j)=Ts_data(fv,j)*100/tot_Ts;
end
Ts_data=Ts_data';
Ts_data(:,2)=roundn(Ts_data(:,2),-1);

str1=num2str(Ts_data(1,2));
str2=num2str(Ts_data(2,2));
str3=num2str(Ts_data(3,2));
str4=num2str(Ts_data(4,2));
str5=num2str(Ts_data(5,2));
str6=num2str(Ts_data(6,2));
str7=num2str(Ts_data(7,2));

figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
hist(dato,3:0.1:3.6)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',3:0.1:3.6,'Fontsize',11,'FontWeight','bold')
title({'Oleaje - Época Seca';'Flujo (Dir NE-E)'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Periodo Significativo (s)','FontName','Arial','FontSize',12,'FontWeight','bold')

text(2.98,Ts_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.08,Ts_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.18,Ts_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.28,Ts_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.38,Ts_data(5,1)+d,str5,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.48,Ts_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.58,Ts_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')

gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')

print('-dpng','-r400','mf_TS_f.png')

clear j d ans fv cv tot_Ts dato str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 Ts_data

% Para histogramas de TMAX
dato=mf_tmax_f;
d=0.3;   %es la distancia para ubicar el texto sobre las barras

Tmax_data=hist(dato,3:0.1:3.6);
[fv,cv]=size(Tmax_data);
tot_Tmax=sum(Tmax_data(fv,:));
for j=1:cv
    Tmax_data(fv+1,j)=Tmax_data(fv,j)*100/tot_Tmax;
end
Tmax_data=Tmax_data';
Tmax_data(:,2)=roundn(Tmax_data(:,2),-1);

str1=num2str(Tmax_data(1,2));
str2=num2str(Tmax_data(2,2));
str3=num2str(Tmax_data(3,2));
str4=num2str(Tmax_data(4,2));
str5=num2str(Tmax_data(5,2));
str6=num2str(Tmax_data(6,2));
str7=num2str(Tmax_data(7,2));

figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
hist(dato,3:0.1:3.6)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',3:0.1:3.6,'Fontsize',11,'FontWeight','bold')
title({'Flujo (Dir NE-E)'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Periodo Máximo (s)','FontName','Arial','FontSize',12,'FontWeight','bold')

text(2.98,Tmax_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.08,Tmax_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.18,Tmax_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.28,Tmax_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.38,Tmax_data(5,1)+d,str5,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.48,Tmax_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.58,Tmax_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')

gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')

print('-dpng','-r400','mf_Tmax_f.png')

clear d ans j fv cv tot_Tmax dato str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 Tmax_data


%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%
%FASE DE REFLUJO
% Para histogramas de HS
dato=olas(t_rf,8);
d=0.6;   %es la distancia para ubicar el texto sobre las barras

hs_data=hist(dato,0:2:14);
[fv,cv]=size(hs_data);
tot_hs=sum(hs_data(fv,:));
for j=1:cv
    hs_data(fv+1,j)=hs_data(fv,j)*100/tot_hs;
end
hs_data=hs_data';
hs_data(:,2)=roundn(hs_data(:,2),-1);

str1=num2str(hs_data(1,2));
str2=num2str(hs_data(2,2));
str3=num2str(hs_data(3,2));
str4=num2str(hs_data(4,2));
str5=num2str(hs_data(5,2));
str6=num2str(hs_data(6,2));
str7=num2str(hs_data(7,2));
str8=num2str(hs_data(8,2));

figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
hist(dato,0:2:14)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',0:2:14,'Fontsize',11,'FontWeight','bold')
title({'Oleaje - Época Seca';'Reflujo'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Altura Significativa (m)','FontName','Arial','FontSize',12,'FontWeight','bold')

text(-0.5,hs_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(1.5,hs_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.5,hs_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(5.5,hs_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(7.5,hs_data(5,1)+d,str5,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(9.5,hs_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(11.5,hs_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(13.5,hs_data(8,1)+d,str8,'Fontsize',12,'Color','magenta','FontWeight','bold')

gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')

print('-dpng','-r400','HS_rf.png')

clear j d fv cv tot_hs dato str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 hs_data


% Para histogramas de HMAX
dato=olas(t_rf,9);
d=0.6;   %es la distancia para ubicar el texto sobre las barras

hmax_data=hist(dato,0:3:18);
[fv,cv]=size(hmax_data);
tot_hmax=sum(hmax_data(fv,:));
for j=1:cv
    hmax_data(fv+1,j)=hmax_data(fv,j)*100/tot_hmax;
end
hmax_data=hmax_data';
hmax_data(:,2)=roundn(hmax_data(:,2),-1);

str1=num2str(hmax_data(1,2));
str2=num2str(hmax_data(2,2));
str3=num2str(hmax_data(3,2));
str4=num2str(hmax_data(4,2));
str5=num2str(hmax_data(5,2));
str6=num2str(hmax_data(6,2));
str7=num2str(hmax_data(7,2));

figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
hist(dato,0:3:18)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',0:3:18,'Fontsize',11,'FontWeight','bold')
title({'Oleaje - Época Seca';'Reflujo'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Altura Máxima (m)','FontName','Arial','FontSize',12,'FontWeight','bold')

text(-1,hmax_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(2,hmax_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(5,hmax_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(8,hmax_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(11,hmax_data(5,1)+d,str5,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(14,hmax_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(17,hmax_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')

gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')

print('-dpng','-r400','Hmax_rf.png')

clear ans j d fv cv tot_hmax dato str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 hmax_data


% Para histogramas de TS
dato=olas(t_rf,10);
d=0.3;   %es la distancia para ubicar el texto sobre las barras

Ts_data=hist(dato,3:0.1:3.6);
[fv,cv]=size(Ts_data);
tot_Ts=sum(Ts_data(fv,:));
for j=1:cv
    Ts_data(fv+1,j)=Ts_data(fv,j)*100/tot_Ts;
end
Ts_data=Ts_data';
Ts_data(:,2)=roundn(Ts_data(:,2),-1);

str1=num2str(Ts_data(1,2));
str2=num2str(Ts_data(2,2));
str3=num2str(Ts_data(3,2));
str4=num2str(Ts_data(4,2));
str5=num2str(Ts_data(5,2));
str6=num2str(Ts_data(6,2));
str7=num2str(Ts_data(7,2));

figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
hist(dato,3:0.1:3.6)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',3:1:8,'Fontsize',11,'FontWeight','bold')
title({'Oleaje - Época Seca','Reflujo'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Periodo Significativo (s)','FontName','Arial','FontSize',12,'FontWeight','bold')

text(2.98,Ts_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.08,Ts_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.18,Ts_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.28,Ts_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.38,Ts_data(5,1)+d,str5,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.48,Ts_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.58,Ts_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')

gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')

print('-dpng','-r400','TS_rf.png')

clear j d ans fv cv tot_Ts dato str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 Ts_data


% Para histogramas de TMAX
dato=olas(t_rf,11);
d=0.4;   %es la distancia para ubicar el texto sobre las barras

Tmax_data=hist(dato,3:0.1:3.6);
[fv,cv]=size(Tmax_data);
tot_Tmax=sum(Tmax_data(fv,:));
for j=1:cv
    Tmax_data(fv+1,j)=Tmax_data(fv,j)*100/tot_Tmax;
end
Tmax_data=Tmax_data';
Tmax_data(:,2)=roundn(Tmax_data(:,2),-1);

str1=num2str(Tmax_data(1,2));
str2=num2str(Tmax_data(2,2));
str3=num2str(Tmax_data(3,2));
str4=num2str(Tmax_data(4,2));
str5=num2str(Tmax_data(5,2));
str6=num2str(Tmax_data(6,2));
str7=num2str(Tmax_data(7,2));

figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
hist(dato,3:0.1:3.6)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',3:0.1:3.6,'Fontsize',11,'FontWeight','bold')
title({'Reflujo'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Periodo Máximo (s)','FontName','Arial','FontSize',12,'FontWeight','bold')

text(2.98,Tmax_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.08,Tmax_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.18,Tmax_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.28,Tmax_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.38,Tmax_data(5,1)+d,str5,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.48,Tmax_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.58,Tmax_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')

gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')

print('-dpng','-r400','Tmax_rf.png')

clear d ans j fv cv tot_Tmax dato str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 Tmax_data

%DIRECCIONES
dato=olas(t_rf,12);
d=0.4;

f_dir=my_angle(dato);
dir_data=hist(f_dir,1:1:8);
[fd,cd]=size(dir_data);
tot_d=sum(dir_data(fd,:));
for j=1:cd
    dir_data(fd+1,j)=dir_data(fd,j)*100/tot_d;
end
dir_data=dir_data';
dir_data(:,2)=roundn(dir_data(:,2),-1);

str1=num2str(dir_data(1,2));
str2=num2str(dir_data(2,2));
str3=num2str(dir_data(3,2));
str4=num2str(dir_data(4,2));
str5=num2str(dir_data(5,2));
str6=num2str(dir_data(6,2));
str7=num2str(dir_data(7,2));
str8=num2str(dir_data(8,2));

figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
hist(f_dir,1:1:8)
grid on
set(gca,'XTickLabel',{'N';'NE';'E';'SE';'S';'SO';'O';'NO'},'FontSize',11,'FontWeight','bold')
title({'Oleaje - Época Seca';'Reflujo'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Dirección','FontName','Arial','FontSize',12,'FontWeight','bold')

text(0.7,dir_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(1.7,dir_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(2.7,dir_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.7,dir_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(4.7,dir_data(5,1)+d,str5,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(5.7,dir_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(6.7,dir_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(7.7,dir_data(8,1)+d,str8,'Fontsize',12,'Color','magenta','FontWeight','bold')

gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')

print('-dpng','-r400','olas_dir_rf.png')
f_dir_rf=f_dir;
clear d f_dir ans j fd cd tot_d Est str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 dir_data dato

mf_dir_rf=find(f_dir_rf==2 | f_dir_rf==3);
mf_hs_rf=olas(mf_dir_rf,8);
mf_hmax_rf=olas(mf_dir_rf,9);
mf_ts_rf=olas(mf_dir_rf,10);
mf_tmax_rf=olas(mf_dir_rf,11);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% REFLUJO
%%%%% Para histogramas de HS
dato=mf_hs_rf;
d=0.5;   %es la distancia para ubicar el texto sobre las barras

hs_data=hist(dato,0:2:14);
[fv,cv]=size(hs_data);
tot_hs=sum(hs_data(fv,:));
for j=1:cv
    hs_data(fv+1,j)=hs_data(fv,j)*100/tot_hs;
end
hs_data=hs_data';
hs_data(:,2)=roundn(hs_data(:,2),-1);

str1=num2str(hs_data(1,2));
str2=num2str(hs_data(2,2));
str3=num2str(hs_data(3,2));
str4=num2str(hs_data(4,2));
str5=num2str(hs_data(5,2));
str6=num2str(hs_data(6,2));
str7=num2str(hs_data(7,2));
str8=num2str(hs_data(8,2));

figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
hist(dato,0:2:14)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',0:2:14,'Fontsize',11,'FontWeight','bold')
title({'Oleaje - Época Seca';'Reflujo (Dir NE-E)'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Altura Significativa (m)','FontName','Arial','FontSize',12,'FontWeight','bold')

text(-0.5,hs_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(1.5,hs_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.5,hs_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(5.5,hs_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(7.5,hs_data(5,1)+d,str5,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(9.5,hs_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(11.5,hs_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(13.5,hs_data(8,1)+d,str8,'Fontsize',12,'Color','magenta','FontWeight','bold')

gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')

print('-dpng','-r400','mf_HS_rf.png')

clear ans j d fv cv tot_hs dato str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 hs_data


% Para histogramas de HMAX
dato=mf_hmax_rf;
d=0.5;   %es la distancia para ubicar el texto sobre las barras

hmax_data=hist(dato,0:3:18);
[fv,cv]=size(hmax_data);
tot_hmax=sum(hmax_data(fv,:));
for j=1:cv
    hmax_data(fv+1,j)=hmax_data(fv,j)*100/tot_hmax;
end
hmax_data=hmax_data';
hmax_data(:,2)=roundn(hmax_data(:,2),-1);

str1=num2str(hmax_data(1,2));
str2=num2str(hmax_data(2,2));
str3=num2str(hmax_data(3,2));
str4=num2str(hmax_data(4,2));
str5=num2str(hmax_data(5,2));
str6=num2str(hmax_data(6,2));
str7=num2str(hmax_data(7,2));

figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
hist(dato,0:3:18)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',0:3:18,'Fontsize',11,'FontWeight','bold')
title({'Oleaje - Época Seca';'Reflujo (Dir NE-E)'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Altura Máxima (m)','FontName','Arial','FontSize',12,'FontWeight','bold')

text(-1,hmax_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(2,hmax_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(5,hmax_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(8,hmax_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(11,hmax_data(5,1)+d,str5,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(14,hmax_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(17,hmax_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')

gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')

print('-dpng','-r400','mf_Hmax_rf.png')

clear ans j d fv cv tot_hmax dato str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 hmax_data

% Para histogramas de TS
dato=mf_ts_rf;
d=0.3;   %es la distancia para ubicar el texto sobre las barras

Ts_data=hist(dato,3:0.1:3.6);
[fv,cv]=size(Ts_data);
tot_Ts=sum(Ts_data(fv,:));
for j=1:cv
    Ts_data(fv+1,j)=Ts_data(fv,j)*100/tot_Ts;
end
Ts_data=Ts_data';
Ts_data(:,2)=roundn(Ts_data(:,2),-1);

str1=num2str(Ts_data(1,2));
str2=num2str(Ts_data(2,2));
str3=num2str(Ts_data(3,2));
str4=num2str(Ts_data(4,2));
str5=num2str(Ts_data(5,2));
str6=num2str(Ts_data(6,2));
str7=num2str(Ts_data(7,2));

figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
hist(dato,3:0.1:3.6)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',3:0.1:3.6,'Fontsize',11,'FontWeight','bold')
title({'Oleaje - Época Seca';'Reflujo (Dir NE-E)'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Periodo Significativo (s)','FontName','Arial','FontSize',12,'FontWeight','bold')

text(2.98,Ts_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.08,Ts_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.18,Ts_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.28,Ts_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.38,Ts_data(5,1)+d,str5,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.48,Ts_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.58,Ts_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')

gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')

print('-dpng','-r400','mf_TS_rf.png')

clear j d ans fv cv tot_Ts dato str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 Ts_data

% Para histogramas de TMAX
dato=mf_tmax_rf;
d=0.3;   %es la distancia para ubicar el texto sobre las barras

Tmax_data=hist(dato,3:0.1:3.6);
[fv,cv]=size(Tmax_data);
tot_Tmax=sum(Tmax_data(fv,:));
for j=1:cv
    Tmax_data(fv+1,j)=Tmax_data(fv,j)*100/tot_Tmax;
end
Tmax_data=Tmax_data';
Tmax_data(:,2)=roundn(Tmax_data(:,2),-1);

str1=num2str(Tmax_data(1,2));
str2=num2str(Tmax_data(2,2));
str3=num2str(Tmax_data(3,2));
str4=num2str(Tmax_data(4,2));
str5=num2str(Tmax_data(5,2));
str6=num2str(Tmax_data(6,2));
str7=num2str(Tmax_data(7,2));

figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
hist(dato,3:0.1:3.6)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',3:0.1:3.6,'Fontsize',11,'FontWeight','bold')
title({'Oleaje - Época Seca';'Reflujo (Dir NE-E)'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Periodo Máximo (s)','FontName','Arial','FontSize',12,'FontWeight','bold')

text(2.975,Tmax_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.075,Tmax_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.175,Tmax_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.275,Tmax_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.375,Tmax_data(5,1)+d,str5,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.475,Tmax_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(3.575,Tmax_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')

gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')
gtext('%','Fontsize',12,'Color','magenta','FontWeight','bold')

print('-dpng','-r400','mf_Tmax_rf.png')

clear d ans j fv cv tot_Tmax dato str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 Tmax_data