remplace NaN*ones par zeros in corTsgLinear.m

 add comments 

tsg_util/corTsgLinear.m

@@ -67,11 +67,13 @@ PROBABLY_BAD  = tsg.qc.hash.PROBABLY_BAD.code;
 VALUE_CHANGED = tsg.qc.hash.VALUE_CHANGED.code;
 
 % Intialisation
-% 01/09/2009 : intialisation to NaN
-% ---------------------------------
+% 01/09/2009 : intialisation to NaN for real and 0 for byte (QC)
+% BE CAREFUL:
+% netcdf toolbox failed with assertion when we write NaN to ncbyte variable
+% -------------------------------------------------------------------------
 if isempty( tsg.([PARA{1} '_ADJUSTED']) )  
   tsg.([PARA{1} '_ADJUSTED'])       = NaN*ones(size(tsg.(PARA{1})));
-  tsg.([PARA{1} '_ADJUSTED_QC'])    = NaN*ones(size(tsg.([PARA{1} '_QC'])));
+  tsg.([PARA{1} '_ADJUSTED_QC'])    = zeros(size(tsg.([PARA{1} '_QC'])));
   tsg.([PARA{1} '_ADJUSTED_ERROR']) = NaN*ones(size(tsg.(PARA{1}))); 
 end
 

tsg_util/corTsgMedian.m

@@ -63,8 +63,10 @@ PROBABLY_BAD  = tsg.qc.hash.PROBABLY_BAD.code;
 BAD           = tsg.qc.hash.BAD.code;
 VALUE_CHANGED = tsg.qc.hash.VALUE_CHANGED.code;
 
-% intialisation
-% -------------
+% intialisation to NaN for real and 0 for byte (QC)
+% BE CAREFUL:
+% netcdf toolbox failed with assertion when we write NaN to ncbyte variable
+% -------------------------------------------------------------------------
 if isempty( tsg.([PARA{1} '_ADJUSTED']) )  
   tsg.([PARA{1} '_ADJUSTED'])       = NaN*ones(size(tsg.(PARA{1})));
   tsg.([PARA{1} '_ADJUSTED_QC'])    = zeros(size(tsg.([PARA{1} '_QC'])));