Skip to content

Commit

Permalink
Browse files Browse the repository at this point in the history
Add radHist template and driver script.
  • Loading branch information
EdwardSafford-NOAA committed Feb 14, 2025
1 parent 92520fa commit 3c0c57f
Show file tree
Hide file tree
Showing 2 changed files with 99 additions and 0 deletions.
26 changes: 26 additions & 0 deletions parm/gfs/test_hist.yaml
Original file line number Diff line number Diff line change
@@ -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
73 changes: 73 additions & 0 deletions parm/templates/radHist.yaml
Original file line number Diff line number Diff line change
@@ -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

0 comments on commit 3c0c57f

Please sign in to comment.