Add comparison to 2002-2021 monthly climatology of MOPITT carbon monoxide

rrbuchholz commented 1 year ago

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). cset4_co_column_mopitt_12

Assistance required?

Yes, I will need some assistance before making a PR

Extra info

Chemistry diagnostic

rrbuchholz commented 1 year ago

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

rrbuchholz commented 1 year ago

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

rrbuchholz commented 4 months ago

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

NCAR / ADF