BBL pressure values of the calculated Bottom Boundary Layer (exact position depends on the criteria)
list_of_BBL_indices corresponding index (eg for interp_pressure or other arrays of the same size)
BBL_range pressure values of "height_above_ground" meters. Follows therefore the batyhmetrie.
list_of_BBL_range_indices corresponding index (eg for interp_pressure or other arrays of the same size)
"""
bathymetrie, list_of_bathymetrie_indices = results[0]
#BBL,list_of_BBL_indices = results[1] #not needed here
BBL_range, list_of_BBL_range_indices = results[2]
eps_N_grid = np.sqrt(eps_N_squared_grid)
#ozmidov scale
ozmidov_scale_grid = np.sqrt(eps_grid / (eps_N_grid**3))
#conversion from pressure coordinates to depth
eps_depth = gsw.z_from_p(
eps_pressure, np.mean(lat)
) #mean lat should be sufficient, because the transect is east-west
bathymetrie_in_m = gsw.z_from_p(bathymetrie, np.mean(lat))
eps_depth_grid = np.reshape(eps_depth,
(1, -1)) * np.ones(np.shape(eps_grid))
distance_from_ground_grid = eps_depth_grid - np.reshape(
bathymetrie_in_m, (-1, 1))
boundary_check_grid = ~(distance_from_ground_grid < ozmidov_scale_grid)
oxygen_flux_Osborn_grid = thesis.get_oxygen_flux_osborn(
eps_Reynolds_bouyancy_grid, eps_grid, eps_N_squared_grid,
eps_oxygen_grid, eps_depth, eps_density_grid)
oxygen_flux_BB_grid = thesis.get_oxygen_flux_BB(
eps_Reynolds_bouyancy_grid, eps_grid, eps_N_squared_grid,
corrected_eps_Reynolds_bouyancy_grid
eps_wiki_Reynolds_bouyancy_grid
corrected_eps_wiki_Reynolds_bouyancy_grid
"""
print("Number of profiles:", number_of_profiles)
results = thesis.find_bottom_and_bottom_currents(
number_of_profiles, eps_pressure, eps_density_grid,
eps_oxygen_grid)
bathymetry, list_of_bathymetrie_indices = results[0]
#conversion from pressure coordinates to depth
eps_depth = gswc.z_from_p(
eps_pressure, np.mean(lat)
) #mean lat should be sufficient, because the transect is east-west
#get oxygen fluxes
oxygen_flux_osborn_grid = thesis.get_oxygen_flux_osborn(
eps_Reynolds_bouyancy_grid, eps_grid, eps_N_squared_grid,
eps_oxygen_grid, eps_depth, eps_density_grid)
oxygen_flux_BB_grid = thesis.get_oxygen_flux_BB(
eps_Reynolds_bouyancy_grid, eps_grid, eps_N_squared_grid,
eps_oxygen_grid, eps_depth, eps_density_grid)
oxygen_flux_Skif_grid = thesis.get_oxygen_flux_skif(
eps_Reynolds_bouyancy_grid, eps_grid, eps_N_squared_grid,
eps_oxygen_grid, eps_depth, eps_density_grid)
example_profile_index = [
np.argmin(np.abs(lon - 20.55)),
