% filtre_surface_0303.m
% créé puis modifié par JG 27/03/04 
% permet de créer les fichiers tetMMYYgoltsX.txt à partir des données tsg30iskaphi.mat
% met des NaN supplementaires sur les 2 premiers bins quand le globalFlag de
% "a_new.mat" est a 0
%+ met des NaN supplementaires quand le module du courant est sup a 120 cm/s et inf à
% 0 cm/s
% xx : numéro du fichier de session (TGS03xxskaphi.mat pour les campagnes
% de 2003)
% exemple : pour TGS0302skaphi.mat xx=02
%           pour TGS0412skaphi.mat xx=12
ici=pwd;

for xx=2:2
    if xx<10
        cc='0';
    elseif xx>=10 & xx<100
        cc='';
    end
cd e:\Gatti2\golts\Valor\Valor080303\golts_aphi2
eval(['load TGS03',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
            u(i,j)=NaN;
            v(i,j)=NaN;
    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(['tet0303golts',int2str(xx),'(:,1:2:119)=u(:,1:60);'])
eval(['tet0303golts',int2str(xx),'(:,2:2:120)=v(:,1:60);'])
eval(['save ../ADCPTethys/transect/tet0303golts',int2str(xx),'.txt tet0303golts',int2str(xx),' -ascii'])
end