%filtre_surface.m
% met dans les fichiers de vitesse issus de "skaphi.mat" des NaN
% supplementaires sur les 5 premiers bins quand le globalFlag de
% "a_new.mat" est a 0
%++ met dans les fichiers de vitesse issus de "skaphi.mat" des NaN
% supplementaires quand le module du courant est sup a 120 cm.s-1 et inf à
% 0 cm.s-1
ici=pwd;

for xx=5:12
    if xx<10
        cc='0';
    elseif xx>=10 & xx<100
        cc='';
    end
cd e:\Gatti2\golts\Survey\golts_2002
eval(['load TGS02',cc,int2str(xx),'aphi_new.mat'])
filtre=globalFlag;
eval(['load TGS02',cc,int2str(xx),'skaphi.mat'])
u=eastCurrentComponent;
v=northCurrentComponent;
cd (ici)
[ampl,angle]=umod(u,v);

n=nEnsembles;
m=nCells;

% suppression des  premiers bins
for j=1:2
    for i=1:n
        if filtre(i,j)==0;
            u(i,j)=NaN;
            v(i,j)=NaN;
        end
    end
end
% suppression des courants trop forts
for i=1:n
    for j=1:m
        if ampl(i,j)<=0 | ampl(i,j)>=120
            u(i,j)=NaN;
            v(i,j)=NaN;
        end
    end
end

eval(['tet1202golts',int2str(xx),'t(:,1:2:119)=u(:,1:60);'])
eval(['tet1202golts',int2str(xx),'t(:,2:2:120)=v(:,1:60);'])
eval(['save ../../Survey/Survey291102/ADCPTethys/tronca/tet1202golts',int2str(xx),'t.txt tet1202golts',int2str(xx),'t -ascii'])

end % boucle xx