Add script to plot QC Flags on map #83
Conversation
kevindougherty-noaa
requested review from
xyzemc,
BenjaminBlake-NOAA and
CoryMartin-NOAA
| import plot_features as features | ||
|
|
||
|
|
||
| def no_data_map(plotmap, domain, metadata): |
There was a problem hiding this comment.
I think this should be moved to a more generic part of emcpy as there may be other uses for this outside of LAMDA
There was a problem hiding this comment.
I am going to move it up one directory because it could be useful in other parts of PyGSI as I plan to add other examples that already exist in pyGSI/plot_diags.py
LAMDA/map_qc_flags.py
Outdated
| plotdir : (str; default='./') path to where figures should be saved | ||
| """ | ||
|
|
||
| diag = Radiance(inputfile) |
There was a problem hiding this comment.
Can we make this generic for all obs types? What if we want to plot QC values of conventional data (say satwinds)? Will it work to just pass the diag object to this function, or do we need to handle all three classes here?
There was a problem hiding this comment.
That is a good question. To my understanding, there are no qc flags for conventional data, only assimilated, rejected and monitored data. So this script is pretty specific to just Radiance data.
There was a problem hiding this comment.
There are QC flags for conventional data for sure, I think there are two Setup_QC_Mark and Prep_QC_Mark but am not totally sure of what should be used (probably Setup)
Creates a 2D map that plots the location of the qc flags. Addresses number 2.), 4th bullet in #84. An example is shown here:
Need to create a yaml and driver script to run.