% Gilles Rougier  Mise à jour Mars 2003
clear all
close all

%
% choix du repertoire de Survey
% et de son fichier de parametres
%

[file_param_Survey,path_Survey]=uigetfile('../../Survey/*',' Fichier Parametre du Survey ');
path_Survey_file=[path_Survey,'CTD/Data_extraites/'];

addpath(path_Survey_file);

%
% Nom du Survey choisi

Survey=file_param_Survey(7:18);

if(Survey=='Survey061201')
	disp(Survey)
	disp(' Pas de profils CTD pendant cette campagne' )
	break
end

fluo=true;
Comment='';
if(Survey=='Survey110602')
	disp(Survey)
	disp(' Pas de fluorescence  pendant cette campagne' )
	Comment=' (pas de fluorescence) ';
	fluo=false;
end

%
% liste des fichiers a tracer contenus dans path_Survey/CTD/Data_extraites
%

d=dir([path_Survey_file,'*.asc']);

% chargement de la position des station CTD (fichier ctd_station)
ctd_station=load('Golts_station') ;

st=1;
% recherche du caractere delimiteur utilise 
delimiteur='_';
file=[path_Survey_file,d(st).name];
fichier=file(findstr('stat',file):findstr('_dc',file)-1);
if isempty(findstr('_',fichier));delimiteur='-';end

impression=input(' Voudrez vous imprimer les profils','s');

for st =1:length(d) ,
        file=[path_Survey_file,d(st).name];
	[p,filename]=fileparts(file);
	ctd=load(file);
	prof_max=ceil(max(ctd(:,6)));
	
	fichier=file(findstr('stat',file):findstr('_dc',file)-1);
	Position_du_delimiteur=findstr(delimiteur,fichier);
	Station=fichier(5:Position_du_delimiteur(1)-1);
	Numero_Station=str2num(Station);
	Date=fichier(Position_du_delimiteur(1)+1:Position_du_delimiteur(2)-1);

	% recuperation de l heure dans le fichier header
	f_hdr=fopen([p,'/',filename,'.hdr'],'r');
	while 1
	  tline=fgetl(f_hdr);
	  if tline(1:20)=='* System UpLoad Time',break,end
	  disp(tline)
	end
	fclose(f_hdr);
	Date=tline(24:length(tline));

	% montee ou descente
	sens=' Up '
	if (ctd(1,6)-ctd(length(ctd)))< 0 sens='Down'; end
	
	set(gcf,'position',[1.8 5  29.7 21],'unit','centimeter',...
		'papertype','A4','PaperOrientation','Landscape',...
		'PaperPosition',[0 0 29.7 21.])

	subplot(2,3,1)
	hold on
	title(' Temperature Potentielle','fontweight','bold','fontsize',14);
	plot(ctd(:,9),ctd(:,6));
	set(gca,'ylim',[0 prof_max],'ydir','reverse');
	grid on;

	subplot(2,3,2)
	hold on
	if fluo
		numero_param=4;
		nom_param=' Fluorescence ';
	else
		numero_param=3;
		nom_param=' Conductivite ';
	end
	title(nom_param,'fontweight','bold','fontsize',14);
	plot(ctd(:,numero_param),ctd(:,6));
	set(gca,'ylim',[0 prof_max],'ydir','reverse');
	grid on;

	subplot(2,3,4)
	hold on
	title(' Salinite ','fontweight','bold','fontsize',14);
	plot(ctd(:,7),ctd(:,6));
	set(gca,'ylim',[0 prof_max],'ydir','reverse');
	grid on;

	subplot(2,3,5)
	hold on
	title(' Densite ','fontweight','bold','fontsize',14);
	plot(ctd(:,10),ctd(:,6));
	set(gca,'ylim',[0 prof_max],'ydir','reverse');
	grid on;

	subplot(2,3,3)
	hold on
	set(gca,'visible','off');
	tS=text(.5,.9,Survey);
	tc=text(.5,.8,Comment);
	ts=text(.5,.6,[' Station ',Station,' / ',sens]);
	td=text(.5,.5,Date);
	set(tS,'Fontsize',22,'Fontweight','Bold','horizontalalignment','center');
	set(tc,'horizontalalignment','center');
	set([ts;td],'Fontsize',18,'Fontweight','Bold','horizontalalignment','center');
	
	h=subplot(2,3,6)
	topo_golts
	hold on
	plot(ctd_station(Numero_Station,3),ctd_station(Numero_Station,2),...
		'ok','markersize',16,'markerfacecolor','black');
	set(h,'position',[0.66 0.10 0.30 0.60])

	set(gcf,'position',[1.8 5  29.7 21],'unit','centimeter',...
		'papertype','A4','PaperOrientation','Landscape',...
		'PaperPosition',[0 0 29.7 21.])

	if (impression(1:1)=='o'|impression(1:1)=='O') print
	else
	ginput(1);
	end
	close 

end
