source: issm/trunk/src/m/solve/solve.m@ 16137

Last change on this file since 16137 was 16137, checked in by Mathieu Morlighem, 12 years ago

merged trunk-jpl and trunk for revision 16135

File size: 3.6 KB
Line 
1function md=solve(md,solutionenum,varargin)
2%SOLVE - apply solution sequence for this model
3%
4% Usage:
5% md=solve(md,solutionenum,varargin)
6% where varargin is a lit of paired arguments of string OR enums
7%
8% solution types available comprise:
9% - StressbalanceSolutionEnum
10% - MasstransportSolutionEnum
11% - ThermalSolutionEnum
12% - SteadystateSolutionEnum
13% - TransientSolutionEnum...
14% - BalancethicknessSolutionEnum
15% - BedSlopeSolutionEnum
16% - SurfaceSlopeSolutionEnum
17% - HydrologySolutionEnum
18% - FlaimSolutionEnum
19%
20% extra options:
21% - loadonly : does not solve. only load results
22%
23% Examples:
24% md=solve(md,StressbalanceSolutionEnum);
25
26%recover and process solve options
27options=pairoptions(varargin{:},'solution_type',solutionenum);
28options=process_solve_options(options);
29
30%recover some fields
31md.private.solution=options.solution_type;
32cluster=md.cluster;
33
34%check model consistency
35disp('checking model consistency');
36if (solutionenum == FlaimSolutionEnum())
37 md.private.isconsistent=true;
38 md=checkconsistency(md.mesh,md,solutionenum);
39 md=checkconsistency(md.flaim,md,solutionenum);
40 if md.private.isconsistent==false,
41 error('Model not consistent, see messages above');
42 end
43else
44 ismodelselfconsistent(md),
45end
46
47%First, build a runtime name that is unique
48c=clock;
49md.private.runtimename=sprintf('%s-%02i-%02i-%04i-%02i-%02i-%02i-%i',md.miscellaneous.name,c(2),c(3),c(1),c(4),c(5),floor(c(6)),feature('GetPid'));
50
51%if running qmu analysis, some preprocessing of dakota files using models
52%fields needs to be carried out.
53if md.qmu.isdakota,
54 md=preqmu(md,options);
55end
56
57%flaim analysis
58if (options.solution_type == FlaimSolutionEnum())
59 md=flaim_sol(md,options);
60 md.private.solution=EnumToString(options.solution_type);
61 return;
62end
63
64%Do we load results only?
65if options.loadonly,
66 md=loadresultsfromcluster(md);
67 return;
68end
69
70%Write all input files
71marshall(md); % bin file
72ToolkitsFile(md.toolkits,[md.miscellaneous.name '.toolkits']); % toolkits file
73BuildQueueScript(cluster,md.private.runtimename,md.miscellaneous.name,md.private.solution,md.settings.io_gather,md.debug.valgrind,md.debug.gprof); % queue file
74
75%Stop here if batch mode
76if strcmpi(options.batch,'yes')
77 disp('batch mode requested: not launching job interactively');
78 disp('launch solution sequence on remote cluster by hand');
79 return;
80end
81
82%Launch job
83modelname = md.miscellaneous.name;
84filelist = {[modelname '.bin '] [modelname '.toolkits ']};
85if ispc,
86 filelist{end+1}=[modelname '.bat '];
87else
88 filelist{end+1}=[modelname '.queue '];
89end
90
91if md.qmu.isdakota,
92 filelist{end+1} = [modelname '.qmu.in'];
93end
94LaunchQueueJob(cluster,md.miscellaneous.name,md.private.runtimename,filelist);
95
96%did we even try to run? if so, wait on lock
97if strcmpi(options.upload,'on'),
98 disp('solve done uploading test decks');
99 return;
100end
101
102%wait on lock
103if isnan(md.settings.waitonlock),
104 %load when user enters 'y'
105 disp('solution launched on remote cluster. log in to detect job completion.');
106 choice=input('Is the job successfully completed? (y/n)','s');
107 if ~strcmp(choice,'y'),
108 disp('Results not loaded... exiting');
109 else
110 md=loadresultsfromcluster(md);
111 end
112elseif md.settings.waitonlock>0,
113 %we wait for the done file
114 done=waitonlock(md);
115 disp('loading results from cluster');
116 md=loadresultsfromcluster(md);
117elseif md.settings.waitonlock==0,
118 disp('Model results must be loaded manually with md=loadresultsfromcluster(md);');
119end
120
121%post processes qmu results if necessary
122if md.qmu.isdakota,
123 if ~strncmpi(options.keep,'y',1)
124 system(['rm -rf qmu' num2str(feature('GetPid'))]);
125 end
126end
Note: See TracBrowser for help on using the repository browser.