gtLoadGrain: can load multiple grains

Signed-off-by: Nicola Vigano <nicola.vigano@esrf.fr>

4_grains/gtLoadGrain.m

@@ -3,11 +3,13 @@ function grain = gtLoadGrain(phase_id, grain_id, varargin)
     conf = parse_pv_pairs(conf, varargin);
 
     phase_dir = fullfile('4_grains', sprintf('phase_%02d', phase_id));
-    grain_file = fullfile(phase_dir, sprintf('grain_%04d.mat', grain_id));
+    for ii = numel(grain_id):-1:1
+        grain_file = fullfile(phase_dir, sprintf('grain_%04d.mat', grain_id(ii)));
 
-    if (isempty(conf.fields))
-        grain = load(grain_file);
-    else
-        grain = load(grain_file, conf.fields{:});
+        if (isempty(conf.fields))
+            grain(ii) = load(grain_file);
+        else
+            grain(ii) = load(grain_file, conf.fields{:});
+        end
     end
 end