%---------------------------------------------------------
% Overflate- og konturlinjeplott av fjelltoppen Berenberg
%---------------------------------------------------------
h0 = 2277;                          % Høyden av toppen av fjellet (m)
R = 4;                              % Mål for radius av fjellet (km)

[x,y] = meshgrid(-10:.1:10);        % Grid for x- og y-verdiene (km)
h = h0./(1+(x.^2+y.^2)./(R.^2));    % Beregn høyden h (m)


v=[0:200:2000];                     %velg hvilke konturlinjer som skal plottes
[C,hh] = contour(x,y,h,v,'b');      % Plott konturlinjene (blå farge)

axis square                         % Sett aksestørrelsene like
clabel(C,hh)                        % Sett verdi på konturlinjene
