Index: /issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp	(revision 26088)
+++ /issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp	(revision 26089)
@@ -927,5 +927,5 @@
 	xDelete<IssmDouble>(bed);
 	xDelete<int>(doflist);
-	if(domaintype!=Domain2DhorizontalEnum){basalelement->DeleteMaterials(); delete basalelement;};
+	if(domaintype!=Domain2DhorizontalEnum && domaintype!=Domain3DsurfaceEnum){basalelement->DeleteMaterials(); delete basalelement;};
 }/*}}}*/
 void           MasstransportAnalysis::UpdateConstraints(FemModel* femmodel){/*{{{*/
Index: /issm/trunk-jpl/src/c/cores/extrudefromtop_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/extrudefromtop_core.cpp	(revision 26088)
+++ /issm/trunk-jpl/src/c/cores/extrudefromtop_core.cpp	(revision 26089)
@@ -23,4 +23,5 @@
 	/*If this is a 2D horizontal domain, return (no need to extrude)*/
 	if(domaintype==Domain2DhorizontalEnum) return;
+	if(domaintype==Domain3DsurfaceEnum) return;
 
 	/*Special method for Penta, otherwise call solution sequence*/