diff --git a/jobs/JGLOBAL_ANALYSIS_STATS b/jobs/JGLOBAL_ANALYSIS_STATS
index e69de29bb2..6543735fb1 100755
--- a/jobs/JGLOBAL_ANALYSIS_STATS
+++ b/jobs/JGLOBAL_ANALYSIS_STATS
@@ -0,0 +1,35 @@
+#! /usr/bin/env bash
+
+source "${HOMEgfs}/ush/preamble.sh"
+source "${HOMEgfs}/ush/jjob_header.sh" -e "anlstat" -c "base anlstat"
+
+##############################################
+# Set variables used in the script
+##############################################
+
+
+##############################################
+# Begin JOB SPECIFIC work
+##############################################
+
+
+###############################################################
+# Run relevant script
+
+EXSCRIPT=${ANLSTATSPY:-${SCRgfs}/exglobal_analysis_stats.py}
+${EXSCRIPT}
+status=$?
+[[ ${status} -ne 0 ]] && exit "${status}"
+
+##############################################
+# End JOB SPECIFIC work
+##############################################
+
+##############################################
+# Final processing
+##############################################
+if [[ -e "${pgmout}" ]] ; then
+  cat "${pgmout}"
+fi
+
+exit 0
\ No newline at end of file
diff --git a/jobs/rocoto/anlstat.sh b/jobs/rocoto/anlstat.sh
index e69de29bb2..52292a582c 100755
--- a/jobs/rocoto/anlstat.sh
+++ b/jobs/rocoto/anlstat.sh
@@ -0,0 +1,18 @@
+#! /usr/bin/env bash
+
+source "${HOMEgfs}/ush/preamble.sh"
+
+###############################################################
+# Source UFSDA workflow modules
+. "${HOMEgfs}/ush/load_ufsda_modules.sh"
+status=$?
+[[ ${status} -ne 0 ]] && exit "${status}"
+
+export job="anlstat"
+export jobid="${job}.$$"
+
+###############################################################
+# Execute the JJOB
+"${HOMEgfs}/jobs/JGLOBAL_ANALYSIS_STATS"
+status=$?
+exit "${status}"
\ No newline at end of file
diff --git a/scripts/exglobal_analysis_stats.py b/scripts/exglobal_analysis_stats.py
index e69de29bb2..045fc914b3 100755
--- a/scripts/exglobal_analysis_stats.py
+++ b/scripts/exglobal_analysis_stats.py
@@ -0,0 +1,20 @@
+#!/usr/bin/env python3
+# exglobal_analysis_stats.py
+# This script will run the OOPS/JEDI code to
+
+import os
+
+from wxflow import Logger, cast_strdict_as_dtypedict
+from pygfs.task.atm_analysis import AtmAnalysis
+
+# Initialize root logger
+logger = Logger(level='DEBUG', colored_log=True)
+
+
+if __name__ == '__main__':
+
+    # Take configuration from environment and cast it as python dictionary
+    config = cast_strdict_as_dtypedict(os.environ)
+
+    # placeholder
+    print('Do nothing right now!')
\ No newline at end of file