Index: /issm/trunk-jpl/src/m/plot/applyoptions.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/applyoptions.m	(revision 15055)
+++ /issm/trunk-jpl/src/m/plot/applyoptions.m	(revision 15056)
@@ -440,37 +440,55 @@
 	%get current plot position
 	cplotpos=get(gca,'pos');
-	%compute inset position
-	PosInset=[cplotpos(1)+insetpos(1)*cplotpos(3),cplotpos(2)+insetpos(2)*cplotpos(4), insetpos(3)*cplotpos(3), insetpos(4)*cplotpos(4)];
-	%show pos
+
 	X1=getfieldvalue(options,'insetx',xlim);
 	Y1=getfieldvalue(options,'insety',ylim);
-	if strcmpi(getfieldvalue(options,'showinset','off'),'on')
-		line(X1([1 2 2 1 1]),Y1([1 1 2 2 1]),zeros(1,5),'Color','k','LineWidth',2);
-	end
-
-	%Get current figure
-	ax1=gca;
-
-	%plot inset
-	axes('pos',PosInset);
-	copyobj(get(ax1,'children'),gca);
-	patch('Faces',[1 2 3 4 1],'Vertices',[X1([1 2 2 1])' Y1([1 1 2 2])'],'FaceColor','None','EdgeColor','k','LineWidth',2);
-
-	%applay options
-	options=removefield(options,'text',0);
-	options=removefield(options,'title',0);
-	options=removefield(options,'xlabel',0);
-	options=removefield(options,'ylabel',0);
-	options=removefield(options,'inset',0);
-	options=removefield(options,'offsetaxispos',0);
-	options=removefield(options,'showregion',0);
-	options=changefieldvalue(options,'colorbar',0);
-	options=changefieldvalue(options,'latlon','off');
-	options=changefieldvalue(options,'axis','equal off');
-	options=changefieldvalue(options,'xlim',getfieldvalue(options,'insetx',xlim));
-	options=changefieldvalue(options,'ylim',getfieldvalue(options,'insety',ylim));
-	applyoptions(md,data,options);
-
-	%back to main gca
-	set(gcf,'CurrentAxes',maingca)
-end
+
+	for i=1:length(getfieldvalue(options,'insetx')),
+		if length(insetpos)==4,
+			insetposi=insetpos;
+		else
+			insetposi=insetpos{i};
+		end
+		PosInseti=[cplotpos(1)+insetposi(1)*cplotpos(3),cplotpos(2)+insetposi(2)*cplotpos(4), insetposi(3)*cplotpos(3), insetposi(4)*cplotpos(4)];
+		%show pos
+		if length(X1==2),
+			X1i=X1;
+		else
+			X1i=X1{i};
+		end
+		if length(Y1==2),
+			Y1i=Y1;
+		else
+			Y1i=Y1{i};
+		end
+		if strcmpi(getfieldvalue(options,'showinset','off'),'on')
+			line(X1i([1 2 2 1 1]),Y1i([1 1 2 2 1]),zeros(1,5),'Color','k','LineWidth',2);
+		end
+
+		%Get current figure
+		ax1=gca;
+
+		%plot inset
+		axes('pos',PosInseti);
+		copyobj(get(ax1,'children'),gca);
+		patch('Faces',[1 2 3 4 1],'Vertices',[X1i([1 2 2 1])' Y1i([1 1 2 2])'],'FaceColor','None','EdgeColor','k','LineWidth',2);
+
+		%apply options
+		options=removefield(options,'text',0);
+		options=removefield(options,'title',0);
+		options=removefield(options,'xlabel',0);
+		options=removefield(options,'ylabel',0);
+		options=removefield(options,'inset',0);
+		options=removefield(options,'offsetaxispos',0);
+		options=removefield(options,'showregion',0);
+		options=changefieldvalue(options,'colorbar',0);
+		options=changefieldvalue(options,'latlon','off');
+		options=changefieldvalue(options,'axis','equal off');
+		options=changefieldvalue(options,'xlim',X1i);
+		options=changefieldvalue(options,'ylim',Y1i);
+		applyoptions(md,data,options);
+
+		%back to main gca
+		set(gcf,'CurrentAxes',maingca)
+	end
+end
