Index: /issm/trunk-jpl/src/m/plot/plot_manager.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_manager.m	(revision 25037)
+++ /issm/trunk-jpl/src/m/plot/plot_manager.m	(revision 25038)
@@ -202,4 +202,9 @@
 end
 
+% clear subplot option
+if exist(options,'subplot')
+    cla;
+end
+
 %plot unit
 plot_unit(x,y,z,elements,data2,is2d,isplanet,datatype,options);
Index: /issm/trunk-jpl/src/m/plot/plotmodel.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plotmodel.m	(revision 25037)
+++ /issm/trunk-jpl/src/m/plot/plotmodel.m	(revision 25038)
@@ -28,9 +28,24 @@
 end
 
+% Add option for subplot, only support one subplot per one plotmodel call
+subindex = 0;
+if (exist(options.list{1},'subplot'))
+    subops = getfieldvalue(options.list{1},'subplot',NaN);
+    nlines = subops(1);
+    ncols = subops(2);
+    subindex = subops(3);
+    % change the width
+    subplotwidth = ncols;
+    % Turn off new figure option for subplot
+    options.list{1}=addfield(options.list{1},'figurestatement','off');
+    % Turn off clf option for subplot
+    options.list{1}=addfield(options.list{1},'clf','off');
+end
+
 %go through subplots
 if numberofplots,
 
 	%Create figure 
-	if strcmpi(getfieldvalue(options.list{1},'figurestatement','on'),'on'),
+	if (strcmpi(getfieldvalue(options.list{1},'figurestatement','on'),'on'))
 		f=figure(figurenumber); 
 	else
@@ -95,6 +110,10 @@
 	try,
 		for i=1:numberofplots,
-			plot_manager(getfieldvalue(options.list{i},'model',md),options.list{i},subplotwidth,nlines,ncols,i);
-			%List all unused options
+            if subindex
+                plot_manager(getfieldvalue(options.list{i},'model',md),options.list{i},subplotwidth,nlines,ncols,subindex);
+            else
+                plot_manager(getfieldvalue(options.list{i},'model',md),options.list{i},subplotwidth,nlines,ncols,i);
+            end
+            %List all unused options
 			if getfieldvalue(options.list{i},'displayunused',1),
 				displayunused(options.list{i})