diff --git a/zUtil_Fit/gtFitFwdSimGrains.m b/zUtil_Fit/gtFitFwdSimGrains.m
index e18e4378a1765e691a4f9abe2ebf817aef791e1d..33c2e566d7a991cfc782602a5ccf69d5cad27f9a 100644
--- a/zUtil_Fit/gtFitFwdSimGrains.m
+++ b/zUtil_Fit/gtFitFwdSimGrains.m
@@ -156,8 +156,8 @@ for ii = 1:length(grain)
         grain{ii}.fsim.plref     = [plref, plref, plref, plref];
         grain{ii}.fsim.refind    = [refind, refind, refind, refind];
         grain{ii}.fsim.hklspind  = grain{ii}.fsim.refind;
-        grain{ii}.fsim.thetatype = [cryst.thetatypesp', cryst.thetatypesp', ...
-            cryst.thetatypesp', cryst.thetatypesp'];
+		thetatype                = cryst.thetatypesp(:)';
+        grain{ii}.fsim.thetatype = [thetatype, thetatype, thetatype, thetatype];
         grain{ii}.fsim.hklind    = grain{ii}.fsim.thetatype;
            
     else
@@ -276,8 +276,8 @@ for ii = 1:length(grain)
         grain{ii}.fsim.plref     = [plref, plref, plref, plref];
         grain{ii}.fsim.refind    = [refind, refind, refind, refind];
         grain{ii}.fsim.hklspind  = grain{ii}.fsim.refind;
-        grain{ii}.fsim.thetatype = [cryst.thetatypesp', cryst.thetatypesp', ...
-            cryst.thetatypesp', cryst.thetatypesp'];
+		thetatype                = cryst.thetatypesp(:)';
+        grain{ii}.fsim.thetatype = [thetatype, thetatype, thetatype, thetatype];
         grain{ii}.fsim.hklind    = grain{ii}.fsim.thetatype;
              
     end