[106] | 1 | % 'struct2nc': create a netcdf file from a Matlab structure |
---|
| 2 | %--------------------------------------------------------------------- |
---|
| 3 | % errormsg=struct2nc(flname,Data) |
---|
| 4 | % |
---|
| 5 | % OUPUT: |
---|
| 6 | % errormsg=error message, =[]: default, no error |
---|
| 7 | % |
---|
| 8 | % INPUT: |
---|
| 9 | % flname: name of the netcdf file to create (must end with the extension '.nc') |
---|
| 10 | % Data: structure containing all the information of the netcdf file (or netcdf object) |
---|
| 11 | % with fields: |
---|
| 12 | % (optional) .ListGlobalAttribute: list (cell array of character strings) of the names of the global attributes Att_1, Att_2... |
---|
| 13 | % .Att_1,Att_2...: values of the global attributes |
---|
| 14 | % (requested) .ListVarName: list of the variable names Var_1, Var_2....(cell array of character strings). |
---|
| 15 | % (requested) .VarDimName: list of dimension names for each element of .ListVarName (cell array of string cells) |
---|
| 16 | % (optional) .VarAttribute: cell array of structures of the form .VarAttribute{ivar}.key=value, defining an attribute key name and value for the variable #ivar |
---|
| 17 | % (requested) .Var1, .Var2....: variables (Matlab arrays) with names listed in .ListVarName |
---|
| 18 | |
---|
| 19 | %AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA |
---|
| 20 | % Copyright Joel Sommeria, 2008, LEGI / CNRS-UJF-INPG, sommeria@coriolis-legi.org. |
---|
| 21 | %AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA |
---|
| 22 | % This file is part of the toolbox UVMAT. |
---|
| 23 | % |
---|
| 24 | % UVMAT is free software; you can redistribute it and/or modify |
---|
| 25 | % it under the terms of the GNU General Public License as published by |
---|
| 26 | % the Free Software Foundation; either version 2 of the License, or |
---|
| 27 | % (at your option) any later version. |
---|
| 28 | % |
---|
| 29 | % UVMAT is distributed in the hope that it will be useful, |
---|
| 30 | % but WITHOUT ANY WARRANTY; without even the implied warranty of |
---|
| 31 | % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
---|
| 32 | % GNU General Public License (file UVMAT/COPYING.txt) for more details. |
---|
| 33 | %AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA |
---|
| 34 | |
---|
| 35 | function errormsg=struct2nc(flname,Data) |
---|
| 36 | if ~ischar(flname) |
---|
[246] | 37 | errormsg='invalid input for the netcf file name'; |
---|
[106] | 38 | return |
---|
| 39 | end |
---|
| 40 | if ~exist('Data','var') |
---|
| 41 | errormsg='no data input for the netcdf file'; |
---|
| 42 | return |
---|
| 43 | end |
---|
[404] | 44 | FilePath=fileparts(flname); |
---|
| 45 | if ~strcmp(FilePath,'') && ~exist(FilePath,'dir') |
---|
| 46 | errormsg=['directory ' FilePath ' needs to be created']; |
---|
| 47 | return |
---|
| 48 | end |
---|
[674] | 49 | %check the validity of the input field structure |
---|
[404] | 50 | [errormsg,ListDimName,DimValue,VarDimIndex]=check_field_structure(Data); |
---|
| 51 | if ~isempty(errormsg) |
---|
| 52 | errormsg=['error in struct2nc:invalid input structure_' errormsg]; |
---|
| 53 | return |
---|
| 54 | end |
---|
| 55 | ListVarName=Data.ListVarName; |
---|
| 56 | nc=netcdf.create(flname,'NC_CLOBBER');%,'clobber'); %create the netcdf file with name flname |
---|
| 57 | %write global constants |
---|
| 58 | if isfield(Data,'ListGlobalAttribute') |
---|
| 59 | keys=Data.ListGlobalAttribute; |
---|
| 60 | for iattr=1:length(keys) |
---|
| 61 | if isfield(Data,keys{iattr}) |
---|
| 62 | testvar=0; |
---|
| 63 | for ivar=1:length(ListVarName)% eliminate possible global attributes with the same name as a variable |
---|
| 64 | if isequal(ListVarName{ivar}, keys{iattr}) |
---|
| 65 | testvar=1; |
---|
| 66 | break |
---|
[106] | 67 | end |
---|
[404] | 68 | end |
---|
| 69 | if ~testvar |
---|
| 70 | eval(['cte=Data.' keys{iattr} ';']) |
---|
| 71 | if (ischar(cte) ||isnumeric(cte)) && ~isempty(cte)%&& ~isequal(cte,'') |
---|
| 72 | %write constant only if it is numeric or char string, and not empty |
---|
| 73 | netcdf.putAtt(nc,netcdf.getConstant('NC_GLOBAL'),keys{iattr},cte) |
---|
[106] | 74 | end |
---|
| 75 | end |
---|
| 76 | end |
---|
| 77 | end |
---|
[404] | 78 | end |
---|
| 79 | %create dimensions |
---|
| 80 | dimid=zeros(1,length(ListDimName)); |
---|
| 81 | for idim=1:length(ListDimName) |
---|
| 82 | dimid(idim) = netcdf.defDim(nc,ListDimName{idim},DimValue(idim)); |
---|
| 83 | end |
---|
| 84 | VarAttribute={}; %default |
---|
| 85 | testattr=0; |
---|
| 86 | if isfield(Data,'VarAttribute') |
---|
| 87 | VarAttribute=Data.VarAttribute; |
---|
| 88 | testattr=1; |
---|
| 89 | end |
---|
| 90 | varid=zeros(1,length(Data.ListVarName)); |
---|
| 91 | for ivar=1:length(ListVarName) |
---|
| 92 | varid(ivar)=netcdf.defVar(nc,ListVarName{ivar},'nc_double',dimid(VarDimIndex{ivar}));%define variable |
---|
| 93 | end |
---|
| 94 | %write variable attributes |
---|
| 95 | if testattr |
---|
| 96 | for ivar=1:min(numel(VarAttribute),numel(ListVarName)) |
---|
| 97 | if isstruct(VarAttribute{ivar}) |
---|
| 98 | attr_names=fields(VarAttribute{ivar}); |
---|
| 99 | for iattr=1:length(attr_names) |
---|
[526] | 100 | attr_val=VarAttribute{ivar}.(attr_names{iattr}); |
---|
| 101 | if ~isempty(attr_names{iattr})&& ~isempty(attr_val)&&~iscell(attr_val) |
---|
[404] | 102 | netcdf.putAtt(nc,varid(ivar),attr_names{iattr},attr_val); |
---|
[106] | 103 | end |
---|
| 104 | end |
---|
| 105 | end |
---|
| 106 | end |
---|
[404] | 107 | end |
---|
| 108 | netcdf.endDef(nc); %put in data mode |
---|
| 109 | for ivar=1:length(ListVarName) |
---|
| 110 | if isfield(Data,ListVarName{ivar}) |
---|
| 111 | VarVal=Data.(ListVarName{ivar}); |
---|
| 112 | %varval=values of the current variable |
---|
| 113 | % VarDimIndex=Data.VarDimIndex{ivar}; %indices of the variable dimensions in the list of dimensions |
---|
| 114 | VarDimName=Data.VarDimName{ivar}; |
---|
| 115 | if ischar(VarDimName) |
---|
| 116 | VarDimName={VarDimName}; |
---|
| 117 | end |
---|
| 118 | siz=size(VarVal); |
---|
| 119 | 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. |
---|
| 120 | testline=isequal(length(siz),2) && isequal(siz(1),1)&& isequal(siz(2), DimValue(VarDimIndex{ivar}));%matlab vector |
---|
| 121 | testcolumn=isequal(length(siz),2) && isequal(siz(1), DimValue(VarDimIndex{ivar}))&& isequal(siz(2),1);%matlab column vector |
---|
| 122 | % if ~testrange && ~testline && ~testcolumn && ~isequal(siz,DimValue(VarDimIndex)) |
---|
[382] | 123 | % errormsg=['wrong dimensions declared for ' ListVarName{ivar} ' in struct2nc.m']; |
---|
| 124 | % break |
---|
| 125 | % end |
---|
[404] | 126 | if testline || testrange |
---|
| 127 | if testrange |
---|
| 128 | VarVal=linspace(VarVal(1),VarVal(2),DimValue(VarDimIndex{ivar}));% restitute the whole array of coordinate values |
---|
| 129 | end |
---|
| 130 | netcdf.putVar(nc,varid(ivar), double(VarVal')); |
---|
| 131 | else |
---|
| 132 | netcdf.putVar(nc,varid(ivar), double(VarVal)); |
---|
| 133 | end |
---|
[106] | 134 | end |
---|
| 135 | end |
---|
[404] | 136 | netcdf.close(nc) |
---|
[106] | 137 | |
---|
[404] | 138 | |
---|
[530] | 139 | %'check_field_structure': check the validity of the field struture representation consistant with the netcdf format |
---|
| 140 | %------------------------------------------------------------------------ |
---|
[674] | 141 | % [errormsg,ListDimName,DimValue,VarDimIndex]=check_field_structure(Data) |
---|
[530] | 142 | % |
---|
| 143 | % OUTPUT: |
---|
| 144 | % errormsg: error message which is not empty when the input structure does not have the right form |
---|
[674] | 145 | % ListDimName: list of dimension names (cell array of cahr strings) |
---|
| 146 | % DimValue: list of dimension values (numerical array with the same dimension as ListDimName) |
---|
| 147 | % VarDimIndex: cell array of dimension index (in the list ListDimName) for each element of Data.ListVarName |
---|
[530] | 148 | % |
---|
| 149 | % INPUT: |
---|
| 150 | % Data: structure containing |
---|
| 151 | % (optional) .ListGlobalAttribute: cell listing the names of the global attributes |
---|
| 152 | % .Att_1,Att_2... : values of the global attributes |
---|
| 153 | % (requested) .ListVarName: list of variable names to select (cell array of char strings {'VarName1', 'VarName2',...} ) |
---|
| 154 | % (requested) .VarDimName: list of dimension names for each element of .ListVarName (cell array of string cells) |
---|
| 155 | % (requested) .Var1, .Var2....: variables (Matlab arrays) with names listed in .ListVarName |
---|
| 156 | |
---|
| 157 | |
---|
| 158 | function [errormsg,ListDimName,DimValue,VarDimIndex]=check_field_structure(Data) |
---|
| 159 | errormsg=''; |
---|
| 160 | ListDimName={}; |
---|
| 161 | DimValue=[]; %default |
---|
| 162 | VarDimIndex={}; |
---|
| 163 | if ~isstruct(Data) |
---|
| 164 | errormsg='input field is not a structure'; |
---|
| 165 | return |
---|
| 166 | end |
---|
| 167 | if isfield(Data,'ListVarName') && iscell(Data.ListVarName) |
---|
| 168 | nbfield=numel(Data.ListVarName); |
---|
| 169 | else |
---|
| 170 | errormsg='input field does not contain the cell array of variable names .ListVarName'; |
---|
| 171 | return |
---|
| 172 | end |
---|
| 173 | %check dimension names |
---|
| 174 | if (isfield(Data,'VarDimName') && iscell(Data.VarDimName)) |
---|
| 175 | if numel(Data.VarDimName)~=nbfield |
---|
| 176 | errormsg=' .ListVarName and .VarDimName have different lengths'; |
---|
| 177 | return |
---|
| 178 | end |
---|
| 179 | else |
---|
| 180 | errormsg='input field does not contain the cell array of dimension names .VarDimName'; |
---|
| 181 | return |
---|
| 182 | end |
---|
| 183 | nbdim=0; |
---|
| 184 | ListDimName={}; |
---|
| 185 | |
---|
| 186 | %% main loop on the list of variables |
---|
[674] | 187 | VarDimIndex=cell(1,nbfield); |
---|
[530] | 188 | for ivar=1:nbfield |
---|
| 189 | VarName=Data.ListVarName{ivar}; |
---|
| 190 | if ~isfield(Data,VarName) |
---|
| 191 | errormsg=['the listed variable ' VarName ' is not found']; |
---|
| 192 | return |
---|
| 193 | end |
---|
| 194 | sizvar=size(Data.(VarName));% sizvar = dimension of variable |
---|
| 195 | DimCell=Data.VarDimName{ivar}; |
---|
| 196 | if ischar(DimCell) |
---|
| 197 | DimCell={DimCell};%case of a single dimension name, defined by a string |
---|
| 198 | elseif ~iscell(DimCell) |
---|
| 199 | errormsg=['wrong format for .VarDimName{' num2str(ivar) ' (must be the cell of dimension names of the variable ' VarName]; |
---|
[674] | 200 | return |
---|
[530] | 201 | end |
---|
| 202 | nbcoord=numel(sizvar);%nbre of coordinates for variable named VarName |
---|
| 203 | testrange=0; |
---|
| 204 | if numel(DimCell)==0 |
---|
| 205 | errormsg=['empty declared dimension .VarDimName{' num2str(ivar) '} for ' VarName]; |
---|
| 206 | return |
---|
| 207 | elseif numel(DimCell)==1% one dimension declared |
---|
| 208 | if nbcoord==2 |
---|
| 209 | if sizvar(1)==1 |
---|
| 210 | sizvar(1)=sizvar(2); |
---|
| 211 | elseif sizvar(2)==1 |
---|
| 212 | else |
---|
| 213 | errormsg=['1 dimension declared in .VarDimName{' num2str(ivar) '} inconsistent with the nbre of dimensions =2 of the variable ' VarName]; |
---|
| 214 | return |
---|
| 215 | end |
---|
| 216 | if sizvar(1)==2 && isequal(VarName,DimCell{1}) |
---|
| 217 | testrange=1;% test for a dimension variable representing a range |
---|
| 218 | end |
---|
| 219 | else |
---|
| 220 | errormsg=['1 dimension declared in .VarDimName{' num2str(ivar) '} inconsistent with the nbre of dimensions =' num2str(nbcoord) ' of the variable ' VarName]; |
---|
| 221 | return |
---|
| 222 | end |
---|
| 223 | else |
---|
| 224 | if numel(DimCell)>nbcoord |
---|
| 225 | sizvar(nbcoord+1:numel(DimCell))=1;% case of singleton dimensions (not seen by the function size) |
---|
| 226 | elseif nbcoord > numel(DimCell) |
---|
| 227 | errormsg=['nbre of declared dimensions in .VarDimName{' num2str(ivar) '} smaller than the nbre of dimensions =' num2str(nbcoord) ' of the variable ' VarName]; |
---|
| 228 | return |
---|
| 229 | end |
---|
| 230 | end |
---|
| 231 | DimIndex=[]; |
---|
| 232 | for idim=1:numel(DimCell) %loop on the coordinates of variable #ivar |
---|
| 233 | DimName=DimCell{idim}; |
---|
| 234 | iprev=find(strcmp(DimName,ListDimName),1);%look for dimension name DimName in the current list |
---|
| 235 | if isempty(iprev)% append the dimension name to the current list |
---|
| 236 | nbdim=nbdim+1; |
---|
| 237 | RangeTest(nbdim)=0; %default |
---|
| 238 | if sizvar(idim)==2 && strcmp(DimName,VarName)%case of a coordinate defined by the two end values (regular spacing) |
---|
| 239 | RangeTest(nbdim)=1; %to be updated for a later variable |
---|
| 240 | end |
---|
| 241 | DimValue(nbdim)=sizvar(idim); |
---|
| 242 | ListDimName{nbdim}=DimName; |
---|
| 243 | DimIndex=[DimIndex nbdim]; |
---|
| 244 | else % DimName is detected in the current list of dimension names |
---|
| 245 | if ~isequal(DimValue(iprev),sizvar(idim)) |
---|
| 246 | if isequal(DimValue(iprev),2)&& RangeTest(iprev) % the dimension has been already detected as a range [min max] |
---|
| 247 | DimValue(iprev)=sizvar(idim); %update with actual value |
---|
| 248 | elseif ~testrange |
---|
| 249 | errormsg=['dimension declaration inconsistent with the size =[' num2str(sizvar) '] for ' VarName]; |
---|
| 250 | return |
---|
| 251 | end |
---|
| 252 | end |
---|
| 253 | DimIndex=[DimIndex iprev]; |
---|
| 254 | end |
---|
| 255 | end |
---|
| 256 | VarDimIndex{ivar}=DimIndex; |
---|
[674] | 257 | end |
---|