Index: /issm/branches/trunk-larour-SLPS2020/src/m/plot/plot_coastlines.m
===================================================================
--- /issm/branches/trunk-larour-SLPS2020/src/m/plot/plot_coastlines.m	(revision 25586)
+++ /issm/branches/trunk-larour-SLPS2020/src/m/plot/plot_coastlines.m	(revision 25587)
@@ -2,5 +2,6 @@
 
 %define coastline: {{{
-coast=[-83.83 -180 ;...
+coast=[-90 -180; ...
+	-83.83 -180 ;...
 -84.33 -178 ;...
 -84.5 -174 ;...
@@ -614,4 +615,7 @@
 -83.5 176 ;...
 -83.83 180 ;...
+-85 180 ;...
+-87 180 ;...
+-90 180 ;...
 NaN NaN ;...
 -77.58 166.33 ;...
@@ -9866,8 +9870,9 @@
 19.77 -156.05 ;...
 19.95 -155.87 ;...
-20.27 -155.88 ];
+20.27 -155.88 ;...
+NaN NaN];
 
-coastlat=coast(:,1);
-coastlon=coast(:,2);
+coastlat=[coast(:,1); coast(:,1)];
+coastlon=[coast(:,2); coast(:,2)-360];
 
 %}}}
@@ -9902,4 +9907,9 @@
 		p=polyshape(coastlon,coastlat);
 		plot(p,'FaceColor','w','FaceAlpha',1);
+		if strcmpi(getfieldvalue(options,'coordcent','atlantic'),'pacific'),
+			xlim(getfieldvalue(options,'xlim',[-360 0]));
+		else 
+			xlim(getfieldvalue(options,'xlim',[-180 180]));
+		end
 		%p=plot(coastlon,coastlat,'k-'); 
 	end
Index: /issm/branches/trunk-larour-SLPS2020/src/m/plot/processdata.m
===================================================================
--- /issm/branches/trunk-larour-SLPS2020/src/m/plot/processdata.m	(revision 25586)
+++ /issm/branches/trunk-larour-SLPS2020/src/m/plot/processdata.m	(revision 25587)
@@ -42,4 +42,10 @@
 end
 
+%special case for coord latlong: 
+if strcmpi(getfieldvalue(options,'coord','xy'),'latlon') | strcmpi(getfieldvalue(options,'coord','xy'),'latlong'),
+	[data datatype] = processdatalatlong(md,data,options);
+	return;
+end
+
 %needed later on
 if isprop(md.mesh,'numberofvertices2d'), 
Index: /issm/branches/trunk-larour-SLPS2020/src/m/plot/processdatalatlong.m
===================================================================
--- /issm/branches/trunk-larour-SLPS2020/src/m/plot/processdatalatlong.m	(revision 25587)
+++ /issm/branches/trunk-larour-SLPS2020/src/m/plot/processdatalatlong.m	(revision 25587)
@@ -0,0 +1,29 @@
+function [data datatype] = processdatalatlong(md,data,options);
+
+	%   datatype = 1 -> elements
+	%   datatype = 2 -> nodes
+	%what is the mesh we are using: 
+	x0=md.mesh.long;
+	y0=md.mesh.lat;
+
+	%add row at lat=90 and lat=-90
+	add=[(-180:.1:-1e-5)';(1e-5:.1:180)'];
+	nadd=length(add);
+	xextra=[add;add];
+	yextra=[90*ones(nadd,1); -90*ones(nadd,1)];
+	x=[x0;xextra];
+	y=[y0;yextra];
+	elements=delaunay(x,y);
+	
+	%with this mesh, interpolate data: 
+	if length(data)==length(md.mesh.long),
+		datatype=2;
+
+		%interpolate data: 
+		extradata=griddata(x0,y0,data,xextra,yextra,'nearest');
+
+		data=[data; extradata];
+	elseif length(data)==length(md.mesh.elements),
+		error('processdatalatlong error message: coord ''latlong'' case not covered for element data ');
+		datatype=1;
+	end
Index: /issm/branches/trunk-larour-SLPS2020/src/m/plot/processmesh.m
===================================================================
--- /issm/branches/trunk-larour-SLPS2020/src/m/plot/processmesh.m	(revision 25586)
+++ /issm/branches/trunk-larour-SLPS2020/src/m/plot/processmesh.m	(revision 25587)
@@ -24,9 +24,20 @@
 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;
+	if strcmpi(getfieldvalue(options,'coord','xy'),'latlon') | strcmpi(getfieldvalue(options,'coord','xy'),'latlong'),
+		x0=md.mesh.long;
+		y0=md.mesh.lat;
+		%add row at lat=90 and lat=-90
+		add=[(-180:.1:-1e-5)';(1e-5:.1:180)'];
+		nadd=length(add);
+		xextra=[add;add];
+		yextra=[90*ones(nadd,1); -90*ones(nadd,1)];
+		x=[x0;xextra];
+		y=[y0;yextra];
+
+		if strcmpi(getfieldvalue(options,'coordcent','atlantic'),'pacific'),
+			pos=find(x>0);  x(pos)=-360+x(pos);
+		end
 		elements=delaunay(x,y);
-		z=md.mesh.lat; z(:)=0;
+		z=x; z(:)=0;
 	end
 	return;
@@ -42,5 +53,5 @@
 else
 	elements=md.mesh.elements;
-	if ~strcmpi(getfieldvalue(options,'coord','xy'),'latlon') ,
+	if ~strcmpi(getfieldvalue(options,'coord','xy'),'latlon') &  ~strcmpi(getfieldvalue(options,'coord','xy'),'latlong') ,
 		x=md.mesh.x;
 		if isprop(md.mesh,'x2d'), x2d=md.mesh.x2d; end
@@ -50,4 +61,7 @@
 		x=md.mesh.long;
 		y=md.mesh.lat;
+		if strcmpi(getfieldvalue(options,'coordcent','atlantic'),'pacific'),
+			pos=find(x>0);  x(pos)-360+x(pos);
+		end
 	end
 end
Index: /issm/branches/trunk-larour-SLPS2020/src/m/plot/subplotmodel.m
===================================================================
--- /issm/branches/trunk-larour-SLPS2020/src/m/plot/subplotmodel.m	(revision 25586)
+++ /issm/branches/trunk-larour-SLPS2020/src/m/plot/subplotmodel.m	(revision 25587)
@@ -25,4 +25,11 @@
 ymin   = 1-vmargin(2)-height; 
 
+box=getfieldvalue(options,'box','off');
+if strcmpi(box,'on'),
+	visible='on';
+else
+	visible='off';
+end
+
 for i = 1:nlines
 	xmin = hmargin(1);
@@ -30,5 +37,6 @@
 		if(((i-1)*ncols+j)==num)
 			ha = axes('Units','normalized', ...
-				'Position',[xmin ymin width height],'XTickLabel','','YTickLabel','','Visible','off');
+				'Position',[xmin ymin width height],'XTickLabel','','YTickLabel','',...
+				'Visible',visible,'box',box);
 			return
 		end