diff --git a/environment.yml b/environment.yml
index c7cbc6e..bdce036 100644
--- a/environment.yml
+++ b/environment.yml
@@ -8,7 +8,7 @@ dependencies:
   - discretize
   - ipywidgets>=0.6.0
   - jupyter
-  - matplotlib<3.5
+  - matplotlib==3.6.*
   - Pillow
   - pip
   - pymatsolver
diff --git a/geoscilabs/dcip/DCWidgetPlate2_5D.py b/geoscilabs/dcip/DCWidgetPlate2_5D.py
index ce06b3e..79d6d1d 100644
--- a/geoscilabs/dcip/DCWidgetPlate2_5D.py
+++ b/geoscilabs/dcip/DCWidgetPlate2_5D.py
@@ -527,7 +527,7 @@ def PLOT(
         label = "Electric Field (V/m)"
         xtype = "F"
         view = "vec"
-        streamOpts = {"color": "w"}
+        streamOpts = {"color": "w", "broken_streamlines": False, "density": 0.3}
         ind = indF
 
         # formatter = LogFormatter(10, labelOnlyBase=False)
diff --git a/geoscilabs/dcip/DCWidgetPlate_2D.py b/geoscilabs/dcip/DCWidgetPlate_2D.py
index 57e84e0..ae2088f 100644
--- a/geoscilabs/dcip/DCWidgetPlate_2D.py
+++ b/geoscilabs/dcip/DCWidgetPlate_2D.py
@@ -528,7 +528,7 @@ def plot_Surface_Potentials(
         label = "Electric Field (V/m)"
         xtype = "F"
         view = "vec"
-        streamOpts = {"color": "w"}
+        streamOpts = {"color": "w", "broken_streamlines": False, "density": 0.3}
         ind = indF
 
         # formatter = LogFormatter(10, labelOnlyBase=False)
@@ -553,7 +553,7 @@ def plot_Surface_Potentials(
         label = "Current density ($A/m^2$)"
         xtype = "F"
         view = "vec"
-        streamOpts = {"color": "w"}
+        streamOpts = {"color": "w", "broken_streamlines": False, "density": 0.3}
         ind = indF
 
         # formatter = LogFormatter(10, labelOnlyBase=False)
diff --git a/geoscilabs/dcip/DCWidgetResLayer2D.py b/geoscilabs/dcip/DCWidgetResLayer2D.py
index 37ce065..705f538 100644
--- a/geoscilabs/dcip/DCWidgetResLayer2D.py
+++ b/geoscilabs/dcip/DCWidgetResLayer2D.py
@@ -583,7 +583,7 @@ def plot_Surface_Potentials(
         label = "Current density ($A/m^2$)"
         xtype = "F"
         view = "vec"
-        streamOpts = {"color": "w"}
+        streamOpts = {"color": "w", "broken_streamlines": False, "density": 0.3}
         ind = indF
 
         # formatter = LogFormatter(10, labelOnlyBase=False)
diff --git a/geoscilabs/dcip/DCWidgetResLayer2_5D.py b/geoscilabs/dcip/DCWidgetResLayer2_5D.py
index b757adc..4aa22b9 100644
--- a/geoscilabs/dcip/DCWidgetResLayer2_5D.py
+++ b/geoscilabs/dcip/DCWidgetResLayer2_5D.py
@@ -563,7 +563,7 @@ def PLOT(
         label = "Electric Field (V/m)"
         xtype = "F"
         view = "vec"
-        streamOpts = {"color": "w"}
+        streamOpts = {"color": "w", "broken_streamlines": False, "density": 0.3}
         ind = indF
 
         # formatter = LogFormatter(10, labelOnlyBase=False)
@@ -588,7 +588,7 @@ def PLOT(
         label = "Current density ($A/m^2$)"
         xtype = "F"
         view = "vec"
-        streamOpts = {"color": "w"}
+        streamOpts = {"color": "w", "broken_streamlines": False, "density": 0.3}
         ind = indF
 
         # formatter = LogFormatter(10, labelOnlyBase=False)
diff --git a/geoscilabs/dcip/DCWidget_Overburden_2_5D.py b/geoscilabs/dcip/DCWidget_Overburden_2_5D.py
index b345d80..240ce16 100644
--- a/geoscilabs/dcip/DCWidget_Overburden_2_5D.py
+++ b/geoscilabs/dcip/DCWidget_Overburden_2_5D.py
@@ -565,7 +565,7 @@ def PLOT(
         label = "Electric Field (V/m)"
         xtype = "F"
         view = "vec"
-        streamOpts = {"color": "w"}
+        streamOpts = {"color": "w", "broken_streamlines": False, "density": 0.3}
         ind = indF
 
         # formatter = LogFormatter(10, labelOnlyBase =False)
@@ -590,7 +590,7 @@ def PLOT(
         label = "Current density ($A/m^2$)"
         xtype = "F"
         view = "vec"
-        streamOpts = {"color": "w"}
+        streamOpts = {"color": "w", "broken_streamlines": False, "density": 0.3}
         ind = indF
 
         # formatter = LogFormatter(10, labelOnlyBase =False)
diff --git a/geoscilabs/dcip/DC_cylinder.py b/geoscilabs/dcip/DC_cylinder.py
index b9ff4ee..867ce74 100644
--- a/geoscilabs/dcip/DC_cylinder.py
+++ b/geoscilabs/dcip/DC_cylinder.py
@@ -408,7 +408,7 @@ def plot_Surface_Potentials(
         label = "Electric Field (V/m)"
         xtype = "F"
         view = "vec"
-        streamOpts = {"color": "w"}
+        streamOpts = {"color": "w", "broken_streamlines": False, "density": 0.3}
         ind = indF
 
         # formatter = LogFormatter(10, labelOnlyBase=False)
@@ -433,7 +433,7 @@ def plot_Surface_Potentials(
         label = "Current density ($A/m^2$)"
         xtype = "F"
         view = "vec"
-        streamOpts = {"color": "w"}
+        streamOpts = {"color": "w", "broken_streamlines": False, "density": 0.3}
         ind = indF
 
         # formatter = LogFormatter(10, labelOnlyBase=False)