% test_ctd.m
% trace les profiles de température et de salinité en fonction de la
% profondeur
clear
ici=pwd;
%path_mat='C:\Documents and Settings\Lob\Bureau\Gatti2\golts\Survey\Survey110602\CTD\Data_extraites\asc';
path_mat='d:\golts\Survey\Survey081203\CTD\Data_extraites\asc';

%choix des profils à tracer
choix=input('combien de profils à représenter? choix= ');
for ii=1:choix
    eval(['[f_mat',int2str(ii),',p_mat',int2str(ii),']=uigetfile([path_mat,''\*.asc'']);'])
    eval(['path_mat_file',int2str(ii),'=[p_mat',int2str(ii),',f_mat',int2str(ii),'];'])
    eval(['fichier',int2str(ii),'=path_mat_file',int2str(ii),...
            '(findstr(''stat'',path_mat_file',int2str(ii),'):findstr(''_dc'',path_mat_file',...
            int2str(ii),')-1);'])
    eval(['ctd',int2str(ii),'=load(path_mat_file',int2str(ii),');'])
    
    % on enlève les bad Flags
    eval(['[m,n]=size(ctd',int2str(ii),');'])
    if n==13
        eval(['indf=find(ctd',int2str(ii),'(:,n)==0);'])
        eval(['ctd',int2str(ii),'=ctd',int2str(ii),'(indf,:);'])
    elseif n==14 % juin 2004
        eval(['ctd',int2str(ii),'=[ctd',int2str(ii),'(:,1:3) ctd',int2str(ii),'(:,5:14)];'])
        n=13;
        eval(['indf=find(ctd',int2str(ii),'(:,n)==0);'])
        eval(['ctd',int2str(ii),'=ctd',int2str(ii),'(indf,:);'])
    elseif n==12
        Flag=0*ones(m,1);
        eval(['ctd',int2str(ii),'=[ctd',int2str(ii),' Flag];'])
        n=13;
    elseif n ==10 % décembre 2004
        Nul=0*ones(m,1);
        eval(['ctd',int2str(ii),'=[ctd',int2str(ii),'(:,1:5) ctd',int2str(ii),'(:,1) ctd',int2str(ii),'(:,6:9) Nul Nul ctd',int2str(ii),'(:,10)];'])
        n=13;
        eval(['indf=find(ctd',int2str(ii),'(:,n)==0);'])
        eval(['ctd',int2str(ii),'=ctd',int2str(ii),'(indf,:);'])
    end  
end

% tracet des profils de température
figure
subplot(121)
hold on
grid on
co=['b' 'r' 'g' 'c' 'm' 'y' 'k' '--b' 'b'];
%co=['r' 'r' 'g' 'y' 'b' 'k' 'k' 'k'];
%co=['k' 'k' 'k' 'k' 'k' 'k' 'r' 'r' 'r' 'r' 'r'];

for ii=1:choix
    eval(['plot(ctd',int2str(ii),'(:,9),-ctd',int2str(ii),...
            '(:,6),co(',int2str(ii),'))'])
    eval(['str(ii)={''',eval(['fichier',int2str(ii),'']),'''}'])
end
set(gca,'ylim',[-300 0],'ytick',[-300:50:-50 -40:10:0],'fontsize',12);
%title('profils de température à la station 4 -Juin 2002');
title('Profils de T potentielle - Décembre 2003','fontsize',12);
ylabel('Profondeur (m)','fontsize',12)
legend(str,0)
hold off

% tracet des profils de salinité
subplot(122)
hold on
grid on
co=['b' 'r' 'g' 'c' 'm' 'y' 'k' '--b' '--r'];
%co=['r' 'r' 'g' 'y' 'b' 'k' 'k' 'k'];

for ii=1:choix
    eval(['plot(ctd',int2str(ii),'(:,7),-ctd',int2str(ii),...
            '(:,6),co(',int2str(ii),'))'])
    eval(['str(ii)={''',eval(['fichier',int2str(ii),'']),'''}'])
end
set(gca,'ylim',[-300 0],'ytick',[-300:50:-50 -40:10:0],'yticklabel',[],'fontsize',12);
%title('profils de salinité à la station    -Décembre 2002');
title('Profils de S - Décembre 2003','fontsize',12);
%ylabel('Profondeur (m)','fontsize',12)
legend(str,0)
hold off