f.savefig(os.path.join(opath, 'map2_ensemble.png')) uwnd = analogs(ens, 'ncep', 'uwnd_850').composite() f = scalar_plot(uwnd, test=0.05, proj='cyl').plot() f.savefig(os.path.join(opath, 'map3_ensemble.png')) # ============================================================================== """ CLIMATE INDICES """ f = indices(ens).plot() f.savefig(os.path.join(opath, 'indices_ensemble.png')) # ============================================================================== """ WEATHER REGIMES """ w = WR(ens, classification='New Zealand') w.probs_anomalies(kind='many') f = w.plot_bar() f.savefig(os.path.join(opath, 'WR_ensemble.png'))
f.savefig(os.path.join(opath, 'indices_proxy.png')) images.append({'id': 'indices_proxy', 'title' : 'Climate Indices', 'filename': 'indices_proxy.png'}) plt.close(f) if verbose: save_progress(opath, 'climate indices', 60) """ Weather Regimes """ w = WR(p, classification='New Zealand') f = w.plot_bar(sig=1) f.savefig(os.path.join(opath, 'NZ_regimes_proxy.png')) images.append({'id': 'NZ_regimes_proxy', 'title' : 'NZ weather regimes (Kidson Types)', 'filename': 'NZ_regimes_proxy.png'}) plt.close(f) if verbose: save_progress(opath, 'NZ weather regimes', 70) w = WR(p, classification='SW Pacific') f = w.plot_bar(sig=1)