Index: /issm/trunk-jpl/src/m/model/plot/applyoptions.m
===================================================================
--- /issm/trunk-jpl/src/m/model/plot/applyoptions.m	(revision 12308)
+++ /issm/trunk-jpl/src/m/model/plot/applyoptions.m	(revision 12309)
@@ -459,5 +459,5 @@
 	data=getfieldvalue(options,'scatter');
 	hold on
-	plot_scatter(data(:,1),data(:,2),getfieldvalue(options,'scattersize',3),data(:,3));
+	plot_scatter(data(:,1),data(:,2),getfieldvalue(options,'scattersize',3),data(:,3),options);
 end
 
Index: /issm/trunk-jpl/src/m/model/plot/plot_overlay.m
===================================================================
--- /issm/trunk-jpl/src/m/model/plot/plot_overlay.m	(revision 12308)
+++ /issm/trunk-jpl/src/m/model/plot/plot_overlay.m	(revision 12309)
@@ -120,5 +120,5 @@
 
 %Select plot area 
-subplot(plotlines,plotcols,i);
+subplotmodel(plotlines,plotcols,i,options);
 
 %Plot: 
@@ -140,2 +140,3 @@
 options=addfielddefault(options,'axis','equal off');           %default axis
 applyoptions(md,data,options);
+drawnow
Index: /issm/trunk-jpl/src/m/model/plot/subplotmodel.m
===================================================================
--- /issm/trunk-jpl/src/m/model/plot/subplotmodel.m	(revision 12309)
+++ /issm/trunk-jpl/src/m/model/plot/subplotmodel.m	(revision 12309)
@@ -0,0 +1,36 @@
+function ha=subplotmodel(nlines,ncols,num,options);
+%SUBPLOTMODEL -  tight subplot that includes margins
+%
+%   Usage:
+%      h=subplotmodel(nlines,ncols,i,options)
+
+%Regular subplot
+if ~exist(options,'tightsubplot')
+	subplot(nlines,ncols,num);
+	return;
+end
+
+gap     = getfieldvalue(options,'gap',[.01 .01]);
+hmargin = getfieldvalue(options,'hmargin',[.01 .01]);
+vmargin = getfieldvalue(options,'vmargin',[.01 .01]);
+
+
+height = (1-sum(vmargin)-(nlines-1)*gap(1))/nlines; 
+width  = (1-sum(hmargin)-(ncols-1)*gap(2))/ncols;
+ymin   = 1-vmargin(2)-height; 
+
+for i = 1:nlines
+	xmin = hmargin(1);
+	for j = 1:ncols
+		if(((i-1)*ncols+j)==num)
+			ha = axes('Units','normalized', ...
+				'Position',[xmin ymin width height],'XTickLabel','','YTickLabel','','Visible','off');
+			return
+		end
+		xmin = xmin+width+gap(2);
+	end
+	ymin = ymin-height-gap(1);
+end
+
+%Activate new axes
+axes(ha);