Index: /issm/trunk-jpl/src/m/classes/inversion.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/inversion.m	(revision 15859)
+++ /issm/trunk-jpl/src/m/classes/inversion.m	(revision 15860)
@@ -127,4 +127,11 @@
 			md = checkfield(md,'inversion.max_parameters','size',[md.mesh.numberofvertices num_controls]);
 
+			%Only SSA, HO and FS are supported right now
+			if solution==StressbalanceSolutionEnum()
+				if ~(md.flowequation.isSSA || md.flowequation.isHO || md.flowequation.isFS),
+					md = checkmessage(md,['inversion can only be performed for SSA, HO or FS ice flow models']);
+				end
+			end
+
 			if solution==BalancethicknessSolutionEnum()
 				md = checkfield(md,'inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
Index: /issm/trunk-jpl/src/m/classes/inversion.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/inversion.py	(revision 15859)
+++ /issm/trunk-jpl/src/m/classes/inversion.py	(revision 15860)
@@ -135,4 +135,9 @@
 		md = checkfield(md,'inversion.max_parameters','size',[md.mesh.numberofvertices,num_controls])
 
+		#Only SSA, HO and FS are supported right now
+		if solution==StressbalanceSolutionEnum():
+			if not (md.flowequation.isSSA or md.flowequation.isHO or md.flowequation.isFS):
+				md.checkmessage("'inversion can only be performed for SSA, HO or FS ice flow models");
+
 		if solution==BalancethicknessSolutionEnum():
 			md = checkfield(md,'inversion.thickness_obs','size',[md.mesh.numberofvertices],'NaN',1)