Changeset 18157

Timestamp:

06/16/14 23:55:55 (11 years ago)

Author:

Eric.Larour

Message:

CHG: handle 2d meshes.

File:

• issm/trunk-jpl/src/m/plot/plot_BC.m (modified) (1 diff)

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

@@ -6 +6 @@
 hold on
 
+[x y z elements is2d isplanet]=processmesh(md,[],options);
+
 %plot dirichlets
 dirichleton=getfieldvalue(options,'dirichlet','on');
 if strcmpi(dirichleton,'on'),
         h1=plot3(...
-                md.mesh.x(find(~isnan(md.stressbalance.spcvx(1:md.mesh.numberofvertices,1)))),...
-                md.mesh.y(find(~isnan(md.stressbalance.spcvx(1:md.mesh.numberofvertices,1)))),...
-                md.mesh.z(find(~isnan(md.stressbalance.spcvx(1:md.mesh.numberofvertices,1)))),...
+                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)))),...
                 'ro','MarkerSize',14,'MarkerFaceColor','r');
         h2=plot3(...
-                md.mesh.x(find(~isnan(md.stressbalance.spcvy(1:md.mesh.numberofvertices,1)))),...
-                md.mesh.y(find(~isnan(md.stressbalance.spcvy(1:md.mesh.numberofvertices,1)))),...
-                md.mesh.z(find(~isnan(md.stressbalance.spcvy(1:md.mesh.numberofvertices,1)))),...
+                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)))),...
                 'bo','MarkerSize',10,'MarkerFaceColor','b');
         h3=plot3(...
-                md.mesh.x(find(~isnan(md.stressbalance.spcvz(1:md.mesh.numberofvertices,1)))),...
-                md.mesh.y(find(~isnan(md.stressbalance.spcvz(1:md.mesh.numberofvertices,1)))),...
-                md.mesh.z(find(~isnan(md.stressbalance.spcvz(1:md.mesh.numberofvertices,1)))),...
+                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)))),...
                 'yo','MarkerSize',6 ,'MarkerFaceColor','y');
 end