Index: /issm/trunk-jpl/src/m/miscellaneous/transientrestart.m
===================================================================
--- /issm/trunk-jpl/src/m/miscellaneous/transientrestart.m	(revision 21768)
+++ /issm/trunk-jpl/src/m/miscellaneous/transientrestart.m	(revision 21769)
@@ -1,10 +1,21 @@
-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)];