# == 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, :, :]
}