1 | function displaydiagnostic(md)
|
---|
2 | %DISPLAYDIAGNOSTIC - display solution parameters of diagnostic
|
---|
3 | %
|
---|
4 | % To avoid clobbering display.m with every field from model md,
|
---|
5 | % we create this routine that displays the solution parameters
|
---|
6 | % from model md, only if requested.
|
---|
7 | %
|
---|
8 | % Usage:
|
---|
9 | % displaydiagnostic(md)
|
---|
10 |
|
---|
11 | disp(sprintf(' Diagnostic solution parameters:'));
|
---|
12 |
|
---|
13 | disp(sprintf('\n %s','Newton convergence criteria:'));
|
---|
14 | fielddisplay(md,'eps_res','mechanical equilibrium residue convergence criterion');
|
---|
15 | fielddisplay(md,'eps_rel','velocity relative convergence criterion, NaN -> not applied');
|
---|
16 | fielddisplay(md,'eps_abs','velocity absolute convergence criterion, NaN -> not applied');
|
---|
17 | fielddisplay(md,'max_nonlinear_iterations','maximum number of nonlinear iterations');
|
---|
18 | fielddisplay(md,'viscosity_overshoot','over-shooting constant new=new+C*(new-old)');
|
---|
19 |
|
---|
20 | disp(sprintf('\n boundary conditions:'));
|
---|
21 | fielddisplay(md,'spcvx','x-axis velocity constraint (NaN means no constraint)');
|
---|
22 | fielddisplay(md,'spcvy','y-axis velocity constraint (NaN means no constraint)');
|
---|
23 | fielddisplay(md,'spcvz','z-axis velocity constraint (NaN means no constraint)');
|
---|
24 | fielddisplay(md,'pressureload','segments on ice front list');
|
---|
25 |
|
---|
26 | disp(sprintf('\n %s','Penalties:'));
|
---|
27 | fielddisplay(md,'penalty_offset','offset used by penalties: penalty = Kmax*10^offset');
|
---|
28 | fielddisplay(md,'min_mechanical_constraints','threshold for instability of mechanical constraints');
|
---|
29 |
|
---|
30 | disp(sprintf('\n %s','Memory management:'));
|
---|
31 | fielddisplay(md,'lowmem','Set to 1 if you are running low on cluster memory');
|
---|
32 | fielddisplay(md,'sparsity','matrix sparsity. Set to .001 for < 1M dof, .0001 for 5M dof, and .00001 for > 10M dof');
|
---|
33 | fielddisplay(md,'connectivity','element connectivity');
|
---|
34 |
|
---|
35 | disp(sprintf('\n %s','Debugging:'));
|
---|
36 | fielddisplay(md,'verbose','level of output statements: none -> 0, yes->1,2. Default is 0');
|
---|
37 | fielddisplay(md,'element_debug','output debug statements for elementswhen possible yes-> 1, no -> 0. Default is 0');
|
---|
38 | fielddisplay(md,'element_debugid','if element_debug on, id of element for which to output messages');
|
---|