diff --git a/zUtil_Boundaries/gtReadBoundaryProperties.m b/zUtil_Boundaries/gtReadBoundaryProperties.m
index 2925bfafdfb333c89b5a208cf6489b60ff537434..12e6afcb3297b49951e14237fe425cc172427299 100644
--- a/zUtil_Boundaries/gtReadBoundaryProperties.m
+++ b/zUtil_Boundaries/gtReadBoundaryProperties.m
@@ -1,9 +1,9 @@
 function gtReadBoundaryProperties(phaseid, boundaries_structure, vol_out, r_vectors)
 % gtReadBoundaryProperties(phaseid, boundaries_structure, vol_out, r_vectors)
 
-disp('have changed this to axis convention matlab1 // instX , matlab2 // instY, matlab3 // instZ')
-disp('AK, 11/8/2010')
-pause(3)
+disp('have changed this to axis convention matlab1 // instX , matlab2 // instY, matlab3 // instZ');
+disp('AK, 11/8/2010');
+pause(3);
 
 %continue to add to boundaries_structure.mat - orientations,
 %disorientation, r-vectors etc.
@@ -23,9 +23,9 @@ symm = gtCrystGetSymmetryOperators([], spacegroup);
 symm = {symm.g3};
 
 % loop through the boundaries
-for ii=1:size(boundaries_structure, 2)
+for ii = 1:size(boundaries_structure, 2)
 
-    if boundaries_structure(ii).count<4 %can still be useful to calc index planes where we have only one of the grains
+    if (boundaries_structure(ii).count < 4) %can still be useful to calc index planes where we have only one of the grains
         continue;
     else
         %Read in the r-vector data for the two grains
@@ -33,23 +33,23 @@ for ii=1:size(boundaries_structure, 2)
         grain2=boundaries_structure(ii).grain2;
 
         if grain1==0
-            g1=[];
-            g1equiv=[];
+            g1 = [];
+            g1equiv = [];
         elseif grain1==-1
             continue;
         else
-            R1= r_vectors(grain1,2:4);
-            boundaries_structure(ii).grain1_R_vector=R1;
+            R1 = r_vectors(grain1, 2:4);
+            boundaries_structure(ii).grain1_R_vector = R1;
             g1 = gtMathsRod2OriMat(R1);
         end
 
         if grain2==0
-            g2=[];
+            g2 = [];
         elseif grain2==-1
-            continue
+            continue;
         else
-            R2= r_vectors(grain2,2:4);
-            boundaries_structure(ii).grain2_R_vector=R2;
+            R2 = r_vectors(grain2,2:4);
+            boundaries_structure(ii).grain2_R_vector = R2;
             g2 = gtMathsRod2OriMat(R2);
         end
 
@@ -58,7 +58,7 @@ for ii=1:size(boundaries_structure, 2)
         if ~isempty(g1) && ~isempty(g2)
             % If we have both grains
             % Need to search symmetry equivalents for the minimum misorientation
-            rot_offset=[];
+            rot_offset = [];
             % warning('g*symm maybe should be symm*g')
             for jj = 1:length(symm)
                 g1equiv = g1*symm{jj};
@@ -171,11 +171,11 @@ for ii=1:size(boundaries_structure, 2)
                 boundaries_structure(ii).misorientation_axis_hex=hkil;
             end
         else
-            disp('crystallography not supported!')
+            disp('crystallography not supported!');
         end
     end % skip small boundaries
 end
 
-save('boundaries_structure.mat', 'boundaries_structure')
+save('boundaries_structure.mat', 'boundaries_structure');
 
 end % end of function