%---------------------------------------------------------
% 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);        % Lager grid for x- og y-verdiene (km)
                                    % fra -10 til +10 med steg 0.1
h = h0./(1+(x.^2+y.^2)./(R.^2));    % Beregn h?yden h (m) i alle gitterpunkter
                                    % Legg merke til symbolet .^2 Sl?yfes
                                    % . mutipliseres ikke vektoren x
                                    % sammen ledd for ledd. Tilsvarende for
                                    % multiplikasjon .* og divisjon ./ 
v=[500:500:2000]                    % Velger konturlinjer fra 500m
                                    % til 2000m med avstand 500m
[C,hh] = contour(x,y,h,v,'b');      % Plott konturlinjene (bl? farge)
clabel(C,hh)                        % Plasser konturlinjene
%surf(x,y,h), shading interp        % Plotter fargekonturer
%colorbar                           % med fargeskala
axis square                         % Sett aksest?rrelsene like
% Onskes konturlinjer ta bort % i surf og colorbar og sett % p?
% de tre foreg?ende kommandoene.
% 3-D plott f?es ved ? g? inn p? tools i plottevinduet.
%NB! NB! Det som st?r bak % tegnet leses som kommentar ikke kommando