Index: /issm/trunk/src/m/classes/@model/model.m
===================================================================
--- /issm/trunk/src/m/classes/@model/model.m	(revision 5518)
+++ /issm/trunk/src/m/classes/@model/model.m	(revision 5519)
@@ -311,4 +311,5 @@
 	md.qmu_mass_flux_profiles=NaN;
 	md.qmu_mass_flux_segments=NaN;
+	md.qmu_relax=0;
 
 	%partitioner:
Index: /issm/trunk/src/m/classes/public/ismodelselfconsistent.m
===================================================================
--- /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 5518)
+++ /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 5519)
@@ -410,6 +410,8 @@
 		end
 	end
-	if md.eps_rel>1.1*10^-5,
-		error(['model not consistent: for qmu analysis, eps_rel should be least than 10^-5, 10^-15 being a better value']);
+	if ~md.qmu_relax,
+		if md.eps_rel>1.1*10^-5,
+			error(['model not consistent: for qmu analysis, eps_rel should be least than 10^-5, 10^-15 being a better value']);
+		end
 	end
 end