Skip to content

Commit

Permalink
Add MinMon templates (NOAA-EMC#12)
Browse files Browse the repository at this point in the history 
Modify `plotObsMon.py `to support minmon plots and add templates to
produce:

- minmon summary plot
- minmon 4 cycle gnorm plot
- minmon single cycle gnorm plot

Please note this PR includes the same changes as NOAA-EMC#11 . I messed up and
started from it rather than from develop. So I'll make this a draft and
convert it to release if/when NOAA-EMC#11 is merged.

Resolves NOAA-EMC#8.
  • Loading branch information
@EdwardSafford-NOAA
EdwardSafford-NOAA authored Dec 26, 2023
1 parent ba89383 commit dafc0de
Show file tree
Hide file tree
Showing 5 changed files with 454 additions and 18 deletions.
14 changes: 14 additions & 0 deletions parm/gfs/gfs_plots.yaml
View file
Original file line number Diff line number Diff line change
Expand Up @@ -50,3 +50,17 @@ satellites:
times: 121
channels: '12,13,14,15,16'
run: ges

minimization:
- net: gfs
plot_list:
- plot: min summary
run: gdas

- plot: min gnorm four cycle
run: gdas
times: 28

- plot: min gnorm one cycle
run: gdas
times: 28
149 changes: 149 additions & 0 deletions parm/gfs/minGnormFourCycle.yaml
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,149 @@
#
# Gnorm 4 cycle plots for MinMon gfs
#
# Generate one plot:
# - gnorm values (all iterations) for last 4 cycles & 7 day avg


# Data read
# ---------
datasets:
- name: gnorm
type: MonDataSpace
control_file:
- {{DATA}}/min_data/glb_allgnorm.ctl
filenames:
- {{DATA}}/min_data/{{ PDATE | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm6 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm12 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm18 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm24 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm30 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm36 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm42 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm48 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm54 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm60 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm66 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm72 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm78 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm84 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm90 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm96 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm102 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm108 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm114 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm120 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm126 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm132 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm138 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm144 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm150 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm156 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm162 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm168 | to_YMDH }}.gnorms.ieee_d

groups:
- name: GsiIeee
variables: &variables ['gnorm']

transforms:
- transform: select time
new name: gnorm::GsiIeee::log_gnorm_cyc1
starting field: gnorm::GsiIeee::log_gnorm
cycle: {{ PDATE | to_YMDH }}
for:
variable: [none]

- transform: select time
new name: gnorm::GsiIeee::log_gnorm_cyc2
starting field: gnorm::GsiIeee::log_gnorm
cycle: {{ PDATEm6 | to_YMDH }}
for:
variable: [none]

- transform: select time
new name: gnorm::GsiIeee::log_gnorm_cyc3
starting field: gnorm::GsiIeee::log_gnorm
cycle: {{ PDATEm12 | to_YMDH }}
for:
variable: [none]

- transform: select time
new name: gnorm::GsiIeee::log_gnorm_cyc4
starting field: gnorm::GsiIeee::log_gnorm
cycle: {{ PDATEm18 | to_YMDH }}
for:
variable: [none]

- transform: select time
new name: gnorm::GsiIeee::log_gnorm_7d
starting field: gnorm::GsiIeee::log_gnorm
start cycle: {{ PDATE | to_YMDH }}
end cycle: {{ PDATEm162 | to_YMDH }}
for:
variable: [none]

graphics:

plotting_backend: Emcpy
figure_list:

# Gnorm single cycle & 7 day hourly mean
# --------------------------------------
- figure:
layout: [1,1]
figure size: [20,18]
tight layout:
title: "Valid: {{ PDATE | to_YMDH }}"
output name: line_plots/minimization/{{NET}}_{{RUN}}.4cycle.gnorms.png
plot logo:
which: 'noaa/nws'
loc: 'upper center'

plots:
- add_xlabel: 'Iteration Number'
add_ylabel: 'log( gnorm )'
add_grid:
axis: 'both'
linestyle: 'dotted'
linewidth: 0.5
color: 'black'
add_legend:
loc: 'upper right'
layers:
- type: LinePlot
x:
variable: gnorm::GsiIeee::iteration
y:
variable: gnorm::GsiIeee::log_gnorm_cyc1
color: 'blue'
label: '{{ PDATE | to_YMDH }}'
- type: LinePlot
x:
variable: gnorm::GsiIeee::iteration
y:
variable: gnorm::GsiIeee::log_gnorm_cyc2
color: 'red'
label: '{{ PDATEm6 | to_YMDH }}'
- type: LinePlot
x:
variable: gnorm::GsiIeee::iteration
y:
variable: gnorm::GsiIeee::log_gnorm_cyc3
color: 'green'
label: '{{ PDATEm12 | to_YMDH }}'
- type: LinePlot
x:
variable: gnorm::GsiIeee::iteration
y:
variable: gnorm::GsiIeee::log_gnorm_cyc4
color: 'orange'
label: '{{ PDATEm18 | to_YMDH }}'
- type: LinePlot
x:
variable: gnorm::GsiIeee::iteration
y:
variable: gnorm::GsiIeee::log_gnorm_7d
color: 'black'
label: '7 day mean'
142 changes: 142 additions & 0 deletions parm/gfs/minGnormOneCycle.yaml
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,142 @@
#
# Gnorm plots for MinMon gfs
#
# Generate two plots:
# - gnorm values (all iterations) for last 4 cycles & 7 day avg
# - single cycle gnorm values (all iterations) & 7 day avg


# Data read
# ---------
datasets:
- name: gnorm
type: MonDataSpace
control_file:
- {{DATA}}/min_data/glb_allgnorm.ctl
filenames:
- {{DATA}}/min_data/{{ PDATE | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm24 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm48 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm72 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm96 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm120 | to_YMDH }}.gnorms.ieee_d
- {{DATA}}/min_data/{{ PDATEm144 | to_YMDH }}.gnorms.ieee_d

groups:
- name: GsiIeee
variables: &variables ['gnorm']

transforms:
- transform: select time
new name: gnorm::GsiIeee::log_gnorm_1cyc
starting field: gnorm::GsiIeee::log_gnorm
cycle: {{ PDATE | to_YMDH }}
for:
variable: [none]

- transform: select time
new name: gnorm::GsiIeee::allgnorm_1cyc
starting field: gnorm::GsiIeee::allgnorm
cycle: {{ PDATE | to_YMDH }}
for:
variable: [none]

- transform: select time
new name: gnorm::GsiIeee::log_gnorm_7d
starting field: gnorm::GsiIeee::log_gnorm
start cycle: {{ PDATE | to_YMDH }}
end cycle: {{ PDATEm144 | to_YMDH }}
for:
variable: [none]

- transform: select time
new name: gnorm::GsiIeee::allgnorm_7d
starting field: gnorm::GsiIeee::allgnorm
start cycle: {{ PDATE | to_YMDH }}
end cycle: {{ PDATEm144 | to_YMDH }}
for:
variable: [none]


graphics:

plotting_backend: Emcpy
figure_list:

# Gnorm single cycle & 7 day hourly mean
# --------------------------------------
- figure:
layout: [1,1]
figure size: [20,18]
tight layout:
title: "Valid: {{ PDATE | to_YMDH }}"
output name: line_plots/minimization/{{NET}}_{{RUN}}.{{ PDATE | to_YMDH }}.gnorms.png
plot logo:
which: 'noaa/nws'
loc: 'upper right'

plots:
- add_xlabel: 'Iteration Number'
add_ylabel: 'log( gnorm )'
add_grid:
axis: 'both'
linestyle: 'dotted'
linewidth: 0.5
color: 'black'
add_legend:
loc: 'center right'
layers:
- type: LinePlot
x:
variable: gnorm::GsiIeee::iteration
y:
variable: gnorm::GsiIeee::log_gnorm_1cyc
color: 'blue'
label: '{{ PDATE | to_YMDH }}'
- type: LinePlot
x:
variable: gnorm::GsiIeee::iteration
y:
variable: gnorm::GsiIeee::log_gnorm_7d
color: 'gray'
label: '7 day mean of 06z cycles'

# Reduction single cycle & 7 day hourly mean
# ------------------------------------------
- figure:
layout: [1,1]
figure size: [20,18]
tight layout:
title: "Valid: {{ PDATE | to_YMDH }}"
output name: line_plots/minimization/{{NET}}_{{RUN}}.{{ PDATE | to_YMDH }}.reduction.png
plot logo:
which: 'noaa/nws'
loc: 'upper right'

plots:
- add_xlabel: 'Iteration Number'
add_ylabel: 'Reduction Value'
add_grid:
axis: 'both'
linestyle: 'dotted'
linewidth: 0.5
color: 'black'
add_legend:
loc: 'center right'

layers:
- type: LinePlot
x:
variable: gnorm::GsiIeee::iteration
y:
variable: gnorm::GsiIeee::allgnorm_1cyc
label: '{{ PDATE | to_YMDH }}'
color: 'blue'
- type: LinePlot
x:
variable: gnorm::GsiIeee::iteration
y:
variable: gnorm::GsiIeee::allgnorm_7d
color: 'gray'
label: '7 day mean of 06z cycles' # NEED {{HH}} value for 06

Loading

0 comments on commit dafc0de

Please sign in to comment.