% jebar_donnee.m
% calcul du transport cross-isobathe à partir de l'énergie potentielle
% ref : Hannah, 1996
% CTD : -juin et dec 2002 n'ont pas de bad flag et ont 12 col
%       -dec 2003 a une 13eme col de bad flag

g=9.81;
omega=7.27e-5;
ici=pwd;
% repertoire dans lequel sont les donnees
path='G:\golts\Survey\Survey081203\CTD\Data_extraites\asc';
cd (path)


% premier point Est (CTD9: 43.1240 ou CTD13: 43.1053 ou CTD33: 43.1051)
%-------------------
lat1=43.1240;
ctd1=load('dstat9_1212_22_dc_ali_ctm_fil_loo_dv_spli_asc.asc');
% on enleve les bad flags pour dec 2003
indflag=find(ctd1(:,13)==0);
ctd1=ctd1(indflag,:);
%
rho1=ctd1(:,10)+1000;
rho01=mean(rho1)
Profz1=-ctd1(:,6); % Prof >0 dans les donnees donc on prend - pour identique à modele
profmin1=min(Profz1)
%% on coupe à prof de la CTD Ouest
% indmin=find(Profz1>=profmin2);
% Profz1=Profz1(indmin);
% rho1=rho1(indmin);
% rho01=mean(rho1);
%%
kmax1=length(Profz1)
f1=2*omega*sin(lat1*pi/180);
qui1=0;

for k=kmax1-1:-1:1
    qui1=qui1+(g/rho01)*Profz1(k).*rho1(k).*abs(Profz1(k)-Profz1(k+1));
end

% deuxième point Ouest (CTD2: 43.07)
%---------------------
lat2=43.07;
ctd2=load('dstat2_1212_25_dc_ali_ctm_fil_loo_dv_spli_asc.asc');
% ctd2=ctd2(1:5343,:); % pour juin 2002
% on enleve les bad flags pour dec 2003
indflag=find(ctd2(:,13)==0);
ctd2=ctd2(indflag,:);
%
rho2=ctd2(:,10)+1000;
rho02=mean(rho2)
Profz2=-ctd2(:,6); % Prof >0 dans les donnees donc on prend - pour identique à modele
profmin2=min(Profz2)
%% on coupe 2 à prof de la CTD Est
% indmin=find(Profz2>=profmin1);
% Profz2=Profz2(indmin);
% rho2=rho2(indmin);
% rho02=mean(rho2);
%%
kmax2=length(Profz2)
f2=2*omega*sin(lat2*pi/180);
qui2=0;

for k=kmax2-1:-1:1
   qui2=qui2+(g/rho02)*Profz2(k).*rho2(k).*abs(Profz2(k)-Profz2(k+1));
end

% transport cross-isobathe
%-------------------------
qui1
qui2
f=(f1+f2)/2;
T=(-1/f)*(qui2-qui1)

cd (ici)


