Changeset 6321 for issm/trunk/src/m/solutions/diagnostic_core.m
- Timestamp:
- 10/15/10 11:57:42 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/solutions/diagnostic_core.m
r5927 r6321 35 35 if ishutter, 36 36 37 displaystring(verbose,'\n%s',['computing hutter velocities...']);37 issmprintf(verbose,'\n%s',['computing hutter velocities...']); 38 38 39 39 %Take the last velocity into account so that the velocity on the MacAyeal domain is not zero … … 49 49 if xor(ismacayealpattyn,isstokes), %if macayealpattyn and stokes, we must use stokescoupling 50 50 51 displaystring(verbose,'\n%s',['computing horizontal velocities...']);51 issmprintf(verbose,'\n%s',['computing horizontal velocities...']); 52 52 femmodel=SetCurrentConfiguration(femmodel,DiagnosticHorizAnalysisEnum); 53 53 femmodel=solver_diagnostic_nonlinear(femmodel,modify_loads); … … 56 56 if ismacayealpattyn & isstokes, 57 57 58 displaystring(verbose,'\n%s',['computing coupling velocities for macayealpattyn and stokes...']);58 issmprintf(verbose,'\n%s',['computing coupling velocities for macayealpattyn and stokes...']); 59 59 femmodel=SetCurrentConfiguration(femmodel,DiagnosticHorizAnalysisEnum); 60 60 femmodel=solver_stokescoupling_nonlinear(femmodel,conserve_loads); … … 63 63 if dim==3 & (ismacayealpattyn | ishutter) & ~isstokes, 64 64 65 displaystring(verbose,'\n%s',['computing vertical velocities...']);65 issmprintf(verbose,'\n%s',['computing vertical velocities...']); 66 66 femmodel=SetCurrentConfiguration(femmodel,DiagnosticVertAnalysisEnum); 67 67 femmodel=solver_linear(femmodel); … … 69 69 70 70 if (solution_type==DiagnosticSolutionEnum & ~control_analysis) 71 displaystring(verbose,'\n%s',['saving results...']);71 issmprintf(verbose,'\n%s',['saving results...']); 72 72 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum); 73 73 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum);
Note:
See TracChangeset
for help on using the changeset viewer.