figtitle = "ECHAM4.6 x CMAP\nCorrel de Prec Acum (mm) - %s %s" % (target_months, hind_period) # figname = "correl_echam46_issue_{0}_{1}.{2}.png".format(fcst_month, target_months, hind_period) # pm.correl_map_shaded(correl, lons_correl_aux-2.8125/2., lats-2.8125/2., fig_name=figname, fig_title=figtitle) # Plot map correlacao globo sem oceano # figname = "correl_echam46_issue_{0}_{1}.{2}_noocean.png".format(fcst_month, target_months, hind_period) # pm.correl_map_shaded(correl, lons_correl_aux-2.8125/2., lats-2.8125/2., fig_name=figname, fig_title=figtitle, ocean_mask=1) # Plot map correlacao Brasil # figname = "correl_echam46_issue_{0}_{1}.{2}_br.png".format(fcst_month, target_months, hind_period) # pm.correl_map_shaded(correl, lons_correl_aux-2.8125/2., lats-2.8125/2., fig_name=figname, fig_title=figtitle, br_country=1) # Plot map correlacao Brasil sem oceano # figname = "correl_echam46_issue_{0}_{1}.{2}_noocean_br.png".format(fcst_month, target_months, hind_period) figname = outdir + '/' + "bra_precip_persistida_hind%s_%s_null_%s_echam46_1dg_cmap_correlacao.png" % (hind_period, target_months, n_fcst_month) pm.correl_map_shaded(correl, lons_correl_aux-2.8125/2., lats-2.8125/2., fig_name=figname, fig_title=figtitle, ocean_mask=1, br_country=1) # Anomalia anomaly, anomaly_pad = cs.compute_anomaly(fsct, hind) # Shift na longitude e nos dados do modelo # Isso é necessário pois a longitude do NetCDF é no formato 0 - 360 # a função aceita valores de longitude no formato -180:+180 # Chefinho quer ver todos os oceanos # anomaly, fsct_lon = shiftgrid(30., anomaly, fsct_lon, start=False) anomaly, lons_anom_aux = shiftgrid(180., anomaly, lons, start=False) # Plot map anomalia globo # a variavel a ser plotada deve ser 2d figtitle = "ECHAM4.6 - Anomalia de Prec Acum (mm)\n{0} {1} - {2} {3}".format(fcst_month, fcst_year, target_months, target_year, hind_period) # figname = "anomalia_echam46_issue_{0}{1}_target_{2}{3}.{4}.png".format(fcst_month, fcst_year, target_months, target_year, hind_period)