Adding a rule for Magnetic intensity map with magnetic and velocity streamlines

plotter.py

@@ -1057,6 +1057,18 @@ class Plotter(Aggregator, BaseProcessor):
                 "Density slice with magnetic field and velocity  overlay",
                 dependencies=["rho", "B_h", "B_v", "speed_h", "speed_v"],
             ),
+            "B_int_B_vel": PlotRule(
+                self,
+                partial(
+                    self._plot_map,
+                    "B_int",
+                    label="B",
+                    unit=cst.T,
+                    overlays=[self._overlay_B, self._overlay_speed],
+                ),
+                "Density slice with magnetic field and velocity  overlay",
+                dependencies=["B_int", "B_h", "B_v", "speed_h", "speed_v"],
+            ),
             "jeans_ratio": PlotRule(
                 self,
                 partial(