diff --git a/cr2/plot_utils.py b/cr2/plot_utils.py
index 113471f..a1a94de 100644
--- a/cr2/plot_utils.py
+++ b/cr2/plot_utils.py
@@ -175,3 +175,14 @@
 
     for ax, run in zip(axis, runs):
         run.thermal_governor.plot_input_power(actor_order, title=run.name, ax=ax)
+
+def plot_output_power(runs, actor_order, width=None, height=None):
+    """Make a multicolumn plot of the output power of each run"""
+    num_runs = len(runs)
+    axis = pre_plot_setup(width=width, height=height, ncols=num_runs)
+
+    if num_runs == 1:
+        axis = [axis]
+
+    for ax, run in zip(axis, runs):
+        run.thermal_governor.plot_output_power(actor_order, title=run.name, ax=ax)