Index: /issm/trunk/src/m/classes/public/extrude.m
===================================================================
--- /issm/trunk/src/m/classes/public/extrude.m	(revision 439)
+++ /issm/trunk/src/m/classes/public/extrude.m	(revision 440)
@@ -134,4 +134,9 @@
 	md.gridonpattyn=project3d(md,md.gridonpattyn,'node');
 	md.gridonstokes=project3d(md,md.gridonstokes,'node');
+
+	%dead grids
+	md.deadgrids=ones(md.numberofgrids,1);
+	md.deadgrids(md.elements(md.elements_type(:,1)~=macayealenum,:))=0;%non macayeal grids are not dead
+	md.deadgrids(find(md.gridonbed))=0;%grids from elements on bed are not dead
 end
 
Index: /issm/trunk/src/m/classes/public/ismodelselfconsistent.m
===================================================================
--- /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 439)
+++ /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 440)
@@ -115,5 +115,5 @@
 
 %SIZE NUMBEROFGRIDS
-fields={'x','y','z','B','drag','gridondirichlet_diag','melting','accumulation','surface','thickness','bed','gridonbed','gridonsurface'};
+fields={'x','y','z','B','drag','gridondirichlet_diag','melting','accumulation','surface','thickness','bed','gridonbed','gridonsurface','deadgrids'};
 for i=1:length(fields),
 	if find(length(eval(['md.' char(fields(i))]))~=md.numberofgrids),