Index: /issm/trunk-jpl/src/m/exp/exp2levelsetfunction.m
===================================================================
--- /issm/trunk-jpl/src/m/exp/exp2levelsetfunction.m	(revision 20178)
+++ /issm/trunk-jpl/src/m/exp/exp2levelsetfunction.m	(revision 20179)
@@ -43,5 +43,7 @@
 		w=[segment.w.x segment.w.y];
 		verts=[x y];
-		t=[x-v(1) y-v(2)]*(w-v)'/norm(w-v)^2; %projection of verts on line defined by v and w
+		dist_vw2=norm(w-v)^2;
+		if(dist_vw2==0.),	t=zeros(size(x)); %cover case where segment has length 0
+		else t=[x-v(1) y-v(2)]*(w-v)'/dist_vw2; end %projection of verts on line defined by v and w
 		dist_vec=(ones(length(x),1)*v+max(0,min(1,t))*(w-v))-verts; %vector of shortest distance between verts and segment v-w
 		dist=sqrt(sum(abs(dist_vec).^2,2));