diff --git a/auxiliary_tools/cdat_regression_testing/892-bottleneck/run_script.cfg b/auxiliary_tools/cdat_regression_testing/892-bottleneck/run_script.cfg
new file mode 100644
index 000000000..3d06adcde
--- /dev/null
+++ b/auxiliary_tools/cdat_regression_testing/892-bottleneck/run_script.cfg
@@ -0,0 +1,13 @@
+[#]
+sets = ["lat_lon"]
+case_id = "ERA5"
+variables = ["U"]
+ref_name = "ERA5"
+reference_name = "ERA5 Reanalysis"
+seasons = ["ANN", "01", "02", "03", "04", "05", "06", "07", "08", "09", "10", "11", "12", "DJF", "MAM", "JJA", "SON"]
+plevs = [850.0]
+test_colormap = "PiYG_r"
+reference_colormap = "PiYG_r"
+contour_levels = [-20, -15, -10, -8, -5, -3, -1, 1, 3, 5, 8, 10, 15, 20]
+diff_levels = [-8, -6, -5, -4, -3, -2, -1, 1, 2, 3, 4, 5, 6, 8]
+regrid_method = "bilinear"
\ No newline at end of file
diff --git a/auxiliary_tools/cdat_regression_testing/892-bottleneck/run_script.py b/auxiliary_tools/cdat_regression_testing/892-bottleneck/run_script.py
new file mode 100644
index 000000000..a26cefb2e
--- /dev/null
+++ b/auxiliary_tools/cdat_regression_testing/892-bottleneck/run_script.py
@@ -0,0 +1,35 @@
+import os
+from e3sm_diags.parameter.core_parameter import CoreParameter
+from e3sm_diags.run import runner
+
+param = CoreParameter()
+
+
+param.reference_data_path = (
+    "/global/cfs/cdirs/e3sm/diagnostics/observations/Atm/time-series"
+)
+param.test_data_path = "/global/cfs/cdirs/e3sm/chengzhu/eamxx/post/data/rgr"
+param.test_name = "eamxx_decadal"
+param.seasons = ["ANN"]
+# param.save_netcdf = True
+
+param.ref_timeseries_input = True
+# Years to slice the ref data, base this off the years in the filenames.
+param.ref_start_yr = "1996"
+param.ref_end_yr = "1996"
+
+prefix = "/global/cfs/cdirs/e3sm/www/cdat-migration-fy24/877-attr-err"
+param.results_dir = os.path.join(prefix, "eamxx_decadal_1996_1107_edv3")
+
+runner.sets_to_run = [
+    "lat_lon",
+    "zonal_mean_xy",
+    "zonal_mean_2d",
+    "zonal_mean_2d_stratosphere",
+    "polar",
+    "cosp_histogram",
+    "meridional_mean_2d",
+    "annual_cycle_zonal_mean",
+]
+
+runner.run_diags([param])