Index: /issm/trunk/src/m/model/extrude.m
===================================================================
--- /issm/trunk/src/m/model/extrude.m	(revision 9159)
+++ /issm/trunk/src/m/model/extrude.m	(revision 9160)
@@ -224,4 +224,5 @@
 md.surface=project3d(md,'vector',md.surface,'type','node');
 md.thickness=project3d(md,'vector',md.thickness,'type','node');
+md.bathymetry=project3d(md,'vector',md.bathymetry,'type','node');
 md.thickness_coeff=project3d(md,'vector',md.thickness_coeff,'type','node');
 md.bed=project3d(md,'vector',md.bed,'type','node');
Index: /issm/trunk/src/m/model/ismodelselfconsistent.m
===================================================================
--- /issm/trunk/src/m/model/ismodelselfconsistent.m	(revision 9159)
+++ /issm/trunk/src/m/model/ismodelselfconsistent.m	(revision 9160)
@@ -372,4 +372,7 @@
 
 			%DIRICHLET IF THICKNESS <= 0
+			%Check the size of spcthickness 
+			fields={'spcthickness'};
+			checksize(md,fields,[md.numberofnodes 1]);
 			if any(md.thickness<=0),
 				pos=find(md.thickness<=0);