Пример #1
0
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")
Пример #2
0
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")