f = scalar_plot(uwnd, test=0.05, proj='cyl').plot() 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 = scalar_plot(uwnd, test=0.05, proj='cyl').plot() 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,'hgt_1000_proxy_NZ_years.png')) images.append({'id': 'hgt_1000_NZ_samples', 'title' : 'Geopotential at 1000 hPa, NZ domain, analog years', 'filename': 'HGT_1000_NZ_sample_proxy.png'}) plt.close(f) if verbose: save_progress(opath, 'HGT 1000 NZ domain analogs', 50) # ============================================================================== """ CLIMATE INDICES """ f = indices(p).plot() 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 """