source: issm/trunk/src/m/solutions/jpl/diagnostic.m@ 4119

Last change on this file since 4119 was 4119, checked in by seroussi, 15 years ago

argument nummodels missing in NewFemModel

File size: 1.3 KB
Line 
1function md=diagnostic(md);
2%DIAGNOSTIC - compute the velocity field of a model
3%
4% Usage:
5% md=diagnostic(md)
6%
7 %timing
8 t1=clock;
9
10 analysis_types=[DiagnosticHorizAnalysisEnum,DiagnosticVertAnalysisEnum,DiagnosticStokesAnalysisEnum,DiagnosticHutterAnalysisEnum,SlopeAnalysisEnum];
11 solution_type=DiagnosticAnalysisEnum;
12
13 displaystring(md.verbose,'%s',['create finite element model']);
14 femmodel=NewFemModel(md,solution_type,analysis_types,5);
15
16 %retrieve parameters
17 verbose=femmodel.parameters.Verbose;
18 qmu_analysis=femmodel.parameters.QmuAnalysis;
19 control_analysis=femmodel.parameters.ControlAnalysis;
20
21 %compute solution
22 if ~qmu_analysis,
23 if ~control_analysis,
24
25 displaystring(verbose,'%s',['call computational core']);
26 femmodel=diagnostic_core(femmodel);
27
28 else,
29
30 displaystring(verbose,'%s',['call computational core']);
31 femmodel=control_core(femmodel);
32
33 end
34
35 displaystring(verbose,'%s',['write results']);
36 OutputResults(femmodel.elements, femmodel.loads, femmodel.nodes, femmodel.vertices, femmodel.materials, femmodel.parameters);
37
38 else
39 %launch dakota driver for diagnostic core solution
40 Qmu(femmodel);
41 end
42
43 %stop timing
44 t2=clock;
45 displaystring(md.verbose,'\n%s\n',['solution converged in ' num2str(etime(t2,t1)) ' seconds']);
Note: See TracBrowser for help on using the repository browser.