Index: /issm/trunk/src/m/model/process_solve_options.m
===================================================================
--- /issm/trunk/src/m/model/process_solve_options.m	(revision 6568)
+++ /issm/trunk/src/m/model/process_solve_options.m	(revision 6569)
@@ -40,4 +40,5 @@
 outoptions.qmufile=getfieldvalue(options,'qmufile','qmu');% qmufile cannot be changed unless ????script.sh is also changed
 outoptions.overwrite=getfieldvalue(options,'overwrite','n');
+outoptions.keep=getfieldvalue(options,'keep','n');
 outoptions.ivar=getfieldvalue(options,'ivar',1);
 outoptions.iresp=getfieldvalue(options,'iresp',1);
Index: /issm/trunk/src/m/model/solve.m
===================================================================
--- /issm/trunk/src/m/model/solve.m	(revision 6568)
+++ /issm/trunk/src/m/model/solve.m	(revision 6569)
@@ -57,5 +57,7 @@
 	md=postqmu(md);
 	cd ..
-	system(['rm -rf qmu' num2str(GetPId)]);
+	if ~strncmpi(options.keep,'y',1)
+		system(['rm -rf qmu' num2str(GetPId)]);
+	end
 end
 
Index: /issm/trunk/src/m/model/solveparallel.m
===================================================================
--- /issm/trunk/src/m/model/solveparallel.m	(revision 6568)
+++ /issm/trunk/src/m/model/solveparallel.m	(revision 6569)
@@ -44,5 +44,7 @@
 	%post processes qmu results if necessary
 	if md.qmu_analysis,
-		system(['rm -rf qmu' num2str(GetPId)]);
+		if ~strncmpi(options.keep,'y',1)
+			system(['rm -rf qmu' num2str(GetPId)]);
+		end
 	end
 else