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)
