From 3c0c57f8fc0cd8288a088655e1fccdffc12c3426 Mon Sep 17 00:00:00 2001
From: "edward.safford" <edward.safford@noaa.gov>
Date: Fri, 14 Feb 2025 18:41:10 +0000
Subject: [PATCH 1/2] Ref #26

Add radHist template and driver script.
---
 parm/gfs/test_hist.yaml     | 26 +++++++++++++
 parm/templates/radHist.yaml | 73 +++++++++++++++++++++++++++++++++++++
 2 files changed, 99 insertions(+)
 create mode 100644 parm/gfs/test_hist.yaml
 create mode 100644 parm/templates/radHist.yaml

diff --git a/parm/gfs/test_hist.yaml b/parm/gfs/test_hist.yaml
new file mode 100644
index 0000000..be4ccc1
--- /dev/null
+++ b/parm/gfs/test_hist.yaml
@@ -0,0 +1,26 @@
+# DA Monitoring source dictionary
+
+# This is a simple script to create a radiance plot using the radHist 
+# (histogram) plot template.
+#
+# It is intended as a proof of concept and isn't an actual part of the 
+# legacy RadMon.
+#
+
+model: gfs
+cycle_interval: 6
+
+# obs-mon hera
+data: '/scratch1/NCEPDEV/da/Edward.Safford/noscrub/test_data/gfs'
+
+satellites:
+
+  - name: g18
+    instruments:
+      - name: abi
+        plot_list:
+          - plot: rad hist
+            times: 1
+            channels: 'all'
+            component: ges
+            run: gdas
diff --git a/parm/templates/radHist.yaml b/parm/templates/radHist.yaml
new file mode 100644
index 0000000..ebb39b1
--- /dev/null
+++ b/parm/templates/radHist.yaml
@@ -0,0 +1,73 @@
+#
+# Generate a histogram plot of omgnbc using legacy RadMon 
+# scan angle data for the specified cycle.
+#
+# This is plot template intended as a proof of concept and isn't 
+# a part of the legacy RadMon plots.  The legacy ConMon makes 
+# histogram # plots but the process of doing so is very cumbersome 
+# and won't be made part of the obs-monitor package.  Instead
+# histograms will be employed when conventional data from JEDI
+# is made available.  Until then the legacy histogram plots will 
+# be made with the legacy GSI-monitor components.  Once JEDI data
+# becomes availble this plot template can be adapted to make the
+# desired histograms with conventional (or any other) data.
+#
+
+
+# Data read
+#-----------
+
+datasets:
+  - name: angle
+    satellite: {{SAT}}
+    sensor: {{SENSOR}}
+    type: MonDataSpace
+    control_file:
+      - {{DATA}}/angle.{{SENSOR}}_{{SAT}}.ctl
+    filenames:
+      - {{ DATA }}/angle.{{ SENSOR }}_{{ SAT }}.{{ PDATE | to_YMDH }}.ieee_d
+
+    channels: {{CHANNELS}}
+    regions: &regions 1
+    groups:
+      - name: GsiIeee
+        variables: &variables ['omgnbc']
+
+graphics:
+
+  plotting_backend: Emcpy
+  figure_list:
+
+    # Obs count plots
+    # ---------------
+    - batch figure:
+        channels: {{CHANNELS}}
+        variables: ['omgbc']
+        sensor: {{SENSOR}}
+
+      figure:
+        layout: [1,1]
+        figure size: [20,18]
+        title: "${variable}, {{SENSOR}}_{{SAT}} channel ${channel}  \n Valid: {{ PDATE | to_YMDH }}"
+        output name: hist_plots/rad/hist/hist.{{SENSOR}}_{{SAT}}.${channel}.${variable}.png
+        plot logo:
+          which: 'noaa/nws'
+          loc: 'upper right'
+
+
+      plots:
+        - add_xlabel: 'omgbc'
+          statistics:
+            fields:
+              - field_name: angle::GsiIeee::omgbc
+                channel: ${channel}
+            statistics_variables:
+            - n
+          layers:
+          - type: Histogram
+            data:
+              variable: angle::GsiIeee::omgbc
+              channel: ${channel}
+            color: 'blue'
+            bins: 100
+            alpha: 0.5        

From dbe6cafcd635073db7b1ba4db637611506081010 Mon Sep 17 00:00:00 2001
From: "edward.safford" <edward.safford@noaa.gov>
Date: Fri, 14 Feb 2025 19:01:39 +0000
Subject: [PATCH 2/2] Ref #26

Clarify comments.
---
 parm/gfs/test_hist.yaml     | 6 +++---
 parm/templates/radHist.yaml | 2 +-
 2 files changed, 4 insertions(+), 4 deletions(-)

diff --git a/parm/gfs/test_hist.yaml b/parm/gfs/test_hist.yaml
index be4ccc1..c0f7691 100644
--- a/parm/gfs/test_hist.yaml
+++ b/parm/gfs/test_hist.yaml
@@ -1,10 +1,10 @@
 # DA Monitoring source dictionary
 
-# This is a simple script to create a radiance plot using the radHist 
-# (histogram) plot template.
+# This is a simple script to create a histogram plot using legacy 
+# GSI-monitor radiance scan angle data.
 #
 # It is intended as a proof of concept and isn't an actual part of the 
-# legacy RadMon.
+# legacy RadMon output.
 #
 
 model: gfs
diff --git a/parm/templates/radHist.yaml b/parm/templates/radHist.yaml
index ebb39b1..fe6486a 100644
--- a/parm/templates/radHist.yaml
+++ b/parm/templates/radHist.yaml
@@ -48,7 +48,7 @@ graphics:
       figure:
         layout: [1,1]
         figure size: [20,18]
-        title: "${variable}, {{SENSOR}}_{{SAT}} channel ${channel}  \n Valid: {{ PDATE | to_YMDH }}"
+        title: "${variable} of scan angles, {{SENSOR}}_{{SAT}} channel ${channel}  \n Valid: {{ PDATE | to_YMDH }}"
         output name: hist_plots/rad/hist/hist.{{SENSOR}}_{{SAT}}.${channel}.${variable}.png
         plot logo:
           which: 'noaa/nws'