Index: /issm/trunk-jpl/src/m/plot/plot_coastlines.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_coastlines.m	(revision 25057)
+++ /issm/trunk-jpl/src/m/plot/plot_coastlines.m	(revision 25058)
@@ -9892,6 +9892,10 @@
 	set(p,'LineWidth',getfieldvalue(options,'coast_linewidth',1));
 else
-	[x,y,z]=AboveGround(coastlat,coastlon,mesh.r(1),1000);
-	hold on, p=plot3(x,y,z,'k-'); 
+	if ~strcmpi(getfieldvalue(options,'coord','xy'),'latlon'),
+		[x,y,z]=AboveGround(coastlat,coastlon,mesh.r(1),1000);
+		hold on, p=plot3(x,y,z,'k-'); 
+	else
+		hold on, p=plot(coastlon,coastlat,'k-'); 
+	end
 	set(p,'Color',getfieldvalue(options,'coast_color','k'));
 	set(p,'LineWidth',getfieldvalue(options,'coast_linewidth',1));
Index: /issm/trunk-jpl/src/m/plot/processmesh.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/processmesh.m	(revision 25057)
+++ /issm/trunk-jpl/src/m/plot/processmesh.m	(revision 25058)
@@ -24,4 +24,10 @@
 if strcmp(domaintype(md.mesh),'3Dsurface'),
 	[x y z elements is2d isplanet] = processmesh(md.mesh,options);
+	if strcmpi(getfieldvalue(options,'coord','xy'),'latlon'),
+		x=md.mesh.long;
+		y=md.mesh.lat;
+		elements=delaunay(x,y);
+		z=md.mesh.lat; z(:)=0;
+	end
 	return;
 end