Changeset 21769

Timestamp:

06/28/17 09:31:22 (8 years ago)

Author:

Mathieu Morlighem

Message:

CHG: optionnaly provide step numer

File:

• issm/trunk-jpl/src/m/miscellaneous/transientrestart.m (modified) (1 diff)

issm/trunk-jpl/src/m/miscellaneous/transientrestart.m

@@ -1 @@
-function md = transientrestart(md)
+function md = transientrestart(md,step)
 %TRANSIENTRESTART - reinitialize model from last transient step
 %
 %   Usage:
 %      md = transientrestart(md)
+%      md = transientrestart(md,step)
+%
+%   By default, transientrestart will use the last step provided in md.results.TransientSolution
 
 %Get result and save it again
-results = md.results.TransientSolution(end);
+if nargin==1,
+        step = numel(md.results.TransientSolution(end));
+end
+if step<1,
+        error('step needs to be >0');
+elseif step>numel(md.results.TransientSolution)
+        error(['md.results.TransientSolution has only ' num2str(numel(md.results.TransientSolution)) ' steps']);
+end
+results = md.results.TransientSolution(step);
 
 newname = ['TransientSolution' num2str(numel(fields(md.results))+1)];