[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_toolbox</title> |
---|
| 6 | <meta name="keywords" content="nc2struct_toolbox"> |
---|
| 7 | <meta name="description" content="'nc2struct_toolbox': transform a netcdf file in a corresponding matlab structure, USE OLD NETCDF LIBRARY"> |
---|
| 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_toolbox.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_toolbox |
---|
| 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_toolbox': transform a netcdf file in a corresponding matlab structure, USE OLD NETCDF LIBRARY</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_toolbox(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_toolbox': transform a netcdf file in a corresponding matlab structure, USE OLD NETCDF LIBRARY |
---|
| 31 | ---------------------------------------------------------------------- |
---|
| 32 | function [Data,var_detect,ichoice]=nc2struct_toolbox(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 | </ul> |
---|
| 78 | This function is called by: |
---|
| 79 | <ul style="list-style-image:url(../matlabicon.gif)"> |
---|
| 80 | <li><a href="nc2struct.html" class="code" title="function [Data,var_detect,ichoice]=nc2struct(nc,ListVarName)">nc2struct</a> 'nc2struct': transform a netcdf file in a corresponding matlab structure</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_toolbox': transform a netcdf file in a corresponding matlab structure, USE OLD NETCDF LIBRARY</span> |
---|
| 86 | 0002 <span class="comment">%----------------------------------------------------------------------</span> |
---|
| 87 | 0003 <span class="comment">% function [Data,var_detect,ichoice]=nc2struct_toolbox(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_toolbox(nc,ListVarName)</a> |
---|
| 129 | 0045 |
---|
| 130 | 0046 <span class="comment">%default output</span> |
---|
| 131 | 0047 Data=[]; |
---|
| 132 | 0048 <span class="comment">% ListIndex=[];</span> |
---|
| 133 | 0049 var_detect=[]; |
---|
| 134 | 0050 ichoice=[];<span class="comment">%default</span> |
---|
| 135 | 0051 <span class="comment">%open the netcdf file for reading</span> |
---|
| 136 | 0052 <span class="keyword">if</span> ischar(nc) |
---|
| 137 | 0053 <span class="keyword">if</span> exist(nc,<span class="string">'file'</span>) |
---|
| 138 | 0054 nc=netcdf(nc,<span class="string">'nowrite'</span>); |
---|
| 139 | 0055 testfile=1; |
---|
| 140 | 0056 <span class="keyword">else</span> |
---|
| 141 | 0057 Data.Txt=[<span class="string">'ERROR:file '</span> nc <span class="string">' does not exist'</span>]; |
---|
| 142 | 0058 <span class="keyword">return</span> |
---|
| 143 | 0059 <span class="keyword">end</span> |
---|
| 144 | 0060 <span class="keyword">else</span> |
---|
| 145 | 0061 testfile=0; |
---|
| 146 | 0062 <span class="keyword">end</span> |
---|
| 147 | 0063 |
---|
| 148 | 0064 <span class="comment">% -------- read global attributes -----------</span> |
---|
| 149 | 0065 att_read=att(nc);<span class="comment">%cell of 'global attributes' (nc objects)</span> |
---|
| 150 | 0066 att_key={};<span class="comment">%default</span> |
---|
| 151 | 0067 iatt_g=0; |
---|
| 152 | 0068 <span class="keyword">for</span> iatt=1:length(att_read) |
---|
| 153 | 0069 aa=att_read{iatt}; |
---|
| 154 | 0070 keystr=name(aa); |
---|
| 155 | 0071 indstr1=regexp(keystr,<span class="string">'\\'</span>);<span class="comment">%replace dots'</span> |
---|
| 156 | 0072 indstr2=regexp(keystr,<span class="string">'\.'</span>);<span class="comment">%replace dots'</span> |
---|
| 157 | 0073 <span class="keyword">if</span> ~isequal(keystr,<span class="string">'title'</span>) <span class="comment">% PROBLEM WITH civx files do not read 'title'</span> |
---|
| 158 | 0074 <span class="keyword">if</span> isempty(indstr1) && isempty(indstr2) <span class="comment">%</span> |
---|
| 159 | 0075 eval([<span class="string">'valuestr=nc.'</span> keystr <span class="string">'(:);'</span>]) |
---|
| 160 | 0076 <span class="keyword">if</span> ischar(valuestr) && length(valuestr)<200 |
---|
| 161 | 0077 iatt_g=iatt_g+1; |
---|
| 162 | 0078 indstr1=regexp(keystr,<span class="string">'\\'</span>);<span class="comment">%replace dots'</span> |
---|
| 163 | 0079 indstr2=regexp(keystr,<span class="string">'\.'</span>);<span class="comment">%replace dots'</span> |
---|
| 164 | 0080 <span class="keyword">if</span> isempty(indstr1) && isempty(indstr2) |
---|
| 165 | 0081 eval([<span class="string">'Data.'</span> keystr <span class="string">'='''</span> valuestr <span class="string">''';'</span>]) |
---|
| 166 | 0082 att_key{iatt_g}=keystr; |
---|
| 167 | 0083 <span class="keyword">end</span> |
---|
| 168 | 0084 <span class="keyword">elseif</span> isempty(valuestr) |
---|
| 169 | 0085 iatt_g=iatt_g+1; |
---|
| 170 | 0086 eval([<span class="string">'Data.'</span> keystr <span class="string">'=[];'</span>]) |
---|
| 171 | 0087 att_key{iatt_g}=keystr; |
---|
| 172 | 0088 <span class="keyword">elseif</span> isnumeric(valuestr) |
---|
| 173 | 0089 iatt_g=iatt_g+1; |
---|
| 174 | 0090 eval([<span class="string">'Data.'</span> keystr <span class="string">'=valuestr;'</span>]) |
---|
| 175 | 0091 att_key{iatt_g}=keystr; |
---|
| 176 | 0092 <span class="keyword">end</span> |
---|
| 177 | 0093 <span class="keyword">end</span> |
---|
| 178 | 0094 <span class="keyword">end</span> |
---|
| 179 | 0095 <span class="keyword">end</span> |
---|
| 180 | 0096 Data.ListGlobalAttribute=att_key; |
---|
| 181 | 0097 nbattr=length(att_key); |
---|
| 182 | 0098 neworder=[nbattr+1 [1:nbattr]]; |
---|
| 183 | 0099 Data=orderfields(Data,neworder); |
---|
| 184 | 0100 |
---|
| 185 | 0101 <span class="comment">% -------- read dimensions -----------</span> |
---|
| 186 | 0102 dim_read=dim(nc);<span class="comment">%cell of variable dimension names (nc objects)</span> |
---|
| 187 | 0103 dim_name={}; |
---|
| 188 | 0104 dim_value=[]; |
---|
| 189 | 0105 <span class="keyword">for</span> idim=1:length(dim_read); |
---|
| 190 | 0106 aa=dim_read{idim}; |
---|
| 191 | 0107 <span class="keyword">if</span> ~isempty(aa) |
---|
| 192 | 0108 dim_name{idim}=name(aa); |
---|
| 193 | 0109 dim_value(idim)=length(aa); |
---|
| 194 | 0110 <span class="keyword">end</span> |
---|
| 195 | 0111 <span class="keyword">end</span> |
---|
| 196 | 0112 <span class="keyword">if</span> ~isempty(dim_name) && ~isempty(dim_value) |
---|
| 197 | 0113 Data.ListDimName=dim_name; |
---|
| 198 | 0114 Data.DimValue=dim_value; |
---|
| 199 | 0115 used=zeros(1,length(dim_value));<span class="comment">%initialize test of used dimensions</span> |
---|
| 200 | 0116 <span class="keyword">end</span> |
---|
| 201 | 0117 |
---|
| 202 | 0118 <span class="comment">% -------- read variables -----------</span> |
---|
| 203 | 0119 var_read={}; <span class="comment">%default</span> |
---|
| 204 | 0120 testmulti=0; |
---|
| 205 | 0121 OutputList=[]; |
---|
| 206 | 0122 <span class="keyword">if</span> isequal(ListVarName,<span class="string">'*'</span>)|| isempty(ListVarName) |
---|
| 207 | 0123 var_read=var(nc);<span class="comment">%cell of all variables</span> |
---|
| 208 | 0124 <span class="keyword">elseif</span> ~isempty(ListVarName) |
---|
| 209 | 0125 sizvar=size(ListVarName); |
---|
| 210 | 0126 testmulti=(sizvar(1)>1); |
---|
| 211 | 0127 <span class="keyword">if</span> testmulti |
---|
| 212 | 0128 OutputList=ListVarName(1,:); |
---|
| 213 | 0129 testend=0; |
---|
| 214 | 0130 <span class="keyword">for</span> iline=1:sizvar(1) |
---|
| 215 | 0131 <span class="keyword">if</span> testend |
---|
| 216 | 0132 <span class="keyword">break</span> |
---|
| 217 | 0133 <span class="keyword">end</span> |
---|
| 218 | 0134 <span class="keyword">for</span> ivar=1:sizvar(2) |
---|
| 219 | 0135 var_read{ivar}=[];<span class="comment">%default</span> |
---|
| 220 | 0136 var_detect(ivar)=0;<span class="comment">%default</span> |
---|
| 221 | 0137 VarName=ListVarName{iline,ivar}; |
---|
| 222 | 0138 <span class="keyword">if</span> ~isempty(VarName) |
---|
| 223 | 0139 var_read{ivar}=nc{VarName};<span class="comment">%select the input variable names</span> |
---|
| 224 | 0140 <span class="keyword">if</span> ivar==1 |
---|
| 225 | 0141 <span class="keyword">if</span> isempty (var_read{ivar}) |
---|
| 226 | 0142 <span class="keyword">break</span><span class="comment">%go to next line if the first nc variable is not found</span> |
---|
| 227 | 0143 <span class="keyword">else</span> |
---|
| 228 | 0144 testend=1; <span class="comment">%this line will be read</span> |
---|
| 229 | 0145 ichoice=iline-1; <span class="comment">%selectedline number in the list of input names of variables</span> |
---|
| 230 | 0146 var_detect(ivar)=1; |
---|
| 231 | 0147 <span class="keyword">end</span> |
---|
| 232 | 0148 <span class="keyword">else</span> |
---|
| 233 | 0149 var_detect(ivar)=~isempty (var_read{ivar}); |
---|
| 234 | 0150 <span class="keyword">end</span> |
---|
| 235 | 0151 <span class="keyword">end</span> |
---|
| 236 | 0152 <span class="keyword">end</span> |
---|
| 237 | 0153 <span class="keyword">end</span> |
---|
| 238 | 0154 <span class="keyword">if</span> ~isempty(find(var_detect)) |
---|
| 239 | 0155 OutputList=OutputList(find(var_detect)); |
---|
| 240 | 0156 <span class="keyword">end</span> |
---|
| 241 | 0157 <span class="keyword">else</span> <span class="comment">%single list of input variables</span> |
---|
| 242 | 0158 var_detect=ones(size(ListVarName)); |
---|
| 243 | 0159 <span class="keyword">for</span> ivar=1:sizvar(2) |
---|
| 244 | 0160 var_read{ivar}=nc{ListVarName{ivar}};<span class="comment">%select the input variable names</span> |
---|
| 245 | 0161 var_detect(ivar)=~isempty(var_read{ivar}); |
---|
| 246 | 0162 <span class="keyword">end</span> |
---|
| 247 | 0163 <span class="keyword">end</span> |
---|
| 248 | 0164 var_read=var_read(find(var_detect)); |
---|
| 249 | 0165 <span class="keyword">end</span> |
---|
| 250 | 0166 <span class="comment">% var_dim_index=[]; %default</span> |
---|
| 251 | 0167 Data.ListVarName={};<span class="comment">%default</span> |
---|
| 252 | 0168 <span class="keyword">for</span> ivar=1:length(var_read) |
---|
| 253 | 0169 vv=var_read{ivar}; |
---|
| 254 | 0170 Data.ListVarName{ivar}=name(vv);<span class="comment">%name of the variable</span> |
---|
| 255 | 0171 <span class="keyword">if</span> testmulti |
---|
| 256 | 0172 Data.ListVarName{ivar}=OutputList{ivar}; |
---|
| 257 | 0173 <span class="keyword">else</span> |
---|
| 258 | 0174 Data.ListVarName{ivar}=name(vv);<span class="comment">%name of the variable</span> |
---|
| 259 | 0175 <span class="keyword">end</span> |
---|
| 260 | 0176 var_dim=dim(vv);<span class="comment">%dimension netcdf object of the variable</span> |
---|
| 261 | 0177 <span class="keyword">for</span> ivardim=1:length(var_dim) |
---|
| 262 | 0178 var_dim_name=name(var_dim{ivardim});<span class="comment">%name of the dimension</span> |
---|
| 263 | 0179 <span class="keyword">for</span> idim=1:length(dim_name)<span class="comment">% find the index of the current dimension in the list of dimensions</span> |
---|
| 264 | 0180 <span class="keyword">if</span> isequal(dim_name{idim},var_dim_name) |
---|
| 265 | 0181 Data.VarDimIndex{ivar}(ivardim)=idim; |
---|
| 266 | 0182 used(idim)=1; |
---|
| 267 | 0183 <span class="keyword">break</span> |
---|
| 268 | 0184 <span class="keyword">end</span> |
---|
| 269 | 0185 <span class="keyword">end</span> |
---|
| 270 | 0186 <span class="keyword">end</span> |
---|
| 271 | 0187 Data.VarDimName{ivar}={}; |
---|
| 272 | 0188 <span class="comment">%variable attributes</span> |
---|
| 273 | 0189 <span class="comment">% if ivar==1</span> |
---|
| 274 | 0190 Data.VarAttribute{ivar}=[];<span class="comment">%initialisation of the list of variable attributes</span> |
---|
| 275 | 0191 <span class="comment">% end</span> |
---|
| 276 | 0192 <span class="comment">%variable attributes</span> |
---|
| 277 | 0193 att_read=att(vv); |
---|
| 278 | 0194 <span class="keyword">for</span> iatt=1:length(att_read) |
---|
| 279 | 0195 aa=att_read{iatt}; |
---|
| 280 | 0196 eval([<span class="string">'valuestr=vv.'</span> name(aa) <span class="string">'(:);'</span>]) |
---|
| 281 | 0197 <span class="keyword">if</span> ischar(valuestr) |
---|
| 282 | 0198 eval([<span class="string">'Data.VarAttribute{ivar}.'</span> name(aa) <span class="string">'='''</span> valuestr <span class="string">''';'</span>]) |
---|
| 283 | 0199 <span class="keyword">elseif</span> isempty(valuestr) |
---|
| 284 | 0200 eval([<span class="string">'Data.VarAttribute{ivar}.'</span> name(aa) <span class="string">'=[];'</span>]) |
---|
| 285 | 0201 <span class="keyword">elseif</span> isnumeric(valuestr) |
---|
| 286 | 0202 eval([<span class="string">'Data.VarAttribute{ivar}.'</span> name(aa) <span class="string">'=valuestr;'</span>]) |
---|
| 287 | 0203 <span class="keyword">end</span> |
---|
| 288 | 0204 <span class="keyword">end</span> |
---|
| 289 | 0205 <span class="keyword">end</span> |
---|
| 290 | 0206 |
---|
| 291 | 0207 <span class="comment">%select the used dimensions</span> |
---|
| 292 | 0208 <span class="keyword">if</span> isempty(var_read) |
---|
| 293 | 0209 <span class="keyword">if</span> isfield(Data,<span class="string">'ListDimName'</span>) && isfield(Data,<span class="string">'DimValue'</span>) |
---|
| 294 | 0210 Data=rmfield(Data,<span class="string">'ListDimName'</span>); |
---|
| 295 | 0211 Data=rmfield(Data,<span class="string">'DimValue'</span>); |
---|
| 296 | 0212 <span class="keyword">end</span> |
---|
| 297 | 0213 <span class="keyword">else</span> |
---|
| 298 | 0214 old_dim_index=find(used); <span class="comment">%dimension indices which are used by the selected variables</span> |
---|
| 299 | 0215 old2new=cumsum(used); |
---|
| 300 | 0216 Data.ListDimName=Data.ListDimName(old_dim_index); |
---|
| 301 | 0217 Data.DimValue=Data.DimValue(old_dim_index); |
---|
| 302 | 0218 <span class="keyword">end</span> |
---|
| 303 | 0219 <span class="keyword">for</span> ivar=1:length(var_read) |
---|
| 304 | 0220 Data.VarDimIndex{ivar}=(old2new(Data.VarDimIndex{ivar})); |
---|
| 305 | 0221 Data.VarDimName{ivar}=(Data.ListDimName(Data.VarDimIndex{ivar})); |
---|
| 306 | 0222 <span class="keyword">end</span> |
---|
| 307 | 0223 <span class="comment">%variable values</span> |
---|
| 308 | 0224 |
---|
| 309 | 0225 <span class="keyword">if</span> ~isempty(ListVarName) |
---|
| 310 | 0226 <span class="keyword">for</span> ivar=1:length(Data.ListVarName) |
---|
| 311 | 0227 vv=var_read{ivar}; |
---|
| 312 | 0228 vdata=vv(:);<span class="comment">%data array of the field variable</span> |
---|
| 313 | 0229 eval([<span class="string">'Data.'</span> Data.ListVarName{ivar} <span class="string">'=vdata;'</span>])<span class="comment">%read the variable data</span> |
---|
| 314 | 0230 <span class="keyword">end</span> |
---|
| 315 | 0231 <span class="keyword">end</span> |
---|
| 316 | 0232 <span class="comment">% -------- close fle-----------</span> |
---|
| 317 | 0233 <span class="keyword">if</span> testfile==1 |
---|
| 318 | 0234 close(nc) |
---|
| 319 | 0235 <span class="keyword">end</span> |
---|
| 320 | 0236</pre></div> |
---|
| 321 | <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> |
---|
| 322 | </body> |
---|
| 323 | </html> |
---|