%lec ctd fichier int format j;raunet
%l.prieur 10 janvier 1998 modifié février 2002
function [ctd, entete,texte]=lecint(fic)
%fic est le nom du fichier (p ex: 'int001')
%ctdint est le tableau du fichier (7 colonne, p,theta s sigmathe,cthem,chal,sommalpha)
%entete retourne les variables indiquee dans texte
%texte='[sta AA juldeb julfin latcent longcent pdeb pfin M_D]';
%station ,annee ,jjulien debut, lat(deg cent), prof de deb, index montee (-1) 
sta=0;
fid =fopen(fic);
ligne=fgetl(fid);lt=length(ligne);nli=1;
for bid=1:3;ligne=fgetl(fid);nli=nli+1; end
%on est a la 4eme ligne du fichier pour lire le no de station
%nli est le compteur de ligne lue

sta=sscanf(ligne(1:5),'%i');
ligne=fgetl(fid);ligne=fgetl(fid);ligne=fgetl(fid);nli=nli+3;
AA=sscanf(ligne(1:5),'%i');
ligne=fgetl(fid);nli=nli+1;MM=sscanf(ligne(1:3),'%i');
ligne=fgetl(fid);nli=nli+1;JJ=sscanf(ligne(1:3),'%i');
ligne=fgetl(fid);nli=nli+1;Hdeb=sscanf(ligne(1:3),'%i');
ligne=fgetl(fid);nli=nli+1;Mdeb=sscanf(ligne(1:3),'%i');
ligne=fgetl(fid);nli=nli+1;Sdeb=sscanf(ligne(1:3),'%i');
ligne=fgetl(fid);nli=nli+1;Hfin=sscanf(ligne(1:3),'%i');
ligne=fgetl(fid);nli=nli+1;Mfin=sscanf(ligne(1:3),'%i');
ligne=fgetl(fid);nli=nli+1;Sfin=sscanf(ligne(1:4),'%i');
ligne=fgetl(fid);ligne=fgetl(fid);nli=nli+2;

ligne=fgetl(fid);nli=nli+1;M_D=sscanf(ligne(1:3),'%i');
for bid=1:24;ligne=fgetl(fid);nli=nli+1;end
ligne=fgetl(fid);nli=nli+1;latdeg=sscanf(ligne(1:20),'%f');
ligne=fgetl(fid);nli=nli+1;longdeg=sscanf(ligne(1:20),'%f');
ligne=fgetl(fid);nli=nli+1;
ligne=fgetl(fid);nli=nli+1;sonde=sscanf(ligne(1:20),'%f');
ligne=fgetl(fid);nli=nli+1;pdeb=sscanf(ligne(1:20),'%f');
ligne=fgetl(fid);nli=nli+1;pfin=sscanf(ligne(1:20),'%f');

while strcmp(ligne,'     P       T        S       C')==0;ligne=fgetl(fid);nli=nli+1;end
ligne=fgetl(fid);nli=nli+1;
%debut du tableau numerique
ctdint=fscanf(fid,'%g',[4 inf]);
ctd=ctdint';
fclose(fid);

if Hfin<Hdeb;Hfin=Hfin+24;end
   juldeb=julien(MM,JJ,Hdeb,Mdeb,Sdeb);
julfin=julien(MM,JJ,Hfin,Mfin,Sfin);
latcent=degcent(latdeg*100);%conversion en degre et centieme
longcent=degcent(longdeg*100);

%je passe donc ctdint,juldeb,julfin,AA,sonde,M_D,pdeb,pfin
texte='[sta AA juldeb julfin latcent longcent pdeb pfin M_D MM JJ]';
 entete=[sta AA juldeb julfin latcent longcent pdeb pfin M_D MM JJ];
