Index: /issm/trunk-jpl/src/c/analyses/DebrisAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/DebrisAnalysis.cpp	(revision 27301)
+++ /issm/trunk-jpl/src/c/analyses/DebrisAnalysis.cpp	(revision 27302)
@@ -604,10 +604,10 @@
 
 	int numnodes = element->GetNumberOfNodes();
-	IssmDouble* thickness    = xNew<IssmDouble>(numnodes);
-	IssmDouble* thicknessold    = xNew<IssmDouble>(numnodes);
-	IssmDouble* newthickness = xNew<IssmDouble>(numnodes);
-	IssmDouble* icethickness    = xNew<IssmDouble>(numnodes);
-	IssmDouble* bedslopex    = xNew<IssmDouble>(numnodes);
-	IssmDouble* surfaceslopex    = xNew<IssmDouble>(numnodes);
+	IssmDouble* thickness     = xNew<IssmDouble>(numnodes);
+	IssmDouble* thicknessold  = xNew<IssmDouble>(numnodes);
+	IssmDouble* newthickness  = xNew<IssmDouble>(numnodes);
+	IssmDouble* icethickness  = xNew<IssmDouble>(numnodes);
+	IssmDouble* bedslopex     = xNew<IssmDouble>(numnodes);
+	IssmDouble* surfaceslopex = xNew<IssmDouble>(numnodes);
 
 	/*Use the dof list to index into the solution vector: */
@@ -701,5 +701,5 @@
 					if(isminthicknessinelement){
 						for(k=0; k<numnodes;k++){
-							if(onsurface[k]){
+							if(onsurface[k]>0.5){
 								slope=fabs(slopex[k]);
 								if(dim==2) slope=pow(pow(slopex[k],2)+pow(slopey[k],2),0.5);