% Innlesing av data
p = load('trykkfelt.dat');
u = load('vindfelt_u.dat');
v = load('vindfelt_v.dat');

% Vi finner st鴕relsen av dataene
disp 'St鴕relsene til p, u og v er:'
size(p)
size(u)
size(v)

% Plott vindfeltet (vektorfeltet) v = ui + vj i gitterpunktene
figure(1)
quiver(u,v)
axis([0 32 0 31])

% Plott vindfeltet med avstand i kilometer langs aksene
figure(2)
x = 0:55:55*30;
y = 0:55:55*29;
quiver(x,y, u,v);
axis([-55 55*31 -55 55*30])

% Beregn vindhastigheten og finn storste verdi
l = sqrt(u.^2 + v.^2);
maks = max(max(l))

% Plott isobarene med 5 hektopascal intervall
figure(3)
isobarer = 980:5:1025;
[C,h] = contour(x,y,p,isobarer);
clabel(C, h,'fontname','Times','fontsize',20);
axis([-55 55*31 -55 55*30])

% Regn ut divergensen og plott konturlinjene
figure(4)
div = divergence(u,v);
[C,h] = contour(x,y,div,-6:1:2);
clabel(C,h,'fontname','Times','fontsize',20);
axis([-55 55*31 -55 55*30])

% Regn ut virvlingen og plott konturlinjene
figure(5)
[curlz,cav] = curl(u,v);
[C,h] = contour(x,y,curlz,-10:5:20);
clabel(C,h,'fontname','Times','fontsize',20);
axis([-55 55*31 -55 55*30])