% ********************************************************************
% FIX DELLA BATIMETRIA
% ********************************************************************
% Esempio di plottaggio della batimetria 
% Nota: la procedura in questo programma consente il plottaggio della
% batimetria utilizzando quella internamente calcolata dal programma
% (sicuri di non sbagliare con i sistemi di riferimento). Ulteriori
% spiegazioni a lezione mercoledì prossimo.

clear

% carica i raggi d'onda 
% (da sostituire con i raggi d'onda calcolati personalmente!)
load A1.R
plot(A1(:,1),A1(:,2))
hold on
%load A2.R
%plot(A2(:,1),A2(:,2))
%load A3.R
%plot(A3(:,1),A3(:,2))
%load A4.R
%plot(A4(:,1),A4(:,2))
%load A5.R
%plot(A5(:,1),A5(:,2))
%load A6.R
%plot(A6(:,1),A6(:,2))
%load A7.R
%plot(A7(:,1),A7(:,2))
%load A8.R
%plot(A8(:,1),A8(:,2))
%load A9.R
%plot(A9(:,1),A9(:,2))
%load A10.R
%plot(A10(:,1),A10(:,2))
%load A11.R
%plot(A11(:,1),A11(:,2))


% carica la profondità calcolata internamente dal programma
load profondi.dat
x=profondi(:,1);
y=profondi(:,2);
h=profondi(:,3);

k=0;
for i=1:40
   for j=1:32
      k=k+1;
      X(i,j)=x(k);
		Y(i,j)=y(k);
      H(i,j)=h(k);
   end
end

V=0.5:0.5:4.5;

% curve di livello della batimetria
[P S]=contour(X,Y,H,V);
clabel(P,S);

% assi da modificare con il proprio dominio calcolato di raggi d'onda

axis([0 200 0 160])


