Index: /issm/trunk/src/m/model/plot/plot_basaldrag.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_basaldrag.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_basaldrag.m	(revision 8472)
@@ -39,10 +39,10 @@
 
 	%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
Index: /issm/trunk/src/m/model/plot/plot_boundaries.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_boundaries.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_boundaries.m	(revision 8472)
@@ -13,5 +13,5 @@
 	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),
Index: /issm/trunk/src/m/model/plot/plot_contour.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_contour.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_contour.m	(revision 8472)
@@ -8,5 +8,5 @@
 
 %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);
 
Index: /issm/trunk/src/m/model/plot/plot_drivingstress.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_drivingstress.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_drivingstress.m	(revision 8472)
@@ -11,5 +11,5 @@
 
 %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,5 +17,5 @@
 %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
Index: /issm/trunk/src/m/model/plot/plot_edges.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_edges.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_edges.m	(revision 8472)
@@ -11,5 +11,5 @@
 
 %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)
Index: /issm/trunk/src/m/model/plot/plot_elementnumbering.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_elementnumbering.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_elementnumbering.m	(revision 8472)
@@ -10,5 +10,5 @@
 
 %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);
 
Index: /issm/trunk/src/m/model/plot/plot_elementstype.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_elementstype.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_elementstype.m	(revision 8472)
@@ -8,5 +8,5 @@
 
 %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);
 
Index: /issm/trunk/src/m/model/plot/plot_gridded.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_gridded.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_gridded.m	(revision 8472)
@@ -8,5 +8,5 @@
 
 %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);
 
Index: /issm/trunk/src/m/model/plot/plot_highlightelements.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_highlightelements.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_highlightelements.m	(revision 8472)
@@ -11,5 +11,5 @@
 
 %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);
 
Index: /issm/trunk/src/m/model/plot/plot_highlightnodes.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_highlightnodes.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_highlightnodes.m	(revision 8472)
@@ -8,5 +8,5 @@
 
 %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);
 
Index: /issm/trunk/src/m/model/plot/plot_importancefactors.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_importancefactors.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_importancefactors.m	(revision 8472)
@@ -64,5 +64,5 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %edgecolor
Index: /issm/trunk/src/m/model/plot/plot_manager.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_manager.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_manager.m	(revision 8472)
@@ -166,5 +166,5 @@
 
 %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,5 +178,5 @@
 
 %plot unit
-plot_unit(x,y,z,elements,data2,is2d,datatype,options);
+plot_unit(x,y,z,elements,data2,is2d,isplanet,datatype,options);
 	
 
Index: /issm/trunk/src/m/model/plot/plot_mesh.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_mesh.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_mesh.m	(revision 8472)
@@ -8,5 +8,5 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %plot mesh
@@ -22,10 +22,16 @@
 	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
 
Index: /issm/trunk/src/m/model/plot/plot_nodenumbering.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_nodenumbering.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_nodenumbering.m	(revision 8472)
@@ -8,5 +8,5 @@
 
 %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);
 
Index: /issm/trunk/src/m/model/plot/plot_overlay.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_overlay.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_overlay.m	(revision 8472)
@@ -8,5 +8,5 @@
 
 %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;
Index: /issm/trunk/src/m/model/plot/plot_parthist.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_parthist.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_parthist.m	(revision 8472)
@@ -8,5 +8,5 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %plot mesh
Index: /issm/trunk/src/m/model/plot/plot_parthistn.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_parthistn.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_parthistn.m	(revision 8472)
@@ -8,5 +8,5 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %plot mesh
Index: /issm/trunk/src/m/model/plot/plot_parthistw.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_parthistw.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_parthistw.m	(revision 8472)
@@ -8,5 +8,5 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %plot mesh
Index: /issm/trunk/src/m/model/plot/plot_penalties.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_penalties.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_penalties.m	(revision 8472)
@@ -8,5 +8,5 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %plot mesh penalties
Index: /issm/trunk/src/m/model/plot/plot_pressureload.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_pressureload.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_pressureload.m	(revision 8472)
@@ -11,5 +11,5 @@
 
 %process mesh and data
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 pressureload=md.pressureload;
 
Index: /issm/trunk/src/m/model/plot/plot_qmu_mass_flux_segments.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_qmu_mass_flux_segments.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_qmu_mass_flux_segments.m	(revision 8472)
@@ -9,5 +9,5 @@
 
 %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;
Index: /issm/trunk/src/m/model/plot/plot_qmumean.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_qmumean.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_qmumean.m	(revision 8472)
@@ -14,5 +14,5 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %find response function
Index: /issm/trunk/src/m/model/plot/plot_qmustddev.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_qmustddev.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_qmustddev.m	(revision 8472)
@@ -14,5 +14,5 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %find response function
Index: /issm/trunk/src/m/model/plot/plot_riftfraction.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_riftfraction.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_riftfraction.m	(revision 8472)
@@ -13,5 +13,5 @@
 
 %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); 
Index: /issm/trunk/src/m/model/plot/plot_riftnumbering.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_riftnumbering.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_riftnumbering.m	(revision 8472)
@@ -8,5 +8,5 @@
 
 %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);
 
Index: /issm/trunk/src/m/model/plot/plot_riftpenetration.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_riftpenetration.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_riftpenetration.m	(revision 8472)
@@ -8,5 +8,5 @@
 
 %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); 
Index: /issm/trunk/src/m/model/plot/plot_riftrelvel.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_riftrelvel.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_riftrelvel.m	(revision 8472)
@@ -30,5 +30,5 @@
 
 %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);
Index: /issm/trunk/src/m/model/plot/plot_rifts.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_rifts.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_rifts.m	(revision 8472)
@@ -8,5 +8,5 @@
 
 %process data and model
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 
 %plot mesh
Index: /issm/trunk/src/m/model/plot/plot_riftvel.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_riftvel.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_riftvel.m	(revision 8472)
@@ -26,5 +26,5 @@
 
 %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);
Index: /issm/trunk/src/m/model/plot/plot_section.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_section.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_section.m	(revision 8472)
@@ -21,5 +21,5 @@
 
 %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,5 +80,5 @@
 			%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])
Index: /issm/trunk/src/m/model/plot/plot_segments.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_segments.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_segments.m	(revision 8472)
@@ -11,5 +11,5 @@
 
 %process mesh and data
-[x y z elements is2d]=processmesh(md,[],options);
+[x y z elements is2d isplanet]=processmesh(md,[],options);
 segments=md.segments;
 
Index: /issm/trunk/src/m/model/plot/plot_streamlines.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_streamlines.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_streamlines.m	(revision 8472)
@@ -6,5 +6,5 @@
 
 %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);
Index: /issm/trunk/src/m/model/plot/plot_tensor_components.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_tensor_components.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_tensor_components.m	(revision 8472)
@@ -24,5 +24,5 @@
 
 %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,30 +36,30 @@
 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
Index: /issm/trunk/src/m/model/plot/plot_tensor_principal.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_tensor_principal.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_tensor_principal.m	(revision 8472)
@@ -37,5 +37,5 @@
 
 %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,18 +46,18 @@
 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
Index: /issm/trunk/src/m/model/plot/plot_tensor_principalaxis.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_tensor_principalaxis.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_tensor_principalaxis.m	(revision 8472)
@@ -11,5 +11,5 @@
 
 %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)),
Index: /issm/trunk/src/m/model/plot/plot_transient_field.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_transient_field.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_transient_field.m	(revision 8472)
@@ -27,5 +27,5 @@
 
 %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,5 +38,5 @@
 	%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);
 
Index: /issm/trunk/src/m/model/plot/plot_transient_movie.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_transient_movie.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_transient_movie.m	(revision 8472)
@@ -21,9 +21,9 @@
 
 		%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);
 		
Index: /issm/trunk/src/m/model/plot/plot_unit.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_unit.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/plot_unit.m	(revision 8472)
@@ -1,7 +1,7 @@
-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,10 +34,15 @@
 			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
 
Index: /issm/trunk/src/m/model/plot/processmesh.m
===================================================================
--- /issm/trunk/src/m/model/plot/processmesh.m	(revision 8471)
+++ /issm/trunk/src/m/model/plot/processmesh.m	(revision 8472)
@@ -1,3 +1,3 @@
-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,2 +91,8 @@
 	z=z*unit;
 end
+
+if md.planet,
+	isplanet=1;
+else
+	isplanet=0;
+end