% Plott histogrammer av snodybde, snotetthet og SWE (snw water equivalent) for Rybekken og Sylsjodammen
% INPUT: .txt-filer med observasjoner fra Rybekken og Sylsjodammen med snodybde,
% snotetthet og SWE i henholdsvis andre, tredje og fjerde kolonne
clc, clear all, close all

%% Last inn data
ry = load('rybekken.txt');
sy = load('sylsjodammen.txt');

ry_depth = ry(:,2);
ry_density = ry(:,3);
ry_SWE = ry(:,4);
sy_depth = sy(:,2);
sy_density = sy(:,3);
sy_SWE = sy(:,4);

%% Plott histogram
figure
subplot(1,2,1)
hist(ry_depth)
xlabel('snodybde(cm)')
ylabel('Antall')
set(gca,'FontSize',20)
title('Histogram for Rybekken','FontSize',16)
subplot(1,2,2)
hist(sy_depth)
xlabel('snodybde(cm)')
ylabel('Antall')
set(gca,'FontSize',20)
title('Histogram for Sylsjodammen','FontSize',16)