1 | <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" |
---|
2 | "http://www.w3.org/TR/REC-html40/loose.dtd"> |
---|
3 | <html> |
---|
4 | <head> |
---|
5 | <title>Description of nc2struct</title> |
---|
6 | <meta name="keywords" content="nc2struct"> |
---|
7 | <meta name="description" content="'nc2struct': transform a netcdf file in a corresponding matlab structure"> |
---|
8 | <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1"> |
---|
9 | <meta name="generator" content="m2html © 2003 Guillaume Flandin"> |
---|
10 | <meta name="robots" content="index, follow"> |
---|
11 | <link type="text/css" rel="stylesheet" href="../m2html.css"> |
---|
12 | </head> |
---|
13 | <body> |
---|
14 | <a name="_top"></a> |
---|
15 | <div><a href="../index.html">Home</a> > <a href="index.html">.</a> > nc2struct.m</div> |
---|
16 | |
---|
17 | <!--<table width="100%"><tr><td align="left"><a href="../index.html"><img alt="<" border="0" src="../left.png"> Master index</a></td> |
---|
18 | <td align="right"><a href="index.html">Index for . <img alt=">" border="0" src="../right.png"></a></td></tr></table>--> |
---|
19 | |
---|
20 | <h1>nc2struct |
---|
21 | </h1> |
---|
22 | |
---|
23 | <h2><a name="_name"></a>PURPOSE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
24 | <div class="box"><strong>'nc2struct': transform a netcdf file in a corresponding matlab structure</strong></div> |
---|
25 | |
---|
26 | <h2><a name="_synopsis"></a>SYNOPSIS <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
27 | <div class="box"><strong>function [Data,var_detect,ichoice]=nc2struct(nc,ListVarName) </strong></div> |
---|
28 | |
---|
29 | <h2><a name="_description"></a>DESCRIPTION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
30 | <div class="fragment"><pre class="comment">'nc2struct': transform a netcdf file in a corresponding matlab structure |
---|
31 | ---------------------------------------------------------------------- |
---|
32 | function [Data,var_detect,ichoice]=nc2struct(nc,ListVarName) |
---|
33 | |
---|
34 | OUTPUT: |
---|
35 | Data: structure containing all the information of the netcdf file (or netcdf object) |
---|
36 | with fields: |
---|
37 | .ListGlobalAttribute: cell listing the names of the global attributes |
---|
38 | .Att_1,Att_2... : values of the global attributes |
---|
39 | .ListDimName: cell listing the names of the array dimensions |
---|
40 | .DimValue: array dimension values (Matlab vector with the same length as .ListDimName |
---|
41 | .ListVarName: cell listing the names of the variables |
---|
42 | .VarDimIndex: cell containing the set of dimension indices (in list .ListDimName) for each variable of .ListVarName |
---|
43 | .VarDimName: cell containing a cell of dimension names (in list .ListDimName) for each variable of .ListVarName |
---|
44 | .VarAttribute: cell of structures s containing names and values of variable attributes (s.name=value) for each variable of .ListVarName |
---|
45 | .Var1, .Var2....: variables (Matlab arrays) with names listed in .ListVarName |
---|
46 | var_detect: vector with same length as ListVarName, with 1 for each detected variable and 0 else. |
---|
47 | ichoice: = line |
---|
48 | |
---|
49 | INPUT: |
---|
50 | nc: name of a netcdf file (char string) or netcdf object |
---|
51 | ListVarName: optional list of variable names to select (cell array of char strings {'VarName1', 'VarName2',...} ) |
---|
52 | if ListVarName=[] or {}, no variables is read (only global attributes and lists of vdimensions, variables and attriburtes) |
---|
53 | if ListVarName is absent, or = '*', ALL the variables are read. |
---|
54 | if ListVarName is a cell array with n lines, the set of variables |
---|
55 | will be sought by order of priority in the list, while output names will be set by the first line |
---|
56 | |
---|
57 | AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA |
---|
58 | Copyright Joel Sommeria, 2008, LEGI / CNRS-UJF-INPG, sommeria@coriolis-legi.org. |
---|
59 | AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA |
---|
60 | This file is part of the toolbox UVMAT. |
---|
61 | |
---|
62 | UVMAT is free software; you can redistribute it and/or modify |
---|
63 | it under the terms of the GNU General Public License as published by |
---|
64 | the Free Software Foundation; either version 2 of the License, or |
---|
65 | (at your option) any later version. |
---|
66 | |
---|
67 | UVMAT is distributed in the hope that it will be useful, |
---|
68 | but WITHOUT ANY WARRANTY; without even the implied warranty of |
---|
69 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
---|
70 | GNU General Public License (file UVMAT/COPYING.txt) for more details. |
---|
71 | AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</pre></div> |
---|
72 | |
---|
73 | <!-- crossreference --> |
---|
74 | <h2><a name="_cross"></a>CROSS-REFERENCE INFORMATION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
75 | This function calls: |
---|
76 | <ul style="list-style-image:url(../matlabicon.gif)"> |
---|
77 | <li><a href="nc2struct_toolbox.html" class="code" title="function [Data,var_detect,ichoice]=nc2struct_toolbox(nc,ListVarName)">nc2struct_toolbox</a> 'nc2struct_toolbox': transform a netcdf file in a corresponding matlab structure, USE OLD NETCDF LIBRARY</li></ul> |
---|
78 | This function is called by: |
---|
79 | <ul style="list-style-image:url(../matlabicon.gif)"> |
---|
80 | <li><a href="civ.html" class="code" title="function varargout = civ(varargin)">civ</a> 'civ': function associated with the interface 'civ.fig' for PIV, spline interpolation and stereo PIV (patch)</li><li><a href="civ_3D.html" class="code" title="function varargout = civ_3D(varargin)">civ_3D</a> 'civ_3D': function associated with the interface 'civ_3D.fig' for PIV in volume</li><li><a href="get_field.html" class="code" title="function varargout = get_field(varargin)">get_field</a> 'get_field': display variables and attributes from a Netcdf file, and plot selected fields</li><li><a href="read_civxdata.html" class="code" title="function [Field,VelTypeOut]=read_civxdata(filename,FieldNames,VelType)">read_civxdata</a> 'read_civxdata': reads civx data from netcdf files</li><li><a href="read_get_field.html" class="code" title="function [SubField,errormsg]=read_get_field(hget_field)">read_get_field</a> 'read_get_field': read the list of selected variables from the GUI get_field (SAME AS THE FIRST</li><li><a href="read_get_field_old.html" class="code" title="function [SubField,errormsg]=read_get_field(hget_field)">read_get_field_old</a> 'read_get_field': read the list of selected variables from the GUI get_field (SAME AS THE FIRST</li><li><a href="series.html" class="code" title="function varargout = series(varargin)">series</a> 'series': master function associated to the GUI series.m for analysis field series</li><li><a href="uvmat.html" class="code" title="function varargout = uvmat(varargin)">uvmat</a> 'uvmat': function associated with the GUI 'uvmat.fig' for images and data field visualization</li></ul> |
---|
81 | <!-- crossreference --> |
---|
82 | |
---|
83 | |
---|
84 | <h2><a name="_source"></a>SOURCE CODE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
85 | <div class="fragment"><pre>0001 <span class="comment">%'nc2struct': transform a netcdf file in a corresponding matlab structure</span> |
---|
86 | 0002 <span class="comment">%----------------------------------------------------------------------</span> |
---|
87 | 0003 <span class="comment">% function [Data,var_detect,ichoice]=nc2struct(nc,ListVarName)</span> |
---|
88 | 0004 <span class="comment">%</span> |
---|
89 | 0005 <span class="comment">% OUTPUT:</span> |
---|
90 | 0006 <span class="comment">% Data: structure containing all the information of the netcdf file (or netcdf object)</span> |
---|
91 | 0007 <span class="comment">% with fields:</span> |
---|
92 | 0008 <span class="comment">% .ListGlobalAttribute: cell listing the names of the global attributes</span> |
---|
93 | 0009 <span class="comment">% .Att_1,Att_2... : values of the global attributes</span> |
---|
94 | 0010 <span class="comment">% .ListDimName: cell listing the names of the array dimensions</span> |
---|
95 | 0011 <span class="comment">% .DimValue: array dimension values (Matlab vector with the same length as .ListDimName</span> |
---|
96 | 0012 <span class="comment">% .ListVarName: cell listing the names of the variables</span> |
---|
97 | 0013 <span class="comment">% .VarDimIndex: cell containing the set of dimension indices (in list .ListDimName) for each variable of .ListVarName</span> |
---|
98 | 0014 <span class="comment">% .VarDimName: cell containing a cell of dimension names (in list .ListDimName) for each variable of .ListVarName</span> |
---|
99 | 0015 <span class="comment">% .VarAttribute: cell of structures s containing names and values of variable attributes (s.name=value) for each variable of .ListVarName</span> |
---|
100 | 0016 <span class="comment">% .Var1, .Var2....: variables (Matlab arrays) with names listed in .ListVarName</span> |
---|
101 | 0017 <span class="comment">% var_detect: vector with same length as ListVarName, with 1 for each detected variable and 0 else.</span> |
---|
102 | 0018 <span class="comment">% ichoice: = line</span> |
---|
103 | 0019 <span class="comment">%</span> |
---|
104 | 0020 <span class="comment">%INPUT:</span> |
---|
105 | 0021 <span class="comment">% nc: name of a netcdf file (char string) or netcdf object</span> |
---|
106 | 0022 <span class="comment">% ListVarName: optional list of variable names to select (cell array of char strings {'VarName1', 'VarName2',...} )</span> |
---|
107 | 0023 <span class="comment">% if ListVarName=[] or {}, no variables is read (only global attributes and lists of vdimensions, variables and attriburtes)</span> |
---|
108 | 0024 <span class="comment">% if ListVarName is absent, or = '*', ALL the variables are read.</span> |
---|
109 | 0025 <span class="comment">% if ListVarName is a cell array with n lines, the set of variables</span> |
---|
110 | 0026 <span class="comment">% will be sought by order of priority in the list, while output names will be set by the first line</span> |
---|
111 | 0027 <span class="comment">%</span> |
---|
112 | 0028 <span class="comment">%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</span> |
---|
113 | 0029 <span class="comment">% Copyright Joel Sommeria, 2008, LEGI / CNRS-UJF-INPG, sommeria@coriolis-legi.org.</span> |
---|
114 | 0030 <span class="comment">%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</span> |
---|
115 | 0031 <span class="comment">% This file is part of the toolbox UVMAT.</span> |
---|
116 | 0032 <span class="comment">%</span> |
---|
117 | 0033 <span class="comment">% UVMAT is free software; you can redistribute it and/or modify</span> |
---|
118 | 0034 <span class="comment">% it under the terms of the GNU General Public License as published by</span> |
---|
119 | 0035 <span class="comment">% the Free Software Foundation; either version 2 of the License, or</span> |
---|
120 | 0036 <span class="comment">% (at your option) any later version.</span> |
---|
121 | 0037 <span class="comment">%</span> |
---|
122 | 0038 <span class="comment">% UVMAT is distributed in the hope that it will be useful,</span> |
---|
123 | 0039 <span class="comment">% but WITHOUT ANY WARRANTY; without even the implied warranty of</span> |
---|
124 | 0040 <span class="comment">% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the</span> |
---|
125 | 0041 <span class="comment">% GNU General Public License (file UVMAT/COPYING.txt) for more details.</span> |
---|
126 | 0042 <span class="comment">%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</span> |
---|
127 | 0043 |
---|
128 | 0044 <a name="_sub0" href="#_subfunctions" class="code">function [Data,var_detect,ichoice]=nc2struct(nc,ListVarName)</a> |
---|
129 | 0045 |
---|
130 | 0046 <span class="keyword">if</span> ~exist(<span class="string">'ListVarName'</span>,<span class="string">'var'</span>) |
---|
131 | 0047 ListVarName=<span class="string">'*'</span>; |
---|
132 | 0048 <span class="keyword">end</span> |
---|
133 | 0049 hhh=which(<span class="string">'netcdf.open'</span>);<span class="comment">% look for built-in matlab netcdf library</span> |
---|
134 | 0050 |
---|
135 | 0051 <span class="keyword">if</span> ~isequal(hhh,<span class="string">''</span>) |
---|
136 | 0052 <span class="comment">%default output</span> |
---|
137 | 0053 Data=[]; |
---|
138 | 0054 <span class="comment">% ListIndex=[];</span> |
---|
139 | 0055 var_detect=[]; |
---|
140 | 0056 ichoice=[];<span class="comment">%default</span> |
---|
141 | 0057 <span class="comment">%open the netcdf file for reading</span> |
---|
142 | 0058 <span class="keyword">if</span> ischar(nc) |
---|
143 | 0059 <span class="keyword">if</span> exist(nc,<span class="string">'file'</span>) |
---|
144 | 0060 nc=netcdf.open(nc,<span class="string">'NC_NOWRITE'</span>); |
---|
145 | 0061 testfile=1; |
---|
146 | 0062 <span class="keyword">else</span> |
---|
147 | 0063 Data.Txt=[<span class="string">'ERROR:file '</span> nc <span class="string">' does not exist'</span>]; |
---|
148 | 0064 <span class="keyword">return</span> |
---|
149 | 0065 <span class="keyword">end</span> |
---|
150 | 0066 <span class="keyword">else</span> |
---|
151 | 0067 testfile=0; |
---|
152 | 0068 <span class="keyword">end</span> |
---|
153 | 0069 [ndims,nvars,ngatts] = netcdf.inq(nc);<span class="comment">%nbre of dimensions, variables, attributes</span> |
---|
154 | 0070 |
---|
155 | 0071 <span class="comment">% -------- read global attributes (constants)-----------</span> |
---|
156 | 0072 att_key={};<span class="comment">%default</span> |
---|
157 | 0073 iatt_g=0; |
---|
158 | 0074 Data.ListGlobalAttribute={};<span class="comment">%default</span> |
---|
159 | 0075 <span class="keyword">for</span> iatt=1:ngatts |
---|
160 | 0076 keystr= netcdf.inqAttName(nc,netcdf.getConstant(<span class="string">'NC_GLOBAL'</span>),iatt-1); |
---|
161 | 0077 indstr1=regexp(keystr,<span class="string">'\\'</span>);<span class="comment">%replace dots'</span> |
---|
162 | 0078 indstr2=regexp(keystr,<span class="string">'\.'</span>);<span class="comment">%replace dots'</span> |
---|
163 | 0079 <span class="keyword">if</span> isempty(indstr1) && isempty(indstr2) |
---|
164 | 0080 valuestr = netcdf.getAtt(nc,netcdf.getConstant(<span class="string">'NC_GLOBAL'</span>),keystr); |
---|
165 | 0081 <span class="keyword">if</span> ischar(valuestr) && length(valuestr)<200 |
---|
166 | 0082 iatt_g=iatt_g+1; |
---|
167 | 0083 indstr1=regexp(keystr,<span class="string">'\\'</span>);<span class="comment">%replace dots'</span> |
---|
168 | 0084 indstr2=regexp(keystr,<span class="string">'\.'</span>);<span class="comment">%replace dots'</span> |
---|
169 | 0085 <span class="keyword">if</span> isempty(indstr1) && isempty(indstr2) |
---|
170 | 0086 eval([<span class="string">'Data.'</span> keystr <span class="string">'='''</span> valuestr <span class="string">''';'</span>]) |
---|
171 | 0087 att_key{iatt_g}=keystr; |
---|
172 | 0088 <span class="keyword">end</span> |
---|
173 | 0089 <span class="keyword">elseif</span> isempty(valuestr) |
---|
174 | 0090 iatt_g=iatt_g+1; |
---|
175 | 0091 eval([<span class="string">'Data.'</span> keystr <span class="string">'=[];'</span>]) |
---|
176 | 0092 att_key{iatt_g}=keystr; |
---|
177 | 0093 <span class="keyword">elseif</span> isnumeric(valuestr) |
---|
178 | 0094 iatt_g=iatt_g+1; |
---|
179 | 0095 eval([<span class="string">'Data.'</span> keystr <span class="string">'=valuestr;'</span>]) |
---|
180 | 0096 att_key{iatt_g}=keystr; |
---|
181 | 0097 <span class="keyword">end</span> |
---|
182 | 0098 <span class="keyword">end</span> |
---|
183 | 0099 <span class="keyword">end</span> |
---|
184 | 0100 Data.ListGlobalAttribute=att_key; |
---|
185 | 0101 |
---|
186 | 0102 <span class="comment">% -------- read dimensions -----------</span> |
---|
187 | 0103 dim_name={}; |
---|
188 | 0104 dim_value=[]; |
---|
189 | 0105 <span class="keyword">for</span> idim=1:ndims<span class="comment">%length(dim_read);</span> |
---|
190 | 0106 [dim_name{idim},dim_value(idim)] = netcdf.inqDim(nc,idim-1); |
---|
191 | 0107 <span class="keyword">end</span> |
---|
192 | 0108 <span class="keyword">if</span> ~isempty(dim_name) && ~isempty(dim_value) |
---|
193 | 0109 Data.ListDimName=dim_name; |
---|
194 | 0110 Data.DimValue=dim_value; |
---|
195 | 0111 <span class="comment">% DimIndices=[1:ndims]; %index of the dimension in the netcdf file</span> |
---|
196 | 0112 dim_used=zeros(1,ndims);<span class="comment">%initialize test of used dimensions</span> |
---|
197 | 0113 <span class="keyword">end</span> |
---|
198 | 0114 |
---|
199 | 0115 <span class="comment">% -------- read variables -----------</span> |
---|
200 | 0116 var_read={}; <span class="comment">%default</span> |
---|
201 | 0117 dimids={}; |
---|
202 | 0118 nbatt=[]; |
---|
203 | 0119 <span class="keyword">for</span> ivar=1:nvars |
---|
204 | 0120 [var_read{ivar},xtype,dimids{ivar},nbatt(ivar)] = netcdf.inqVar(nc,ivar-1); |
---|
205 | 0121 <span class="keyword">end</span> |
---|
206 | 0122 var_index=1:nvars; <span class="comment">%default set of variable indices in the netcdf file</span> |
---|
207 | 0123 testmulti=0; |
---|
208 | 0124 OutputList=[]; |
---|
209 | 0125 |
---|
210 | 0126 <span class="comment">%select input variables, if requested by the input ListVarName</span> |
---|
211 | 0127 <span class="keyword">if</span> isequal(ListVarName,<span class="string">'*'</span>)||isemptyl(ListVarName) |
---|
212 | 0128 sizvar=size(ListVarName); |
---|
213 | 0129 testmulti=(sizvar(1)>1); |
---|
214 | 0130 var_index=zeros(1,sizvar(2));<span class="comment">%default</span> |
---|
215 | 0131 <span class="keyword">if</span> testmulti |
---|
216 | 0132 OutputList=ListVarName(1,:); |
---|
217 | 0133 testend=0; |
---|
218 | 0134 <span class="keyword">for</span> iline=1:sizvar(1) |
---|
219 | 0135 <span class="keyword">if</span> testend |
---|
220 | 0136 <span class="keyword">break</span> |
---|
221 | 0137 <span class="keyword">end</span> |
---|
222 | 0138 <span class="comment">% var_index=zeros(size(ListVarName));%default</span> |
---|
223 | 0139 <span class="keyword">for</span> ivar=1:sizvar(2) |
---|
224 | 0140 <span class="keyword">if</span> ~isempty(ListVarName{iline,ivar}) |
---|
225 | 0141 <span class="keyword">for</span> ilist=1:nvars |
---|
226 | 0142 <span class="keyword">if</span> isequal(var_read{ilist},ListVarName{iline,ivar}) |
---|
227 | 0143 var_index(ivar)=ilist; |
---|
228 | 0144 <span class="comment">% var_detect(ivar)=1;</span> |
---|
229 | 0145 <span class="keyword">break</span> |
---|
230 | 0146 <span class="keyword">end</span> |
---|
231 | 0147 <span class="keyword">end</span> |
---|
232 | 0148 <span class="keyword">if</span> ivar==1 |
---|
233 | 0149 <span class="keyword">if</span> var_index(ivar)==0 |
---|
234 | 0150 <span class="keyword">break</span><span class="comment">%go to next line if the first nc variable is not found</span> |
---|
235 | 0151 <span class="keyword">else</span> |
---|
236 | 0152 testend=1; <span class="comment">%this line will be read</span> |
---|
237 | 0153 ichoice=iline-1; <span class="comment">%selectedline number in the list of input names of variables</span> |
---|
238 | 0154 <span class="keyword">end</span> |
---|
239 | 0155 <span class="keyword">end</span> |
---|
240 | 0156 <span class="keyword">end</span> |
---|
241 | 0157 <span class="keyword">end</span> |
---|
242 | 0158 <span class="keyword">end</span> |
---|
243 | 0159 <span class="keyword">else</span> <span class="comment">%single list of input variables</span> |
---|
244 | 0160 <span class="keyword">for</span> ivar=1:sizvar(2) |
---|
245 | 0161 <span class="keyword">for</span> ilist=1:nvars |
---|
246 | 0162 <span class="keyword">if</span> isequal(var_read{ilist},ListVarName{ivar}) |
---|
247 | 0163 var_index(ivar)=ilist; |
---|
248 | 0164 var_detect(ivar)=1; |
---|
249 | 0165 <span class="keyword">break</span> |
---|
250 | 0166 <span class="keyword">end</span> |
---|
251 | 0167 <span class="keyword">end</span> |
---|
252 | 0168 <span class="keyword">end</span> |
---|
253 | 0169 <span class="keyword">end</span> |
---|
254 | 0170 list_index=find(var_index); |
---|
255 | 0171 <span class="keyword">if</span> ~isempty(list_index) |
---|
256 | 0172 <span class="keyword">if</span> testmulti |
---|
257 | 0173 OutputList=OutputList(list_index); |
---|
258 | 0174 <span class="keyword">end</span> |
---|
259 | 0175 var_index=var_index(list_index); |
---|
260 | 0176 var_detect=(var_index~=0); |
---|
261 | 0177 var_read=var_read(var_index); |
---|
262 | 0178 <span class="keyword">end</span> |
---|
263 | 0179 <span class="keyword">end</span> |
---|
264 | 0180 |
---|
265 | 0181 |
---|
266 | 0182 <span class="comment">%select variable attributes and associate dimensions</span> |
---|
267 | 0183 <span class="comment">% var_dim_index=[]; %default</span> |
---|
268 | 0184 Data.ListVarName={};<span class="comment">%default</span> |
---|
269 | 0185 VarDimIndex={};<span class="comment">%default</span> |
---|
270 | 0186 <span class="keyword">for</span> ivar=1:length(var_read) |
---|
271 | 0187 <span class="keyword">if</span> testmulti |
---|
272 | 0188 Data.ListVarName{ivar}=OutputList{ivar};<span class="comment">%new name given by ListVarName(1,:)</span> |
---|
273 | 0189 <span class="keyword">else</span> |
---|
274 | 0190 Data.ListVarName{ivar}=var_read{ivar};<span class="comment">%name of the variable</span> |
---|
275 | 0191 <span class="keyword">end</span> |
---|
276 | 0192 var_dim=dimids{var_index(ivar)}+1; <span class="comment">%dimension indices used by the variable</span> |
---|
277 | 0193 dim_used(var_dim)=ones(size(var_dim)); |
---|
278 | 0194 VarDimIndex{ivar}=var_dim; |
---|
279 | 0195 |
---|
280 | 0196 <span class="comment">%variable attributes</span> |
---|
281 | 0197 <span class="keyword">if</span> ivar==1 |
---|
282 | 0198 Data.VarAttribute={};<span class="comment">%initialisation of the list of variable attributes</span> |
---|
283 | 0199 <span class="keyword">end</span> |
---|
284 | 0200 <span class="comment">%variable attributes</span> |
---|
285 | 0201 <span class="keyword">for</span> iatt=1:nbatt(var_index(ivar)) |
---|
286 | 0202 attname = netcdf.inqAttName(nc,var_index(ivar)-1,iatt-1); |
---|
287 | 0203 valuestr= netcdf.getAtt(nc,var_index(ivar)-1,attname); |
---|
288 | 0204 <span class="keyword">if</span> ischar(valuestr) |
---|
289 | 0205 eval([<span class="string">'Data.VarAttribute{ivar}.'</span> attname <span class="string">'='''</span> valuestr <span class="string">''';'</span>]) |
---|
290 | 0206 <span class="keyword">elseif</span> isempty(valuestr) |
---|
291 | 0207 eval([<span class="string">'Data.VarAttribute{ivar}.'</span> attname <span class="string">'=[];'</span>]) |
---|
292 | 0208 <span class="keyword">elseif</span> isnumeric(valuestr) |
---|
293 | 0209 eval([<span class="string">'Data.VarAttribute{ivar}.'</span> attname <span class="string">'=valuestr;'</span>]) |
---|
294 | 0210 <span class="keyword">end</span> |
---|
295 | 0211 <span class="keyword">end</span> |
---|
296 | 0212 <span class="keyword">end</span> |
---|
297 | 0213 |
---|
298 | 0214 <span class="comment">%select the used dimensions</span> |
---|
299 | 0215 <span class="keyword">if</span> isempty(var_read) |
---|
300 | 0216 <span class="keyword">if</span> isfield(Data,<span class="string">'ListDimName'</span>) && isfield(Data,<span class="string">'DimValue'</span>) |
---|
301 | 0217 Data=rmfield(Data,<span class="string">'ListDimName'</span>); |
---|
302 | 0218 Data=rmfield(Data,<span class="string">'DimValue'</span>); |
---|
303 | 0219 <span class="keyword">end</span> |
---|
304 | 0220 <span class="keyword">else</span> |
---|
305 | 0221 <span class="comment">% list_dim=1:ndims;</span> |
---|
306 | 0222 dim_index=find(dim_used); |
---|
307 | 0223 <span class="comment">% list_dim=list_dim(dim_index);</span> |
---|
308 | 0224 old2new=cumsum(dim_used); |
---|
309 | 0225 Data.ListDimName=Data.ListDimName(dim_index); |
---|
310 | 0226 Data.DimValue=Data.DimValue(dim_index); |
---|
311 | 0227 <span class="keyword">end</span> |
---|
312 | 0228 <span class="keyword">for</span> ivar=1:length(var_read) |
---|
313 | 0229 Data.VarDimIndex{ivar}=old2new(VarDimIndex{ivar});<span class="comment">% ENLEVER Data.VarDimIndex ulterieurement</span> |
---|
314 | 0230 Data.VarDimName{ivar}=Data.ListDimName(Data.VarDimIndex{ivar}); |
---|
315 | 0231 <span class="keyword">end</span> |
---|
316 | 0232 <span class="comment">%variable values</span> |
---|
317 | 0233 <span class="keyword">if</span> ~isempty(ListVarName) |
---|
318 | 0234 <span class="keyword">for</span> ivar=1:length(Data.ListVarName) |
---|
319 | 0235 eval([<span class="string">'Data.'</span> Data.ListVarName{ivar} <span class="string">'=netcdf.getVar(nc,var_index(ivar)-1);'</span>])<span class="comment">%read the variable data</span> |
---|
320 | 0236 eval([<span class="string">'Data.'</span> Data.ListVarName{ivar} <span class="string">'=Data.'</span> Data.ListVarName{ivar} <span class="string">''';'</span>])<span class="comment">%read the variable data</span> |
---|
321 | 0237 <span class="keyword">end</span> |
---|
322 | 0238 <span class="keyword">end</span> |
---|
323 | 0239 <span class="comment">% -------- close fle-----------</span> |
---|
324 | 0240 <span class="keyword">if</span> testfile==1 |
---|
325 | 0241 netcdf.close(nc) |
---|
326 | 0242 <span class="keyword">end</span> |
---|
327 | 0243 <span class="keyword">else</span> |
---|
328 | 0244 [Data,var_detect,ichoice]=<a href="nc2struct_toolbox.html" class="code" title="function [Data,var_detect,ichoice]=nc2struct_toolbox(nc,ListVarName)">nc2struct_toolbox</a>(nc,ListVarName); |
---|
329 | 0245 <span class="keyword">end</span></pre></div> |
---|
330 | <hr><address>Generated on Fri 13-Nov-2009 11:17:03 by <strong><a href="http://www.artefact.tk/software/matlab/m2html/">m2html</a></strong> © 2003</address> |
---|
331 | </body> |
---|
332 | </html> |
---|