add comments

if no data for <var>_WS and <var>_EXT, dimension DAYD_WS or DAYD_EXT are not defined

tsg_io/writeTSGDataNetCDF.m

@@ -13,7 +13,9 @@ function [error] = writeTSGDataNetCDF( hTsgGUI, filename)
 %
 % The data are store using setappdata - Variable name : 'tsg_data'
 %
-% Caution : replace the fill-value with NaN
+% Actually, with DAYD, DAYD_WS and DAT_EXT dimensions defined: 61 variables
+%                DAYD, DAYD_WS dimensions defined:             57 variables
+%                DAYD only dimension defined:                  53 variables
 %
 % $Id$
 
@@ -52,8 +54,8 @@ tsg = getappdata( hTsgGUI, 'tsg_data');
 nca             = tsg_nc('tsg_ncattr.csv');
 nca_keys        = keys(nca);
 
-% Get GF3 variables list codes from class tsg_nc with file 'tsg_ncvar.csv'
-% -----------------------------------------------------------------------
+% Get 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);
@@ -68,10 +70,31 @@ fprintf('...writing %s : ', filename);
 % Variable dimensions
 % -------------------
 nc('DAYD')      = numel(tsg.DAYD); 
-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
+
+% if dimension is empty (no data), don't create, because create empty
+% dimension quive it unlimited dimension, but you can only have one
+% unlimited dimension in a file
+% --------------------------------------------------------------------
+
+% number of water samples
+% -----------------------
+if ~isempty(tsg.DAYD_WS)
+  nc('DAYD_WS')   = numel(tsg.DAYD_WS);  
+end
+
+% number of external SSPS or SSTP comparison
+% ------------------------------------------
+if ~isempty(tsg.DAYD_EXT)
+  nc('DAYD_EXT')  = numel(tsg.DAYD_EXT); 
+end
+
+% number of calibration coefficients
+% ----------------------------------
+nc('NCOEF_CAL') = 5;       
+
+% number of linearisation coefficients
+% ------------------------------------
+nc('NCOEF_LIN') = 2;                   
 
 % Fixed dimensions
 % ----------------
@@ -99,8 +122,8 @@ for i=1:numel(ncv_keys)
   dimension = get(ncv, variable, ncv_fieldNames{1});
   nctype    = get(ncv, variable, ncv_fieldNames{2});
   
-  % define dimension, value must be a cell
-  % --------------------------------------
+  % create variable dimension, value must be a cell
+  % -----------------------------------------------
   switch nctype
     case 'double'
       nc{variable} = ncdouble(dimension);
@@ -118,13 +141,19 @@ for i=1:numel(ncv_keys)
       nc{variable} = ncfloat(dimension); 
   end
   
-  % jump dimension and nctype field
-  % -------------------------------
+  % create variable attribute, jump dimension and nctype fields
+  % -----------------------------------------------------------
   for j=3:numel(ncv_fieldNames)
     fieldName = ncv_fieldNames{j};
     value = get(ncv, variable, fieldName);
     
+    % if field is empy, it's not a attribute for this netcdf variable
+    % ---------------------------------------------------------------
     if ~isempty(value)
+      
+      % only two type for variable attribute: char or float, could be
+      % modified in the futur ...
+      % -------------------------------------------------------------
       if ischar(value)
         nc{variable}.(fieldName) = ncchar(value);
       else  
@@ -158,30 +187,22 @@ end
 % Write measurements (variables)
 % -------------------------------
 for i=1:numel(ncv_keys)
-  
+
   % get netcdf variable name
   % ------------------------
   variable = ncv_keys{i};
-  
+
   % display waitbar
   % ---------------
   waitbar( i/numel(ncv_keys), wb, ['writing ' variable ' variable']);
-  
-  % DATE as two dimension
-  % don't work under windows
-  % nc{variable}(:) = tsg.(variable);
-  %if strmatch('DATE', variable)
-  %  nc{variable}(:,:) = tsg.(variable);  
-  %else
-    % other one dimension
-    % under Linux, when you write an with empty matrix to an unlimited
-    % variable dimension, nc was corrupted to double object
-    % works well on windows
-    if ~isempty(tsg.(variable))
-      %nc{variable}(1:length(tsg.(variable))) = tsg.(variable); 
-      nc{variable}(:) = tsg.(variable); 
-    end
-  %end
+
+  % don't write value if empty data, netcdf use FillValue_ attribute
+  % in case of unlimited dimension, for example only, use:
+  % nc{variable}(1:length(tsg.(variable))) = tsg.(variable);
+  % ----------------------------------------------------------------
+  if ~isempty(tsg.(variable))
+    nc{variable}(:) = tsg.(variable);
+  end
 end
 
 % Close waitbar