except:
values.append(data)
ocean_mask_varname = 'oceanmask'
if ocean_mask_varname in nc.variables.keys():
ocean_mask.append(np.squeeze(ppt.permute(nc.variables['mask'])))
else:
ocean_mask.append(np.zeros_like(data))
nc.close()
# set the print mode
if print_mode in 'height':
ntn = nt
if (ntn == 2):
lw, pad_inches = ppt.set_mode(print_mode, aspect_ratio=.75)
if (ntn == 3):
lw, pad_inches = ppt.set_mode(print_mode, aspect_ratio=.55)
elif (ntn == 4):
lw, pad_inches = ppt.set_mode(print_mode, aspect_ratio=.35)
elif (ntn == 5):
lw, pad_inches = ppt.set_mode(print_mode, aspect_ratio=.25)
else:
lw, pad_inches = ppt.set_mode(print_mode)
else:
lw, pad_inches = ppt.set_mode(print_mode)
# make a separate colorbar (if requested)
if colorbar:
fig = plt.figure()
if obs_file:
obs = ObservationsDataset(obs_file, varname)
for flux_gate in flux_gates:
flux_gate.add_observations(obs)
# Add experiments to flux gates
for k, filename in enumerate(args):
# id = re.search("id_(\b0*([1-9][0-9]*|0)\b)", filename).group(1)
# pid = int(filename.split("id_")[1].split("_")[0])
id = k
experiment = ExperimentDataset(id, filename, varname)
for flux_gate in flux_gates:
flux_gate.add_experiment(experiment)
# set the print mode
lw, pad_inches = ppt.set_mode(print_mode, aspect_ratio=aspect_ratio)
ne = len(flux_gates[0].experiments)
ng = len(flux_gates)
if table_file and obs_file:
export_latex_table_flux(table_file, flux_gates, label_params)
# make figure for each flux gate
for gate in flux_gates:
if make_figures:
gate.make_line_plot(label_param_list=label_params)
else:
if not gate.has_fluxes:
gate.calculate_fluxes()
if gate.has_observations:
values.append(np.ma.array(data, mask=mask))
except:
values.append(data)
ocean_mask_varname = "oceanmask"
if ocean_mask_varname in list(nc.variables.keys()):
ocean_mask.append(np.squeeze(ppt.permute(nc.variables["mask"])))
else:
ocean_mask.append(np.zeros_like(data))
nc.close()
# set the print mode
if print_mode in "height":
ntn = nt
if ntn == 2:
lw, pad_inches = ppt.set_mode(print_mode, aspect_ratio=0.75)
if ntn == 3:
lw, pad_inches = ppt.set_mode(print_mode, aspect_ratio=0.55)
elif ntn == 4:
lw, pad_inches = ppt.set_mode(print_mode, aspect_ratio=0.35)
elif ntn == 5:
lw, pad_inches = ppt.set_mode(print_mode, aspect_ratio=0.25)
else:
lw, pad_inches = ppt.set_mode(print_mode)
else:
lw, pad_inches = ppt.set_mode(print_mode)
# make a separate colorbar (if requested)
if colorbar:
fig = plt.figure()
# import and convert colormap
cdict = ppt.gmtColormap(colormap)
cmap = colors.LinearSegmentedColormap('my_colormap', cdict)
# Init Unit system
sys = System()
# Plotting styles
axisbg = '0.9'
shadow_color = '0.25'
numpoints = 1
aspect_ratio = golden_mean
# set the print mode
lw, pad_inches = ppt.set_mode(print_mode, aspect_ratio=aspect_ratio)
plt.rcParams['legend.fancybox'] = True
ne = len(args)
cNorm = colors.Normalize(vmin=0, vmax=ne)
scalarMap = cmx.ScalarMappable(norm=cNorm, cmap=cmap)
profile_axis_o_units = 'm'
filename = args[0]
try:
nc0 = NC(filename, 'r')
except: