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
"""