| 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,'spcvelocity','constraints flag list (first 3 columns) and values [m/yr] (last 3 columns)');
|
|---|
| 22 | fielddisplay(md,'pressureload','segments on ice front list');
|
|---|
| 23 |
|
|---|
| 24 | disp(sprintf('\n %s','Penalties:'));
|
|---|
| 25 | fielddisplay(md,'penalty_offset','offset used by penalties: penalty = Kmax*10^offset');
|
|---|
| 26 | fielddisplay(md,'min_mechanical_constraints','threshold for instability of mechanical constraints');
|
|---|
| 27 |
|
|---|
| 28 | disp(sprintf('\n %s','Memory management:'));
|
|---|
| 29 | fielddisplay(md,'lowmem','Set to 1 if you are running low on cluster memory');
|
|---|
| 30 | fielddisplay(md,'sparsity','matrix sparsity. Set to .001 for < 1M dof, .0001 for 5M dof, and .00001 for > 10M dof');
|
|---|
| 31 | fielddisplay(md,'connectivity','element connectivity');
|
|---|
| 32 |
|
|---|
| 33 | disp(sprintf('\n %s','Debugging:'));
|
|---|
| 34 | fielddisplay(md,'verbose','level of output statements: none -> 0, yes->1,2. Default is 0');
|
|---|
| 35 | fielddisplay(md,'element_debug','output debug statements for elementswhen possible yes-> 1, no -> 0. Default is 0');
|
|---|
| 36 | fielddisplay(md,'element_debugid','if element_debug on, id of element for which to output messages');
|
|---|