ADF

A number of plots are provided from ADF. The full output from the stand-alone ADF configuration is in the link below.

Note that in standalone format (eg, CUPiD run not through CESM workflow), ADF is currently run by users via the following process:

1. Install ADF and activate cupid-analysis

2. Use the CUPiD/helper_scripts/generate_adf_config_file.py script to generate an ADF config file based on a CUPiD configuration file.

   • cd CUPiD/examples/external_diag_packages

   • ../../helper_scripts/generate_adf_config_file.py --cupid-config-loc . --adf-template ../../externals/ADF/config_amwg_default_plots.yaml --out-file ADF_config.yaml

3. Run ADF with the newly created configuration file.

   • ../../externals/ADF/run_adf_diag ADF_config.yaml

Show code cell content

Hide code cell content

import os

from IPython.core.display import HTML, Image
from IPython.display import display
import pandas as pd

Show code cell content

Hide code cell content

adf_root = "."
case_name = None
base_case_name = None
start_date = ""
end_date = ""
base_start_date = None
base_end_date = None
key_plots = None
compare_obs: False
# adf_root will be external_diag_packages/computed_notebooks/ADF/

# Parameters
case_name = "b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.313"
base_case_name = "b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.308"
case_nickname = "BLT1850_313"
base_case_nickname = "BLT1850_308"
CESM_output_dir = "/glade/campaign/cesm/development/cross-wg/diagnostic_framework/CESM_output_for_testing"
start_date = "0001-01-01"
end_date = "0044-01-01"
climo_start_year = 24
climo_end_year = 44
base_start_date = "0001-01-01"
base_end_date = "0105-01-01"
base_climo_start_year = 85
base_climo_end_year = 105
obs_data_dir = (
    "/glade/campaign/cesm/development/cross-wg/diagnostic_framework/CUPiD_obs_data"
)
ts_dir = None
lc_kwargs = {"threads_per_worker": 1}
serial = False
adf_root = "../../examples/313v308/ADF_output/"
key_plots = [
    "Surface_Wind_Stress_ANN_LatLon_Vector_Mean.png",
    "PRECT_ANN_LatLon_Mean.png",
    "PS_DJF_SHPolar_Mean.png",
    "TaylorDiag_ANN_Special_Mean.png",
]
compare_obs = False
subset_kwargs = {}
product = "/glade/u/home/dsanders/cesm_dev/CUPiD/examples/313v308/computed_notebooks//atm/ADF.ipynb"

Show code cell content

Hide code cell content

# Want some base case parameter defaults to equal control case values
if base_case_name is not None:
    if base_start_date is None:
        base_start_date = start_date

    if base_end_date is None:
        base_end_date = end_date
if compare_obs:
    base_case_name = "Obs"
    base_start_date = None
    base_end_date = None

Show code cell content

Hide code cell content

# convert start-date and end-date to year range
case_year_range = [int(start_date.split("-")[0]), int(end_date.split("-")[0])]

if (not base_start_date) and (not base_end_date):
    base_case_year_range = None
else:
    base_case_year_range = [
        int(base_start_date.split("-")[0]),
        int(base_end_date.split("-")[0]),
    ]

if base_case_year_range:
    base_case_yr_range_str = f"_{base_case_year_range[0]}_{base_case_year_range[1]}"
    alt_base_case_yr_range_str = (
        f"_{base_case_year_range[0]}_{str(int(base_case_year_range[1])-1)}"
    )
else:
    base_case_yr_range_str = ""
    alt_base_case_yr_range_str = ""

possible_adf_comparison_names = [
    f"{case_name}_{case_year_range[0]}_{case_year_range[1]}_vs_{base_case_name}{base_case_yr_range_str}"
]
possible_adf_comparison_names.append(
    f"{case_name}_{case_year_range[0]}_{str(int(case_year_range[1])-1)}_vs_{base_case_name}{base_case_yr_range_str}"
)
possible_adf_comparison_names.append(
    f"{case_name}_{case_year_range[0]}_{str(int(case_year_range[1])-1)}_vs_{base_case_name}{alt_base_case_yr_range_str}"
)
possible_adf_comparison_names.append(
    f"{case_name}_{case_year_range[0]}_{case_year_range[1]}_vs_{base_case_name}{alt_base_case_yr_range_str}"
)

Show code cell source

Hide code cell source

adf_found = False
for adf_comparison_name in possible_adf_comparison_names:
    adf_root_candidate = os.path.join(adf_root, adf_comparison_name)
    if os.path.exists(adf_root_candidate):
        adf_found = True
        adf_root = adf_root_candidate
        display(
            HTML(
                f'<a href="../ADF/{adf_comparison_name}/website/index.html" target="_blank" style="font-size: 30px">Full ADF output</a>'
            )
        )
        break

if not adf_found:
    print("No ADF output found for the specified case and date range.")

Full ADF output

Key Metrics from ADF

Some important things to look at from ADF include a comparison table and a few maps:

Show code cell source

Hide code cell source

comparison_table = os.path.join(adf_root, "amwg_table_comp.csv")
if os.path.isfile(comparison_table):
    table = pd.read_csv(comparison_table)
    display(HTML(table.to_html(index=False, float_format="{:6g}".format)))

Show code cell source

Hide code cell source

for path_to_key_plot in key_plots:
    full_path = os.path.join(adf_root, path_to_key_plot)
    if os.path.isfile(full_path):
        display(Image(full_path))