Enable contact overpotential display (#5247)

Co-authored-by: jiaxiu han <[email protected]>
Co-authored-by: Valentin Sulzer <[email protected]>
Co-authored-by: Robert Timms <[email protected]>

Author: 4 people

src/pybamm/plotting/plot_voltage_components.py

@@ -110,6 +110,12 @@ def plot_voltage_components(
             "Ohmic negative electrode overpotential",
             "Ohmic positive electrode overpotential",
         ]
+    # Only add the contact overpotential label if its values are not (numerically) all zero
+    if not np.allclose(
+        solution["Contact overpotential [V]"].entries, 0, atol=1e-12, equal_nan=True
+    ):
+        overpotentials.append("Contact overpotential [V]")
+        labels.append("Contact overpotential")
 
     # Plot
     # Initialise
@@ -152,8 +158,9 @@ def plot_voltage_components(
     # Plot components
     for overpotential, label in zip(overpotentials, labels, strict=False):
         # negative overpotentials are positive for a discharge and negative for a charge
+        # Contact overpotential is positive for a discharge and negative for a charge
         # so we have to multiply by -1 to show them correctly
-        sgn = -1 if "egative" in overpotential else 1
+        sgn = -1 if ("egative" in overpotential or "Contact" in overpotential) else 1
         multiplier = sgn if "attery" in overpotential else sgn * num_cells
         bottom = top + multiplier * solution[overpotential].entries
         ax.fill_between(time, bottom, top, **kwargs_fill, label=label)