source: issm/trunk/src/m/model/display/displaydiagnostic.m@ 5901

Last change on this file since 5901 was 3201, checked in by Eric.Larour, 15 years ago

Added freezing of constraints for rifts that zigzag.
Also added max_nonlinear_iterations, to control how many iterations can be carried out in Newton iterations.
Add min_mechanical_constraints, acts like min_thermal_constraints, but for rift constraints.

File size: 1.9 KB
Line 
1function 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
11disp(sprintf(' Diagnostic solution parameters:'));
12
13disp(sprintf('\n %s','Newton convergence criteria:'));
14fielddisplay(md,'eps_res','mechanical equilibrium residue convergence criterion');
15fielddisplay(md,'eps_rel','velocity relative convergence criterion, NaN -> not applied');
16fielddisplay(md,'eps_abs','velocity absolute convergence criterion, NaN -> not applied');
17fielddisplay(md,'max_nonlinear_iterations','maximum number of nonlinear iterations');
18fielddisplay(md,'viscosity_overshoot','over-shooting constant new=new+C*(new-old)');
19
20disp(sprintf('\n boundary conditions:'));
21fielddisplay(md,'spcvelocity','constraints flag list (first 3 columns) and values [m/yr] (last 3 columns)');
22fielddisplay(md,'pressureload','segments on ice front list');
23
24disp(sprintf('\n %s','Penalties:'));
25fielddisplay(md,'penalty_offset','offset used by penalties: penalty = Kmax*10^offset');
26fielddisplay(md,'min_mechanical_constraints','threshold for instability of mechanical constraints');
27
28disp(sprintf('\n %s','Memory management:'));
29fielddisplay(md,'lowmem','Set to 1 if you are running low on cluster memory');
30fielddisplay(md,'sparsity','matrix sparsity. Set to .001 for < 1M dof, .0001 for 5M dof, and .00001 for > 10M dof');
31fielddisplay(md,'connectivity','element connectivity');
32
33disp(sprintf('\n %s','Debugging:'));
34fielddisplay(md,'verbose','level of output statements: none -> 0, yes->1,2. Default is 0');
35fielddisplay(md,'element_debug','output debug statements for elementswhen possible yes-> 1, no -> 0. Default is 0');
36fielddisplay(md,'element_debugid','if element_debug on, id of element for which to output messages');
Note: See TracBrowser for help on using the repository browser.