Newer
Older
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
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
function [error] = tsg_writeTSGDataNetCDF( hTsgGUI, filename)
% Function to write TSG data in NetCDF file.
%
% Input
% -----
% hTsgGUI ............ Handel to the main user interface
% filename ........... Data filename
%
% Output
% ------
% error ............. 1 : OK
% ............. -1 : an error occured
%
% The data are store using setappdata - Variable name : 'tsg_data'
%
% Caution : replace the fill-value with NaN
% $Id$
% Open the file
% -------------
nc = netcdf(filename, 'clobber');
if isempty(nc)
msg_error = ['TSG_GOSUD file_lecture : Open file error : ' filename];
warndlg( msg_error, 'NetCDF write dialog');
error = -1;
return;
end
% Initialize loading NetCDF file waitbar
% --------------------------------------
wb = waitbar(0, 'Initializing NetCDF file. Please wait...');
% We cannot change title property to 'none' (default set as 'tex') inside
% waitbar function (line 81), see
% http://www.mathworks.co.uk/matlabcentral/newsreader/view_thread/115725
% -----------------------------------------------------------------------
hchild = get(wb,'children');
htitle = get(hchild,'title');
set(htitle,'Interpreter','None');
% Get the data from the application GUI
% -------------------------------------
tsg = getappdata( hTsgGUI, 'tsg_data');
% Get global attributes list from class tsg_nc with file 'tsg_ncattr.csv'
% -----------------------------------------------------------------------
nca = tsg_nc('tsg_ncattr.csv');
nca_keys = keys(nca);
% Get GF3 variables list codes from class tsg_nc with file 'tsg_ncvar.csv'
% -----------------------------------------------------------------------
ncv = tsg_nc('tsg_ncvar.csv');
ncv_keys = keys(ncv);
ncv_fieldNames = get_fieldnames(ncv);
% Create NetCDF file:
% -------------------
% Variable dimensions
% -------------------
nc('DAYD') = numel(tsg.DAYD); % number of TSG measurements
nc('DAYD_WS') = numel(tsg.DAYD_WS); % number of water samples
nc('DAYD_EXT') = numel(tsg.DAYD_EXT); % number of external SSPS or SSTP comparison
nc('NCOEF_CAL') = 5; % number of calibration coefficients
nc('NCOEF_LIN') = 2; % number of linearisation coefficients
% Fixed dimensions
% ----------------
nc('STRING256') = 256;
nc('STRING14') = 14;
nc('STRING4') = 4;
nc('N1') = 1;
% display filename after 'Interpreter','None' initialization to prevent
% display console warning
% ---------------------------------------------------------------------
waitbar(1/50,wb,['Writing file: ' filename ' Please wait...']);
% Create NetCDF global attributes
% -------------------------------
for i=1:numel(nca_keys)
nc.(nca_keys{i}) = tsg.(nca_keys{i});
end
% Create NetCDF variables and attributes
% ---------------------------------------------------
for i=1:numel(ncv_keys)
for j=1:numel(ncv_fieldNames)
fn = ncv_fieldNames{j};
% define dimension
if strcmp(fn,'dimension')
nc{ncv_keys{i}} = ncfloat(get(ncv, ncv_keys{i}, fn));
% define variable atributes
elseif ~isempty(get(ncv, ncv_keys{i}, fn))
if ischar(get(ncv, ncv_keys{i}, fn))
nc{ncv_keys{i}}.(fn) = ncchar(get(ncv, ncv_keys{i}, fn));
else
nc{ncv_keys{i}}.(fn) = ncfloat(get(ncv, ncv_keys{i}, fn));
end
end
end
end
% Create data to NetCDF file:
% ---------------------------
% Write global attributes
% -----------------------
for i=1:numel(nca_keys)
nc.(nca_keys{i}) = tsg.(nca_keys{i});
end
% Convert Matlab julian days (datenum) to 1950 reference
% ------------------------------------------------------
d = strmatch('DAYD',ncv_keys);
for i=1:numel(d)
tsg.(ncv_keys{d(i)}) = datenumToJulian(tsg.(ncv_keys{d(i)}));
end
% Write measurements (variables)
% -------------------------------
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
for i=1:numel(ncv_keys)
% display waitbar
waitbar( i/numel(ncv_keys), wb, ['writing ' ncv_keys{i} ' variable']);
% DATE as two dimension
if strmatch('DATE',ncv_keys(i))
nc{ncv_keys{i}}(:,:) = tsg.(ncv_keys{i});
else
% other one dimension
nc{ncv_keys{i}}(:) = tsg.(ncv_keys{i});
end
end
% Close waitbar
% -------------
close(wb)
% Close NetCDF file
% -----------------
endef(nc)
close(nc)
% Everything OK
% -------------
error = 1;
end