diff --git a/engine/performance_stats.py b/engine/performance_stats.py
new file mode 100644
index 00000000..536bf616
--- /dev/null
+++ b/engine/performance_stats.py
@@ -0,0 +1,85 @@
+"""
+CLI for performance checks
+
+This module provides a command-line interface (CLI) tool for evaluating the
+performance of a current experiment by comparing its timing data with a
+reference.
+It assesses whether the current experiment's runtime is within acceptable limits
+based on specified parameters and provides feedback on whether the performance
+check has passed or failed.
+"""
+
+import sys
+
+import click
+import numpy as np
+import xarray as xr
+
+from util.click_util import cli_help
+from util.log_handler import logger
+from util.tree import TimingTree
+
+
+@click.command()
+@click.option(
+    "--timing-current",
+    default="database",
+    help=cli_help["timing_current"],
+)
+@click.option(
+    "--timing-stats",
+    default='timer_stats.nc',
+    help='name of the timing stats file',
+)
+@click.option("--i-table", type=int, help=cli_help["i_table"], default=-1)
+def performance_stats(
+    timing_current,
+    timing_stats,
+    i_table,
+):  # pylint: disable=too-many-positional-arguments
+
+    timing_tree = TimingTree.from_json(timing_current)
+
+    timer_names = ['model_init', 'total', 'integrate_nh', 'nh_solve', 'nh_hdiff', 'transport', 'physics']
+
+    aggregate_timer_names = {'dycore':('nh_solve', 'nh_hdiff', 'transport')}
+
+    all_timer_names = timer_names + list(aggregate_timer_names.keys())
+
+    dims = ('name', 'metric')
+
+    coords = {'name': all_timer_names, 'metric': ['mean', 'std']}
+
+    print(all_timer_names)
+
+    # Create an empty DataArray with np.nan
+    timer_stats = xr.DataArray(
+        data=np.full((len(coords['name']), len(coords['metric'])), np.nan),  # Fill with NaN
+        coords=coords,
+        dims=dims
+    )
+
+
+    for timer_name in timer_names:
+      times = np.asarray(
+          timing_tree.data[i_table].loc[(timer_name, slice(None)), "total max (s)"].values
+      )
+      timer_stats.loc[timer_name, 'mean'] = np.mean(times)
+      timer_stats.loc[timer_name, 'std'] = np.std(times)
+
+
+    for aggregated_timer_name in aggregate_timer_names.keys():
+
+      times = np.zeros_like(np.asarray(
+          timing_tree.data[i_table].loc[(aggregate_timer_names[aggregated_timer_name][0], slice(None)), "total max (s)"].values
+      ))
+      for timer_name in aggregate_timer_names[aggregated_timer_name]:
+        times += np.asarray(
+            timing_tree.data[i_table].loc[(timer_name, slice(None)), "total max (s)"].values
+        )
+      timer_stats.loc[aggregated_timer_name, 'mean'] = np.mean(times)
+      timer_stats.loc[aggregated_timer_name, 'std'] = np.std(times)
+
+
+    # Save the dataset to a NetCDF file
+    timer_stats.to_netcdf(timing_stats)
\ No newline at end of file
diff --git a/probtest.py b/probtest.py
index 2dc0f7e2..f5312a00 100755
--- a/probtest.py
+++ b/probtest.py
@@ -14,6 +14,7 @@
 from engine.init import init
 from engine.performance import performance
 from engine.performance_check import performance_check
+from engine.performance_stats import performance_stats
 from engine.perturb import perturb
 from engine.run_ensemble import run_ensemble
 from engine.select_members import select_members
@@ -57,6 +58,7 @@ def cli(ctx, log_level, log_file):
 cli.add_command(performance_plot)
 cli.add_command(performance_meta_data)
 cli.add_command(performance_check)
+cli.add_command(performance_stats)
 cli.add_command(cdo_table)
 cli.add_command(cdo_table_reader)