% maketr_0603.m
% permet de créer les fichiers *.mat correspondant aux transects
% de la campagne Golts de Jun03 à partir des fichiers tet0603aXX.txt et
% tet0603goltsXX.txt
% from AP le 10/09/02
% modifié par JG

zc=(-244:4:-8)'; % d'apres ConfigADCP du fichier .mat

for xx=3:4

cd e:\Gatti2\golts\Survey\Survey110603\ADCPTethys\tronca
eval(['load tet0603a',int2str(xx),'t.txt'])
eval(['xytold=tet0603a',int2str(xx),'t;'])
eval(['load tet0603golts',int2str(xx),'t.txt'])
eval(['uvold=tet0603golts',int2str(xx),'t;'])

if xx==3
   
   xyt=xytold(:,1:97)';
   uv=uvold(1:97,:)'/100;
   save g0603tr0t2at.mat xyt uv zc
      
   xyt=xytold(:,115:162)';
   uv=uvold(115:162,:)'/100; % pour avoir en m/s
   save g0603tr2tBat.mat xyt uv zc
   
   xyt=xytold(:,264:396)';
   uv=uvold(264:396,:)'/100; % pour avoir en m/s
   save g0603trBtIat.mat xyt uv zc
   
   xyt=xytold(:,699:742)';
   uv=uvold(699:742,:)'/100; % pour avoir en m/s
   save g0603trIt1at.mat xyt uv zc
   
   %xyt=xytold(:,742:end)';
   %uv=uvold(742:end,:)'/100; % pour avoir en m/s
   xyt=xytold(:,699:end)';
   uv=uvold(699:end,:)'/100; % pour avoir en m/s
end

if xx==4

%    xyt(71:148,:)=xytold(:,1:78)';
%    uv(:,71:148)=uvold(1:78,:)'/100; % pour avoir en m/s
   xyt(114:191,:)=xytold(:,1:78)';
   uv(:,114:191)=uvold(1:78,:)'/100; % pour avoir en m/s
   %save g0603tr1t9at.mat xyt uv zc
   save g0603trIt9at.mat xyt uv zc
   
   clear xyt uv
   
   xyt=xytold(:,78:280)';
   uv=uvold(78:280,:)'/100; % pour avoir en m/s
   uv(107:108,100)=NaN;
   save g0603tr9t4at.mat xyt uv zc
   
   xyt=xytold(:,280:415)';
   uv=uvold(280:415,:)'/100; % pour avoir en m/s
   uv(107:108,58)=NaN;
   uv(109:114,30)=NaN;
   save g0603tr4t6at.mat xyt uv zc
   
   xyt=xytold(:,416:754)';
   uv=uvold(416:754,:)'/100; % pour avoir en m/s
   save g0603tr6t7at.mat xyt uv zc
   
   clear xyt uv xytold uvold
end
end
for xx=10:11

cd e:\Gatti2\golts\Survey\Survey110603\ADCPTethys\tronca
eval(['load tet0603a',int2str(xx),'t.txt'])
eval(['xytold=tet0603a',int2str(xx),'t;'])
eval(['load tet0603golts',int2str(xx),'t.txt'])
eval(['uvold=tet0603golts',int2str(xx),'t;'])

if xx==10
   
   xyt=xytold(:,1:end)';
   uv=uvold(1:end,:)'/100;
   save g0603tr7t8at.mat xyt uv zc
   
   xyt=xytold(:,1:49)';
   uv=uvold(1:end,:)'/100; % pour avoir en m/s
end
if xx==11

   xyt(50:107,:)=xytold(:,1:58)';
   uv(:,50:107)=uvold(1:end,:)'/100; % pour avoir en m/s
   save g0603tr7t9at.mat xyt uv zc
   
   clear xyt uv
   
   xyt=xytold(:,1:end)';
   uv=uvold(1:end,:)'/100; % pour avoir en m/s
   save g0603tr8t9at.mat xyt uv zc
    
   clear xyt uv xytold uvold
end
end