Add comparison to 2002-2021 monthly climatology of MOPITT carbon monoxide
New feature type
New plot and/or plot enhancement
What is this new feature?
Compare model output against satellite-based measurements of carbon monoxide. Requires horizontal and vertical regridding to the observation space and applying averaging kernels.
Attached file is an example comparison for December from the old Diagnostics (using MOPITT CO climatology 2003-2012).

Assistance required?
Yes, I will need some assistance before making a PR
Extra info
Chemistry diagnostic
- Dataset details: MOPITT V9J 2002-2021
- Variables to load: CO column, a priori, averaging kernel
- Location: /glade/campaign/acom/acom-da/buchholz/MOPITT_v9_climo
Applied for data archiving through NCAR GDEX, Data Reference: https://gdex.ucar.edu/dataset/369_buchholz.html
Update: we have a working Jupyter Notebook comparison code. Some notes:
- Currently plots all 4 options (MOPITT, MOPITT SD, Smoothed model, Difference) for one month. You can choose the month at the start of the plotting cell, will turn need a loop (or parallel option) to compare all months.
- Compares one year of model output - we will need to create climatology of CO from the CESM/CAM-chem output.
- Only works with the finite volume output currently.
- Horizontal regridding is xarray linear interpolation. There might be a better/more accurate regridding through xesmf.
- Consider alternate colormaps.
Image below compares to /glade/scratch/tilmes/archive/FCnudged_MAM4_f09.carma_trop_strat09.aqchem.2001_2020

Upated code is available on the UCAR HPC: /glade/work/shawnh/derecho/ADF_O3DIAG/scripts/plotting/MOPITT.py