1 | %'sub_field': combines two input fields
|
---|
2 | %
|
---|
3 | % the two fields are subtstracted when of the same nature (scalar or
|
---|
4 | % vector), if the coordinates do not coincide, the second field is
|
---|
5 | % interpolated on the cooridintes of the first one
|
---|
6 | %
|
---|
7 | % when scalar and vectors are combined, the fields are just merged in a single matlab structure for common visualisation
|
---|
8 | %-----------------------------------------------------------------------
|
---|
9 | % function SubData=sub_field(Field,XmlData,Field_1)
|
---|
10 | %
|
---|
11 | % OUTPUT:
|
---|
12 | % SubData: structure representing the resulting field
|
---|
13 | %
|
---|
14 | % INPUT:
|
---|
15 | % Field: matlab structure representing the first field
|
---|
16 | % Field_1:matlab structure representing the second field
|
---|
17 |
|
---|
18 | %=======================================================================
|
---|
19 | % Copyright 2008-2021, LEGI UMR 5519 / CNRS UGA G-INP, Grenoble, France
|
---|
20 | % http://www.legi.grenoble-inp.fr
|
---|
21 | % Joel.Sommeria - Joel.Sommeria (A) legi.cnrs.fr
|
---|
22 | %
|
---|
23 | % This file is part of the toolbox UVMAT.
|
---|
24 | %
|
---|
25 | % UVMAT is free software; you can redistribute it and/or modify
|
---|
26 | % it under the terms of the GNU General Public License as published
|
---|
27 | % by the Free Software Foundation; either version 2 of the license,
|
---|
28 | % or (at your option) any later version.
|
---|
29 | %
|
---|
30 | % UVMAT is distributed in the hope that it will be useful,
|
---|
31 | % but WITHOUT ANY WARRANTY; without even the implied warranty of
|
---|
32 | % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
---|
33 | % GNU General Public License (see LICENSE.txt) for more details.
|
---|
34 | %=======================================================================
|
---|
35 |
|
---|
36 | function SubData=sub_field(Field,XmlData,Field_1)
|
---|
37 |
|
---|
38 | SubData=[];
|
---|
39 | if strcmp(Field,'*')
|
---|
40 | return
|
---|
41 | end
|
---|
42 | if nargin<3
|
---|
43 | SubData=Field;
|
---|
44 | return
|
---|
45 | end
|
---|
46 | if ~isfield(Field_1,'VarAttribute')
|
---|
47 | Field_1.VarAttribute={};
|
---|
48 | end
|
---|
49 |
|
---|
50 | %% global attributes
|
---|
51 | SubData.ListGlobalAttribute={};%default
|
---|
52 | %transfer global attributes of Field
|
---|
53 | if isfield(Field,'ListGlobalAttribute')
|
---|
54 | SubData.ListGlobalAttribute=Field.ListGlobalAttribute;
|
---|
55 | for ilist=1:numel(Field.ListGlobalAttribute)
|
---|
56 | AttrName=Field.ListGlobalAttribute{ilist};
|
---|
57 | SubData.(AttrName)=Field.(AttrName);
|
---|
58 | end
|
---|
59 | end
|
---|
60 | %transfer global attributes of Field_1
|
---|
61 | if isfield(Field_1,'ListGlobalAttribute')
|
---|
62 | for ilist=1:numel(Field_1.ListGlobalAttribute)
|
---|
63 | AttrName=Field_1.ListGlobalAttribute{ilist};
|
---|
64 | AttrNameNew=AttrName;
|
---|
65 | while ~isempty(find(strcmp(AttrNameNew,SubData.ListGlobalAttribute)))&&~isequal(Field_1.(AttrNameNew),Field.(AttrNameNew))
|
---|
66 | AttrNameNew=[AttrNameNew '_1'];
|
---|
67 | end
|
---|
68 | if ~isfield(Field,AttrName) || ~isequal(Field_1.(AttrName),Field.(AttrName))
|
---|
69 | SubData.ListGlobalAttribute=[SubData.ListGlobalAttribute {AttrNameNew}];
|
---|
70 | SubData.(AttrNameNew)=Field_1.(AttrName);
|
---|
71 | end
|
---|
72 | end
|
---|
73 | end
|
---|
74 |
|
---|
75 | %% variables
|
---|
76 | %reproduce variables of the first field and list its dimensions
|
---|
77 | SubData.ListVarName=Field.ListVarName;
|
---|
78 | SubData.VarDimName=Field.VarDimName;
|
---|
79 | if isfield(Field,'VarAttribute')
|
---|
80 | SubData.VarAttribute=Field.VarAttribute;
|
---|
81 | end
|
---|
82 | ListDimName={};
|
---|
83 | for ilist=1:numel(Field.ListVarName)
|
---|
84 | VarName=Field.ListVarName{ilist};
|
---|
85 | SubData.(VarName)=Field.(VarName);
|
---|
86 | SubData.VarAttribute{ilist}.CheckSub=0;
|
---|
87 | DimCell=Field.VarDimName{ilist};
|
---|
88 | if ischar(DimCell)
|
---|
89 | DimCell={DimCell};
|
---|
90 | end
|
---|
91 | for idim=1:numel(DimCell)
|
---|
92 | if isempty(find(strcmp(DimCell{idim},ListDimName)))
|
---|
93 | ListDimName=[ListDimName DimCell(idim)];
|
---|
94 | end
|
---|
95 | end
|
---|
96 | end
|
---|
97 |
|
---|
98 | %% field request
|
---|
99 | ProjModeRequest=cell(size(Field.ListVarName));
|
---|
100 | if isfield(Field,'VarAttribute')
|
---|
101 | for ilist=1:numel(Field.VarAttribute)
|
---|
102 | if isfield(Field.VarAttribute{ilist},'ProjModeRequest')
|
---|
103 | ProjModeRequest{ilist}=Field.VarAttribute{ilist}.ProjModeRequest;
|
---|
104 | end
|
---|
105 | end
|
---|
106 | end
|
---|
107 | ProjModeRequest_1=cell(size(Field_1.ListVarName));
|
---|
108 | if isfield(Field_1,'VarAttribute')
|
---|
109 | for ilist=1:numel(Field_1.VarAttribute)
|
---|
110 | if isfield(Field_1.VarAttribute{ilist},'ProjModeRequest')
|
---|
111 | ProjModeRequest_1{ilist}=Field_1.VarAttribute{ilist}.ProjModeRequest;
|
---|
112 | end
|
---|
113 | end
|
---|
114 | end
|
---|
115 |
|
---|
116 | %% rename the dimensions of the second field if identical to those of the first
|
---|
117 | for ilist=1:numel(Field_1.VarDimName)
|
---|
118 | DimCell=Field_1.VarDimName{ilist};
|
---|
119 | if ischar(DimCell)
|
---|
120 | DimCell={DimCell};
|
---|
121 | end
|
---|
122 | for idim=1:numel(DimCell)
|
---|
123 | ind_dim=find(strcmp(DimCell{idim},ListDimName));
|
---|
124 | if ~isempty(ind_dim)
|
---|
125 | if ischar(Field_1.VarDimName{ilist})
|
---|
126 | Field_1.VarDimName{ilist}=Field_1.VarDimName(ilist);
|
---|
127 | end
|
---|
128 | Field_1.VarDimName{ilist}{idim}=[ListDimName{ind_dim} '_1'];
|
---|
129 | end
|
---|
130 | end
|
---|
131 | end
|
---|
132 |
|
---|
133 | %% look for coordinates common to Field in Field_1
|
---|
134 | ind_remove=false(size(Field_1.ListVarName));
|
---|
135 | % loop on the variables of the second field Field_1
|
---|
136 | for ilist=1:numel(Field_1.VarAttribute)
|
---|
137 | % case of variable with a single dimension
|
---|
138 | if ~isempty(Field_1.VarAttribute{ilist}) && isfield(Field_1.VarAttribute{ilist},'Role')&&~isempty(regexp(Field_1.VarAttribute{ilist}.Role,'^coord'))% if variable with Role coord... is found.
|
---|
139 | OldDimName=Field_1.VarDimName{ilist};
|
---|
140 | if ischar(OldDimName), OldDimName={OldDimName}; end% transform char string to cell if relevant
|
---|
141 | if numel(OldDimName)==1
|
---|
142 | OldDim=Field_1.(Field_1.ListVarName{ilist});% get variable
|
---|
143 | %look for the existence of the variable OldDim in the first field Field
|
---|
144 | for i1=1:numel(Field.ListVarName)
|
---|
145 | if isequal(Field.(Field.ListVarName{i1}),OldDim) &&...
|
---|
146 | ((isempty(ProjModeRequest{i1}) && isempty(ProjModeRequest_1{ilist})) || strcmp(ProjModeRequest{i1},ProjModeRequest_1{ilist}))
|
---|
147 | ind_remove(ilist)=1;
|
---|
148 | NewDimName=Field.VarDimName{i1};
|
---|
149 | if ischar(NewDimName), NewDimName={NewDimName}; end %transform char chain to cell if needed
|
---|
150 | Field_1.VarDimName=regexprep_r(Field_1.VarDimName,['^' OldDimName{1} '$'],NewDimName{1});% change the var name of Field_1 to the corresponding var name of Field
|
---|
151 | end
|
---|
152 | end
|
---|
153 | end
|
---|
154 | end
|
---|
155 | end
|
---|
156 | if ~isempty(find(ind_remove, 1))
|
---|
157 | Field_1.ListVarName(ind_remove)=[];%removes the redondent coordinate
|
---|
158 | Field_1.VarDimName(ind_remove)=[];
|
---|
159 | Field_1.VarAttribute(ind_remove)=[];
|
---|
160 | end
|
---|
161 |
|
---|
162 | %% append the other variables of the second field, modifying their name if needed
|
---|
163 | ListVarNameSub=Field_1.ListVarName;
|
---|
164 | ListVarNameNew=ListVarNameSub;
|
---|
165 | check_rename=zeros(size(ListVarNameSub));
|
---|
166 | check_remove=zeros(size(ListVarNameSub));
|
---|
167 | VarDimNameSub=Field_1.VarDimName;
|
---|
168 | VarAttributeSub={};
|
---|
169 | if isfield(Field_1,'VarAttribute')&&~isempty(Field_1.VarAttribute)
|
---|
170 | for ilist=1:numel(ListVarNameSub)
|
---|
171 | ind_prev=find(strcmp(ListVarNameSub{ilist},Field.ListVarName),1);% look for duplicated variable name
|
---|
172 | if ~isempty(ind_prev)% variable name exists in Field
|
---|
173 | if isfield(Field_1.VarAttribute{ilist},'Role')&&...
|
---|
174 | ismember(Field_1.VarAttribute{ilist}.Role,{'coord_x','coord_y','scalar','vector_x','vector_y','errorflag'})
|
---|
175 | ListVarNameNew{ilist}=[ListVarNameSub{ilist} '_1'];%modify the name of the second variable
|
---|
176 | check_rename(ilist)=1;
|
---|
177 | else
|
---|
178 | check_remove(ilist)=1;% variable will be removed
|
---|
179 | end
|
---|
180 | end
|
---|
181 | end
|
---|
182 | ListVarNameSub=ListVarNameSub(~check_remove); %eliminate removed variables from the list of the second field
|
---|
183 | ListVarNameNew=ListVarNameNew(~check_remove); % %list of renaimed varaibles corresponding to ListVarNameSub
|
---|
184 | VarDimNameSub=Field_1.VarDimName(~check_remove);
|
---|
185 | if numel(Field_1.VarAttribute)<max(find(~check_remove))
|
---|
186 | for ilist=numel(Field_1.VarAttribute)+1:max(find(~check_remove))
|
---|
187 | Field_1.VarAttribute{ilist}={};
|
---|
188 | end
|
---|
189 | end
|
---|
190 | VarAttributeSub=Field_1.VarAttribute(~check_remove);
|
---|
191 | check_rename=check_rename(~check_remove);
|
---|
192 | end
|
---|
193 | % apply the variable renaming and mark the second field variables with the attribute .CheckSub
|
---|
194 | for ilist=1:numel(ListVarNameSub)
|
---|
195 | SubData.(ListVarNameNew{ilist})=Field_1.(ListVarNameSub{ilist});% copy the variable content to the new name
|
---|
196 | if check_rename(ilist)
|
---|
197 | % replace name in field expression FieldName, e.g. 'norm(U,V)'-> 'norm(U_1,V_1)'
|
---|
198 | if isfield(VarAttributeSub{ilist},'FieldName')
|
---|
199 | for ivar=1:numel(find(check_rename))
|
---|
200 | VarAttributeSub{ilist}.FieldName=regexprep_r(VarAttributeSub{ilist}.FieldName,...
|
---|
201 | ListVarNameSub{ivar},ListVarNameNew{ivar});
|
---|
202 | end
|
---|
203 | end
|
---|
204 | end
|
---|
205 | VarAttributeSub{ilist}.CheckSub=1;% mark that the field needs to be substracted as an attribute
|
---|
206 | end
|
---|
207 |
|
---|
208 | SubData.ListVarName=[SubData.ListVarName ListVarNameNew];
|
---|
209 | SubData.VarDimName=[SubData.VarDimName VarDimNameSub];
|
---|
210 | SubData.VarAttribute=[SubData.VarAttribute VarAttributeSub];
|
---|
211 |
|
---|
212 | %% substrat fields when possible
|
---|
213 | [CellInfo,NbDim,errormsg]=find_field_cells(SubData);
|
---|
214 | ind_remove=false(size(SubData.ListVarName));
|
---|
215 | ivar=[];
|
---|
216 | ivar_1=[];
|
---|
217 | for icell=1:numel(CellInfo)
|
---|
218 | if ~isempty(CellInfo{icell})
|
---|
219 | % if two scalar are in the same cell
|
---|
220 | if isfield(CellInfo{icell},'VarIndex_scalar') && numel(CellInfo{icell}.VarIndex_scalar)==2 && SubData.VarAttribute{CellInfo{icell}.VarIndex_scalar(2)}.CheckSub;
|
---|
221 | ivar=[ivar CellInfo{icell}.VarIndex_scalar(1)];
|
---|
222 | ivar_1=[ivar_1 CellInfo{icell}.VarIndex_scalar(2)];
|
---|
223 | end
|
---|
224 | % if two vector u components are in the same cell
|
---|
225 | if isfield(CellInfo{icell},'VarIndex_vector_x') && numel(CellInfo{icell}.VarIndex_vector_x)==2 && SubData.VarAttribute{CellInfo{icell}.VarIndex_vector_x(2)}.CheckSub;
|
---|
226 | ivar=[ivar CellInfo{icell}.VarIndex_vector_x(1)];
|
---|
227 | ivar_1=[ivar_1 CellInfo{icell}.VarIndex_vector_x(2)];
|
---|
228 | end
|
---|
229 | % if two vector v components are in the same cell
|
---|
230 | if isfield(CellInfo{icell},'VarIndex_vector_y') && numel(CellInfo{icell}.VarIndex_vector_y)==2 && SubData.VarAttribute{CellInfo{icell}.VarIndex_vector_y(2)}.CheckSub;
|
---|
231 | ivar=[ivar CellInfo{icell}.VarIndex_vector_y(1)];
|
---|
232 | ivar_1=[ivar_1 CellInfo{icell}.VarIndex_vector_y(2)];
|
---|
233 | end
|
---|
234 | % merge the error flags if needed
|
---|
235 | if isfield(CellInfo{icell},'VarIndex_errorflag') && numel(CellInfo{icell}.VarIndex_errorflag)==2 && SubData.VarAttribute{CellInfo{icell}.VarIndex_vector_y(2)}.CheckSub;
|
---|
236 | ivar_flag=CellInfo{icell}.VarIndex_errorflag(1);
|
---|
237 | ivar_flag_1=CellInfo{icell}.VarIndex_errorflag(2);
|
---|
238 | VarName=SubData.ListVarName{ivar_flag};
|
---|
239 | VarName_1=SubData.ListVarName{ivar_flag_1};
|
---|
240 | SubData.(VarName)=SubData.(VarName)~=0 | SubData.(VarName_1)~=0;% combine the error flags of the two fields
|
---|
241 | ind_remove(ivar_flag_1)=1;
|
---|
242 | end
|
---|
243 | end
|
---|
244 | end
|
---|
245 | % subtract fields if relevant
|
---|
246 | for imod=1:numel(ivar)
|
---|
247 | VarName=SubData.ListVarName{ivar(imod)};
|
---|
248 | VarName_1=SubData.ListVarName{ivar_1(imod)};
|
---|
249 | SubData.(VarName)=double(SubData.(VarName))-double(SubData.(VarName_1));
|
---|
250 | ind_remove(ivar_1(imod))=1;
|
---|
251 | end
|
---|
252 | SubData.ListVarName(ind_remove)=[];
|
---|
253 | SubData.VarDimName(ind_remove)=[];
|
---|
254 | SubData.VarAttribute(ind_remove)=[];
|
---|
255 |
|
---|
256 | function OutputCell=regexprep_r(InputCell,search_string,new_string)
|
---|
257 | if ischar(InputCell); InputCell={InputCell}; end
|
---|
258 | OutputCell=InputCell;%default
|
---|
259 | for icell=1:numel(InputCell)
|
---|
260 | OutputCell{icell}=regexprep(InputCell{icell},search_string,new_string);
|
---|
261 | end
|
---|
262 |
|
---|
263 |
|
---|