source: issm/trunk/src/m/solutions/jpl/bedslope.m@ 4113

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

cleanup of serial solutions, to look similar to parallel solutions.

File size: 958 bytes
RevLine 
[4113]1function md=bedslopecompute(md);
2%SLOPECOMPUTE - compute the bed slope of a model
[4102]3%
4% Usage:
[4113]5% md=bedslope(md)
[4102]6%
7 %timing
8 t1=clock;
9
10 analysis_types=SlopeAnalysisEnum;
11 solution_type=SlopeAnalysisEnum;
12
13 displaystring(md.verbose,'%s',['create finite element model']);
14 femmodel=NewFemModel(md,solution_type,analysis_types);
15
16 %retrieve parameters
17 verbose=femmodel.parameters.Verbose;
18 qmu_analysis=femmodel.parameters.QmuAnalysis;
19
20 %compute solution
21 if ~qmu_analysis,
22 displaystring(verbose,'%s',['call computational core']);
[4113]23 femmodel=bedslope_core(femmodel);
24
25 displaystring(verbose,'%s',['write results']);
26 OutputResults(femmodel.elements, femmodel.loads, femmodel.nodes, femmodel.vertices, femmodel.materials, femmodel.parameters);
27
[4102]28 else
29 %launch dakota driver for diagnostic core solution
30 Qmu(femmodel);
31 end
32
33 %stop timing
34 t2=clock;
35 displaystring(md.verbose,'\n%s\n',['solution converged in ' num2str(etime(t2,t1)) ' seconds']);
Note: See TracBrowser for help on using the repository browser.