Index: /issm/trunk/src/m/classes/public/MeltingGroundingLines.m
===================================================================
--- /issm/trunk/src/m/classes/public/MeltingGroundingLines.m	(revision 2499)
+++ /issm/trunk/src/m/classes/public/MeltingGroundingLines.m	(revision 2499)
@@ -0,0 +1,26 @@
+function md=MeltingGroundingLines(md,distance,value)
+%MELTINGGROUNDINGLINES - set melting near grounding lines to a constant value
+%
+%   Usage:
+%      md=MeltingGroundingLines(md,distance,value)
+%
+
+%get nodes on ice sheet and on ice shelf
+pos_shelf=find(~md.gridonicesheet);
+pos_GL=intersect(unique(md.elements(find(md.elementonicesheet),:)),unique(md.elements(find(md.elementoniceshelf),:)));
+
+for i=1:length(pos_shelf)
+
+	if (mod(i,100)==0),
+		fprintf('\b\b\b\b\b\b\b%5.2f%s',i/length(pos_shelf)*100,' %');
+	end
+
+	%search the grid on ice sheet the closest to i
+	[d posd]=min(sqrt((md.x(pos_shelf(i))-md.x(pos_GL)).^2+(md.y(pos_shelf(i))-md.y(pos_GL)).^2));
+
+	if d<distance,
+
+		md.melting(pos_shelf(i))=value;
+
+	end
+end
Index: /issm/trunk/src/m/classes/public/tres.m
===================================================================
--- /issm/trunk/src/m/classes/public/tres.m	(revision 2498)
+++ /issm/trunk/src/m/classes/public/tres.m	(revision 2499)
@@ -15,5 +15,7 @@
 	end
 	md.vel=md.results.diagnostic.vel;
-	md.pressure=md.results.diagnostic.pressure;
+	if isfield(md.results.diagnostic,'pressure'),
+		md.pressure=md.results.diagnostic.pressure;
+	end
 	if md.numrifts,
 		md.riftproperties=md.results.diagnostic.riftproperties;