Changeset 21705

Timestamp:

05/04/17 14:24:09 (8 years ago)

Author:

Mathieu Morlighem

Message:

CHG: allow for layer plot BC

Location:

issm/trunk-jpl/src/m/plot

Files:

• plot_BC.m (modified) (2 diffs)
• plot_icefront.m (modified) (1 diff)
• plot_manager.m (modified) (1 diff)

issm/trunk-jpl/src/m/plot/plot_BC.m ¶

@@ -7 +7 @@
 
 [x y z elements is2d isplanet]=processmesh(md,[],options);
+spcvx=processdata(md,md.stressbalance.spcvx,options);
+spcvy=processdata(md,md.stressbalance.spcvy,options);
+spcvz=processdata(md,md.stressbalance.spcvz,options);
+nbv = numel(x);
 
 %plot dirichlets
@@ -12 +16 @@
 if strcmpi(dirichleton,'on'),
         h1=plot3(...
-                x(find(~isnan(md.stressbalance.spcvx(1:md.mesh.numberofvertices,1)))),...
-                y(find(~isnan(md.stressbalance.spcvx(1:md.mesh.numberofvertices,1)))),...
-                z(find(~isnan(md.stressbalance.spcvx(1:md.mesh.numberofvertices,1)))),...
+                x(find(~isnan(spcvx(1:nbv,1)))),...
+                y(find(~isnan(spcvx(1:nbv,1)))),...
+                z(find(~isnan(spcvx(1:nbv,1)))),...
                 'ro','MarkerSize',14,'MarkerFaceColor','r');
         h2=plot3(...
-                x(find(~isnan(md.stressbalance.spcvy(1:md.mesh.numberofvertices,1)))),...
-                y(find(~isnan(md.stressbalance.spcvy(1:md.mesh.numberofvertices,1)))),...
-                z(find(~isnan(md.stressbalance.spcvy(1:md.mesh.numberofvertices,1)))),...
+                x(find(~isnan(spcvy(1:nbv,1)))),...
+                y(find(~isnan(spcvy(1:nbv,1)))),...
+                z(find(~isnan(spcvy(1:nbv,1)))),...
                 'bo','MarkerSize',10,'MarkerFaceColor','b');
         h3=plot3(...
-                x(find(~isnan(md.stressbalance.spcvz(1:md.mesh.numberofvertices,1)))),...
-                y(find(~isnan(md.stressbalance.spcvz(1:md.mesh.numberofvertices,1)))),...
-                z(find(~isnan(md.stressbalance.spcvz(1:md.mesh.numberofvertices,1)))),...
+                x(find(~isnan(spcvz(1:nbv,1)))),...
+                y(find(~isnan(spcvz(1:nbv,1)))),...
+                z(find(~isnan(spcvz(1:nbv,1)))),...
                 'yo','MarkerSize',6 ,'MarkerFaceColor','y');
 end

issm/trunk-jpl/src/m/plot/plot_icefront.m ¶

@@ -19 +19 @@
 elementzeroice=sum(zeroice(md.mesh.elements),2);
 
-if dimension(md.mesh)==2,
+if is2d,
         icefront=(elementice & elementnoice) & ~(elementice==2 & elementzeroice);
 

issm/trunk-jpl/src/m/plot/plot_manager.m ¶

@@ -21 +21 @@
                         plot_boundaries(md,options,subplotwidth,i);
                         return;
-                case 'BC',
+                case {'BC','bc'},
                         plot_BC(md,options,subplotwidth,i,data);
                         return;