% levhar.m
% -----------------------------------------------------------------------------
% Analyse Harmonique des Marées
% Référence : Pawlowicz, R., B. Beardsley and S. Lentz. 2002.
%             Classical harmonic analysis including error estimates in MATLAB,
%             using T_TIDE. Computers and geosciences.

% Cours : Analyse et modélisation du milieu marin - Hiver 2002.

% V. G. Koutitonsky
% ----------------------------------------------------------------------------

load level.dat;                         % charger la série de niveau
xin=level(:,1);                         % extraire, si nécéssaire la bonne colonne
start=datenum('06-24-1998 14:08:00');   % préciser la date de début (mm-jj-aa hh:mi:ss) 
lat=(30+20/60);                         % préciser la lattitude
N=length(xin);                          % calculer la longueur de la série
dt=.1667;                                   % préciser le pas temporel (en heure) 
t=[0:N-1];    
infername=['P1';'K2'];
inferfrom=['K1';'S2'];
infamp=[.33093;.27215];
infphase=[-7.07;-22.40];
 
[tidestruc,pout]=t_tide(xin,...
'start_time',start,...
'interval',dt,...
'latitude',lat,...
'error','cboot',...
'inference',infername,inferfrom,infamp,infphase,...
'output','levhares.dat');               % donner un nom au fichier de résultats.

res=xin-pout;
subplot(3,1,1)
plot(t,xin-mean(xin));
title('Observed')
subplot(3,1,2)
plot(t,pout-mean(pout));
title('Low-passed')
subplot(3,1,3)
plot(t,res-mean(res));
title('High-passed')
xlabel('temps (heure)')
ylabel('niveau (m)')