diff --git a/boundaries/boundaryhandling.py b/boundaries/boundaryhandling.py
index 2f4b4c82b3a3fe12615d48ea4f92ee375e7dc4ab..95dd77d825ee21f699ae608ad3cc46fd6f29a890 100644
--- a/boundaries/boundaryhandling.py
+++ b/boundaries/boundaryhandling.py
@@ -76,11 +76,6 @@ class BoundaryHandling:
     def flag_array_name(self):
         return self.flag_interface.flag_field_name
 
-    def get_boundary_name_to_flag_dict(self):
-        result = {b_obj.name: b_info.flag for b_obj, b_info in self._boundary_object_to_boundary_info.items()}
-        result['domain'] = self.flag_interface.domain_flag
-        return result
-
     def get_mask(self, slice_obj, boundary_obj, inverse=False):
         if isinstance(boundary_obj, str) and boundary_obj.lower() == 'domain':
             flag = self.flag_interface.domain_flag
diff --git a/plot2d.py b/plot2d.py
index 8f7faa2ed9d63983e8c6cdda073a78a1b79241da..bb0b18feb69f0a166f67964946742c57e3ccf775 100644
--- a/plot2d.py
+++ b/plot2d.py
@@ -196,8 +196,9 @@ def vector_field_animation(run_function, step=2, rescale=True, plot_setup_functi
     if rescale:
         max_norm = np.max(norm(field, axis=2, ord=2))
         field = field / max_norm
-        if 'scale' not in kwargs:
-            kwargs['scale'] = 1.0
+        kwargs.setdefault('scale', 1 / step)
+        kwargs.setdefault('angles', 'xy')
+        kwargs.setdefault('scale_units', 'xy')
 
     quiver_plot = vector_field(field, step=step, **kwargs)
     plot_setup_function(quiver_plot)