Index: ../trunk-jpl/src/m/classes/stressbalance.py
===================================================================
--- ../trunk-jpl/src/m/classes/stressbalance.py	(revision 15615)
+++ ../trunk-jpl/src/m/classes/stressbalance.py	(revision 15616)
@@ -119,11 +119,6 @@
 		md = checkfield(md,'diagnostic.isnewton','numel',[1],'values',[0,1,2])
 		md = checkfield(md,'diagnostic.FSreconditioning','size',[1],'NaN',1)
 		md = checkfield(md,'diagnostic.viscosity_overshoot','size',[1],'NaN',1)
-		if md.mesh.dimension==2:
-			md = checkfield(md,'diagnostic.icefront','size',[float('NaN'),4],'NaN',1)
-		else:
-			md = checkfield(md,'diagnostic.icefront','size',[float('NaN'),6],'NaN',1)
-		md = checkfield(md,'diagnostic.icefront[:,-1]','values',[0,1,2])
 		md = checkfield(md,'diagnostic.maxiter','size',[1],'>=',1)
 		md = checkfield(md,'diagnostic.referential','size',[md.mesh.numberofvertices,6])
 		md = checkfield(md,'diagnostic.loadingforce','size',[md.mesh.numberofvertices,3])