%SCRIPTS PARA EL ANÁLISIS DE OLAS ÉPOCA HÚMEDA
olas_h=xlsread('olas_hum.xlsx');
t=datenum(olas_h(:,3),olas_h(:,1),olas_h(:,2),olas_h(:,4),olas_h(:,5),olas_h(:,6));

scrsz = get(0,'ScreenSize');

% Para histogramas de HS
dato=olas_h(:,7)*100;
d=0.1;   %es la distancia para ubicar el texto sobre las barras

hs_data=hist(dato,0:2:12);
[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:12)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',0:2:12,'Fontsize',11,'FontWeight','bold')
title({'Oleaje - Época Húmeda'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Número de datos','FontName','Arial','FontSize',12,'FontWeight','bold')
xlabel('Altura Significativa (cm)','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')
 
print('-dpng','-r400','HS_hum.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 TS
dato=olas_h(:,9);
d=0.1;   %es la distancia para ubicar el texto sobre las barras

Ts_data=hist(dato,6.5:0.05:6.7);
[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,6.5:0.05:6.7)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',6.5:0.05:6.7,'Fontsize',11,'FontWeight','bold')
title({'Oleaje - Época Húmeda'},'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(6.49,Ts_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(6.54,Ts_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(6.59,Ts_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(6.64,Ts_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(6.67,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','TS_hum.png')

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


%DIRECCIONES
dato=olas_h(:,12);
d=0.1;

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 Húmeda'},'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_hum.png')

clear d ans j fd cd tot_d Est str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 dir_data dato

mf30_dir=find(f_dir==6);
mf30_hs=olas_h(mf30_dir,7)*100;
mf30_ts=olas_h(mf30_dir,9);

mf20_dir=find(f_dir==2 | f_dir==8);
mf20_hs=olas_h(mf20_dir,7)*100;
mf20_ts=olas_h(mf20_dir,9);


% Para histogramas de HS
dato=mf30_hs;
d=0.05;   %es la distancia para ubicar el texto sobre las barras

hs_data=hist(dato,6:0.4:8);
[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));

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

text(5.9,hs_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(6.3,hs_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(6.7,hs_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(7.1,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(7.9,hs_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','HS_hum.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 TS
dato=olas_h(:,9);
d=0.1;   %es la distancia para ubicar el texto sobre las barras

Ts_data=hist(dato,6.5:0.05:6.7);
[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,6.5:0.05:6.7)
grid on
%SE DEBE CAMBIAR MC = Q ARRIBA
set(gca,'XTick',6.5:0.05:6.7,'Fontsize',11,'FontWeight','bold')
title({'Oleaje - Época Húmeda'},'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(6.49,Ts_data(1,1)+d,str1,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(6.54,Ts_data(2,1)+d,str2,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(6.59,Ts_data(3,1)+d,str3,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(6.64,Ts_data(4,1)+d,str4,'Fontsize',12,'Color','magenta','FontWeight','bold')
text(6.67,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','TS_hum.png')

clear j d ans fv cv tot_Ts dato str1 str2 str3 str4 str5 str6 str7 str8 str9 str10 str11 Ts_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_h(:,7)*100,'LineWidth',2)
grid on
hold on
%title({'Oleaje - Época Húmeda'},'FontSize',14,'FontName','Arial','FontWeight','bold')
ylabel('Altura de olas (cm)','FontName','Arial','FontSize',12,'FontWeight','bold')
datetick('x',13)
xlabel('25-26 Enero 2012','FontName','Arial','FontSize',12,'FontWeight','bold')

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