Changeset 7142

Timestamp:

01/23/11 11:39:42 (14 years ago)

Author:

Eric.Larour

Message:

Added new showbasins and modified call to basinzoom

File:

• issm/trunk/src/m/model/plot/applyoptions.m (modified) (3 diffs)

issm/trunk/src/m/model/plot/applyoptions.m

@@ -89 +89 @@
 %Basinzoom
 if exist(options,'basin');
-        basinzoom(getfieldvalue(options,'basin'),getfieldvalue(options,'unit',1));
-end
+        if strcmpi(md.hemisphere,'n'), hemisphere=+1; else hemisphere=-1; end
+        basinzoom('basin',getfieldvalue(options,'basin'),'hemisphere',getfieldvalue(options,'hemisphere',hemisphere),'basindeltax',getfieldvalue(options,'basindeltax',300000),'basindeltay',getfieldvalue(options,'basindeltay',300000),'unit',getfieldvalue(options,'unit',1));
+end
+
+%ShowBasins
+if strcmpi(getfieldvalue(options,'showbasins','off'),'on')
+        if strcmpi(md.hemisphere,'n'), hemisphere=+1; else hemisphere=-1; end
+        showbasins('hemisphere',getfieldvalue(options,'hemisphere',hemisphere),'fontsize',fontsize);
+end
+
 
 %Caxis
@@ -202 +210 @@
         end
         if exist(options,'colorbartitle'),
-                set(get(c,'title'),'String',getfieldvalue(options,'colorbartitle'),'Color',getfieldvalue(options,'FontColor','k'));
+                set(get(c,'title'),'FontSize',getfieldvalue(options,'colorbarfontsize',fontsize),'String',getfieldvalue(options,'colorbartitle'),'Color',getfieldvalue(options,'FontColor','k'));
         end
         if exist(options,'colorbarwidth'),
@@ -356 +364 @@
         %if we are zooming on a basin, don't take the mesh for the boundaries!
         if exist(options,'basin'),
-                [mdx mdy]=basinzoom(getfieldvalue(options,'basin'));
-        else
-                mdx=[min(md.x) max(md.x)];
-                mdy=[min(md.y) max(md.y)];
-        end
-        line(A.x,A.y,'color','b');
-        patch([Ax(1)  Ax(2)  Ax(2)  Ax(1) Ax(1)],[Ay(1)  Ay(1)  Ay(2)  Ay(2) Ay(1)],[1 1 1],'EdgeColor',[0 0 0],'LineWidth',1,'FaceLighting','none')
-        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','r','FaceAlpha',0.5)
-        colorbar('off');
+                if strcmpi(md.hemisphere,'n'), hemisphere=+1; else hemisphere=-1; end
+                [mdx mdy]=basinzoom('basin',getfieldvalue(options,'basin'),'hemisphere',getfieldvalue(options,'hemisphere',hemisphere),'basindeltax',getfieldvalue(options,'basindeltax',300000),'basindeltay',getfieldvalue(options,'basindeltay',300000),'unit',getfieldvalue(options,'unit',1));
+else
+        mdx=[min(md.x) max(md.x)];
+        mdy=[min(md.y) max(md.y)];
+end
+line(A.x,A.y,'color','b');
+patch([Ax(1)  Ax(2)  Ax(2)  Ax(1) Ax(1)],[Ay(1)  Ay(1)  Ay(2)  Ay(2) Ay(1)],[1 1 1],'EdgeColor',[0 0 0],'LineWidth',1,'FaceLighting','none')
+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','r','FaceAlpha',0.5)
+colorbar('off');
 end
 
 %flag edges of a partition
 if exist(options,'partitionedges')
-        [xsegments ysegments]=flagedges(md.elements,md.x,md.y,md.part);
-        xsegments=xsegments*getfieldvalue(options,'unit',1);
-        ysegments=ysegments*getfieldvalue(options,'unit',1);
-        color=getfieldvalue(options,'partitionedgescolor','r-');
-        linewidth=getfieldvalue(options,'linewidth',1);
-        hold on;
-        for i=1:length(xsegments),
-                plot(xsegments(i,:),ysegments(i,:),color,'LineWidth',linewidth);
-        end
+[xsegments ysegments]=flagedges(md.elements,md.x,md.y,md.part);
+xsegments=xsegments*getfieldvalue(options,'unit',1);
+ysegments=ysegments*getfieldvalue(options,'unit',1);
+color=getfieldvalue(options,'partitionedgescolor','r-');
+linewidth=getfieldvalue(options,'linewidth',1);
+hold on;
+for i=1:length(xsegments),
+        plot(xsegments(i,:),ysegments(i,:),color,'LineWidth',linewidth);
+end
 end