Index: /issm/trunk/src/m/model/plot/applyoptions.m
===================================================================
--- /issm/trunk/src/m/model/plot/applyoptions.m	(revision 7650)
+++ /issm/trunk/src/m/model/plot/applyoptions.m	(revision 7651)
@@ -69,4 +69,9 @@
 	end
 end
+
+if exist(options,'z'),
+	axis square;
+end
+
 
 %xlim
Index: /issm/trunk/src/m/model/plot/plot_manager.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_manager.m	(revision 7650)
+++ /issm/trunk/src/m/model/plot/plot_manager.m	(revision 7651)
@@ -173,4 +173,5 @@
 %plot unit
 h=plot_unit(x,y,z,elements,data2,is2d,datatype,options);
+	
 
 %apply all options
Index: /issm/trunk/src/m/model/plot/plot_unit.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_unit.m	(revision 7650)
+++ /issm/trunk/src/m/model/plot/plot_unit.m	(revision 7651)
@@ -21,6 +21,5 @@
 		pos=find(~isnan(data)); %needed for element on water
 		if is2d,
-			A=elements(pos,1); B=elements(pos,2); C=elements(pos,3); 
-			h=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', data(pos),'FaceColor','flat','EdgeColor',edgecolor);
+			A=elements(pos,1); B=elements(pos,2); C=elements(pos,3); h=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', data(pos),'FaceColor','flat','EdgeColor',edgecolor);
 		else
 			A=elements(pos,1); B=elements(pos,2); C=elements(pos,3); D=elements(pos,4); E=elements(pos,5); F=elements(pos,6);
Index: /issm/trunk/src/m/model/plot/processmesh.m
===================================================================
--- /issm/trunk/src/m/model/plot/processmesh.m	(revision 7650)
+++ /issm/trunk/src/m/model/plot/processmesh.m	(revision 7651)
@@ -26,9 +26,19 @@
 		y2d=md.y2d;
 	end
+
+	z_field=getfieldvalue(options,'z',md.z);
+	if ischar(z_field),
+		z=md.(z_field);
+	elseif isnumeric(z_field),
+		z=z_field;
+	else
 		z=md.z;
-		z2d=md.z2d;
-		elements2d=md.elements2d;
-		elements=md.elements;
-		elements_type2d=md.elements_type2d;
+	end
+
+
+	z2d=md.z2d;
+	elements2d=md.elements2d;
+	elements=md.elements;
+	elements_type2d=md.elements_type2d;
 
 	%is it a 2d plot?