Make time of emergence PDF bar plots from historical vs. historicalNat Choose which variable and domain to work on Method1: read lat/rho ToE and average in the chosen domain Method2: read ToE computed from the already averaged signal and noise in the chosen domain """ import numpy as np import matplotlib.pyplot as plt from netCDF4 import Dataset as open_ncfile from maps_matplot_lib import defVarmme, averageDom from modelsDef import defModels from libToE import findToE, ToEdomainhistvshistNat # ----- Work ------ models = defModels() varname = defVarmme('salinity'); v = 'S' #varname = defVarmme('temp'); v = 'T' #varname = defVarmme('depth'); v = 'Z' # -- Choose method for computing ToE #method = 'average_ToE' # Determine 2D lat/rho ToE then average in the box method = 'average_signal' # Average signal and noise in the box, then compute ToE # Method 2 only has ToE for salinity for now # -- Choose which 'noise' to use for the ToE calculation (only with method 2 : average_signal) # method_noise = 'average_std' # Average the standard deviation of PiC in the specified domains method_noise = 'average_histNat' # Average histNat in the specified domains then determine the std of this averaged value
import matplotlib.pyplot as plt from maps_matplot_lib import defVarmme, averageDom from modelsDef import defModels from libToE import ToEdomainhistvshistNat from matplotlib.ticker import AutoMinorLocator, MultipleLocator import glob # ----- Workspace ------ indir_histrcp85 = '/home/ysilvy/Density_bining/Yona_analysis/data/hist_rcp85/' indir_histNat = '/data/ericglod/Density_binning/Prod_density_april15/mme_histNat/' indir_piC = '/data/ericglod/Density_binning/Prod_density_april15/mme_piControl/' indir_noise = '/home/ysilvy/Density_bining/Yona_analysis/data/noise_estimate/' models = defModels() # ----- Work ------ varname = defVarmme('salinity'); v = 'S' method = 'average_signal' # Average signal and noise in the box, then compute ToE # -- Choose which 'noise' to use # method_noise = 'average_std' # Average the standard deviation of histNat or PiC in the specified domains method_noise = 'average_histNat' # Average histNat or PiC in the specified domains then determine the std # of this averaged value domains = ['Southern ST', 'SO', 'Northern ST', 'North Atlantic', 'North Pacific'] domain_name = 'Southern ST' idomain = 0