diff --git a/@netcdf_native/tsgqc_netcdf.xls b/@netcdf_native/tsgqc_netcdf.xls
index 5789d1ef1d2dd9b532bf5983045c44b4595d0dfc..836a1cc96f32113e12968a39039dab9b18c6d0e2 100644
Binary files a/@netcdf_native/tsgqc_netcdf.xls and b/@netcdf_native/tsgqc_netcdf.xls differ
diff --git a/@tsg_nc/tsg_platform.csv b/@tsg_nc/tsg_platform.csv
index 06d599ec2c9e695c480e4fe881926b78c8dd39f1..7d6c4dbde55ff82ac64e972ebdd25407c80e84e2 100644
--- a/@tsg_nc/tsg_platform.csv
+++ b/@tsg_nc/tsg_platform.csv
@@ -1,5 +1,6 @@
-% tsg_platform.csv: 14 members & 0 lines with data; VERSION = 1.4850
+% tsg_platform.csv: 15 members & 1 lines with data; VERSION = 1.5000
%
%
id;PLATFORM_NAME;SHIP_CALL_SIGN;SHIP_MMSI;PROJECT_NAME;TYPE_TSG;NUMBER_TSG;TYPE_TINT;NUMBER_TINT;SAMPLING_PERIOD;SSPS_DEPH;SSPS_DEPH_MIN;SSPS_DEPH_MAX;SSTP_DEPH;SSTP_DEPH_MIN;SSTP_DEPH_MAX;DATE_TSG;endl
-char;char;char;char;char;char;char;char;char;integer;integer;float;float;float;float;float;float;float
+char;char;char;char;char;char;char;char;char;integer;integer;float;float;float;float;float;float;float
+#;COLIBRI;;;;;;;;;5;4.000000;6.000000;NaN;NaN;NaN;;#
diff --git a/tsg_tools/CO2/SSJT.png b/tsg_tools/CO2/SSJT.png
new file mode 100644
index 0000000000000000000000000000000000000000..cf08c4c932205c41b20cce2af524720fb6a6b6ab
Binary files /dev/null and b/tsg_tools/CO2/SSJT.png differ
diff --git a/tsg_tools/CO2/SSJT_QC.png b/tsg_tools/CO2/SSJT_QC.png
new file mode 100644
index 0000000000000000000000000000000000000000..7c6cd573d952d8bb6007093e2ca704e692a18f11
Binary files /dev/null and b/tsg_tools/CO2/SSJT_QC.png differ
diff --git a/tsg_tools/CO2/SSPS_QC.png b/tsg_tools/CO2/SSPS_QC.png
new file mode 100644
index 0000000000000000000000000000000000000000..56d30b9372a8540b2fb7d98bc44ff0f328b55906
Binary files /dev/null and b/tsg_tools/CO2/SSPS_QC.png differ
diff --git a/tsg_tools/CO2/SSPS_et_co2_SSPS.png b/tsg_tools/CO2/SSPS_et_co2_SSPS.png
new file mode 100644
index 0000000000000000000000000000000000000000..477c7d729f51de3714517bbc30b574ea83626d4d
Binary files /dev/null and b/tsg_tools/CO2/SSPS_et_co2_SSPS.png differ
diff --git a/tsg_tools/CO2/TD280 Oxygen Optode Calculations.xls b/tsg_tools/CO2/TD280 Oxygen Optode Calculations.xls
new file mode 100644
index 0000000000000000000000000000000000000000..e8dc4a27f95ab3b8934a0254d574858f990d3308
Binary files /dev/null and b/tsg_tools/CO2/TD280 Oxygen Optode Calculations.xls differ
diff --git a/tsg_tools/CO2/concatCO2.m b/tsg_tools/CO2/concatCO2.m
new file mode 100644
index 0000000000000000000000000000000000000000..14a20a746da12e1a36bde9a257de9ef43655102d
--- /dev/null
+++ b/tsg_tools/CO2/concatCO2.m
@@ -0,0 +1,378 @@
+function concatCO2
+%
+% Cette fonction permet de concaténer les fichiers CO2 de type *data.txt
+% Elle crée un fichier unique au format Excel, *.CSV, dont les colonnes
+% sont séparées par des points virgules.
+% La lecture des fichiers peut se faire dans n'importe quel ordre, les
+% enregistrement sont triés, avant écriture, par ordre croissannt de date.
+%
+% ATTENTION
+% 1 - On considère que le nombre de colonnes est fixe : 34. Il est difficile
+% de tester le nombre de colonnes des fichiers car les noms d'entêtes
+% ne sont pas tous des chaines de caractères liées. Ils peuvent être
+% séparés par des blancs (exemple : "GPS TIME"). Autre difficulté :
+% l'absence de point GPS se traduit par une heure au format '::' au
+% lieu du format habituel hh:mmm:ss
+% 2 - Le format des fichiers Craig Neil diffère du format des
+% fichiers G.O.
+% Les fichiers G.O. peuvent contenir des colonnes supplémentaires :
+% celles-ci ne sont pas lues
+% La variable "cond temp" est placée en colonne 32 dans les fichiers GO
+% et en colonne 28 dans les fichiers Craig Neil. Le fichier final est
+% au format Craig Neil
+% 3 - On considère que les entêtes des colonnes ne varient pas. Un test
+% contrôle la chaine de caractères de l'entête qui est lue avec
+% la chaine de caractères "HeaderIn".
+% 4 - Le séparateur numérique des nombres est le '.'. Selon la
+% configuration des paramètres internationaux de l'ordinateur, le
+% séparateur peut être une ','. Pourqu'Excel considère que ce sont
+% bien des nombres il faudra remplacer tous le '.' par des ','
+% 5 - Les dates au format jj/mm/aaaa hh:mm:ss dans un fichier *.CSV sont
+% affichées dans Excel avec le format jj/mm/aaaa hh:mm. Il manque
+% les secondes.
+% Si le fichier devait être modifié et sauvegardé au format *.CSV
+% il convient de modifier le format de la colonne des date pour que
+% celle-ci affiche bien les secondes.
+%
+% Le programme
+% 1 - permet de regrouper les colonnes "Date" et "Time" (du PC) dans une
+% seule colonne appelée : "Date Time".
+% 2 - permet de sélectionner plusieurs fichiers.
+% Le repertoire où se trouvent les fichiers peut être modifié à la ligne
+% XXXX - Variable repIn
+% 3 - permet de donner le nom du fichier concaténé
+% Le repertoire où sera sauvegardé le fichier peut être modifié à la ligne
+% XXXX - Variable repOut
+%
+% ERREUR
+% Le programme teste si le fichier n'est pas complet ou si le nombre de
+% colonnes n'est pas régulier. En particulier si la première colonne
+% contient plusieurs chaines de caractères comme : "GO TO SLEEP",
+% "SHUT DOWN", "WAKE UP".
+% Le programme ne lit pas ces lignes.
+%
+
+clear all
+close all
+clc
+
+% Constantes
+% ----------
+% repIn = 'D:\0-US191\UR-US Clientes\UR182 - LOCEAN\UR182 - CO2 et navires marchands\Données';
+% repOut = 'D:\0-US191\UR-US Clientes\UR182 - LOCEAN\UR182 - CO2 et navires marchands\Données\';
+FileType = '*dat.txt';
+
+% eval( ['cd ''' repIn ' '' ']);
+% eval( ['cd ''' repOut ' '' ']);
+
+% Nombre de colonnes des fichiers à concaténer
+% Appareil type Craig Neil : 34
+% Appareil type GO : 36
+% --------------------------------------------
+NbColCN = 34;
+
+% Entête des fichiers d'origine
+% Appareil Craig Neil
+% -----------------------------
+HeaderCN = ['Type error Date Time GPS time latitude longitude '...
+ 'equ temp Std CO2 raw CO2 um/m H2O raw H2O mm/m '...
+ 'cell temp cell press atm press equ press H2O flow '...
+ 'licor flow equ pump vent flow cond temp atm cond '...
+ 'equ cond drip 1 drip 2 PC temp DB temp'];
+
+% Appareil GO - N'inclut pas les colonnes au-delà de la 34ème
+% ------------------------------------------------------------
+HeaderGO = ['Type error Date PC Time gps time latitude longitude '...
+ 'equ temp std val CO2 mv CO2 um/m H2O mv H2O mm/m ',...
+ 'licor temp licor press atm press equ press H2O flow '...
+ 'licor flow equ pump vent flow atm cond '...
+ 'equ cond drip 1 drip 2 cond temp dry box temp deck box temp'];
+
+% Appareil GO -type fichier 2 - N'inclut pas les colonnes au-delà de la 34ème
+% ------------------------------------------------------------
+HeaderGO2 = ['Type error PC Date PC Time GPS date gps time latitude '...
+ 'longitude equ temp std val CO2a W CO2b W CO2 um/m H2Oa W H2Ob W '...
+ 'H2O mm/m licor temp licor press atm press equ press H2O flow '...
+ 'licor flow equ pump vent flow atm cond equ cond drip 1 na '...
+ 'cond temp dry box temp deck box temp'];
+
+% Format de lecture. Lit les entiers avec %f car la conversion
+% d'un tableau de cellule en tableau numérique via "cell2mat" n'est
+% possible que si les données ont le même format.
+% Le format NoGPS rempalce la lecture de l'heure GPS "%f:%f:%f" par "::"
+% car si l'acquistion GPS n'est pas possible aucun digit n'est enregistré
+% Le format In2 considère les fichiers GO de type 2, qui possèdent plus de
+% paramètres
+% Le format NoGPS2 rempalce la lecture de l'heure GPS "%f:%f:%f" par "::"
+% pour les fichiers de type 2 car si l'acquistion GPS n'est pas possible
+% aucun digit n'est enregistré
+% -----------------------------------------------------------------------
+formatIn = ['%s %f %f/%f/%f %f:%f:%f %f:%f:%f '...
+ '%f %f %f %f %f %f %f %f %f %f %f '...
+ '%f %f %f %f %f %f %f %f %f %f %f %f'];
+formatIn2 = ['%s %f %f/%f/%f %f:%f:%f %f/%f/%f '...
+ '%f:%f:%f %f %f %f %f %f %f %f %f %f %f %f '...
+ '%f %f %f %f %f %f %f %f %f %f %f %f %f %f'];
+formatNoGPS = ['%s %f %f/%f/%f %f:%f:%f :: '...
+ '%f %f %f %f %f %f %f %f %f %f %f '...
+ '%f %f %f %f %f %f %f %f %f %f %f %f'];
+formatNoGPS2 = ['%s %f %f/%f/%f %f:%f:%f '...
+ '%f %f %f %f %f %f %f %f %f %f %f '...
+ '%f %f %f %f %f %f %f %f %f %f %f %f'];
+% Entête en écriture
+% Séparation des colonnes par un ";" compatible avec le format Excel *.CSV
+% les 2 colonnes Date et Time sont rassemblées dans la même colonne
+% pour faciliter le traitement sous excel
+% ----------------------------------------------------------------
+HeaderOut = ['Type;error;Date Time; GPS time; latitude; longitude;'...
+ 'equ temp;Std; CO2 raw; CO2 um/m; H2O raw; H2O mm/m;'...
+ 'licor temp; licor press; atm press; equ press; H2O flow;'...
+ 'licor flow; equ pump;vent flow; cond temp; atm cond;'...
+ 'equ cond; drip 1; drip 2; dry box temp; deck box temp'];
+
+HeaderOut = ['DATE_TIME;GPS_TIME;TYPE;ERROR;LATX;LONX;'...
+ 'EQU_T;STD;CO2_RAW;CO2_PHYS;H2O_RAW;H2O_PHYS;'...
+ 'LICOR_T;LICOR_P;ATM_P;EQU_P;H2O_FLOW;LICOR_FLOW;'...
+ 'EQU_PUMP;VENT_FLOW;COND_T;COND_ATM;COND_EQU;'...
+ 'DRIP_1;DRIP_2;DRY_BOX_T;DECK_BOX_T'];
+
+% Format en écriture. Colonnes séparées par des ";" pour une
+% lecture sous Excel
+% Les colonnes Date et Time sont concaténées (placées entre 2 ";")
+% ---------------------------------------------------------------
+formatOut = ['%02d/%02d/%4d %02d:%02d:%02d;%02d:%02d:%02d;'...
+ '%s;%d;%.4f;%.4f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;'...
+ '%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f\n'];
+
+% Sélection des fichiers à concaténer via une fenêtre utilisateur
+% Le programme ne fonctionne pas si un seul fichier est sélectionné
+% La multiselection de garantit par l'ordre des fichiers
+% -----------------------------------------------------------------
+[FileIn, PathIn] = uigetfile( FileType, 'Fichiers en lecture', 'MultiSelect','on');
+
+if ~isempty( PathIn )
+
+ % Initialisation du tableau de cellules qui contient les données concaténées
+ % --------------------------------------------------------------------------
+ data = {};
+
+ % Boucle sur les fichiers à concaténer
+ % ------------------------------------
+ length = size(FileIn);
+ for ifile = 1:length(2)
+
+ % Ouverture d'un fichier
+ % ----------------------
+ fileIn = char([PathIn FileIn{ifile}]);
+ fid = fopen( fileIn, 'r' );
+
+ if fid ~= -1
+
+ % Variable qui permet de tester la date au format 2011 ou 11
+ % fmtDate = 0 si format 11
+ % fmtDate = 2000 si format 2011
+ % ----------------------------------------------------------
+ fmtDate = 0;
+
+ disp(' ');
+ disp(char(FileIn{ifile}))
+
+ % Lecture de l'entete et comparaison pour déterminer l'origine du
+ % fichier
+ % L'entete détermline le type de fichier Craig Neil ou GO
+ % Not case sensitive
+ % ------------------------------------------------------
+ header = fgetl( fid );
+ if strncmpi( header, HeaderCN, size(HeaderCN,2) )
+ disp('Lecture d''un fichier type Craig Neil' );
+ TypeApp = 'CN';
+ elseif strncmpi( header, HeaderGO, size(HeaderGO,2) )
+ disp('Lecture d''un fichier type G.O.' );
+ TypeApp = 'GO';
+ elseif strncmpi( header, HeaderGO2, size(HeaderGO2,2) )
+ disp('Lecture d''un fichier type G.O. (type2)' );
+ TypeApp = 'GO2';
+ formatIn = formatIn2; % Changement foramt de fichiers car plus de paramètres
+ else
+ disp('Type de fichier inconnu - Pb d''entête');
+ end
+
+ % Lecture ligne à ligne pour détecter les lignes incomplètes
+ % ----------------------------------------------------------
+ nligne = 0;
+ while ~feof(fid)
+ try
+ nligne = nligne + 1;
+ ligne = fgetl( fid );
+
+ % Ne lit que les ligne commençant par les 3 caractères :
+ % EQU - STD - ATM - Test not case sensitive
+ % ------------------------------------------------------
+ temp = {};
+ if strcmpi(ligne(1:3),'EQU') || strcmpi(ligne(1:3),'STD') || strcmpi(ligne(1:3),'ATM')
+ temp = textscan( ligne, formatIn, 1, 'ReturnOnError', 0 );
+ end
+
+ catch
+
+ % Gestion des erreurs de lecture
+ % ------------------------------
+ if strfind(ligne, '::')
+
+ % Si le GPS ne fonctionne pas, l'heure n'est pas enregistrée.
+ % Uniquement la chaine de caractère '::', utilise le second
+ % format d electure (NoGPS) et remplace l'heure GPS par des
+ % NaN.
+ % ------------------------------------------------------------
+ try
+ temp = textscan( ligne, formatNoGPS, 1, 'ReturnOnError', 0 );
+ temp2 = {temp{1:8} NaN NaN NaN NaN NaN NaN temp{9:end}};
+ temp = temp2;
+ catch
+
+ % Si un autre problème de lecture. Emission d'un message
+ % La cellule temp sera vide
+ % ------------------------------------------------------
+ disp(['Problème à la Lecture de la ligne n° : ' num2str(nligne + 1)] );
+ disp( ligne );
+ end
+
+ elseif strcmp( TypeApp, 'GO2' )
+
+ % Pour les fichiers GO de type 2, si le GPS ne fonctionne pas,
+ % l'heure n'est pas enregistrée.
+ % Pas de chaîne de caractère, utilise le second
+ % format d electure (NoGPS2) et remplace l'heure GPS et la postition
+ % par des NaN.
+ % ------------------------------------------------------------
+ try
+ temp = textscan( ligne, formatNoGPS2, 1, 'ReturnOnError', 0 );
+ temp2 = {temp{1:2} NaN NaN NaN NaN NaN NaN temp{3:5} NaN NaN NaN NaN NaN temp{9:end}};
+ temp = temp2;
+ catch
+
+ % Si un autre problème de lecture. Emission d'un message
+ % La cellule temp sera vide
+ % ------------------------------------------------------
+ disp(['Problème à la Lecture de la ligne n° : ' num2str(nligne + 1)] );
+ disp( ligne );
+ end
+
+ else
+ % Si un autre problème de lecture. Emission d'un message
+ % La cellule temp sera vide
+ % ------------------------------------------------------
+ disp(['Problème à la Lecture de la ligne n° : ' num2str(nligne + 1)] );
+ disp( ligne );
+ end
+ end
+
+ if ~isempty( temp )
+ % La variable 'cond temp' ne se trouve pas dans la même colonne
+ % pour l'appareil GO. On la change de position
+ % Pour les fichiers de type GO2, réarrangement des colonnes
+ % -------------------------------------------------------------
+ if strcmp( TypeApp, 'GO' )
+ temp = {temp{1:27} temp{32} temp{28:31} temp{33:34}};
+ elseif strcmp( TypeApp, 'GO2' )
+ temp = {temp{1:2} temp{9:14} temp{6:8} temp{15:19} temp{21:22} temp{24:31} temp{37} temp{32:36} temp{38:39}};
+ end
+
+ % Remplace les cellules vides par des Nan
+ % ---------------------------------------
+ for i = 2:NbColCN
+ if isempty( temp{i} )
+ temp{i} = NaN;
+ end
+ end
+
+ % Ajoute 2000 à la date
+ % ---------------------
+ if (2000-temp{5}) < 0
+ fmtDate = 2000;
+ else
+ temp{5} = temp{5} + 2000;
+ end
+
+ % Concatène les tableaux de cellules
+ % ----------------------------------
+ if isempty(data)
+ data = temp;
+ else
+ for i = 1:NbColCN
+ data{i} = [data{i}; temp{i}];
+ end
+ end
+ end
+ end
+ end
+
+ % Test sur le format de la date pour affichage d'un avertissement
+ % ---------------------------------------------------------------
+ if fmtDate == 2000
+ disp('Attention date au format 2000' );
+ end
+
+ fclose( fid );
+ end
+end
+
+% Préparation des tableaux pour l'écriture par fprintf.
+% Je n'ai pas réussi à écrire un tableau de cellules contenant
+% des chaines de caractères et des flottants.
+% Sépare le tableau de celllules en un vecteur de chaine de
+% caractères "A" et un tableau de nombres "B"
+% ------------------------------------------------------------
+% Convertit la cellule de chaine de caractères en un vecteur
+A(:,:)= char(data{1}(:));
+% Remplace la cellule 1 (chaine de caractères) par la cellule 2
+% (flottants)
+data{1} = data{2};
+% conversion du tableau de cellules en un tableau de doubles.
+B = cell2mat(data);
+% élimine la première colonne dans laquelle se trouvait à l'origine
+% les chaines de caratères.
+B(:,1) = [];
+
+% Teste que les dates sont bien croissantes
+% Conversion des dates heures en une matrice de nombres :
+% datenum accepts only double array
+% --------------------------------------------------------------
+date = double( datenum(B(:,4), B(:,3), B(:,2), B(:,5), B(:,6), B(:,7)) );
+
+if any( diff( date) <= 0 )
+
+ [sortDate, IX] = sort(date);
+ A(:,:) = A(IX,:);
+ B = B(IX,:);
+
+ disp(' ');
+ disp('Les dates ne sont pas strictement croissantes');
+ disp('Le fichier final a été trié');
+end
+
+% Nom du fichier en écriture et ouverture
+% ---------------------------------------
+[FileOut,PathOut] = uiputfile('*.csv','Fichier en écriture');
+fidOut = fopen( [PathOut FileOut], 'w' );
+
+if fidOut ~= -1
+
+ % Ecriture entete
+ % ---------------
+ fprintf( fidOut,'%s\n', HeaderOut);
+
+ % Ecriture données
+ % ----------------
+ [m,n] = size(B);
+ for i=1:m
+ fprintf( fidOut, formatOut, B(i,[2:10]),A(i,:),B(i,[1,11:end]));
+ end
+
+ fclose( fidOut );
+
+ disp(' ');
+ disp ('Ecriture fichier ok' );
+end
+end
+
diff --git "a/tsg_tools/CO2/donn\303\251es oxyg\303\250ne \303\240 int\303\251grer aux donn\303\251es CO2.txt" "b/tsg_tools/CO2/donn\303\251es oxyg\303\250ne \303\240 int\303\251grer aux donn\303\251es CO2.txt"
new file mode 100644
index 0000000000000000000000000000000000000000..2c2abb5ba3823b7f75417e264634bb422248b4a2
--- /dev/null
+++ "b/tsg_tools/CO2/donn\303\251es oxyg\303\250ne \303\240 int\303\251grer aux donn\303\251es CO2.txt"
@@ -0,0 +1,667 @@
+
+Sujet :
+données oxygène à intégrer aux données CO2
+De :
+Denis DIVERRES <denis.diverres@ird.fr>
+Date :
+13/06/2019 Ã 16:27
+Pour :
+Jacques Grelet <Jacques.Grelet@ird.fr>
+
+Salut Jacques
+
+
+Si tu as un peu de temps libre, pourrais-tu regarder comment concaténer les fichiers O2 et CO2? Cette manip concerne le Cap San Lorenzo uniquement: mesures O2 + TSG + CO2.
+
+Le but est d'avoir les données O2 dans le même fichier que celles du CO2.
+
+Les fichiers O2 sont acquis indépendamment par le logiciel oxyview qui tourne sur un PC portable en même temps que SODA.
+
+
+Jusqu'à présent, grâce à un programme matlab écrit par Yves (interpTSG_CO2 voir ci joint), le fichier TSG au format .tsgqc et le fichier CO2 au format csv étaient transformés en un seul fichier csv.
+
+L'idée est donc d'ajouter les données O2 à ce fichier, sans doute en faisant une interpolation. Il y a une difficulté supplémentaire: les données O2 doivent être corrigées de la salinité.
+
+J'imagine qu'on doit d'abord interpoler les données pour avoir un fichier unique (TSG+CO2+O2), puis corriger la valeur O2 brute en utilisant la salinité correspondante.
+
+
+J'ai mis des exemples de fichiers sur ftp://ftp.ifremer.fr/ifremer/ird/diverres/O2/
+
+cslo1902.tsgqc est le fichier TSG initial.
+
+cslo1902i.csv est le fichier CO2 initial.
+
+cslo1902.csv est le fichier TSG + CO2 final obtenu par interp TSG_CO2 Ã partir de CO2 cslo1902i.csv et de cslo1902.tsgqc
+
+cslo1902.oxy est le fichier O2 initial
+
+2019-04-01Â Â Â 07:07:36Â Â Â MEASUREMENTÂ Â Â 3835Â Â Â 1939 Oxygen:Â Â Â 310.33Â Â Â Saturation:Â Â Â 88.12Â Â Â Temperature: 10.01Â Â Â DPhase:Â Â Â 36.21Â Â Â BPhase:Â Â Â 32.95Â Â Â RPhase: 0.00Â Â Â BAmp:Â Â Â 241.06Â Â Â BPot:Â Â Â 78.00Â Â Â RAmp:Â Â Â 0.00 RawTem.:Â Â Â 480.68
+2019-04-01Â Â Â 07:07:37Â Â Â MEASUREMENTÂ Â Â 3835Â Â Â 1939 Oxygen:Â Â Â 310.36Â Â Â Saturation:Â Â Â 88.13Â Â Â Temperature: 10.00Â Â Â DPhase:Â Â Â 36.21Â Â Â BPhase:Â Â Â 32.95Â Â Â RPhase: 0.00Â Â Â BAmp:Â Â Â 241.27Â Â Â BPot:Â Â Â 78.00Â Â Â RAmp:Â Â Â 0.00 RawTem.:Â Â Â 480.76
+2019-04-01Â Â Â 07:08:33Â Â Â MEASUREMENTÂ Â Â 3835Â Â Â 1939 Oxygen:Â Â Â 310.65Â Â Â Saturation:Â Â Â 88.25Â Â Â Temperature: 10.02Â Â Â DPhase:Â Â Â 36.18Â Â Â BPhase:Â Â Â 32.92Â Â Â RPhase: 0.00Â Â Â BAmp:Â Â Â 241.01Â Â Â BPot:Â Â Â 78.00Â Â Â RAmp:Â Â Â 0.00 RawTem.:Â Â Â 480.19
+2019-04-01Â Â Â 07:12:42Â Â Â MEASUREMENTÂ Â Â 3835Â Â Â 1939 Oxygen:Â Â Â 309.90Â Â Â Saturation:Â Â Â 87.99Â Â Â Temperature: 10.00Â Â Â DPhase:Â Â Â 36.23Â Â Â BPhase:Â Â Â 32.97Â Â Â RPhase: 0.00Â Â Â BAmp:Â Â Â 238.13Â Â Â BPot:Â Â Â 78.00Â Â Â RAmp:Â Â Â 0.00 RawTem.:Â Â Â 480.78
+
+Les champs qui nous intéressent sont la date, l'heure, Oxygen (310.33 par exemple) et la temperature (10.01).
+
+
+Les calculs de correction de l'O2 brut sont détaillés page 31-33 de la doc Aanderaa "Oxygen Optode Manual".
+
+Dans notre cas, seule la salinité et la température influent sur la correction, la pression n'a pas d'incidence (pression nulle). Ci-joint un programme excel pour le calcul de la correction si tu veux faire des tests.
+
+
+Voilà , j'espère que ce n'est pas trop compliqué. Dis-moi s'il manque quelque chose.
+
+Bonne mission!
+
+Denis
+
+
+
+
+
+--Â
+Denis Diverres
+
+Unité de Service 191 IMAGO
+"Instrumentation, Moyens Analytiques, Observatoires en Geophysique et
+Oceanographie"
+
+IRD Centre de Bretagne
+Technopole Ifremer
+BP 70
+29 280 Plouzane
+FRANCE
+
+Email: Denis.Diverres@ird.fr
+Tel: +33 (0)2 98 22 46 55
+
+
+
+
+Pour nous remonter une erreur de filtrage, veuillez vous rendre ici : http://f.security-mail.net/3wJGIQZ7
+
+interpTSG_CO2.m
+
+function interpTSG_CO2
+%
+% Interpolation des mesures TSG Ã la date des mesures CO2.
+% Par rapport au fichier *.csv qui est lu, les variables suivantes sont
+% ajoutées au fichier final :
+% SSPS salinité interpolée
+% SSPS_QC code qualité de la salinité
+% SSJT température de cuve du TSG
+% SSJT_QC code qualité de la température de cuve
+% SSJT_COR température de cuve corrigée (-9 si pas de correction)
+% EQU_T_COR température équilibrateur corrigée (-9 si pas de correction)
+%
+% En entrée :
+% 1 - Fichier au format *.csv concaténé par le programme "concatCO2"
+% Le fichier est ouvert via la fonction "readConcatCO2"
+% 2 - Fichier ASCII de mesures TSG
+% Le fichier est ouvert via la fonction "readAsciiTsgCO2"
+%
+% Fonctions externes appelées :
+% readConcatCO2, readAsciiTsgCO2
+%
+% Fonctions internes en fin de ce fichier :
+% selectTSG, interp
+
+clear all
+close all
+clc
+
+% Lecture du fichier CO2
+% ----------------------
+[co2, co2File, error] = readConcatCO2;
+
+if error
+
+ % Lecture du fichier .ini
+ % -----------------------
+ [ini, error1] = readIniCo2(co2File);
+
+ if error1
+
+ % Lecture du fichier TSG
+ % ----------------------
+ [tsg, error2] = readAsciiTsgCO2;
+
+ if error2
+
+ figure( 'Name', 'PlotYY' )
+ [AX,H1,H2] = plotyy(tsg.DAYD, tsg.SSPS, co2.DAYD, co2.CO2_PHYS );
+ datetick(AX(1), 'x', 'mmm-dd', 'keeplimits', 'keepticks')
+ datetick(AX(2), 'x', 'mmm-dd', 'keeplimits', 'keepticks')
+ set(H1, 'LineStyle', 'None', 'Marker', '.', 'MarkerEdgeColor', 'k');
+ set(H2, 'LineStyle', 'None', 'Marker', '.', 'MarkerEdgeColor', 'r');
+ legend('SSPS','CO2 PHYS', 'Location','SouthWest')
+
+ % Selection des mesures TSG. On sélectionne dans l'ordre :
+ % 1 - les données corrigées (ADJUSTED) - Code qualité 5
+ % 2 - les données non corrigées mais ayant un code qualité de 1 (GOOD) ou
+ % 2 (PROBABLY GOOD)
+ % 3 - les autres mesures de la série TSG sont remplacées par 35 pour la
+ % salinité (SSPS) et -9 pour la température de cuve (SSJT). Dans ce cas
+ % le code qualité est positionné à 9 (MISSING VALUE).
+ %
+ % La fonction selectTSG se trouve en fin de fichier
+ % -------------------------------------------------
+ [tsg] = selectTS(tsg);
+
+ % -------------------------------------------------------------------
+ % INTERPOLATION
+ % Attention : l'ordre dans lequel des fonctions d'interpolation est
+ % important
+
+ % Interpolation des données TSG et traitement des cas particuliers
+ % ----------------------------------------------------------------
+ [co2] = interp(co2, tsg);
+
+ % Interpolation des positions
+ % ---------------------------
+ [co2] = interp_POS(co2, tsg);
+
+ % si un fichier .ini existe, effecuter les corrections et modifier le
+ % format d'écriture
+ % -------------------------------------------------------------------
+ co2.SSJT_COR = -999 * ones(size(co2.SSJT));
+ co2.EQU_T_COR = -999 * ones(size(co2.EQU_T));
+% if ~isempty( ini )
+%
+% % Correction temperature SSJT
+% % ---------------------------
+% if ini.dT_SSJT ~= 0
+% co2.SSJT_COR = co2.SSJT;
+% ind = find( co2.SSJT > -900 );
+% co2.SSJT_COR(ind) = co2.SSJT(ind) + ini.dT_SSJT;
+% % HeaderOut = [HeaderOut ';SSJT_COR'];
+% % formatOut = [formatOut ';%.2f'];
+% end
+%
+% % Correction temperature équilibrateur
+% % ------------------------------------
+% if ini.A_EQU_T ~= 1 || ini.B_EQU_T ~= 0
+% co2.EQU_T_COR = ini.A_EQU_T * co2.EQU_T + ini.B_EQU_T;
+% % HeaderOut = [HeaderOut ';EQU_T_COR'];
+% % formatOut = [formatOut ';%.2f'];
+% end
+% end
+
+ % Compléter le format d'écriture par un retour chariot
+ % ----------------------------------------------------
+% formatOut = [formatOut '\n'];
+
+ % -------------------------------------------------------------------
+ % TRACE
+ % -------------------------------------------------------------------
+ figure( 'Name', 'Latitude');
+ hold on;
+ plot( co2.DAYD, co2.LATX, 'b.');
+ plot( co2.DAYD, co2.LATX_INT, 'r.');
+ datetick('x', 'mmm-dd', 'keeplimits', 'keepticks');
+ legend('LATX','LATX_INT Interp', 'Location','SouthWest');
+
+ figure( 'Name', 'Longitude');
+ hold on;
+ plot( co2.DAYD, co2.LONX, 'b.');
+ plot( co2.DAYD, co2.LONX_INT, 'r.');
+ datetick('x', 'mmm-dd', 'keeplimits', 'keepticks');
+ legend('LONX','LONX_INT Interp', 'Location','SouthWest');
+
+ figure( 'Name', 'SSPS_et_co2_SSPS');
+ % plot( tsg.DAYD, tsg.SSPS_SEL, 'k.');
+ hold on;
+ plot( tsg.DAYD, tsg.SSPS_SEL, 'b.');
+ plot( co2.DAYD, co2.SSPS, 'r.');
+ datetick('x', 'mmm-dd', 'keeplimits', 'keepticks');
+ legend('SSPS','SSPS Interp', 'Location','SouthWest');
+
+ % Sauvegarde la figure au format *.png
+ eval( ['print -dpng ' get(gcf, 'Name')]);
+
+ figure( 'Name', 'SSPS_QC');
+ plot( co2.DAYD, co2.SSPS_QC, 'k.');
+ datetick('x', 'mmm-dd', 'keeplimits', 'keepticks');
+ legend('SSPS QC', 'Location','SouthWest');
+
+ % Sauvegarde la figure au format *.png
+ eval( ['print -dpng ' get(gcf, 'Name')]);
+
+ figure( 'Name', 'SSJT');
+ % plot( tsg.DAYD, tsg.SSJT_SEL, 'k.');
+ hold on;
+ plot( tsg.DAYD, tsg.SSJT_SEL, 'b.');
+ plot( co2.DAYD, co2.SSJT, 'r.');
+ datetick('x', 'mmm-dd', 'keeplimits', 'keepticks');
+ legend('SSJT','SSJT Interp', 'Location','SouthWest');
+
+ % Sauvegarde la figure au format *.png
+ eval( ['print -dpng ' get(gcf, 'Name')]);
+
+ figure( 'Name', 'SSJT_QC');
+ plot( co2.DAYD, co2.SSJT_QC, 'k.');
+ datetick('x', 'mmm-dd', 'keeplimits', 'keepticks');
+ legend('SSJT QC', 'Location','SouthWest');
+
+ % Sauvegarde la figure au format *.png
+ eval( ['print -dpng ' get(gcf, 'Name')]);
+
+ % Ecriture
+ % --------
+ HeaderOut = ['DATE_TIME;GPS_TIME;TYPE;ERROR;LATX;LONX;'...
+ 'LATX_INT;LONX_INT;'...
+ 'EQU_T;STD;CO2_RAW;CO2_PHYS;H2O_RAW;H2O_PHYS;'...
+ 'LICOR_T;LICOR_P;ATM_P;EQU_P;H2O_FLOW;LICOR_FLOW;'...
+ 'EQU_PUMP;VENT_FLOW;COND_T;COND_ATM;COND_EQU;'...
+ 'DRIP_1;DRIP_2;DRY_BOX_T;DECK_BOX_T;'...
+ 'SSPS;SSPS_QC;SSJT;SSJT_QC;SSJT_COR;EQU_T_COR;'...
+ ];
+
+ % Ordre des champs de la structure co2. Les champs seront réordonnés
+ % dans l'ordre ci-dessous pour la préparation de l'écriture dans
+ % un fichier
+ % -------------------------------------------------------------------
+ StrucOrder = {'DAYD';'DATE';'GPS_TIME';'TYPE';'ERROR';'LATX';'LONX';...
+ 'LATX_INT';'LONX_INT';...
+ 'EQU_T';'STD';'CO2_RAW';'CO2_PHYS';'H2O_RAW';'H2O_PHYS';...
+ 'LICOR_T';'LICOR_P';'ATM_P';'EQU_P';'H2O_FLOW';'LICOR_FLOW';...
+ 'EQU_PUMP';'VENT_FLOW';'COND_T';'COND_ATM';'COND_EQU';...
+ 'DRIP_1';'DRIP_2';'DRY_BOX_T';'DECK_BOX_T';...
+ 'SSPS';'SSPS_QC';'SSJT';'SSJT_QC';'SSJT_COR';'EQU_T_COR';...
+ };
+
+ % Format en écriture. Colonnes séparées par des ";" pour une
+ % lecture sous Excel
+ % Les colonnes Date et Time sont concaténées (placées entre 2 ";")
+ % ---------------------------------------------------------------
+ formatOut = ['%02d/%02d/%4d %02d:%02d:%02d;'...
+ '%s; %s; %d;%.4f;%.4f;'...
+ '%.4f;%.4f;'...
+ '%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;'...
+ '%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;%.2f;'...
+ '%.3f;%.0f;%.3f;%.0f;%.2f;%.2f;'...
+ '\n'];
+
+ % Nom du fichier en écriture et ouverture
+ % ---------------------------------------
+ [FileOut,PathOut] = uiputfile('*.csv','Fichier en écriture');
+ fidOut = fopen( [PathOut FileOut], 'w' );
+
+ % Ecriture entete
+ % ---------------
+ if fidOut == -1
+ msg_error = ['Open file error : ' FileOut];
+ warndlg( msg_error, 'ASCII error dialog');
+ else
+
+ % Réordonne les champs
+ % --------------------
+ co2b = orderfields(co2, StrucOrder);
+
+ % conversion de la structure en cellule de cellule
+ % ------------------------------------------------
+ c = struct2cell(co2b);
+
+ % Nombre de cellules dans c : correspond au nombre de colonnes
+ % ------------------------------------------------------------
+ n = size(c,1);
+
+ % Nombre de ligne d'une cellule
+ % ----------------------------
+ m = size(c{1},1);
+
+ % Extraction des entiers et reformatage de la matrice
+ % ---------------------------------------------------
+ x = cell2mat(c(5:n));
+ y = reshape(x,m,n-4);
+
+ % Conversion de la cellule contenant la date et l'heure
+ % -----------------------------------------------------
+ [Y, M, D, H, MN, S] = datevec(c{1});
+ date = [D,M,Y,H,MN,S];
+
+ % ecriture
+ % --------
+ fprintf( fidOut,'%s\n', HeaderOut);
+
+ for i=1:m
+ fprintf( fidOut, formatOut, date(i,:), co2.GPS_TIME(i,:), co2.TYPE(i,:),y(i,:));
+ end
+ fclose(fidOut);
+ end
+
+ end % if error2
+ end % if error1
+end % if error
+end
+
+%% Fonction selectTS
+%
+function [tsg] = selectTS(tsg)
+
+% Le code suivant permet :
+% 1 - d'initialiser la série temporelle de salinité (SSPS) à 35 et le code
+% qualité à 9 (missing value)
+% 2 - de sélectionner la bonne variable.
+% Soit tsg.SSPS avec un code qualité 1 (GOOD) ou 2 (PROBABLY GOOD),
+% Soit tsg.SSPS_ADJUSTED avec un code 5 (VALUE CHANGED)
+% ---------------------------------------------------------------------
+
+
+probablygoodQC = 2;
+nocontrolQC = 0;
+adjustedQC = 5;
+harbourQC = 6;
+nomeasureQC = 9;
+
+PARA = {'SSPS', 'SSJT'};
+
+for par = PARA
+
+ para = char(par);
+
+ % intialisation des valeurs par défaut
+ % ------------------------------------
+ switch para
+ case 'SSPS'
+ defaultPARA = 35;
+ case 'SSJT'
+ defaultPARA = -999;
+ otherwise
+ disp('Nom de variable inconnu');
+ return;
+ end
+
+ % initialisation des variables intermédiaires
+ % -------------------------------------------
+ A = defaultPARA * ones(size(tsg.(para)));
+ A_QC = nomeasureQC * ones(size(tsg.(para)));
+
+ % Vérifie si la variable non ADJUSTED est vide.
+ % On sélectionne les données brutes ayant un code qualité 1 ou 2
+ % Par commodité on positionne le code qualité à 2 (PROBABLY GOOD)
+ % --------------------------------------------------------------
+ ind = find( isnan( tsg.(para) ) == 0 & (tsg.([para '_QC']) < probablygoodQC+1 | tsg.([para '_QC']) == harbourQC));
+ if ~isempty( ind )
+ A( ind ) = tsg.(para)( ind );
+ A_QC( ind ) = probablygoodQC;
+
+ % Si les mesures n'ont pas été controlées/validées, leur QC (0)est
+ % conservé
+ % ----------------------------------------------------------------
+ ind = find( tsg.([para '_QC']) == nocontrolQC);
+ if ~isempty(ind)
+ A_QC( ind ) = nocontrolQC;
+ end
+
+ end
+
+ % Vérifie si la variable corrigée (ADJUSTED) n'est pas vide
+ % Le code QC est forcé à 5
+ % Si des données ADJUSTED existent, elles écrasent les mesures brutes
+ % -------------------------------------------------------------------
+ ind = find( isnan( tsg.([para '_ADJUSTED'])) == 0);
+ if ~isempty( ind )
+ A( ind ) = tsg.([para '_ADJUSTED'])( ind );
+ A_QC( ind ) = adjustedQC ;
+ end
+
+ tsg.([para '_SEL']) = A;
+ tsg.([para '_QC_SEL']) = A_QC;
+
+end
+
+% Test de tracé
+% -------------
+figure( 'Name', 'selectTSG' )
+%plot( tsg.DAYD,tsg.SSPS, 'k.')
+hold on
+%plot( tsg.DAYD,tsg.SSPS_SEL, 'm*')
+legendTxt = '';
+ind = find(tsg.SSPS_QC_SEL < 10);
+if ~isempty( ind )
+ plot( tsg.DAYD(ind), tsg.SSPS_SEL(ind), 'r*');
+ legendTxt{1,1} = 'SSPS NoMeasure';
+end
+ind = find(tsg.SSPS_QC_SEL < 6);
+if ~isempty( ind )
+ plot( tsg.DAYD(ind), tsg.SSPS_SEL(ind), 'g*');
+ legendTxt{1,2} = 'SSPS ADJUSTED';
+end
+ind = find(tsg.SSPS_QC_SEL < 3);
+if ~isempty( ind )
+ plot( tsg.DAYD(ind), tsg.SSPS_SEL(ind), 'b*');
+ legendTxt{1,3} = 'SSPS QC = 1,2 ou 6';
+end
+ind = find(tsg.SSPS_QC_SEL == 0);
+if ~isempty( ind )
+ plot( tsg.DAYD(ind), tsg.SSPS_SEL(ind), 'b*');
+ legendTxt{1,4} = 'SSPS NoControl';
+end
+datetick('x', 'mmm-dd', 'keeplimits', 'keepticks')
+title( 'Mesures SSPS selectionnees');
+legend( legendTxt, 'Location','SouthWest');
+% legend('SSPS NoMeasure', 'SSPS ADJUSTED', 'SSPS QC = 1 ou 2' , 'SSPS NoControl', 'Location','SouthWest');
+
+% Sauvegarde la figure au format *.png
+eval( ['print -dpng ' get(gcf, 'name')]);
+
+end
+
+%% Function interpolation
+%
+function [co2] = interp(co2, tsg)
+%
+% Interpolation et traitement des cas particuliers
+% -------------------------------------------------
+
+% Initialisation
+% --------------
+dTmax = 2;
+nocontrolQC = 0;
+probablygoodQC = 2;
+adjustedQC = 5;
+nomeasureQC = 9;
+
+% L'interpolation ne devrait pas être possible lorsqu'il y a des
+% données manquantes sur des périodes supérieures à dT heures.
+% --------------------------------------------------------------
+dT = datenum(0,0,0,dTmax,0,0);
+
+PARA = {'SSPS', 'SSJT'};
+
+for par = PARA
+
+ para = char(par);
+
+ % intialisation des valeurs par défaut
+ % ------------------------------------
+ switch para
+ case 'SSPS'
+ defaultPARA = 35;
+ case 'SSJT'
+ defaultPARA = -999;
+ otherwise
+ disp('Nom de variable inconnu');
+ return;
+ end
+
+ co2.(para) = defaultPARA * ones(size(co2.DAYD));
+ co2.([para '_QC']) = nomeasureQC * ones(size(co2.DAYD));
+
+ % Test les limites temporelles min et max pour lesquelles les mesures
+ % du TSG peuvent être interpolées à la position des mesures CO2.
+ % On n'utilise que les bonnes mesures TSG (code qualité < 9)
+ % -------------------------------------------------------------------
+ ind = find( tsg.([para '_QC_SEL']) < nomeasureQC );
+
+ indmin = 1;
+ if co2.DAYD(1) < tsg.DAYD(ind(1))
+ A = find( co2.DAYD > tsg.DAYD(ind(1)) );
+ if ~isempty(A)
+ indmin = A(1);
+ end
+ end
+
+ indmax = size(co2.DAYD,1);
+ if co2.DAYD(end) > tsg.DAYD(ind(end))
+ A = find( co2.DAYD < tsg.DAYD(ind(end)) );
+ if ~isempty(A)
+ indmax = A(end);
+ end
+ end
+
+ % Interpolation des mesures TSG Ã la position des mesures CO2
+ % On n'utilise que les bonnes mesures (code qualité < 9)
+ % ------------------------------------------------------------
+ co2.(para)(indmin:indmax) = interp1(tsg.DAYD(ind), tsg.([para '_SEL'])(ind), co2.DAYD(indmin:indmax));
+ co2.([para '_QC'])(indmin:indmax) = interp1(tsg.DAYD(ind), tsg.([para '_QC_SEL'])(ind), co2.DAYD(indmin:indmax));
+
+ % Code pour remplacer les valeurs interpolées sur des périodes
+ % supérieures à 2h, par les valeurs par défaut
+ % ----------------------------------------------------------
+ if ~isempty( ind )
+
+ % Calcul la différence de temps entre 2 mesures TSG correctes
+ % -----------------------------------------------------------
+ dTime = diff(tsg.DAYD(ind));
+
+ % Recherche les indices pour lesquels la différence de temps
+ % entre 2 mesures TSG correctes est supérieure à dT heures
+ % ---------------------------------------------------------------
+ ind2 = find(dTime > dT);
+
+ if ~isempty(ind2)
+
+ % Boucle sur les périodes où les différences de temps sont
+ % supérieures à dT
+ % --------------------------------------------------------
+ for i = 1:size(ind2,1)
+
+ % Indices des heures de début et fin où les mesures ont été
+ % interpolées sur plus de dT heures
+ % ---------------------------------------------------------
+ i1 = ind(ind2(i));
+ i2 = ind(ind2(i)+1);
+
+ % Recherche les dates des mesures CO2 correspondant Ã
+ % des interpolations TSG supérieures à dT heures.
+ % Remplace les mesures TSG interpolées par les valeurs par défaut
+ % ---------------------------------------------------------------
+ ind3 = find( co2.DAYD > tsg.DAYD(i1) & co2.DAYD < tsg.DAYD(i2));
+ co2.(para)(ind3) = defaultPARA;
+ co2.([para '_QC'])(ind3) = nomeasureQC;
+ end
+ end
+ end
+
+ % Traitement des QC des mesures pour lesquelles il y a eu interpolation
+ % entre des mesures avec QC = probablyGOOD et des mesures corrigées
+ % avec QC = measureQC. On positionne les QC de ces mesures interpolées
+ % Ã measureQC
+ % ---------------------------------------------------------------------
+ ind = find( co2.([para '_QC']) > probablygoodQC & co2.([para '_QC']) <= adjustedQC);
+ if ~isempty(ind)
+ co2.([para '_QC'])(ind) = adjustedQC;
+ end
+
+ % Traitement des QC des mesures pour lesquelles il y a eu interpolation
+ % entre des mesures avec QC = probablyGOOD et des mesures non validées
+ % avec QC = nocontrolQC. On positionne les QC de ces mesures interpolées
+ % Ã nocontrolQC
+ % ---------------------------------------------------------------------
+ ind = find( co2.([para '_QC']) >= nocontrolQC & co2.([para '_QC']) < probablygoodQC);
+ if ~isempty(ind)
+ co2.([para '_QC'])(ind) = nocontrolQC;
+ end
+
+ end
+ end
+%% Function interpolation Position
+%
+function [co2] = interp_POS(co2, tsg)
+%
+% Interpolation et traitement des cas particuliers
+% -------------------------------------------------
+
+PARA = {'LATX', 'LONX'};
+
+for par = PARA
+
+ para = char(par);
+
+ % intialisation des valeurs par défaut
+ % ------------------------------------
+ defaultPARA = NaN;
+
+ co2.([para '_INT']) = defaultPARA * ones(size(co2.DAYD));
+
+ % Test les limites temporelles min et max pour lesquelles les mesures
+ % du TSG peuvent être interpolées à la position des mesures CO2.
+ % -------------------------------------------------------------------
+
+ indmin = 1;
+ if co2.DAYD(1) < tsg.DAYD(1)
+ A = find( co2.DAYD > tsg.DAYD(1) );
+ if ~isempty(A)
+ indmin = A(1);
+ end
+ end
+
+ indmax = size(co2.DAYD,1);
+ if co2.DAYD(end) > tsg.DAYD(end)
+ A = find( co2.DAYD < tsg.DAYD(end) );
+ if ~isempty(A)
+ indmax = A(end);
+ end
+ end
+
+ % Interpolation des positions à la datedes mesures CO2
+ % ------------------------------------------------------------
+ co2.([para '_INT'])(indmin:indmax) = interp1(tsg.DAYD, tsg.(para), co2.DAYD(indmin:indmax));
+
+% % Code pour remplacer les valeurs interpolées sur des périodes
+% % supérieures à 2h, par les valeurs par défaut
+% % ----------------------------------------------------------
+% if ~isempty( ind )
+%
+% % Calcul la différence de temps entre 2 mesures TSG correctes
+% % -----------------------------------------------------------
+% dTime = diff(tsg.DAYD(ind));
+%
+% % Recherche les indices pour lesquels la différence de temps
+% % entre 2 mesures TSG correctes est supérieure à dT heures
+% % ---------------------------------------------------------------
+% ind2 = find(dTime > dT);
+%
+% if ~isempty(ind2)
+%
+% % Boucle sur les périodes où les différences de temps sont
+% % supérieures à dT
+% % --------------------------------------------------------
+% for i = 1:size(ind2,1)
+%
+% % Indices des heures de début et fin où les mesures ont été
+% % interpolées sur plus de dT heures
+% % ---------------------------------------------------------
+% i1 = ind(ind2(i));
+% i2 = ind(ind2(i)+1);
+%
+% % Recherche les dates des mesures CO2 correspondant Ã
+% % des interpolations TSG supérieures à dT heures.
+% % Remplace les mesures TSG interpolées par les valeurs par défaut
+% % ---------------------------------------------------------------
+% ind3 = find( co2.DAYD > tsg.DAYD(i1) & co2.DAYD < tsg.DAYD(i2));
+% co2.(para)(ind3) = defaultPARA;
+% co2.([para '_QC'])(ind3) = nomeasureQC;
+% end
+% end
+% end
+
+
+end
+end
+
+
+Pièces jointesÂ