def parallel_ops_budg(s, bmk_mon):
ref = get_filepath(gcc_vs_gcc_refdir, col, bmk_mon)
dev = get_filepath(gcc_vs_gcc_devdir, col, bmk_mon)
plot_dir = join(gcc_vs_gcc_budgetdir, bmk_mon_yr_strs[s])
bmk.make_benchmark_operations_budget(gcc_ref_version,
ref,
gcc_dev_version,
dev,
bmk_sec_per_month[s],
benchmark_type=bmk_type,
label=bmk_mon_yr_strs[s],
dst=gcc_vs_gcc_tablesdir)
def parallel_mass_table(s, bmk_mon):
ref = get_filepath(gcc_vs_gcc_refrstdir, col, bmk_mon)
dev = get_filepath(gcc_vs_gcc_devrstdir, col, bmk_mon)
bmk.make_benchmark_mass_tables(ref,
gcc_vs_gcc_refstr,
dev,
gcc_vs_gcc_devstr,
dst=gcc_vs_gcc_tablesdir,
subdst=bmk_mon_yr_strs[s],
label="at 01{}".format(
bmk_mon_yr_strs[s]),
overwrite=True,
spcdb_dir=spcdb_dir)
def parallel_ops_budg(s, bmk_mon):
# Create budget table for each benchmark month (ewl??)
ref = get_filepath(gchp_vs_gcc_refdir, col, bmk_mon)
dev = get_filepath(gchp_vs_gcc_devdir,
col,
bmk_mons_mid[s],
is_gchp=True)
plot_dir = join(gchp_vs_gcc_budgetdir, bmk_mon_yr_strs[s])
bmk.make_benchmark_operations_budget(gcc_ref_version,
ref,
gchp_dev_version,
dev,
bmk_sec_per_month[s],
benchmark_type=bmk_type,
label=bmk_mon_yr_strs[s],
operations=[
"Chemistry", "Convection",
"EmisDryDep", "Mixing",
"WetDep"
],
compute_accum=False,
dst=gcc_vs_gcc_tablesdir)
# Diagnostic collections to read
col = "SpeciesConc"
colmet = "StateMet"
colmet_gchp = "StateMet_avg" # Use this for benchmarks prior to 13.0
# Create concentration plots for each benchmark month
for t in range(bmk_n_months):
# Time & date quantities
reftime = bmk_mons_ref[t]
devtime = bmk_mons_dev[t]
datestr = bmk_mon_yr_strs_dev[t]
# Seasonal diagnostic collection files to read
ref = get_filepath(gcc_vs_gcc_refdir, col, reftime)
dev = get_filepath(gcc_vs_gcc_devdir, col, devtime)
refmet = get_filepath(gcc_vs_gcc_refdir, colmet, reftime)
devmet = get_filepath(gcc_vs_gcc_devdir, colmet, devtime)
bmk.make_benchmark_conc_plots(ref,
gcc_vs_gcc_refstr,
dev,
gcc_vs_gcc_devstr,
refmet=refmet,
devmet=devmet,
dst=gcc_vs_gcc_resultsdir,
subdst=datestr,
weightsdir=weightsdir,
benchmark_type=bmk_type,
restrict_cats=restrict_cats,
# Create GCC vs GCC benchmark plots and tables
# ======================================================================
if gcc_vs_gcc:
#---------------------------------------------------------------
# GCC vs GCC Concentration plots
#
# Includes lumped species and separates by category if plot_by_spc_cat
# is true; otherwise excludes lumped species and writes to one file
#---------------------------------------------------------------
if plot_conc:
title = "\n%%% Creating GCC vs. GCC concentration plots %%%"
# Diagnostic collection files to read
col = "SpeciesConc"
ref = get_filepath(gcc_vs_gcc_refdir, col, gcc_date)
dev = get_filepath(gcc_vs_gcc_devdir, col, gcc_date)
# Meteorology data needed for calculations
colmet = "StateMet"
refmet = get_filepath(gcc_vs_gcc_refdir, colmet, gcc_date)
devmet = get_filepath(gcc_vs_gcc_devdir, colmet, gcc_date)
# Make concentration plots
bmk.make_benchmark_conc_plots(ref,
gcc_vs_gcc_refstr,
dev,
gcc_vs_gcc_devstr,
refmet=refmet,
devmet=devmet,
dst=gcc_vs_gcc_resultsdir,
# GCC vs GCC Concentration plots
#
# Includes lumped species and separates by category if plot_by_spc_cat
# is true; otherwise excludes lumped species and writes to one file.
# --------------------------------------------------------------
if plot_conc:
print("\n%%% Creating GCC vs. GCC concentration plots %%%")
# Diagnostic collections to read
col = "SpeciesConc"
colmet = "StateMet"
# Create concentration plots for each benchmark month
for s, bmk_mon in enumerate(bmk_mons):
ref = get_filepath(gcc_vs_gcc_refdir, col, bmk_mon)
dev = get_filepath(gcc_vs_gcc_devdir, col, bmk_mon)
refmet = get_filepath(gcc_vs_gcc_refdir, colmet, bmk_mon)
devmet = get_filepath(gcc_vs_gcc_devdir, colmet, bmk_mon)
bmk.make_benchmark_conc_plots(ref,
gcc_vs_gcc_refstr,
dev,
gcc_vs_gcc_devstr,
refmet=refmet,
devmet=devmet,
dst=gcc_vs_gcc_resultsdir,
subdst=bmk_mon_yr_strs[s],
weightsdir=weightsdir,
benchmark_type=bmk_type,
plot_by_spc_cat=plot_by_spc_cat,
overwrite=True,
