k=1;
for i = 1:length(HD)
    HD(i,6)=k;
    k=k+1;
end
HD_Mon=[];
k=1;

%Number of month
for i = 1:length(HD)
    if HD(i,2) == 12
        HD_Mon(k,:) = HD(i,:);
        k=k+1;
    end
end
HD_Mon_hs=[];
k=1;

%Value of Hs
for i = 1:length(HD_Mon)
    if HD_Mon(i,5)<=1
         HD_Mon_hs(k,:) = HD_Mon(i,:);
         k=k+1;
    end
end
for i=1:minus(length(HD_Mon_hs),1)
    if minus(HD_Mon_hs(i+1,6), HD_Mon_hs(i,6))>1
        HD_Mon_hs(i,7)=1;
    end
end

Num_Intervals=0;
for i = 1:length(HD_Mon_hs)
    if HD_Mon_hs(i,7)==1
        Num_Intervals=Num_Intervals+1;
    end
end
Num_Intervals=Num_Intervals+1
Mean_Duration=length(HD_Mon_hs)*3/Num_Intervals
Percantage_Time=(length(HD_Mon_hs)/length(HD_Mon))*100