clc; clear; clf; close all
%% Failide struktuur
ncdisp('energy_timmean_map_1993_2023.nc')
ncdisp('energy_gor_daily_ts_1993_2023.nc')
ncdisp('energy_gor_yeardaymean.nc')
ncdisp('energy_gor_yeardaystd.nc')
ncdisp('energy_gor_yearday_mean_std.nc')

%% time mean map (Ülesanne 1. osa)
E = ncread('energy_timmean_map_1993_2023.nc','energy');
lon = ncread('energy_timmean_map_1993_2023.nc','longitude');
lat = ncread('energy_timmean_map_1993_2023.nc','latitude');
% Kontroll
size(E)
length(lon)
length(lat)
min(E(:))
mean(E(:))
max(E(:))
% Kaart
figure(1)
imagesc(lon, lat, E')
set(gca,'YDir','normal')
colorbar
xlabel('Pikkuskraad')
ylabel('Laiuskraad')
title('Tuuleenergia tiheduse ajakeskmine(1993-2023)')

%% Gulf of Riga time series
E_ts = ncread('energy_gor_daily_ts_1993_2023.nc','energy');
E_ts = squeeze(E_ts);
time = ncread('energy_gor_daily_ts_1993_2023.nc','time');
% Suuruse kontroll
size(E_ts)
% Graafik:
figure(2)
plot(E_ts)
axis tight
xlabel('Päev')
ylabel('Energia (W m^{-2})')
title('Tuuleenergia tiheduse päevakeskmine aegrida (GoR)')
grid on

%% yearday mean + yearday std (Ülesanne 2. osa)
E_ymean = ncread('energy_gor_yeardaymean.nc','energy');
E_ystd  = ncread('energy_gor_yeardaystd.nc','energy');
E_ymean = squeeze(E_ymean); E_ystd = squeeze(E_ystd);
% Suuruste kontroll
length(E_ymean)
length(E_ystd)

%% Põhigraafik
d = 1:length(E_ymean);
figure(3)
plot(d, E_ymean, 'LineWidth', 2)
axis tight
hold on
plot(d, E_ymean + E_ystd, 'r-')
plot(d, E_ymean - E_ystd,'m-')
hold off
xlabel('Aastapäev')
ylabel('Energia (W m^{-2})')
title('Tuuleenergia tiheduse keskmine ja standardhälve päeavade kaupa aastas (yearday mean +- std, GoR)')
legend('Keskmine','Keskminde + standarthälve','Keskminde - standarthälve',Location='best')
grid on

%% Jooniste salvestamine
outdir = 'wind_joonised';
if ~exist(outdir,'dir')
    mkdir(outdir);
end

figs = findall(0,'Type','figure');
for k = 1:length(figs)
    exportgraphics(figs(k), fullfile(outdir, sprintf('figure_%02d.png', k)),'Resolution', 300);
end