Index: /issm/trunk/src/m/utils/Mesh/zerothickness_icesheetfront.m
===================================================================
--- /issm/trunk/src/m/utils/Mesh/zerothickness_icesheetfront.m	(revision 1247)
+++ /issm/trunk/src/m/utils/Mesh/zerothickness_icesheetfront.m	(revision 1247)
@@ -0,0 +1,26 @@
+function [thickness,surface,bed]=zerothickness_icesheetfront(md,default_value)
+%ZEROTHICKNESS_ICESHEETFRONT   zero the thickness, surface and bed at a grounded icefront 
+%
+% Usage: [thickness,surface,bed]=zerothickness_icesheetfront(md,default)
+% 
+% See also: parameterize.m
+
+%retrieve thickness, surface and bed
+thickness=md.thickness;
+bed=md.bed;
+surface=md.surface;
+
+
+%retrieve elements on ice sheet front
+segments=md.segmentonneumann_diag;
+els=segments(:,3);
+pos=find(md.elementonicesheet(els));
+segments_icesheet=segments(pos,:);
+
+%retrieve grids on ice sheet front
+grid_icesheetfront=intersect(segments_icesheet(:,1),segments_icesheet(:,2));
+
+%modify thickness, surface and bed to be a default value
+thickness(grid_icesheetfront)=default_value;
+bed(grid_icesheetfront)=surface(grid_icesheetfront)-thickness(grid_icesheetfront);
+