Index: /issm/trunk-jpl/src/m/classes/geometry.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/geometry.m	(revision 14810)
+++ /issm/trunk-jpl/src/m/classes/geometry.m	(revision 14811)
@@ -25,13 +25,17 @@
 		end % }}}
 		function md = checkconsistency(obj,md,solution,analyses) % {{{
-
-			md = checkfield(md,'geometry.surface'  ,'NaN',1,'size',[md.mesh.numberofvertices 1]);
-			md = checkfield(md,'geometry.bed'      ,'NaN',1,'size',[md.mesh.numberofvertices 1]);
-			md = checkfield(md,'geometry.thickness','NaN',1,'size',[md.mesh.numberofvertices 1],'>',0);
-			if any((obj.thickness-obj.surface+obj.bed)>10^-9),
-				md = checkmessage(md,['equality thickness=surface-bed violated']);
-			end 
-			if solution==TransientSolutionEnum() & md.transient.isgroundingline,
-				md = checkfield(md,'geometry.bathymetry','NaN',1,'size',[md.mesh.numberofvertices 1]);
+				
+			if (solution==TransientSolutionEnum() & md.transient.isgia) | (solution==GiaSolutionEnum()),
+				md = checkfield(md,'geometry.thickness','forcing',1,'NaN',1,'>=',0);
+			else
+				md = checkfield(md,'geometry.surface'  ,'NaN',1,'size',[md.mesh.numberofvertices 1]);
+				md = checkfield(md,'geometry.bed'      ,'NaN',1,'size',[md.mesh.numberofvertices 1]);
+				md = checkfield(md,'geometry.thickness','NaN',1,'size',[md.mesh.numberofvertices 1],'>',0);
+				if any((obj.thickness-obj.surface+obj.bed)>10^-9),
+					md = checkmessage(md,['equality thickness=surface-bed violated']);
+				end 
+				if solution==TransientSolutionEnum() & md.transient.isgroundingline,
+					md = checkfield(md,'geometry.bathymetry','NaN',1,'size',[md.mesh.numberofvertices 1]);
+				end
 			end
 		end % }}}