From 4e8a9dc19d895ce1d0af49f13c7f17164497f1ab Mon Sep 17 00:00:00 2001
From: Yves Gouriou <yves.gouriou@ird.fr>
Date: Tue, 27 Nov 2007 17:26:33 +0000
Subject: [PATCH] ajout de dev_mergesample.m. Concaten les fichiers bouteilles
---
tsg_util/dev_mergesample.m | 73 ++++++++++++++++++++++++++++++++++++++
1 file changed, 73 insertions(+)
create mode 100644 tsg_util/dev_mergesample.m
diff --git a/tsg_util/dev_mergesample.m b/tsg_util/dev_mergesample.m
new file mode 100644
index 0000000..ba9869a
--- /dev/null
+++ b/tsg_util/dev_mergesample.m
@@ -0,0 +1,73 @@
+function [sample] = dev_mergesample( bucketNETCDF, bucketASCII)
+%
+% function [sample] = dev_mergesample( bucketNETCDF, bucketASCII)
+%
+% Merge the 2 types of 'bucket' records, if they exists
+%
+% Input
+% bucketNETCDF ........ Struct. bucket data included in the TSG NetCDF file
+% bucketASCII ......... Struct. bucket data store in an ASCII file
+%
+% Ouptut
+% sample .............. Struct. gathering all bucket measurements
+%
+
+% Create an empty structure with no fields
+% ----------------------------------------
+sample = struct('TIME',{} ,'LATITUDE',{} ,'LONGITUDE',{} ,'PSAL',{} , ...
+ 'PSAL_QC',{} ,'PSAL_DIF',{} ,'PSAL_SMOOTH',{} ,'TYPE', {});
+
+% Fill the structure sample with the bucketNETCDF data.
+% This measurements are of TYPE 1
+% ----------------------------------------------------
+if ~isempty(bucketNETCDF)
+
+ % Complete the bucket structure
+ % -----------------------------
+ bucketNETCDF.PSAL_DIF = zeros(size(bucketNETCDF.PSAL));
+ bucketNETCDF.PSAL_SMOOTH = zeros(size(bucketNETCDF.PSAL));
+ bucketNETCDF.TYPE = ones(size(bucketNETCDF.PSAL));
+
+ sample = bucketNETCDF;
+
+end
+
+% Fill the structure sample with the bucketNETCDF data.
+% This measurements are of TYPE 2
+% ----------------------------------------------------
+if ~isempty(bucketASCII)
+
+ % Complete the bucketNETCDF structure
+ % -----------------------------------
+ bucketASCII.PSAL_DIF = zeros(size(bucketASCII.PSAL));
+ bucketASCII.PSAL_SMOOTH = zeros(size(bucketASCII.PSAL));
+ bucketASCII.TYPE = 2 * ones(size(bucketASCII.PSAL));
+
+ % Fill in the sample structure
+ % ----------------------------
+ if isempty( sample )
+ sample = bucketASCII;
+ else
+ sample.TIME = [sample.TIME; bucketASCII.TIME];
+ sample.LATITUDE = [sample.LATITUDE; bucketASCII.LATITUDE];
+ sample.LONGITUDE = [sample.LONGITUDE; bucketASCII.LONGITUDE];
+ sample.PSAL = [sample.PSAL; bucketASCII.PSAL];
+ sample.PSAL_QC = [sample.PSAL_QC; bucketASCII.PSAL_QC];
+ sample.PSAL_DIF = [sample.PSAL_DIF; bucketASCII.PSAL_DIF];
+ sample.PSAL_SMOOTH = [sample.PSAL_SMOOTH; bucketASCII.PSAL_SMOOTH];
+ sample.TYPE = [sample.TYPE; bucketASCII.TYPE];
+ end
+end
+
+% Sort the struct sample - increasing time.
+% -----------------------------------------
+if ~isempty(sample)
+ [sample.TIME, iOrder] = sort(sample.TIME);
+ sample.LATITUDE = sample.LATITUDE(iOrder);
+ sample.LONGITUDE = sample.LONGITUDE(iOrder);
+ sample.PSAL = sample.PSAL(iOrder);
+ sample.PSAL_QC = sample.PSAL_QC(iOrder);
+ sample.PSAL_DIF = sample.PSAL_DIF(iOrder);
+ sample.PSAL_SMOOTH = sample.PSAL_SMOOTH(iOrder);
+ sample.TYPE = sample.TYPE(iOrder);
+end
--
GitLab