Index: /issm/trunk-jpl/src/c/analyses/ExtrapolationAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/ExtrapolationAnalysis.cpp	(revision 17427)
+++ /issm/trunk-jpl/src/c/analyses/ExtrapolationAnalysis.cpp	(revision 17428)
@@ -28,4 +28,8 @@
 			counter++;
 		}
+	}
+	
+	if(iomodel->meshtype==Mesh3DEnum){
+		iomodel->FetchDataToInput(elements,MeshElementonbedEnum);
 	}
 }
@@ -199,8 +203,4 @@
 	}/*}}}*/
 
-	//TESTING
-	_printf_("in element: " << basalelement->Id() << "\n");
-	Ke->Echo();
-
 	/*Clean up and return*/
 	xDelete<IssmDouble>(xyz_list);
@@ -209,4 +209,5 @@
 	xDelete<IssmDouble>(dlsf);
 	xDelete<IssmDouble>(normal);
+	if(meshtype!=Mesh2DhorizontalEnum){basalelement->DeleteMaterials(); delete basalelement;};
 	delete gauss;
 	return Ke;
@@ -219,8 +220,9 @@
 
 	/*Intermediaries */
-	int i;
+	int i, meshtype;
 	
 	/*Fetch number of nodes */
 	int numnodes = basalelement->GetNumberOfNodes();
+	basalelement->FindParam(&meshtype,MeshTypeEnum);
 
 	/*Initialize Element vector*/
@@ -228,4 +230,5 @@
 	for(i=0;i<numnodes;i++) 
 		pe->values[i]=0.; 
+	if(meshtype!=Mesh2DhorizontalEnum){basalelement->DeleteMaterials(); delete basalelement;};
 	return pe;
 }/*}}}*/