Index: ../trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.m
===================================================================
--- ../trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.m	(revision 17474)
+++ ../trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.m	(revision 17475)
@@ -51,7 +51,7 @@
 md.mask.ice_levelset(find(vertexonicefront))=0;
 
 %First find segments that are not completely on the front
-if(md.mesh.dimension==2),
+if ~strcmp(md.mesh.meshtype(),'3D'),
 	numbernodesfront=2;
 else 
 	numbernodesfront=4;
Index: ../trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py
===================================================================
--- ../trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py	(revision 17474)
+++ ../trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py	(revision 17475)
@@ -51,7 +51,7 @@
 	md.mask.ice_levelset[pos]=0
 
 	#First find segments that are not completely on the front
-	if md.mesh.dimension==2:
+	if not strcmp(md.mesh.meshtype(),'3D')
 		numbernodesfront=2
 	else:
 		numbernodesfront=4