% maketr_0604.m
% permet de créer les fichiers *.mat correspondant aux transects
% de la campagne Golts de Jun04 à partir des fichiers tet0604aXX.txt et
% tet0604goltsXX.txt
% from AP le 10/09/02
% modifié par JG
ici=pwd;

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

for xx=14:16

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

if xx==14
   xyt=xytold(:,1:132)';
   uv=uvold(1:132,:)'/100;   % pour avoir en m/s 
   save g1204tr0t2at.mat xyt uv zc
      
   xyt=xytold(:,175:247)';
   uv=uvold(175:247,:)'/100;
   uv(71:74,73)=NaN;
   save g1204tr2t3at.mat xyt uv zc
   
   xyt=xytold(:,1:247)';
   uv=uvold(1:247,:)'/100;
   %uv(71:74,73)=NaN;
   save g1204tr0t3at.mat xyt uv zc
   
   xyt=xytold(:,247:386)';
   uv=uvold(247:386,:)'/100; % pour avoir en m/s
   save g1204tr3t10at.mat xyt uv zc
   
   xyt=xytold(:,386:482)';
   uv=uvold(386:482,:)'/100; % pour avoir en m/s
   save g1204tr10t0at.mat xyt uv zc
   
   clear xyt uv xytold uvold  
end

if xx==15
   xyt=xytold(:,1:130)';
   uv=uvold(1:130,:)'/100; % pour avoir en m/s
   uv(61:62,130)=NaN;
   uv(21:28,70)=NaN;
   save g1204tr0t2bt.mat xyt uv zc
   
   xyt=xytold(:,164:300)';
   uv=uvold(164:300,:)'/100; % pour avoir en m/s
   uv(93:106,109)=NaN;
   save g1204tr2t9at.mat xyt uv zc
   
   xyt=xytold(:,300:436)';
   uv=uvold(300:436,:)'/100; % pour avoir en m/s
   uv(:,29:30)=NaN;
   uv(83:88,28)=NaN;
   uv(43:46,69)=NaN;
   save g1204tr9t2at.mat xyt uv zc
   
   xyt=xytold(:,436:696)';
   uv=uvold(436:696,:)'/100; % pour avoir en m/s
   uv(83:88,167)=NaN;
   save g1204tr2t6at.mat xyt uv zc
   
   xyt=xytold(:,696:952)';
   uv=uvold(696:952,:)'/100; % pour avoir en m/s
   save g1204tr6t2at.mat xyt uv zc
   
   xyt=xytold(:,952:1082)';
   uv=uvold(952:1082,:)'/100; % pour avoir en m/s
   uv(95:102,2)=NaN;
   save g1204tr2t0at.mat xyt uv zc   
   
   clear xyt uv xytold uvold
end

if xx==16  
   xyt=xytold(:,1:121)';
   uv=uvold(1:121,:)'/100; % pour avoir en m/s
   save g1204tr0t1at.mat xyt uv zc
   
   xyt=xytold(:,121:632)';
   uv=uvold(121:632,:)'/100; % pour avoir en m/s
   uv(91:92,511)=NaN;
   save g1204tr1t68at.mat xyt uv zc
   
   xyt=xytold(:,653:1075)';
   uv=uvold(653:1075,:)'/100; % pour avoir en m/s
   
   uv(35:36,423)=NaN;
   save g1204tr68t7at.mat xyt uv zc
   
   xyt=xytold(:,1075:1452)';
   uv=uvold(1075:1452,:)'/100; % pour avoir en m/s
   uv(35:36,1:2)=NaN;
   uv(95:100,378)=NaN;
   uv(91:92,378)=NaN;
   uv(85:88,378)=NaN;
   uv(87:92,316)=NaN;
   uv(15:18,178)=NaN;
   uv(78:88,293)=NaN;
   uv(81:82,318)=NaN;
   uv(31:36,183)=NaN;
   save g1204tr7t6at.mat xyt uv zc
   
   xyt=xytold(:,1452:1715)';
   uv=uvold(1452:1715,:)'/100; % pour avoir en m/s
   uv(83:100,1)=NaN;
   save g1204tr6t2bt.mat xyt uv zc
   
   xyt=xytold(:,1715:1851)';
   uv=uvold(1715:1851,:)'/100; % pour avoir en m/s
   uv(82:120,109:111)=NaN;
   uv(37:38,137)=NaN;
   uv(43:46,70)=NaN;
   uv(83:90,2)=NaN;
   uv(81,110)=NaN;
   save g1204tr2t9bt.mat xyt uv zc
   
   xyt=xytold(:,1851:2108)';
   uv=uvold(1851:2108,:)'/100; % pour avoir en m/s
   uv(65:66,258)=NaN;
   uv(91:92,258)=NaN;
   uv(37:38,1)=NaN;
   uv(85:86,167)=NaN;
   save g1204tr9t4at.mat xyt uv zc
   
   xyt=xytold(:,2108:2316)';
   uv=uvold(2108:2316,:)'/100; % pour avoir en m/s
   uv(:,1)=NaN;
   save g1204tr4t9at.mat xyt uv zc
   
   xyt=xytold(:,2316:2451)';
   uv=uvold(2316:2451,:)'/100; % pour avoir en m/s
   uv(:,29)=NaN;
   uv(87:98,27)=NaN;
   uv(43:44,69)=NaN;
   uv(65:66,135)=NaN;
   uv(37:38,1:2)=NaN;
   save g1204tr9t2bt.mat xyt uv zc   
   
   xyt=xytold(:,2451:2564)';
   uv=uvold(2451:2564,:)'/100; % pour avoir en m/s
   uv(31:32,56)=NaN;
   uv(15:20,64)=NaN;
   uv(5:6,112)=NaN;
   save g1204tr2t0bt.mat xyt uv zc
   
   clear xyt uv xytold uvold    
end

end

cd (ici)