def oc_plots(oc_raw, oc_ast, year):
oc_raw_conc = h.map_ij(h.ppbv_to_conc, oc_raw)
oc_ast_conc = h.map_ij(h.ppbv_to_conc, oc_ast)
oc_log_ratio = h.map_ij(np.log, h.map_ij(h.divider(default=1), oc_ast_conc, oc_raw_conc))
h.plot_ij(oc_raw_conc,
filename=os.path.join(plt_dir,
"oc_raw_{}.png".format(year)),
title="GEOS-Chem Raw OC {}".format(year), cblabel="ug/m^3",
vmin=0, vmax=3.5)
h.plot_ij(oc_ast_conc,
filename=os.path.join(plt_dir,
"oc_ast_{}.png".format(year)),
title="GEOS-Chem Assimilated OC {}".format(year), cblabel="ug/m^3",
vmin=0, vmax=3.5)
h.plot_ij(oc_log_ratio,
filename=os.path.join(plt_dir,
"oc_ratio_{}.png".format(year)),
title="OC Ln(opt/prior) {}".format(year),
cblabel='',
vmin=-_oc_lim, vmax=_oc_lim,
cmap="RdBu")
def bc_plots(bc_raw, bc_ast, year):
bc_raw_conc = h.map_ij(h.ppbv_to_conc, bc_raw)
bc_ast_conc = h.map_ij(h.ppbv_to_conc, bc_ast)
bc_log_ratio = h.map_ij(np.log, h.map_ij(h.divider(default=1), bc_ast_conc, bc_raw_conc))
h.plot_ij(bc_raw_conc,
filename=os.path.join(plt_dir,
"bc_raw_{}".format(year)),
title="GEOS-Chem Raw BC {}".format(year), cblabel="ug/m^3",
vmin=0, vmax=0.6)
h.plot_ij(bc_ast_conc,
filename=os.path.join(plt_dir,
"bc_ast_{}.png".format(year)),
title="GEOS-Chem Assimilated BC {}".format(year), cblabel="ug/m^3",
vmin=0, vmax=0.6)
h.plot_ij(bc_log_ratio,
filename=os.path.join(plt_dir,
"bc_ratio_{}.png".format(year)),
title="BC Ln(opt/prior) {}".format(year),
cblabel='', norm=h.SymNorm(),
vmin=-_bc_lim, vmax=_bc_lim,
cmap="RdBu")