Formatting and commenting...

Signed-off-by: Laura Nervo <laura.nervo@esrf.fr>

git-svn-id: https://svn.code.sf.net/p/dct/code/trunk@987 4c865b51-4357-4376-afb4-474e03ccb993

6_rendering/gtShowSampleSurface.m

@@ -58,7 +58,7 @@ else %more difficult...
     end
     
     
-    for x=[1:size(section, 2)] %along the rows
+    for x=1:size(section, 2) %along the rows
 
         %estimate how thick a layer we should look at to get near surface
         %grains
@@ -99,7 +99,7 @@ set(get(get(section_id, 'CurrentAxes'), 'Title'), 'String',sprintf('Section at a
 set(get(get(section_id, 'CurrentAxes'), 'yLabel'),'String', 'z axis')
 set(get(get(section_id, 'CurrentAxes'), 'xLabel'),'String', 'y axis')
 axis equal
-
+axis square
 
 
 if 0
@@ -152,10 +152,10 @@ for i=1:length(grains)
     %basal trace
     x_tilt=atand(-basal_trace(3)/basal_trace(2));
     %display for user
-    disp(sprintf('c_axis of grain %d - in 3D coordinates, after sample rotation, is:', grain))
-    disp(sprintf('%0.2f %0.2f %0.2f', c_axis(1), c_axis(2), c_axis(3)))
-    disp(sprintf('x_tilt for grain %d should be %0.2f degrees', grain, x_tilt))
-    disp(sprintf('y_tilt for grain %d should be %0.2f degrees', grain, y_tilt))
+    fprintf('c_axis of grain %d - in 3D coordinates, after sample rotation, is:\n', grain)
+    fprintf('%0.2f %0.2f %0.2f\n', c_axis(1), c_axis(2), c_axis(3))
+    fprintf('x_tilt for grain %d should be %0.2f degrees\n', grain, x_tilt)
+    fprintf('y_tilt for grain %d should be %0.2f degrees\n', grain, y_tilt)
     disp('~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~')
     
     else
@@ -170,3 +170,4 @@ grains(find(dum))=[];
 %change gtPlotHexagon to take the two coordinate systems into account
 gtPlotHexagon(grains, r_vectors, angle, cmap)
 
+end % end of function