[37] | 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> |
---|