1 | function ismodelselfconsistent(md),
|
---|
2 | %ISMODELSELFCONSISTENT - check that model forms a closed form solvable problem.
|
---|
3 | %
|
---|
4 | % Usage:
|
---|
5 | % ismodelselfconsistent(md),
|
---|
6 |
|
---|
7 | %tolerance we use in litmus checks for the consistency of the model
|
---|
8 | tolerance=10^-12;
|
---|
9 | %check usage {{{1
|
---|
10 | if nargin~=1,
|
---|
11 | help ismodelselfconsistent
|
---|
12 | error('ismodelselfconsistent error message: wrong usage');
|
---|
13 | end
|
---|
14 | %}}}
|
---|
15 |
|
---|
16 | %recursive call for TRANSIENT {{{1
|
---|
17 | if (md.analysis_type==Transient2DSolutionEnum | md.analysis_type==Transient3DSolutionEnum),
|
---|
18 | if md.dt<=0,
|
---|
19 | error('model not consistent: field dt must be positive for a transient run')
|
---|
20 | end
|
---|
21 |
|
---|
22 | %recursive call to ismodelselfconsistent
|
---|
23 | if (md.dim==2),
|
---|
24 | analysis=[DiagnosticSolutionEnum PrognosticSolutionEnum];
|
---|
25 | else
|
---|
26 | analysis=[DiagnosticSolutionEnum PrognosticSolutionEnum ThermalSolutionEnum];
|
---|
27 | end
|
---|
28 |
|
---|
29 | for i=1:length(analysis),
|
---|
30 | md.analysis_type=analysis(i);
|
---|
31 | ismodelselfconsistent(md);
|
---|
32 | end
|
---|
33 | end
|
---|
34 | %}}}
|
---|
35 |
|
---|
36 | % Common checks
|
---|
37 | %COUNTER {{{1
|
---|
38 | if md.counter<3,
|
---|
39 | error(['model ' md.name ' is not correctly configured. You forgot one step in the following sequence (mesh, geography, parameterize,setelementstype)!']);
|
---|
40 | end
|
---|
41 | %}}}
|
---|
42 | %NAME{{{1
|
---|
43 | if isempty(md.name),
|
---|
44 | error(['model is not correctly configured: missing name!']);
|
---|
45 | end
|
---|
46 | %}}}
|
---|
47 | %ELEMENTS{{{1
|
---|
48 | fields={'elements'};
|
---|
49 | if (md.dim==2),
|
---|
50 | checksize(md,fields,[md.numberofelements 3]);
|
---|
51 | else
|
---|
52 | checksize(md,fields,[md.numberofelements 6]);
|
---|
53 | end
|
---|
54 | %}}}
|
---|
55 | %ELEMENTSTYPE{{{1
|
---|
56 | %Check the size of elements_type
|
---|
57 | fields={'elements_type'};
|
---|
58 | checksize(md,fields,[md.numberofelements 2]);
|
---|
59 | %Check the values of elements_type(1)
|
---|
60 | checkvalues(md,{'elements_type(:,1)'},[MacAyealFormulationEnum() HutterFormulationEnum() PattynFormulationEnum()]);
|
---|
61 | %Check the values of elements_type(2)
|
---|
62 | checkvalues(md,{'elements_type(:,2)'},[StokesFormulationEnum() NoneFormulationEnum()]);
|
---|
63 | if (md.dim==2),
|
---|
64 | checkvalues(md,{'elements_type(:,1)'},[MacAyealFormulationEnum() HutterFormulationEnum()]);
|
---|
65 | end
|
---|
66 | if (md.ismacayealpattyn==0 && md.ishutter==0 && md.isstokes==0),
|
---|
67 | error(['model not consistent: no elements type set for this model. at least one of ismacayealpattyn, ishutter and isstokes need to be =1']);
|
---|
68 | end
|
---|
69 | %}}}
|
---|
70 | %VERTICESTYPE{{{1
|
---|
71 | %Check the size of verticess_type
|
---|
72 | fields={'vertices_type'};
|
---|
73 | checksize(md,fields,[md.numberofgrids 2]);
|
---|
74 | %Check the values of elements_type(1)
|
---|
75 | checkvalues(md,{'elements_type(:,1)'},[MacAyealFormulationEnum() HutterFormulationEnum() PattynFormulationEnum() MacAyealPattynFormulationEnum()]);
|
---|
76 | %Check the values of elements_type(2)
|
---|
77 | checkvalues(md,{'elements_type(:,2)'},[StokesFormulationEnum() NoneFormulationEnum()]);
|
---|
78 | if (md.dim==2),
|
---|
79 | checkvalues(md,{'elements_type(:,1)'},[MacAyealFormulationEnum() HutterFormulationEnum()]);
|
---|
80 | end
|
---|
81 | if (md.ismacayealpattyn==0 && md.ishutter==0 && md.isstokes==0),
|
---|
82 | error(['model not consistent: no elements type set for this model. at least one of ismacayealpattyn, ishutter and isstokes need to be =1']);
|
---|
83 | end
|
---|
84 | %}}}
|
---|
85 | %DG {{{1
|
---|
86 | if md.prognostic_DG==1;
|
---|
87 | %CHECK THE COMPATIBILITY OF THE EDGES
|
---|
88 | fields={'edges'};
|
---|
89 | checksize(md,fields,[NaN 4]);
|
---|
90 | checknan(md,fields);
|
---|
91 | end
|
---|
92 | %}}}
|
---|
93 | %PRESSURELOAD{{{1
|
---|
94 | if (md.dim==2),
|
---|
95 | fields={'pressureload'};
|
---|
96 | checksize(md,fields,[NaN 4]);
|
---|
97 | elseif md.dim==3,
|
---|
98 | fields={'pressureload'};
|
---|
99 | checksize(md,fields,[NaN 6]);
|
---|
100 | else
|
---|
101 | error('dim should be either 2 3');
|
---|
102 | end
|
---|
103 | checkvalues(md,{'pressureload(:,end)'},[WaterEnum() AirEnum()]);
|
---|
104 | %}}}
|
---|
105 | %NO NAN {{{1
|
---|
106 | fields={'numberofelements','numberofgrids','x','y','z','drag_coefficient','drag_type','drag_p','drag_q',...
|
---|
107 | 'rho_ice','rho_water','rheology_B','elementoniceshelf','surface','thickness','bed','g','lowmem','sparsity','nsteps','maxiter',...
|
---|
108 | 'tolx','np','eps_res','max_nonlinear_iterations','exclusive','rheology_n','gridonbed','gridonsurface','elementonbed','elementonsurface','deltaH','DeltaH','timeacc','timedec'};
|
---|
109 | checknan(md,fields);
|
---|
110 | %}}}}
|
---|
111 | %FIELDS >= 0 {{{1
|
---|
112 | fields={'numberofelements','numberofgrids','elements','drag_coefficient','drag_type','drag_p','drag_q',...
|
---|
113 | 'rho_ice','rho_water','rheology_B','elementoniceshelf','thickness','g','eps_res','max_nonlinear_iterations','eps_rel','eps_abs','nsteps','maxiter','tolx','exclusive',...
|
---|
114 | 'sparsity','lowmem','rheology_n','gridonbed','gridonsurface','elementonbed','elementonsurface','deltaH','DeltaH','timeacc','timedec'};
|
---|
115 | checkgreater(md,fields,0);
|
---|
116 | %}}}
|
---|
117 | %FIELDS > 0 {{{1
|
---|
118 | fields={'numberofelements','numberofgrids','elements','drag_type','drag_p',...
|
---|
119 | 'rho_ice','rho_water','rheology_B','thickness','g','max_nonlinear_iterations','eps_res','eps_rel','eps_abs','maxiter','tolx',...
|
---|
120 | 'sparsity','deltaH','DeltaH','timeacc','timedec'};
|
---|
121 | checkgreaterstrict(md,fields,0);
|
---|
122 | %}}}
|
---|
123 | %SIZE NUMBEROFELEMENTS {{{1
|
---|
124 | fields={'drag_p','drag_q','elementoniceshelf','rheology_n','elementonbed'};
|
---|
125 | checksize(md,fields,[md.numberofelements 1]);
|
---|
126 | %}}}
|
---|
127 | %SIZE NUMBEROFGRIDS {{{1
|
---|
128 | fields={'x','y','z','rheology_B','drag_coefficient','melting_rate','accumulation_rate','surface','thickness','bed','gridonbed','gridonsurface'};
|
---|
129 | checksize(md,fields,[md.numberofgrids 1]);
|
---|
130 | %}}}
|
---|
131 | %OTHER SIZES {{{1
|
---|
132 | fields={'spcvelocity'};
|
---|
133 | checksize(md,fields,[md.numberofgrids 6]);
|
---|
134 | %}}}
|
---|
135 | %THICKNESS = SURFACE - BED {{{1
|
---|
136 | if any((md.thickness-md.surface+md.bed)>tolerance),
|
---|
137 | error(['model not consistent: model ' md.name ' violates the equality thickness=surface-bed!']);
|
---|
138 | end
|
---|
139 | %}}}
|
---|
140 | %RIFTS{{{1
|
---|
141 | if md.numrifts,
|
---|
142 | if ~(md.dim==2),
|
---|
143 | error(['model not consistent: models with rifts are only supported in 2d for now!']);
|
---|
144 | end
|
---|
145 | if ~isstruct(md.rifts),
|
---|
146 | error(['model not consistent: md.rifts should be a structure!']);
|
---|
147 | end
|
---|
148 | if ~isempty(find(md.segmentmarkers>=2)),
|
---|
149 | %We have segments with rift markers, but no rift structure!
|
---|
150 | error(['model not consistent: model ' md.name ' should be processed for rifts (run meshprocessrifts)!']);
|
---|
151 | end
|
---|
152 | %Check that rifts are filled with proper material
|
---|
153 | checkvalues(md,{'rifts.fill'},[WaterEnum() AirEnum() IceEnum() MelangeEnum()]);
|
---|
154 | else
|
---|
155 | if ~isnans(md.rifts),
|
---|
156 | error(['model not consistent: md.rifts shoud be NaN since md.numrifts is 0!']);
|
---|
157 | end
|
---|
158 | end
|
---|
159 | %}}}
|
---|
160 | %FLAGS (0 or 1){{{1
|
---|
161 | if ~ismember(md.artificial_diffusivity,[0 1]),
|
---|
162 | error('model not consistent: artificial_diffusivity should be a scalar (1 or 0)');
|
---|
163 | end
|
---|
164 | if ~ismember(md.prognostic_DG,[0 1]),
|
---|
165 | error('model not consistent: prognostic_DG should be a scalar (1 or 0)');
|
---|
166 | end
|
---|
167 | if ~ismember(md.lowmem,[0 1]),
|
---|
168 | error(['model not consistent: model ' md.name ' lowmem field should be 0 or 1']);
|
---|
169 | end
|
---|
170 | %}}}
|
---|
171 | %PARAMETEROUTPUT {{{1
|
---|
172 | if md.numoutput~=length(md.parameteroutput),
|
---|
173 | error('model not consistent: numoutput should be the same size as parameteroutput');
|
---|
174 | end
|
---|
175 | %}}}
|
---|
176 | %CONNECTIVITY {{{1
|
---|
177 | if (md.dim==2),
|
---|
178 | if md.connectivity<9,
|
---|
179 | error('model not consistent: connectivity should be at least 9 for 2d models');
|
---|
180 | end
|
---|
181 | end
|
---|
182 | if md.dim==3,
|
---|
183 | if md.connectivity<24,
|
---|
184 | error('model not consistent: connectivity should be at least 24 for 3d models');
|
---|
185 | end
|
---|
186 | end
|
---|
187 | %}}}
|
---|
188 | %PARALLEL{{{1
|
---|
189 | if ~strcmpi(md.cluster,'none'),
|
---|
190 |
|
---|
191 | %NAN VALUES
|
---|
192 | fields={'time','np'};
|
---|
193 | checknan(md,fields);
|
---|
194 |
|
---|
195 | %FIELD > 0
|
---|
196 | fields={'time','np'};
|
---|
197 | checkgreaterstrict(md,fields,0);
|
---|
198 |
|
---|
199 | end
|
---|
200 | %}}}
|
---|
201 |
|
---|
202 | % Solution checks
|
---|
203 | %DIAGNOSTIC{{{1
|
---|
204 | if md.analysis_type==DiagnosticSolutionEnum,
|
---|
205 |
|
---|
206 | %HUTTER ON ICESHELF WARNING
|
---|
207 | if any(md.elements_type(:,1)==HutterFormulationEnum & md.elementoniceshelf),
|
---|
208 | disp(sprintf('\n !!! Warning: Hutter''s model is not consistent on ice shelves !!!\n'));
|
---|
209 | end
|
---|
210 |
|
---|
211 | %SINGULAR
|
---|
212 | if ~any(sum(md.spcvelocity(:,1:2),2)==2),
|
---|
213 | error(['model not consistent: model ' md.name ' is not well posed (singular). You need at least one grid with fixed velocity!'])
|
---|
214 | end
|
---|
215 |
|
---|
216 | %DIRICHLET IF THICKNESS <= 0
|
---|
217 | if any(md.thickness<=0),
|
---|
218 | pos=find(md.thickness<=0);
|
---|
219 | if any(find(md.spcthickness(pos,1)==0)),
|
---|
220 | error(['model not consistent: model ' md.name ' has some grids with 0 thickness']);
|
---|
221 | end
|
---|
222 | end
|
---|
223 |
|
---|
224 | %INITIAL VELOCITY
|
---|
225 | if length(md.vx)==md.numberofgrids & length(md.vy)==md.numberofgrids,
|
---|
226 | fields={'vx','vy'};
|
---|
227 | checknan(md,fields);
|
---|
228 | end
|
---|
229 | end
|
---|
230 | %}}}
|
---|
231 | %PROGNOSTIC{{{1
|
---|
232 | if md.analysis_type==PrognosticSolutionEnum,
|
---|
233 |
|
---|
234 | %INITIAL VELOCITIES
|
---|
235 | fields={'vx','vy'};
|
---|
236 | checksize(md,fields,[md.numberofgrids 1]);
|
---|
237 | checknan(md,fields);
|
---|
238 |
|
---|
239 | %CHECK THAT WE ARE NOT FULLY CONSTRAINED
|
---|
240 | if (md.dim==2),
|
---|
241 | if isempty(find(~md.spcthickness(:,1))),
|
---|
242 | error(['model not consistent: model ' md.name ' is totally constrained for prognostic, no need to solve!']);
|
---|
243 | end
|
---|
244 | end
|
---|
245 |
|
---|
246 | end
|
---|
247 | %}}}
|
---|
248 | %STEADYSTATE{{{1
|
---|
249 | if md.analysis_type==SteadystateSolutionEnum,
|
---|
250 |
|
---|
251 | %NDT
|
---|
252 | if md.dt~=0,
|
---|
253 | error(['model not consistent: for a steadystate computation, dt must be zero.']);
|
---|
254 | end
|
---|
255 |
|
---|
256 | %PRESSURE
|
---|
257 | if isnans(md.pressure),
|
---|
258 | error(['model not consistent: for a steadystate computation, the model must have an initial pressure, even lithostatic will do.']);
|
---|
259 | end
|
---|
260 |
|
---|
261 | %eps:
|
---|
262 | if isnan(md.eps_rel),
|
---|
263 | error(['model not consistent: for a steadystate computation, eps_rel (relative convergence criterion) must be defined!']);
|
---|
264 | end
|
---|
265 |
|
---|
266 | %dim:
|
---|
267 | if (md.dim==2),
|
---|
268 | error(['model not consistent: for a steadystate computation, model needs to be 3d']);
|
---|
269 | end
|
---|
270 | end
|
---|
271 | %}}}
|
---|
272 | %THERMAL {{{1
|
---|
273 | %THERMAL STEADY AND THERMAL TRANSIENT
|
---|
274 | if md.analysis_type==ThermalSolutionEnum,
|
---|
275 |
|
---|
276 | %EXTRUSION
|
---|
277 | if (md.dim==2),
|
---|
278 | error(['model not consistent: for a ' md.analysis_type ' computation, the model must be 3d, extrude it first!'])
|
---|
279 | end
|
---|
280 |
|
---|
281 | %CHECK THAT WE ARE NOT FULLY CONSTRAINED
|
---|
282 | if isempty(find(~md.spctemperature(:,1))),
|
---|
283 | error(['model not consistent: model ' md.name ' is totally constrained for temperature, no need to solve!']);
|
---|
284 | end
|
---|
285 |
|
---|
286 | %VELOCITIES AND PRESSURE
|
---|
287 | fields={'vx','vy','vz','pressure'};
|
---|
288 | checksize(md,fields,[md.numberofgrids 1]);
|
---|
289 | checknan(md,fields);
|
---|
290 |
|
---|
291 | end
|
---|
292 |
|
---|
293 | %THERMAL TRANSIENT
|
---|
294 | if md.analysis_type==ThermalSolutionEnum & md.dt~=0,
|
---|
295 |
|
---|
296 | %DT and NDT
|
---|
297 | fields={'dt','ndt'};
|
---|
298 | checkgreaterstrict(md,fields,0);
|
---|
299 |
|
---|
300 | %INITIAL TEMPERATURE, MELTING AND ACCUMULATION
|
---|
301 | fields={'temperature','accumulation_rate','melting_rate'};
|
---|
302 | checksize(md,fields,[md.numberofgrids 1]);
|
---|
303 | checknan(md,fields);
|
---|
304 |
|
---|
305 | %INITIAL TEMPERATURE
|
---|
306 | fields={'temperature','spctemperature(:,2)','observed_temperature'};
|
---|
307 | checkgreater(md,fields,0)
|
---|
308 |
|
---|
309 | end
|
---|
310 | %}}}
|
---|
311 | %BALANCEDTHICKNESS{{{1
|
---|
312 | if md.analysis_type==BalancedthicknessSolutionEnum
|
---|
313 |
|
---|
314 | %VELOCITIES MELTING AND ACCUMULATION
|
---|
315 | fields={'vx','vy','accumulation_rate','melting_rate','dhdt'};
|
---|
316 | checksize(md,fields,[md.numberofgrids 1]);
|
---|
317 | checknan(md,fields);
|
---|
318 |
|
---|
319 | %SPC
|
---|
320 | if ~md.prognostic_DG,
|
---|
321 | if any(md.spcthickness(find(md.gridonboundary))~=1),
|
---|
322 | error(['model not consistent: model ' md.name ' should have all the nodes on boundary constrained in field spcthickness']);
|
---|
323 | end
|
---|
324 | end
|
---|
325 | end
|
---|
326 | %}}}
|
---|
327 | %BALANCEDVELOCITIES{{{1
|
---|
328 | if md.analysis_type==BalancedvelocitiesSolutionEnum
|
---|
329 |
|
---|
330 | %VELOCITIES MELTING AND ACCUMULATION
|
---|
331 | fields={'vx','vy','accumulation_rate','melting_rate'};
|
---|
332 | checksize(md,fields,[md.numberofgrids 1]);
|
---|
333 | checknan(md,fields);
|
---|
334 |
|
---|
335 | %SPC
|
---|
336 | if any(md.spcvelocity(find(md.gridonboundary),[1:2])~=1),
|
---|
337 | error(['model not consistent: model ' md.name ' should have all the nodes on boundary constrained in field spcvelocity']);
|
---|
338 | end
|
---|
339 | end
|
---|
340 | %}}}
|
---|
341 | %CONTROL{{{1
|
---|
342 | if md.control_analysis,
|
---|
343 |
|
---|
344 | %CONTROL TYPE
|
---|
345 | if ~ischar(md.control_type),
|
---|
346 | error('model not consistent: control_type should be a string');
|
---|
347 | end
|
---|
348 | if ~(strcmpi(md.control_type,'rheology_B') | strcmpi(md.control_type,'drag_coefficient')),
|
---|
349 | error('model not consistent: control_type should be rheology_B or drag_coefficient');
|
---|
350 | end
|
---|
351 |
|
---|
352 | %LENGTH CONTROL FIELDS
|
---|
353 | fields={'maxiter','optscal','fit','cm_jump'};
|
---|
354 | checksize(md,fields,[md.nsteps 1]);
|
---|
355 |
|
---|
356 | %FIT
|
---|
357 | checkvalues(md,{'fit'},[0 1 2 3 4]);
|
---|
358 |
|
---|
359 | %WEIGHTS
|
---|
360 | fields={'weights'};
|
---|
361 | checksize(md,fields,[md.numberofgrids 1]);
|
---|
362 | checkgreater(md,fields,0);
|
---|
363 |
|
---|
364 | %OBSERVED VELOCITIES
|
---|
365 | fields={'vx_obs','vy_obs'};
|
---|
366 | checksize(md,fields,[md.numberofgrids 1]);
|
---|
367 | checknan(md,fields);
|
---|
368 |
|
---|
369 | %DIRICHLET IF THICKNESS <= 0
|
---|
370 | if any(md.thickness<=0),
|
---|
371 | pos=find(md.thickness<=0);
|
---|
372 | if any(find(md.spcthickness(pos,1)==0)),
|
---|
373 | error(['model not consistent: model ' md.name ' has some grids with 0 thickness']);
|
---|
374 | end
|
---|
375 | end
|
---|
376 |
|
---|
377 | %parameters
|
---|
378 | fields={'cm_noisedmp'};
|
---|
379 | checknan(md,fields);
|
---|
380 | end
|
---|
381 | %}}}
|
---|
382 | %QMU {{{1
|
---|
383 | if md.qmu_analysis,
|
---|
384 | if md.qmu_params.evaluation_concurrency~=1,
|
---|
385 | error(['model not consistent: concurrency should be set to 1 when running dakota in library mode']);
|
---|
386 | end
|
---|
387 | if ~isempty(md.part),
|
---|
388 | if numel(md.part)~=md.numberofgrids,
|
---|
389 | error(['model not consistent: user supplied partition for qmu analysis should have size md.numberofgrids x 1 ']);
|
---|
390 | end
|
---|
391 | if find(md.part)>=md.numberofgrids,
|
---|
392 | error(['model not consistent: user supplied partition should be indexed from 0 (c-convention)']);
|
---|
393 | end
|
---|
394 | if min(md.part)~=0,
|
---|
395 | error(['model not consistent: partition vector not indexed from 0 on']);
|
---|
396 | end
|
---|
397 | if max(md.part)>=md.numberofgrids,
|
---|
398 | error(['model not consistent: partition vector cannot have maximum index larger than number of grids']);
|
---|
399 | end
|
---|
400 | if ~isempty(find(md.part<0)),
|
---|
401 | error(['model not consistent: partition vector cannot have values less than 0']);
|
---|
402 | end
|
---|
403 | if ~isempty(find(md.part>=md.npart)),
|
---|
404 | error(['model not consistent: partition vector cannot have values more than md.npart-1']);
|
---|
405 | end
|
---|
406 | if max(md.part)>=md.npart,
|
---|
407 | error(['model not consistent: for qmu analysis, partitioning vector cannot go over npart, number of partition areas']);
|
---|
408 | end
|
---|
409 | end
|
---|
410 | if md.eps_rel>1.1*10^-5,
|
---|
411 | error(['model not consistent: for qmu analysis, eps_rel should be least than 10^-5, 10^-15 being a better value']);
|
---|
412 | end
|
---|
413 | end
|
---|
414 |
|
---|
415 | if strcmpi(md.analysis_type,'qmu'),
|
---|
416 | if ~strcmpi(md.cluster,'none'),
|
---|
417 | if md.waitonlock==0,
|
---|
418 | error(['model is not correctly configured: waitonlock should be activated when running qmu in parallel mode!']);
|
---|
419 | end
|
---|
420 | end
|
---|
421 | end
|
---|
422 | %}}}
|
---|
423 |
|
---|
424 | end
|
---|
425 |
|
---|
426 | %checks additional functions
|
---|
427 | %checklength {{{1
|
---|
428 | function checklength(md,fields,fieldlength)
|
---|
429 | %CHECKSIZE - check length of a field
|
---|
430 | for i=1:length(fields),
|
---|
431 | if length(eval(['md.' fields{i}]))~=fieldlength,
|
---|
432 | error(['model not consistent: field ' fields{i} ' length should be ' num2str(fieldlength)]);
|
---|
433 | end
|
---|
434 | end
|
---|
435 | end
|
---|
436 | %}}}
|
---|
437 | %checksize {{{1
|
---|
438 | function checksize(md,fields,fieldsize)
|
---|
439 | %CHECKSIZE - check size of a field
|
---|
440 | for i=1:length(fields),
|
---|
441 | if isnan(fieldsize(1)),
|
---|
442 | if (size(eval(['md.' fields{i}]),2)~=fieldsize(2)),
|
---|
443 | %LOG warnings:
|
---|
444 | if strcmpi(fields{i},'pressureload'),
|
---|
445 | disp(' ');
|
---|
446 | disp(' If pressureload does not have the right size (one column missing), this is due to a change in md.pressure.');
|
---|
447 | disp(' To update your model, use the following line');
|
---|
448 | disp(' md.pressureload=[md.pressureload WaterEnum*md.elementoniceshelf(md.pressureload(:,end))+AirEnum*md.elementonicesheet(md.pressureload(:,end))];');
|
---|
449 | disp(' ');
|
---|
450 | end
|
---|
451 | error(['model not consistent: field ' fields{i} ' should have ' num2str(fieldsize(2)) ' columns']);
|
---|
452 | end
|
---|
453 | elseif isnan(fieldsize(2)),
|
---|
454 | if (size(eval(['md.' fields{i}]),1)~=fieldsize(1)),
|
---|
455 | error(['model not consistent: field ' fields{i} ' should have ' num2str(fieldsize(1)) ' rows']);
|
---|
456 | end
|
---|
457 | else
|
---|
458 | if ((size(eval(['md.' fields{i}]),1)~=fieldsize(1)) | (size(eval(['md.' fields{i}]),2)~=fieldsize(2)))
|
---|
459 | %LOG warnings:
|
---|
460 | if strcmpi(fields{i},'weights'),
|
---|
461 | disp(' ');
|
---|
462 | disp(' ''weights'' is a new field that should be of length numberofgrids.');
|
---|
463 | disp(' To update your model, use the following line');
|
---|
464 | disp(' md.weights=ones(md.numberofgrids,1)');
|
---|
465 | disp(' ');
|
---|
466 | end
|
---|
467 | error(['model not consistent: field ' fields{i} ' size should be ' num2str(fieldsize(1)) ' x ' num2str(fieldsize(2))]);
|
---|
468 | end
|
---|
469 | end
|
---|
470 | end
|
---|
471 | end
|
---|
472 | %}}}
|
---|
473 | %checknan {{{1
|
---|
474 | function checknan(md,fields)
|
---|
475 | %CHECKNAN - check nan values of a field
|
---|
476 | for i=1:length(fields),
|
---|
477 | if any(isnan(eval(['md.' fields{i}]))),
|
---|
478 | error(['model not consistent: NaN values in field ' fields{i}]);
|
---|
479 | end
|
---|
480 | end
|
---|
481 | end
|
---|
482 | %}}}
|
---|
483 | %checkreal{{{1
|
---|
484 | function checkreal(md,fields)
|
---|
485 | %CHECKREAL - check real values of a field
|
---|
486 | for i=1:length(fields),
|
---|
487 | if any(eval(['~isreal(md.' fields{i} ')'])),
|
---|
488 | error(['model not consistent: complex values in field ' fields{i}]);
|
---|
489 | end
|
---|
490 | end
|
---|
491 | end
|
---|
492 | %}}}
|
---|
493 | %checkgreaterstrict{{{1
|
---|
494 | function checkgreaterstrict(md,fields,lowerbound)
|
---|
495 | %CHECKGREATERSTRICT - check values of a field
|
---|
496 | for i=1:length(fields),
|
---|
497 | if any(eval(['md.' fields{i} '<=' num2str(lowerbound) ])),
|
---|
498 | error(['model not consistent: field ' fields{i} ' should have values stricly above ' num2str(lowerbound)]);
|
---|
499 | end
|
---|
500 | end
|
---|
501 | end
|
---|
502 | %}}}
|
---|
503 | %checkgreater{{{1
|
---|
504 | function checkgreater(md,fields,lowerbound)
|
---|
505 | %CHECKGREATER - check values of a field
|
---|
506 | for i=1:length(fields),
|
---|
507 | if any(eval(['md.' fields{i} '<' num2str(lowerbound) ])),
|
---|
508 | error(['model not consistent: field ' fields{i} ' should have values above ' num2str(lowerbound)]);
|
---|
509 | end
|
---|
510 | end
|
---|
511 | end
|
---|
512 | %}}}
|
---|
513 | %checklessstrict{{{1
|
---|
514 | function checklessstrict(md,fields,upperbound)
|
---|
515 | %CHECKLESSSTRICT - check values of a field
|
---|
516 | for i=1:length(fields),
|
---|
517 | if any(eval(['md.' fields{i} '>=' num2str(upperbound) ])),
|
---|
518 | error(['model not consistent: field ' fields{i} ' should have values stricly below ' num2str(upperbound)]);
|
---|
519 | end
|
---|
520 | end
|
---|
521 | end
|
---|
522 | %}}}
|
---|
523 | %checkless{{{1
|
---|
524 | function checkless(md,fields,upperbound)
|
---|
525 | %CHECKLESS - check values of a field
|
---|
526 | for i=1:length(fields),
|
---|
527 | if any(eval(['md.' fields{i} '>' num2str(upperbound) ])),
|
---|
528 | error(['model not consistent: field ' fields{i} ' should have values below ' num2str(upperbound)]);
|
---|
529 | end
|
---|
530 | end
|
---|
531 | end
|
---|
532 | %}}}
|
---|
533 | %checkvalues {{{1
|
---|
534 | function checkvalues(md,fields,values)
|
---|
535 | %CHECKVALUE - check that a field has specified values
|
---|
536 | for i=1:length(fields),
|
---|
537 | if eval(['any(~ismember( md.' fields{i} ',values))']),
|
---|
538 | error(['model not consistent: field ' fields{i} ' should have values in ' num2str(values)]);
|
---|
539 | end
|
---|
540 | end
|
---|
541 | end
|
---|
542 | %}}}
|
---|