# == Target grid for remapping == # WOA13 grid targetz = [ 0, 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 100, 125, 150, 175, 200, 225, 250, 275, 300, 325, 350, 375, 400, 425, 450, 475, 500, 550, 600, 650, 700, 750, 800, 850, 900, 950, 1000, 1050, 1100, 1150, 1200, 1250, 1300, 1350, 1400, 1450, 1500, 1550, 1600, 1650, 1700, 1750, 1800, 1850, 1900, 1950, 2000, 2100, 2200, 2300, 2400, 2500, 2600, 2700, 2800, 2900, 3000, 3100, 3200, 3300, 3400, 3500, 3600, 3700, 3800, 3900, 4000, 4100, 4200, 4300, 4400, 4500, 4600, 4700, 4800, 4900, 5000, 5100, 5200, 5300, 5400, 5500 ] # == Remap == fieldz, zbowl = remaptoz(vardiffr, pseudo_depth, targetz) print('field remapped') density_z, zbowl2 = remaptoz(density3d, pseudo_depth, targetz) print('density remapped') if name == 'mme_hist' and modelAgree: var_agreez, zbowl3 = remaptoz(var_agreer, pseudo_depth, targetz) # import xarray as xr # density_z = xr.DataArray(density_z,dims=['basin','pseudo-depth','latitude'], # coords=[[0,1,2,3],targetz,lat]) # density_z = density_z.where(xr.DataArray(targetz,dims=['pseudo-depth'],coords=[targetz])>=zbowl) # -- Make variable bundles for each basin varAtl = { 'name': 'Atlantic', 'var_change': fieldz[1, :, :],
pseudo_depth = pickle.load(open(indir_z + file_z, "rb")) # -- Target grid for remapping - WOA13 grid -- targetz = [ 0, 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 100, 125, 150, 175, 200, 225, 250, 275, 300, 325, 350, 375, 400, 425, 450, 475, 500, 550, 600, 650, 700, 750, 800, 850, 900, 950, 1000, 1050, 1100, 1150, 1200, 1250, 1300, 1350, 1400, 1450, 1500, 1550, 1600, 1650, 1700, 1750, 1800, 1850, 1900, 1950, 2000, 2100, 2200, 2300, 2400, 2500, 2600, 2700, 2800, 2900, 3000, 3100, 3200, 3300, 3400, 3500, 3600, 3700, 3800, 3900, 4000, 4100, 4200, 4300, 4400, 4500, 4600, 4700, 4800, 4900, 5000, 5100, 5200, 5300, 5400, 5500 ] # -- Remap -- medianToEz, zbowl = remaptoz(medianToEr, pseudo_depth, targetz) medianToEz = np.ma.around(medianToEz) print('Median remapping done') rangeToEz, zbowlbis = remaptoz(rangeToEr, pseudo_depth, targetz) rangeToEz = np.ma.around(rangeToEz) print('Range remapping done') noagreez, zbowlbis = remaptoz(noagreer, pseudo_depth, targetz) print('noagree remapping done') norangez, zbowlbis = remaptoz(noranger, pseudo_depth, targetz) print('norange remapping done') density_z, zbowlbis = remaptoz(density3d, pseudo_depth, targetz) print('density remapping done') # -- Make variable bundles for each basin varAtlmedian = { 'name': 'Atlantic',
# == Target grid for remapping == # WOA13 grid targetz = [ 0, 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 100, 125, 150, 175, 200, 225, 250, 275, 300, 325, 350, 375, 400, 425, 450, 475, 500, 550, 600, 650, 700, 750, 800, 850, 900, 950, 1000, 1050, 1100, 1150, 1200, 1250, 1300, 1350, 1400, 1450, 1500, 1550, 1600, 1650, 1700, 1750, 1800, 1850, 1900, 1950, 2000, 2100, 2200, 2300, 2400, 2500, 2600, 2700, 2800, 2900, 3000, 3100, 3200, 3300, 3400, 3500, 3600, 3700, 3800, 3900, 4000, 4100, 4200, 4300, 4400, 4500, 4600, 4700, 4800, 4900, 5000, 5100, 5200, 5300, 5400, 5500 ] # == Remap == fieldz = remaptoz(vardiffr, pseudo_depth, targetz) print('field remapped') density_z = remaptoz(density3d, pseudo_depth, targetz) print('density remapped') if name == 'mme_hist' and modelAgree: var_agreez = remaptoz(var_agreer, pseudo_depth, targetz) # -- Make variable bundles for each basin varAtl = { 'name': 'Atlantic', 'var_change': fieldz[1, :, :], 'bowl1': None, 'bowl2': None, 'labBowl': None, 'density': density_z[1, :, :] }