epoch check for arrays

pymagglobal/core.py

@@ -121,8 +121,10 @@ class Model(object):
         float
             The epoch in years.
         '''
-        if epoch < self.t_min or self.t_max < epoch:
-            warn(f'Epoch {epoch} is outside of the models range [{self.t_min},'
+        epoch = np.asarray(epoch)
+        out_of_range = epoch[(epoch < self.t_min) + (self.t_max < epoch)]
+        if np.any(out_of_range):
+            warn(f'Epochs {out_of_range} are outside of the models range [{self.t_min},'
                  f' {self.t_max}]. The result will be extrapolated and not '
                  f'robust.',
                  category=OutOfRangeWarning)