path = '/Users/iljagolikov/Documents/ssh/study_Ilja/baka/period/fsum_icevol_'

pikkused = [];
aastad = 1993:1:2020;
algus= [];

lopp=[];

andmed = [];

for aasta = string(1993:1:2020)

ncf = append(path,aasta,'.nc');

data = ncread(ncf, 'sithick');
data = data(:);

st = find(data,1,'first');
ed = find(data,1,'last');

if sum(data)~=0
    data(1:st)= NaN;
    data(end:-1:ed)=NaN; 
    algus(end+1)=st;
    lopp(end+1)=ed;
    pikkused(end+1) = ed - st + 1;
    andmed = cat(2,andmed,data); 
else
    pikkused(end+1) = NaN;     %kui polnud jaad siis pikkus = 0
    data(1:243)=0;           % andmed  = NaN
    algus(end+1)=NaN;         % algus = NaN
    lopp(end+1)=NaN;           % lopp = NaN
    andmed = cat(2,andmed,data);
end



k = plot(data)
hold on
k.LineWidth = 0.8
title('Kogu jää pindala 1993/94-2020/21')
xlabel('Päevade arv alates 1-st oktoobrist') 
ylabel('Jää pindala [km^2]') 
disp(aasta);
end

pikkused = pikkused(:);
aastad = aastad(:);
algus=algus(:);
lopp=lopp(:);

tabel = table(aastad,algus,lopp,pikkused);
%%

% lm1 = fitlm(aastad,algus)
% figure(1)
% subplot(2,2,1)
% plot(lm1)
% title('aastad/algus')
% 
% lm2 = fitlm(aastad,lopp)
% subplot(2,2,2)
% plot(lm2)
% title('aastad/lopp')
% 
% lm3 = fitlm(aastad,pikkused)
% subplot(2,2,3)
% plot(lm3)
% title('aastad/pikkused')
% 
% % Kogu jaapindala aastas
% sum_vol = nansum(andmed)
% lm4 = fitlm(aastad,sum_vol)
% subplot(2,2,4)
% plot(lm4)
% title('aastad/vol')

%% Statistika tabel
% stats = ["min";"max";"mean";"std";"trend";"p-value"];
% 
% algus_stat = [min(algus),max(algus),nanmean(algus),nanstd(algus),table2array(lm1.Coefficients(2,1)),table2array(lm1.Coefficients(2,4))];
% algus_stat=algus_stat(:);
% lopp_stat = [min(lopp),max(lopp),nanmean(lopp),nanstd(lopp),table2array(lm2.Coefficients(2,1)),table2array(lm2.Coefficients(2,4))];
% lopp_stat=lopp_stat(:);
% pikkused_stat = [nanmin(pikkused),nanmax(pikkused),nanmean(pikkused),nanstd(pikkused),table2array(lm3.Coefficients(2,1)),table2array(lm3.Coefficients(2,4))];
% pikkused_stat=pikkused_stat(:);
% 
% % mean paevade kaupa ja std keskmistest paevade vaartustest
% ice_vol = [min(min(andmed)),max(max(andmed)),nanmean(nanmean(andmed,2)),nanstd(nanmean(andmed,2)),table2array(lm4.Coefficients(2,1)),table2array(lm4.Coefficients(2,4))]
% ice_vol=ice_vol(:)
% 
% stat = table(stats,algus_stat,lopp_stat,pikkused_stat,ice_vol)
%%
data = andmed;
data(end:-1:220,:) = NaN;
%% koik pindalad 3D

paevad = 1:243;
f = figure(1);
f.Position = [100 100 1000 600];

fake = zeros(243,1)
%data = cat(2,data,fake);  % Kui runnid esimest korda votta komeaat maha
years = 1993:2021;
years = years(:);
data(data == 0) = NaN

time_num = [];
for aeg = 1:243;
        time_num(end+1)=datenum(1993,10,aeg) ;
end
test = datestr(time_num());
time = datetime(test)

pcolor(time,years,transpose(data));

shading flat
cb=colorbar;
grid on
cb.Label.String = 'Jää ruumala [km^3]'
cb.Label.FontWeight = 'bold'
cb.FontSize = 20
set(gca,'ytick',aastad(1:end))

ax = gca;
ax.XAxis.FontSize = 16;
ax.YAxis.FontSize = 16;
ax.LineWidth = 1;

datetick('x','mmm')
xlabel('Kuud','fontweight','bold','FontSize', 20) 
ylabel('Aastad','fontweight','bold','FontSize', 20)




