Index: /issm/trunk/src/m/classes/public/ismodelselfconsistent.m
===================================================================
--- /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 2051)
+++ /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 2052)
@@ -194,4 +194,8 @@
 		end
 
+	end
+	if md.eps_rel>10^-3,
+		disp(['for qmu analysis, eps_rel should be least than 10^-5, 10^-15 being a better value']);
+		bool=0;return;
 	end
 end
Index: /issm/trunk/test/Verification/IceSheetIceFrontM2d_Dakota_25/Square.par
===================================================================
--- /issm/trunk/test/Verification/IceSheetIceFrontM2d_Dakota_25/Square.par	(revision 2051)
+++ /issm/trunk/test/Verification/IceSheetIceFrontM2d_Dakota_25/Square.par	(revision 2052)
@@ -58,2 +58,4 @@
 md.npart=md.numberofgrids;
 md.qmu_analysis=1;
+
+md.eps_rel=10^-15; %tighten for qmu analysese