source: issm/trunk/src/m/solutions/jpl/prognostic2.m@ 3534

Last change on this file since 3534 was 3534, checked in by Mathieu Morlighem, 15 years ago

Some bug fix in prognostic2, still not working

File size: 1.6 KB
Line 
1function md=prognostic2(md)
2%PROGNOSITC2 - prognostic2 solution sequence.
3%
4% Usage:
5% md=prognostic2(md)
6
7 %timing
8 t1=clock;
9
10 %Build all models requested for diagnostic simulation
11 models.analysis_type=Prognostic2AnalysisEnum; %needed for processresults
12
13 displaystring(md.verbose,'%s',['reading prognostic2 model data']);
14 md.analysis_type=Prognostic2AnalysisEnum; models.p=CreateFemModel(md);
15
16 % figure out number of dof: just for information purposes.
17 md.dof=modelsize(models);
18
19 %initialize inputs
20 displaystring(md.verbose,'\n%s',['setup inputs...']);
21 inputs=inputlist;
22 inputs=add(inputs,'vx',models.p.parameters.vx_g,'doublevec',1,models.p.parameters.numberofvertices);
23 inputs=add(inputs,'vy',models.p.parameters.vy_g,'doublevec',1,models.p.parameters.numberofvertices);
24 inputs=add(inputs,'thickness',models.p.parameters.h_g,'doublevec',1,models.p.parameters.numberofvertices);
25 inputs=add(inputs,'melting',models.p.parameters.m_g,'doublevec',1,models.p.parameters.numberofvertices);
26 inputs=add(inputs,'accumulation',models.p.parameters.a_g,'doublevec',1,models.p.parameters.numberofvertices);
27 inputs=add(inputs,'dt',models.p.parameters.dt*models.p.parameters.yts,'double');
28
29 displaystring(md.verbose,'\n%s',['call computational core:']);
30 results=prognostic2_core(models,inputs,Prognostic2AnalysisEnum(),NoneAnalysisEnum());
31
32 displaystring(md.verbose,'\n%s',['load results...']);
33 if ~isstruct(md.results), md.results=struct(); end
34 md.results.prognostic2=processresults(models, results);
35
36 %stop timing
37 t2=clock;
38 displaystring(md.verbose,'\n%s\n',['solution converged in ' num2str(etime(t2,t1)) ' seconds']);
Note: See TracBrowser for help on using the repository browser.