Index: /issm/trunk-jpl/src/m/parameterization/reinitializelevelset.m
===================================================================
--- /issm/trunk-jpl/src/m/parameterization/reinitializelevelset.m	(revision 22452)
+++ /issm/trunk-jpl/src/m/parameterization/reinitializelevelset.m	(revision 22452)
@@ -0,0 +1,16 @@
+function levelsetnew = reinitializelevelset(md,levelset)
+%REINITIALIZELEVELSET - reinitialize levelset as a signed distance function
+%
+%   Usage:
+%      levelsetnew = reinitializelevelset(md,levelset)
+
+
+%First: extract segments
+contours=contourlevelzero(md,levelset,0);
+
+%Now, make this a distance field (might not be closed)
+levelsetnew=abs(ExpToLevelSet(md.mesh.x,md.mesh.y,contours));
+
+%Finally, change sign
+pos = find(levelset<0);
+levelsetnew(pos) = -levelsetnew(pos);