diff --git a/4_grains/gtGrainAllblobsFilterOrder.m b/4_grains/gtGrainAllblobsFilterOrder.m
index 15828f20c7950440eefa36cad7880af804bcc638..b820c728280670f83d5ca7692828363684d9a307 100644
--- a/4_grains/gtGrainAllblobsFilterOrder.m
+++ b/4_grains/gtGrainAllblobsFilterOrder.m
@@ -3,11 +3,13 @@ function gr = gtGrainAllblobsFilterOrder(gr, refgr_omind)
 
     gr_omind = gr.allblobs.omind;
     gr_omind_ind = [ ...
-        find(gr_omind == 1); find(gr_omind == 2); ...
-        find(gr_omind == 3); find(gr_omind == 4) ];
+        find(gr_omind == 1)'; find(gr_omind == 2)'; ...
+        find(gr_omind == 3)'; find(gr_omind == 4)' ];
+    gr_omind_ind = reshape(gr_omind_ind, [], 1);
     refgr_omind_ind = [ ...
-        find(refgr_omind == 1); find(refgr_omind == 2); ...
-        find(refgr_omind == 3); find(refgr_omind == 4) ];
+        find(refgr_omind == 1)'; find(refgr_omind == 2)'; ...
+        find(refgr_omind == 3)'; find(refgr_omind == 4)' ];
+    refgr_omind_ind = reshape(refgr_omind_ind, [], 1);
 
     fnames = fieldnames(gr.allblobs);
     for ii_f = 1:numel(fnames)