%SCRIPTS PARA EL ANÁLISIS DE OLAS ÉPOCA SECA
olas=xlsread('EstacionW(1).xlsx','olas');
t=datenum(olas(:,1:6));

scrsz = get(0,'ScreenSize');

% Para histogramas de HS
dato=olas(:,8);
d=2;   %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));
%str9=num2str(hs_data(9,2));
%str10=num2str(hs_data(10,2));
%str11=num2str(hs_data(11,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'},'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')
%text(78,hs_data(9,1)+d,str9,'Fontsize',12,'Color','magenta','FontWeight','bold')
%text(88,hs_data(10,1)+d,str10,'Fontsize',12,'Color','magenta','FontWeight','bold')
%text(98,hs_data(11,1)+d,str11,'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')
  pause, %quitar losporcentajes 0
print('-dpng','-r400','HS.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(:,9);
d=2;   %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'},'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')
%text(20,hmax_data(8,1)+d,str8,'Fontsize',12,'Color','magenta','FontWeight','bold')
%text(23,hmax_data(9,1)+d,str9,'Fontsize',12,'Color','magenta','FontWeight','bold')
%text(26,hmax_data(10,1)+d,str10,'Fontsize',12,'Color','magenta','FontWeight','bold')
%text(29,hmax_data(11,1)+d,str11,'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.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(:,10);
d=1.5;   %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'},'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')
%text(8.7,Ts_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')
%text(9.7,Ts_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')

pause, %quitar losporcentajes 0
print('-dpng','-r400','TS.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(:,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));
%str6=num2str(Tmax_data(6,2));
%str7=num2str(Tmax_data(7,2));
%str8=num2str(Tmax_data(8,2));
%str9=num2str(Tmax_data(9,2));
%str10=num2str(Tmax_data(10,2));
%str11=num2str(Tmax_data(11,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'},'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')
%text(7.7,Tmax_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')
%text(8.7,Tmax_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')
%text(9.7,Tmax_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')

pause, %quitar losporcentajes 0
print('-dpng','-r400','Tmax.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(:,12);
d=2;

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'},'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.png')
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=find(f_dir==2 | f_dir==3);
mf_hs=olas(mf_dir,8);
mf_hmax=olas(mf_dir,9);
mf_ts=olas(mf_dir,10);
mf_tmax=olas(mf_dir,11);




% Para histogramas de HS
dato=mf_hs;
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));
%str9=num2str(hs_data(9,2));
%str10=num2str(hs_data(10,2));
%str11=num2str(hs_data(11,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 (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')
%text(78,hs_data(9,1)+d,str9,'Fontsize',12,'Color','magenta','FontWeight','bold')
%text(88,hs_data(10,1)+d,str10,'Fontsize',12,'Color','magenta','FontWeight','bold')
%text(98,hs_data(11,1)+d,str11,'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')
  pause, %quitar losporcentajes 0
print('-dpng','-r400','mf_HS.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;
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));
%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 (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.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;
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 (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.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')
%text(8.7,Ts_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')
%text(9.7,Ts_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')

print('-dpng','-r400','mf_TS.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;
d=1;   %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));
%str6=num2str(Tmax_data(6,2));
%str7=num2str(Tmax_data(7,2));
%str8=num2str(Tmax_data(8,2));
%str9=num2str(Tmax_data(9,2));
%str10=num2str(Tmax_data(10,2));
%str11=num2str(Tmax_data(11,2));


subplot(5,1,1:3)
hist(dato,3:1:7)
grid onfigure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',3:1:7,'Fontsize',11,'FontWeight','bold')
title({'Oleaje - Época Seca (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.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')
%text(7.7,Tmax_data(6,1)+d,str6,'Fontsize',12,'Color','magenta','FontWeight','bold')
%text(8.7,Tmax_data(7,1)+d,str7,'Fontsize',12,'Color','magenta','FontWeight','bold')
%text(9.7,Tmax_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')

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

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

%%%%SERIE DE TIEMPO DE ALTURA DE OLAS
figure('Position',[1 1 (3.*scrsz(3))/5 (5.*scrsz(4))/6])
subplot(5,1,1:3)
plot(t,olas(:,8),'LineWidth',2)
grid on
hold on
%title({'Oleaje - Época Seca (Dir NE-E)'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Altura de olas (m)','FontName','Arial','FontSize',12,'FontWeight','bold')
datetick('x',13)
xlabel('Diciembre 1 - 2011','FontName','Arial','FontSize',12,'FontWeight','bold')

print('-dpng','-r400','ST_olas_ep seca.png')