diff --git a/CHANGELOG.rst b/CHANGELOG.rst
index 1cbeb498..4e3bcbf6 100644
--- a/CHANGELOG.rst
+++ b/CHANGELOG.rst
@@ -5,6 +5,7 @@ Change Log
 -[ADDED] heat_consumer plotting
 -[CHANGED] switched from setup.py to pyproject.toml
 - [CHANGED] variable "alpha_w_per_m2k" to "u_w_per_m2k"
+- [FIXED] Pressure plot not working for circ pump
 - [ADDED] variable "u_w_per_m2k" to std_type pipe
 - [ADDED] standard district heating pipe types
 
diff --git a/src/pandapipes/plotting/pipeflow_results.py b/src/pandapipes/plotting/pipeflow_results.py
index bc32fc81..a7a64f3b 100644
--- a/src/pandapipes/plotting/pipeflow_results.py
+++ b/src/pandapipes/plotting/pipeflow_results.py
@@ -54,7 +54,7 @@ def plot_pressure_profile(net, ax=None, x0_junctions=None, plot_pressure_control
         axis to plot to, by default None
     x0_junctions : Any[list[int], pd.Index[int]], optional
         list of junction indices which should be at position x0. If None, all in service slack junctions are considered,
-        by default None
+        by default None. For circ pumps, if no x0_junctions are given the flow junctions are chosen automatically.
     plot_pressure_controller : bool, optional
         Whether vertical lines should be plotted to display the pressure drop of the pressure controller,
         by default True
@@ -92,9 +92,9 @@ def plot_pressure_profile(net, ax=None, x0_junctions=None, plot_pressure_control
     if x0_junctions is None:
         x0_junctions = set(net.ext_grid[net.ext_grid.in_service].junction.values)
         if hasattr(net, "circ_pump_pressure"):
-            x0_junctions |= set(net.circ_pump_pressure[net.circ_pump_pressure.in_service].junction.values)
+            x0_junctions |= set(net.circ_pump_pressure[net.circ_pump_pressure.in_service].flow_junction.values)
         if hasattr(net, "circ_pump_mass"):
-            x0_junctions |= set(net.circ_pump_mass[net.circ_pump_mass.in_service].junction.values)
+            x0_junctions |= set(net.circ_pump_mass[net.circ_pump_mass.in_service].flow_junction.values)
         x0_junctions = list(x0_junctions)
 
     unsupplied_junctions = list(top.unsupplied_junctions(net, slacks=set(x0_junctions)))