def time_serie_Arctic(path, var):
yr, xr, time, depth = loading.extracting_coord(path)
TorS = loading.extracting_var(path, var)
make_plot.points_on_map(xr, yr, var, basin)
S = loading.extracting_var(path, 'sal')
TorS[np.where(S == 0.)] = np.nan
i_zmax = np.max(np.where(depth < simu.zmax))
print(np.shape(TorS))
mean_Arctic = np.mean(np.nanmean(TorS[:, 0:i_zmax, :], axis=0), axis=0)
return time, depth, mean_Arctic
def time_serie_Arctic(path,var,lat_min,basin):
yr, xr, time, depth = loading.extracting_coord(path)
TorS = loading.extracting_var(path, var)
S = loading.extracting_var(path,'sal')
TorS[np.where( S==0. )] = np.nan
if basin =='undefined':
mask = loading.latitudinal_band_mask(yr,lat_min,90)
else:
mask = loading.Ocean_mask(xr,yr,basin)
make_plot.points_on_map(xr[mask],yr[mask],var,basin)
arctic = TorS[mask,:,:]
mean_Arctic = np.nanmean(arctic,axis=0)
return time, depth, mean_Arctic
self.vmin = 25
self.vmax = 28.1
self.output_file = str(var) + '_' + str(month) + str(
self.y1) + '_' + str(sufix)
self.output_file2 = str(var) + '_' + str(month) + str(
self.y1) + '_' + str(self.y2) + '_' + str(sufix)
#//
return #//
#//
#//////////////////////////////////////////////////////////////////////////////////////////
simu = yearly_LongPeriod(option, var, y1, y2, basin, lat_min, month)
yr, xr, time, depth = loading.extracting_coord(simu.path)
VarArray_simuRho = loading.extracting_var(simu.path,
'rho') # Practical Salinity
index_y1 = np.min(
np.where(simu.first_year + time[:] / (3600 * 24 * 364.5) > simu.y1))
index_y2 = np.min(
np.where(simu.first_year + time[:] / (3600 * 24 * 364.5) > y2))
if basin == 'undefined':
mask = loading.latitudinal_band_mask(yr, lat_min, 90)
else:
mask = loading.Ocean_mask(xr, yr, basin)
make_plot.points_on_map(xr[mask], yr[mask], var, basin)
VarArray_simuRho[np.where(VarArray_simuRho == 0)] = np.nan