Changeset 919 for trunk/src/struct2nc.m

Timestamp:

Jun 26, 2015, 5:14:45 PM (8 years ago)

Author:

sommeria

Message:

netcdf files stored in 32 bits (not double)

File:

• trunk/src/struct2nc.m (modified) (6 diffs)

trunk/src/struct2nc.m

@@ -49 +49 @@
     return
 end
-%check the validity of the input field structure
+
+%% check the validity of the input field structure
 [errormsg,ListDimName,DimValue,VarDimIndex]=check_field_structure(Data);
 if ~isempty(errormsg)
@@ -56 +57 @@
 end
 ListVarName=Data.ListVarName;
-%nc=netcdf.create(flname,'NC_CLOBBER');%,'clobber'); %create the netcdf file with name flname 
+
+%% create the netcdf file with name flname in format NETCDF4
 cmode = netcdf.getConstant('NETCDF4');
 cmode = bitor(cmode, netcdf.getConstant('CLASSIC_MODEL'));
 cmode = bitor(cmode, netcdf.getConstant('CLOBBER'));
-nc = netcdf.create(flname, cmode) 
-
-%write global constants
+nc = netcdf.create(flname, cmode); 
+
+%% write global constants
 if isfield(Data,'ListGlobalAttribute')
     keys=Data.ListGlobalAttribute;
@@ -84 +86 @@
     end
 end
-%create dimensions
+
+%% create the dimensions
 dimid=zeros(1,length(ListDimName));
 for idim=1:length(ListDimName)
@@ -95 +98 @@
     testattr=1;
 end
-varid=zeros(1,length(Data.ListVarName));
+
+
+%% create the variables
+varid=nan(1,length(Data.ListVarName));
 for ivar=1:length(ListVarName)
-    varid(ivar)=netcdf.defVar(nc,ListVarName{ivar},'nc_double',dimid(VarDimIndex{ivar}));%define variable  
-end
- %write variable attributes
+    if isfield(Data,ListVarName{ivar})
+        VarClass=class(Data.(ListVarName{ivar}));
+        VarType='';
+        switch VarClass
+            case {'single','double'}
+                VarType='nc_float'; % store all floating reals as single
+            case {'uint8','int16','uint16','int32','uint32','int64','uint64'}
+                VarType='nc_int';
+            case  'logical'
+                VarType='nc_byte';
+        end
+        if ~isempty(VarType)
+            varid(ivar)=netcdf.defVar(nc,ListVarName{ivar},VarType,dimid(VarDimIndex{ivar}));%define variable
+        end
+    end
+end
+
+%% write variable attributes
 if testattr
     for ivar=1:min(numel(VarAttribute),numel(ListVarName))  
-        if isstruct(VarAttribute{ivar})
+        if isstruct(VarAttribute{ivar}) && ~isnan(varid(ivar))
             attr_names=fields(VarAttribute{ivar});
             for iattr=1:length(attr_names)
@@ -114 +135 @@
 end
 netcdf.endDef(nc); %put in data mode
+
+%% fill the variables with input data
 for ivar=1:length(ListVarName)
-    if isfield(Data,ListVarName{ivar})
+    if ~isnan(varid(ivar))
         VarVal=Data.(ListVarName{ivar}); 
         %varval=values of the current variable 
-%        VarDimIndex=Data.VarDimIndex{ivar}; %indices of the variable dimensions in the list of dimensions
         VarDimName=Data.VarDimName{ivar};
         if ischar(VarDimName)
@@ -126 +148 @@
         testrange=(numel(VarDimName)==1 && strcmp(VarDimName{1},ListVarName{ivar}) && numel(VarVal)==2);% case of a coordinate defined on a regular mesh by the first and last values.
         testline=isequal(length(siz),2) && isequal(siz(1),1)&& isequal(siz(2), DimValue(VarDimIndex{ivar}));%matlab vector
-        testcolumn=isequal(length(siz),2) && isequal(siz(1), DimValue(VarDimIndex{ivar}))&& isequal(siz(2),1);%matlab column vector
-%             if ~testrange && ~testline && ~testcolumn && ~isequal(siz,DimValue(VarDimIndex))
-%                 errormsg=['wrong dimensions declared for ' ListVarName{ivar} ' in struct2nc.m'];
-%                 break
-%             end 
+        %testcolumn=isequal(length(siz),2) && isequal(siz(1), DimValue(VarDimIndex{ivar}))&& isequal(siz(2),1);%matlab column vector
         if testline || testrange
             if testrange