Changeset 3061

Timestamp:

02/20/10 15:07:42 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added show region option

Location:

issm/trunk/src/m/classes/public/plot

Files:

• applyoptions.m (modified) (2 diffs)
• plot_overlay.m (modified) (1 diff)
• plotdoc.m (modified) (1 diff)

issm/trunk/src/m/classes/public/plot/applyoptions.m

@@ -6 +6 @@
 %
 %   See also: PLOTMODEL, PARSE_OPTIONS
+
+global ISSM_DIR
 
 %fontsize
@@ -251 +253 @@
 end
 
+%showregion
+if strcmpi(getfieldvalue(options,'showregion','off'),'on'),
+        %get inset relative position (x,y,width,height)
+        insetpos=getfieldvalue(options,'insetpos',[0.02 0.70 0.18 0.18]);
+        %get current plos 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)];
+        axes('pos',PosInset);
+        axis equal off
+        box off
+        A=expread([ ISSM_DIR '/../models/Antarctica/Exp_Par/DomainOutline.exp']);
+        Ax=[min(A.x) max(A.x)];
+        Ay=[min(A.y) max(A.y)];
+        mdx=[min(md.x) max(md.x)];
+        mdy=[min(md.y) max(md.y)];
+        patch([Ax(1)  Ax(2)  Ax(2)  Ax(1)],[Ay(1)  Ay(1)  Ay(2)  Ay(2)],[1 1 1],'edgecolor','none')
+        line(A.x,A.y,'color','b');
+        patch( [mdx(1) mdx(2) mdx(2) mdx(1)],[mdy(1) mdy(1) mdy(2) mdy(2)],ones(4,1),'EdgeColor',[0 0 0],'FaceColor',[0.5 0.5 0.5],'FaceAlpha',0.5)
+        colorbar('off');
+end
 %flag edges of a partition
 if exist(options,'partitionedges')

issm/trunk/src/m/classes/public/plot/plot_overlay.m

@@ -94 +94 @@
 
 %Select plot area 
-%subplot(plotlines,plotcols,i);
+subplot(plotlines,plotcols,i);
+%uncomment following lines to have more space
 %P=get(gca,'pos');
 %P(3)=P(3)+0.05;

issm/trunk/src/m/classes/public/plot/plotdoc.m

@@ -114 +114 @@
 disp('       ''offset'': mesh offset used by ''rifts'', default is 500');
 disp('       ''scaleruler'': add a scale ruler, ''on'' for default value or [x0 y0 length width numberofticks] where (x0,y0) are the coordinates of the lower left corner');
+disp('       ''showregion'': show domain in Antarctica on an inset, use ''insetpos'' properties');
 disp('       ''wrapping'': repeat ''n'' times the colormap (''n'' must be an integer)');
 disp('       ''unit'': by default, in m, otherwise, ''km'' is available');