Index: /issm/branches/trunk-larour-NatGeoScience2016/src/m/plot/applyoptions.m
===================================================================
--- /issm/branches/trunk-larour-NatGeoScience2016/src/m/plot/applyoptions.m	(revision 22217)
+++ /issm/branches/trunk-larour-NatGeoScience2016/src/m/plot/applyoptions.m	(revision 22218)
@@ -312,4 +312,9 @@
 if exist(options,'contourlevels'),
 	plot_contour(md,data,options);
+end
+
+%coastlines
+if exist(options,'coastline'),
+	plot_coastlines(md.mesh,options);
 end
 
Index: /issm/branches/trunk-larour-NatGeoScience2016/src/m/plot/plot_coastlines.m
===================================================================
--- /issm/branches/trunk-larour-NatGeoScience2016/src/m/plot/plot_coastlines.m	(revision 22218)
+++ /issm/branches/trunk-larour-NatGeoScience2016/src/m/plot/plot_coastlines.m	(revision 22218)
@@ -0,0 +1,18 @@
+function plot_coastlines(mesh,varargin);
+
+if nargin==2, 
+	if strcmpi(class(varargin{1}),'pairoptions'),
+		options=varargin{1};
+	end
+else
+	%process options: 
+	options=pairoptions(varargin{:});
+end
+
+%retrieve some variables: 
+color=getfieldvalue(options,'coast_color','w');
+
+load ./Data/coastlines
+[x,y,z]=AboveGround(coastlat,coastlon,mesh.r(1),1000);
+hold on, p=plot3(x,y,z,'k-'); set(p,'Color',color);
+
Index: /issm/branches/trunk-larour-NatGeoScience2016/src/m/plot/plot_contour.m
===================================================================
--- /issm/branches/trunk-larour-NatGeoScience2016/src/m/plot/plot_contour.m	(revision 22217)
+++ /issm/branches/trunk-larour-NatGeoScience2016/src/m/plot/plot_contour.m	(revision 22218)
@@ -21,5 +21,5 @@
 if datatype==1,
 	%elements -> take average
-	data=averaging(md,data,0);
+	data=averaging(md,data,0)
 elseif datatype==2,
 	%nodes -> do nothing
@@ -246,5 +246,5 @@
 			if isplanet,
 				if (max(xc)-min(xc)+max(yc)-min(yc)+max(zc)-min(zc))<dist, continue; end
-				h=[h;patch('Xdata',[xc;NaN],'Ydata',[yc;NaN],'Zdata',[zc;NaN],'facecolor','none','linewidth',linewidth)];
+				h=[h;patch('Xdata',[xc;NaN],'Ydata',[yc;NaN],'Zdata',[zc;NaN],'facecolor','none','edgecolor',color,'linewidth',linewidth)];
 				c = horzcat([level, xc'; length(xc), yc'; length(xc), zc']);
 			else