Changeset 21759 for issm/branches/trunk-larour-NatGeoScience2016/src/m/plot/plot_BC.m

Timestamp:

06/07/17 10:50:54 (8 years ago)

Author:

Eric.Larour

Message:

CHG: merged branch back to trunk-jpl 21754.

File:

• issm/branches/trunk-larour-NatGeoScience2016/src/m/plot/plot_BC.m (modified) (2 diffs)

issm/branches/trunk-larour-NatGeoScience2016/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