Newer
Older
Yves Gouriou
committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
function [sample] = tsg_mergesample( hMainFig )
%
% 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
%
% Get bucket data in ASCII or NETCDF format. If the data does not exist
% getappdata return an empty matrix
% ---------------------------------------------------------------------
bucketASCII = getappdata( hMainFig, 'bucketASCII');
bucketNETCDF = getappdata( hMainFig, 'bucketNETCDF');
% 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