close all;

clear;

T = readtable('efhk_weather_data_14_15.csv');

startDate = datenum(table2array(T(1,1)));

endDate = datenum(table2array(T(end,1)));

temps = table2array(T(1:end,3));

dewpts = table2array(T(1:end,6));

alldates = linspace(startDate,endDate,height(T));

axisdates = linspace(startDate,endDate,floor(height(T)/30.4));



figure('Position', [300 300 1000 500])

subplot(2,1,1);

plot(alldates,temps,'r',alldates,dewpts,'g');

grid on;

title('Hki-Vantaa keskilämpötila ja -kastepiste');

legend('Keskilämpötila','Keskikastepiste','Location','SouthWest');

ax = gca;

ax.XTick = axisdates;

datetick('x','mm/yyyy','keepticks');

ax.YAxis.TickLabelFormat = '%g^\\circC';

xlabel('Kuukausi');

ylabel('Lämpötila','rot',0);

ylab = get(gca,'ylabel');

set(ylab,'Position',get(ylab,'Position')-[-15 -28 0]);



diff = temps - dewpts;



subplot(2,1,2);

plot(alldates,diff,'b');

grid on;

title('Hki-Vantaa keskilämpötilan ja -kastepisteen erotus');

legend('Keskilämpötilan ja -kastepisteen erotus','Location','SouthWest');

ax = gca;

ax.XTick = axisdates;

datetick('x','mm/yyyy','keepticks');

ax.YAxis.TickLabelFormat = '%g^\\circC';

xlabel('Kuukausi');

ylabel('Lämpötila','rot',0);

ylab = get(gca,'ylabel');

set(ylab,'Position',get(ylab,'Position')-[-15 -14 0]);



fprintf('Keskilämpötilojen keskiarvo tarkasteluvälillä: %.2f\nKeskimääräisten kastepisteiden keskiarvo tarkasteluvälillä: %.2f\n', mean(temps), mean(dewpts));