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)),