Changeset 8472

Timestamp:

06/01/11 22:53:48 (14 years ago)

Author:

Eric.Larour

Message:

can now plot a planet shell

Location:

issm/trunk/src/m/model/plot

Files:

• plot_basaldrag.m (modified) (1 diff)
• plot_boundaries.m (modified) (1 diff)
• plot_contour.m (modified) (1 diff)
• plot_drivingstress.m (modified) (2 diffs)
• plot_edges.m (modified) (1 diff)
• plot_elementnumbering.m (modified) (1 diff)
• plot_elementstype.m (modified) (1 diff)
• plot_gridded.m (modified) (1 diff)
• plot_highlightelements.m (modified) (1 diff)
• plot_highlightnodes.m (modified) (1 diff)
• plot_importancefactors.m (modified) (1 diff)
• plot_manager.m (modified) (2 diffs)
• plot_mesh.m (modified) (2 diffs)
• plot_nodenumbering.m (modified) (1 diff)
• plot_overlay.m (modified) (1 diff)
• plot_parthist.m (modified) (1 diff)
• plot_parthistn.m (modified) (1 diff)
• plot_parthistw.m (modified) (1 diff)
• plot_penalties.m (modified) (1 diff)
• plot_pressureload.m (modified) (1 diff)
• plot_qmu_mass_flux_segments.m (modified) (1 diff)
• plot_qmumean.m (modified) (1 diff)
• plot_qmustddev.m (modified) (1 diff)
• plot_riftfraction.m (modified) (1 diff)
• plot_riftnumbering.m (modified) (1 diff)
• plot_riftpenetration.m (modified) (1 diff)
• plot_riftrelvel.m (modified) (1 diff)
• plot_rifts.m (modified) (1 diff)
• plot_riftvel.m (modified) (1 diff)
• plot_section.m (modified) (2 diffs)
• plot_segments.m (modified) (1 diff)
• plot_streamlines.m (modified) (1 diff)
• plot_tensor_components.m (modified) (2 diffs)
• plot_tensor_principal.m (modified) (2 diffs)
• plot_tensor_principalaxis.m (modified) (1 diff)
• plot_transient_field.m (modified) (2 diffs)
• plot_transient_movie.m (modified) (1 diff)
• plot_unit.m (modified) (2 diffs)
• processmesh.m (modified) (2 diffs)

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

@@ -39 +39 @@
 
         %process data and model
-        [x y z elements is2d]=processmesh(md,[],options);
+        [x y z elements is2d isplanet]=processmesh(md,[],options);
         [basal_drag datatype]=processdata(md,drag,options);
 
         %plot basaldrag
         subplot(width,width,i); 
-        plot_unit(x,y,z,elements,basal_drag,is2d,datatype,options);
+        plot_unit(x,y,z,elements,basal_drag,is2d,isplanet,datatype,options);
 
         %apply options

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

@@ -13 +13 @@
         options=removefield(options,'layer',1);
 end
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 for i=1:size(md.segments,1),

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

@@ -8 +8 @@
 
 %process data and model
-[x y z index is2d]=processmesh(md,[],options);
+[x y z index is2d isplanet]=processmesh(md,[],options);
 [data datatype]=processdata(md,datain,options);
 

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

@@ -11 +11 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 [dstress datatype]=processdata(md,s,options);
 dstress=dstress/1000;
@@ -17 +17 @@
 %plot mesh quivervel
 subplot(width,width,i); 
-plot_unit(x,y,z,elements,dstress,is2d,datatype,options)
+plot_unit(x,y,z,elements,dstress,is2d,isplanet,datatype,options)
 
 %apply options

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

@@ -11 +11 @@
 
 %process mesh and data
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 edges=md.edges;
 if isnan(edges)

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

@@ -10 +10 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 [elementnumbers datatype]=processdata(md,[1:md.numberofelements]',options);
 

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

@@ -8 +8 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 [data datatype]=processdata(md,md.elements_type,options);
 

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

@@ -8 +8 @@
 
 %process mesh and data
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 [data datatype]=processdata(md,data,options);
 

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

@@ -11 +11 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 [elementnumbers datatype]=processdata(md,[1:md.numberofelements]',options);
 

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

@@ -8 +8 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 [nodenumbers datatype]=processdata(md,[1:md.numberofnodes]',options);
 

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

@@ -64 +64 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %edgecolor

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

@@ -166 +166 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,data,options);
+[x y z elements is2d isplanet]=processmesh(md,data,options);
 [data2 datatype]=processdata(md,data,options);
 
@@ -178 +178 @@
 
 %plot unit
-plot_unit(x,y,z,elements,data2,is2d,datatype,options);
+plot_unit(x,y,z,elements,data2,is2d,isplanet,datatype,options);
         
 

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

@@ -8 +8 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %plot mesh
@@ -22 +22 @@
         patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','none','EdgeColor',edgecolor,'linewidth',linewidth);
 else
-        A=elements(:,1); B=elements(:,2); C=elements(:,3); D=elements(:,4); E=elements(:,5); F=elements(:,6);
-        patch( 'Faces', [A B C],  'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','none','EdgeColor',edgecolor,'linewidth',linewidth);
-        patch( 'Faces', [D E F],  'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','none','EdgeColor',edgecolor,'linewidth',linewidth);
-        patch( 'Faces', [A B E D],'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','none','EdgeColor',edgecolor,'linewidth',linewidth);
-        patch( 'Faces', [B E F C],'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','none','EdgeColor',edgecolor,'linewidth',linewidth);
-        patch( 'Faces', [C A D F],'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','none','EdgeColor',edgecolor,'linewidth',linewidth);
+        if ~md.planet,
+                A=elements(:,1); B=elements(:,2); C=elements(:,3); D=elements(:,4); E=elements(:,5); F=elements(:,6);
+                patch( 'Faces', [A B C],  'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','none','EdgeColor',edgecolor,'linewidth',linewidth);
+                patch( 'Faces', [D E F],  'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','none','EdgeColor',edgecolor,'linewidth',linewidth);
+                patch( 'Faces', [A B E D],'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','none','EdgeColor',edgecolor,'linewidth',linewidth);
+                patch( 'Faces', [B E F C],'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','none','EdgeColor',edgecolor,'linewidth',linewidth);
+                patch( 'Faces', [C A D F],'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','none','EdgeColor',edgecolor,'linewidth',linewidth);
+        else
+                A=elements(:,1); B=elements(:,2); C=elements(:,3); D=elements(:,4); 
+                patch( 'Faces', [A B C D],  'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','none','EdgeColor',edgecolor,'linewidth',linewidth);
+
+        end
 end
 

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

@@ -8 +8 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 [nodenumbers datatype]=processdata(md,[1:md.numberofnodes]',options);
 

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

@@ -8 +8 @@
 
 %process mesh and data
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 if strcmpi(data,'none'),
         radaronly=1;

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

@@ -8 +8 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %plot mesh

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

@@ -8 +8 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %plot mesh

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

@@ -8 +8 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %plot mesh

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

@@ -8 +8 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %plot mesh penalties

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

@@ -11 +11 @@
 
 %process mesh and data
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 pressureload=md.pressureload;
 

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

@@ -9 +9 @@
 
 %process mesh and data
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 allsegments=md.qmu_mass_flux_segments;

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

@@ -14 +14 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %find response function

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

@@ -14 +14 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %find response function

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

@@ -13 +13 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 subplot(nlines,ncols,index); 

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

@@ -8 +8 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 fontsize=getfieldvalue(options,'FontSize',8);
 

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

@@ -8 +8 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 subplot(nlines,ncols,index); 

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

@@ -30 +30 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 [vel datatype]=processdata(md,[u v],options);
 [quivers,palette]=quiver_process(x,y,vel(:,1),vel(:,2),options);

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

@@ -8 +8 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %plot mesh

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

@@ -26 +26 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 [vel datatype]=processdata(md,[u v],options);
 [quivers,palette]=quiver_process(x,y,vel(:,1),vel(:,2),options);

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

@@ -21 +21 @@
 
 %process model
-[x_m y_m z_m elements_m is2d]=processmesh(md,[],options);
+[x_m y_m z_m elements_m is2d isplanet]=processmesh(md,[],options);
 
 %Get number of curves and generate random colors
@@ -80 +80 @@
                         %plot section only with labels
                         subplot(nlines,ncols,index2)
-                        plot_unit(x_m,y_m,z_m,elements_m,data(:,i),is2d,datatype,options)
+                        plot_unit(x_m,y_m,z_m,elements_m,data(:,i),is2d,isplanet,datatype,options)
                         hold on
                         text(x(1),y(1),'1','backgroundcolor',[0.8 0.9 0.8])

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

@@ -11 +11 @@
 
 %process mesh and data
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 segments=md.segments;
 

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

@@ -6 +6 @@
 
 %process data and model
-[x y z index is2d]=processmesh(md,[],options);
+[x y z index is2d isplanet]=processmesh(md,[],options);
 [u datatype]=processdata(md,md.vx,options);
 [v datatype]=processdata(md,md.vy,options);

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

@@ -24 +24 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 [tensor.xx datatype]=processdata(md,tensor.xx,options);
 [tensor.yy datatype]=processdata(md,tensor.yy,options);
@@ -36 +36 @@
 if ((md.dim==2)),
         subplot(2*width,2*width,index1),
-        plot_unit(x,y,z,elements,tensor.xx,is2d,datatype,options)
+        plot_unit(x,y,z,elements,tensor.xx,is2d,isplanet,datatype,options)
         Apply_options_tensor(md,options,type,'xx')
         subplot(2*width,2*width,index2),
-        plot_unit(x,y,z,elements,tensor.yy,is2d,datatype,options)
+        plot_unit(x,y,z,elements,tensor.yy,is2d,isplanet,datatype,options)
         Apply_options_tensor(md,options,type,'yy')
         subplot(2*width,2*width,index3),
-        plot_unit(x,y,z,elements,tensor.xy,is2d,datatype,options)
+        plot_unit(x,y,z,elements,tensor.xy,is2d,isplanet,datatype,options)
         Apply_options_tensor(md,options,type,'xy')
 else
         subplot(3*width,3*width,index1),
-        plot_unit(x,y,z,elements,tensor.xx,is2d,datatype,options)
+        plot_unit(x,y,z,elements,tensor.xx,is2d,isplanet,datatype,options)
         Apply_options_tensor(md,options,type,'xx')
         subplot(3*width,3*width,index2),
-        plot_unit(x,y,z,elements,tensor.yy,is2d,datatype,options)
+        plot_unit(x,y,z,elements,tensor.yy,is2d,isplanet,datatype,options)
         Apply_options_tensor(md,options,type,'yy')
         subplot(3*width,3*width,index3),
-        plot_unit(x,y,z,elements,tensor.zz,is2d,datatype,options)
+        plot_unit(x,y,z,elements,tensor.zz,is2d,isplanet,datatype,options)
         Apply_options_tensor(md,options,type,'zz')
         subplot(3*width,3*width,index4),
-        plot_unit(x,y,z,elements,tensor.xy,is2d,datatype,options)
+        plot_unit(x,y,z,elements,tensor.xy,is2d,isplanet,datatype,options)
         Apply_options_tensor(md,options,type,'xy')
         subplot(3*width,3*width,index5),
-        plot_unit(x,y,z,elements,tensor.xz,is2d,datatype,options)
+        plot_unit(x,y,z,elements,tensor.xz,is2d,isplanet,datatype,options)
         Apply_options_tensor(md,options,type,'xz')
         subplot(3*width,3*width,index6),
-        plot_unit(x,y,z,elements,tensor.yz,is2d,datatype,options)
+        plot_unit(x,y,z,elements,tensor.yz,is2d,isplanet,datatype,options)
         Apply_options_tensor(md,options,type,'yz')
 end

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

@@ -37 +37 @@
 
 %now plot principal values
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 [tensor.principalvalue1 datatype]=processdata(md,tensor.principalvalue1,options);
 [tensor.principalvalue2 datatype]=processdata(md,tensor.principalvalue2,options);
@@ -46 +46 @@
 if ((md.dim==2)),
         subplot(2*width,2*width,index3)
-        plot_unit(x,y,z,elements,tensor.principalvalue1,is2d,datatype,options)
+        plot_unit(x,y,z,elements,tensor.principalvalue1,is2d,isplanet,datatype,options)
         Apply_options_tensor(md,options,type,'principal value 1')
         subplot(2*width,2*width,index4)
-        plot_unit(x,y,z,elements,tensor.principalvalue2,is2d,datatype,options)
+        plot_unit(x,y,z,elements,tensor.principalvalue2,is2d,isplanet,datatype,options)
         Apply_options_tensor(md,options,type,'principal value 2')
 else
         subplot(3*width,3*width,index4)
-        plot_unit(x,y,z,elements,tensor.principalvalue1,is2d,datatype,options)
+        plot_unit(x,y,z,elements,tensor.principalvalue1,is2d,isplanet,datatype,options)
         Apply_options_tensor(md,options,type,'principal value 1')
         subplot(3*width,3*width,index5)
-        plot_unit(x,y,z,elements,tensor.principalvalue2,is2d,datatype,options)
+        plot_unit(x,y,z,elements,tensor.principalvalue2,is2d,isplanet,datatype,options)
         Apply_options_tensor(md,options,type,'principal value 2')
         subplot(3*width,3*width,index6)
-        plot_unit(x,y,z,elements,tensor.principalvalue3,is2d,datatype,options)
+        plot_unit(x,y,z,elements,tensor.principalvalue3,is2d,isplanet,datatype,options)
         Apply_options_tensor(md,options,type,'principal value 3')
 end

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

@@ -11 +11 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 if ((md.dim==2)),

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

@@ -27 +27 @@
 
 %process mes only once
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %plot data for all steps
@@ -38 +38 @@
         %create plot of step i
         subplot(subplotwidth,subplotwidth,i);
-        plot_unit(x,y,z,elements,data,is2d,datatype,options)
+        plot_unit(x,y,z,elements,data,is2d,isplanet,datatype,options)
         applyoptions(md,data,options);
 

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

@@ -21 +21 @@
 
                 %process data
-                [x y z elements is2d]=processmesh(md,results(i).(field),options);
+                [x y z elements is2d isplanet]=processmesh(md,results(i).(field),options);
                 [data datatype]=processdata(md,results(i).(field),options);
 
                 titlestring=[field ' at time ' num2str(results(i).time/md.yts) ' year'];
-                plot_unit(x,y,z,elements,data,is2d,datatype,options)
+                plot_unit(x,y,z,elements,data,is2d,isplanet,datatype,options)
                 apply_options_movie(md,options,titlestring);
                 

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

@@ -1 @@
-function plot_unit(x,y,z,elements,data,is2d,datatype,options)
+function plot_unit(x,y,z,elements,data,is2d,isplanet,datatype,options)
 %PLOT_UNIT - unit plot, display data
 %
 %   Usage:
-%      plot_unit(x,y,z,elements,data,is2d,datatype,options);
+%      plot_unit(x,y,z,elements,data,is2d,isplanet,datatype,options);
 %
 %   See also: PLOTMODEL, PLOT_MANAGER
@@ -34 +34 @@
                         patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', data(:),'FaceColor','interp','EdgeColor',edgecolor);
                 else
-                        A=elements(:,1); B=elements(:,2); C=elements(:,3); D=elements(:,4); E=elements(:,5); F=elements(:,6);
-                        patch( 'Faces', [A B C],  'Vertices', [x y z],'FaceVertexCData', data(:),'FaceColor','interp','EdgeColor',edgecolor);
-                        patch( 'Faces', [D E F],  'Vertices', [x y z],'FaceVertexCData', data(:),'FaceColor','interp','EdgeColor',edgecolor);
-                        patch( 'Faces', [A B E D],'Vertices', [x y z],'FaceVertexCData', data(:),'FaceColor','interp','EdgeColor',edgecolor);
-                        patch( 'Faces', [B E F C],'Vertices', [x y z],'FaceVertexCData', data(:),'FaceColor','interp','EdgeColor',edgecolor);
-                        patch( 'Faces', [C A D F],'Vertices', [x y z],'FaceVertexCData', data(:),'FaceColor','interp','EdgeColor',edgecolor);
+                        if ~isplanet,
+                                A=elements(:,1); B=elements(:,2); C=elements(:,3); D=elements(:,4); E=elements(:,5); F=elements(:,6);
+                                patch( 'Faces', [A B C],  'Vertices', [x y z],'FaceVertexCData', data(:),'FaceColor','interp','EdgeColor',edgecolor);
+                                patch( 'Faces', [D E F],  'Vertices', [x y z],'FaceVertexCData', data(:),'FaceColor','interp','EdgeColor',edgecolor);
+                                patch( 'Faces', [A B E D],'Vertices', [x y z],'FaceVertexCData', data(:),'FaceColor','interp','EdgeColor',edgecolor);
+                                patch( 'Faces', [B E F C],'Vertices', [x y z],'FaceVertexCData', data(:),'FaceColor','interp','EdgeColor',edgecolor);
+                                patch( 'Faces', [C A D F],'Vertices', [x y z],'FaceVertexCData', data(:),'FaceColor','interp','EdgeColor',edgecolor);
+                        else
+                                A=elements(:,1); B=elements(:,2); C=elements(:,3); D=elements(:,4); 
+                                patch( 'Faces', [A B C D],  'Vertices', [x y z],'FaceVertexCData', data(:),'FaceColor','interp','EdgeColor',edgecolor);
+                        end
                 end
 

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

@@ -1 @@
-function [x y z elements is2d]=processmesh(md,data,options);
+function [x y z elements is2d isplanet]=processmesh(md,data,options);
 %PROCESSMESH - process mesh to be plotted
 %
@@ -91 +91 @@
         z=z*unit;
 end
+
+if md.planet,
+        isplanet=1;
+else
+        isplanet=0;
+end