Index: /issm/trunk-jpl/src/c/analyses/DebrisAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/DebrisAnalysis.cpp	(revision 27303)
+++ /issm/trunk-jpl/src/c/analyses/DebrisAnalysis.cpp	(revision 27304)
@@ -170,5 +170,5 @@
 		}else{
 			for(int i=0;i<numvertices;i++){
-				//if(onsurface[i]){
+				//if(onsurface[i]>.5){
 				slope=fabs(slopex[i]);
 				if((atan(slope)*rad2deg)>25.){
@@ -736,5 +736,5 @@
 						int kk=0;
 						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);
@@ -837,5 +837,5 @@
 						test=0;
 						for(k=0; k<numnodes;k++){
-							if(onsurface[k]){
+							if(onsurface[k]>.5){
 								slope=pow(slopex[k]*slopex[k],0.5);
 								if(dim==2) slope=pow(pow(slopex[k],2)+pow(slopey[k],2),0.5);
@@ -859,5 +859,5 @@
 						IssmDouble min=1e14;
 						for(k=0; k<numnodes;k++){
-							if(onsurface[k]){
+							if(onsurface[k]>.5){
 								if(surface[k]<min){
 									index=k;
@@ -867,5 +867,5 @@
 						}
 						for(k=0; k<numnodes;k++){
-							if(onsurface[k]){
+							if(onsurface[k]>.5){
 								if(k==index){
 									debristhickness[k]=debristhickness[k]+debrissum;